High computational complexity due to rapidly increasing numerical stiffness is a difficult problem for simulating a supersonic reactive flow by using the uncoupled method. On the basis of our previous work, this paper proposes a dual adaptive method to ensure high calculation efficiency and good robustness in simulating stiff cases. The principle of this method is to realize adaptive coordination for the advection and reaction time steps in accordance with the non‐uniform feature of stiffness in the space and time dimensions. The proposed method can advance by a small time step in strong stiffness while with a large one in weak stiffness through the “prediction‐correction‐recovery” strategy. Some classical problems are chosen to verify the performance of the proposed method. The proposed method improved the computation efficiency by at least comparing with the previous method [1] and widened the error tolerance of the initial time step.