Convergence of a finite difference method for combustion model problems (Q2386574)

The author studies the convergence of a numerical scheme intended to solve 1D combustion models. This scheme combines a projection method near the combustion wave and an upwind finite difference scheme far from this combustion wave. The equations of the 1D combustion model are: \(\partial \left( u+qz\right) /\partial t+\partial \left( f\left( u\right) \right) /\partial x=0\) and \(\partial z/\partial t=-K\phi \left( u\right) z\), where \(u \) denotes the density, or the velocity, or the temperature of the gas, \(z\) represents the fraction of the unburnt gas (\(z\in \left[ 0,1\right] \)), \(q\) and \(K\) are positive constants. \(u\) and \(z\) start at \(t=0\) from initial data \(u_{0}\) and \(z_{0}\) respectively. The author refers to the literature for the existence of a solution of this problem, assuming that \(f\) and \(\phi \) satisfy appropriate hypotheses. Starting from the classical explicit upwind scheme, the author considers two other schemes introducing intermediate points. The author proves that these schemes work under the classical CFL condition and hypotheses on the initial data \(u_{0}\) and \(z_{0}\). The main part of the paper is devoted to the proof of the convergence of this numerical scheme. The author also proves the existence of a curve \(\Gamma \) which separates the zones \(z=0\) and \(z=1\) and which is a combustion wave satisfying a Rankine-Hugoniot condition.