  Solving the reactive low-Mach Navier-Stokes equations with high-order
adaptive methods in time is still a challenging problem, in particular due to
the handling of the algebraic variables involved in the mass constraint. We
focus on the one-dimensional configuration, where this challenge has long
existed in the combustion community. We consider a model of solid propellant
combustion, which possesses the characteristic difficulties encountered in the
homogeneous or spray combustion cases, with the added complication of an active
interface. The system obtained after semi-discretisation in space is shown to
be differential-algebraic of index 1. A numerical strategy relying on stiffly
accurate Runge-Kutta methods is introduced, with a specific discretisation of
the algebraic constraints and time adaptation. High order is shown to be
reached on all variables, while handling the constraints properly. Three
challenging test cases are investigated: ignition, limit cycle, and unsteady
response with detailed gas-phase kinetics. We show that the time integration
method can greatly affect the ability to predict the dynamics of the system.
The proposed numerical strategy exhibits high efficiency and accuracy for all
cases compared to traditional schemes used in the combustion literature.
