LAUSR

Abstract A numerical investigation was conducted to study ignition and the initial phases of methane combustion in a two-dimensional setup consisting of a pre- (PC) and main (MC) chamber. The effect of the wall thermal boundary condition (isothermal T w = 500  K or adiabatic), initial mixture temperature ( T u = 300 and 800 K) and equivalence ratio in the main chamber ( ϕ M C = 0.5 and 1.0) was studied. A stoichiometric mixture was used in the PC and the mixtures were initially quiescent in all cases. Flame propagation in the prechamber is affected mainly by the initial temperature, while the thermal state of the wall plays a minor role. The transient jet that is generated at the exit of the orifice connecting the two chambers creates an intense flow field with vortical structures that, depending on the initial temperature, persist for a long time or dissipate quickly affecting combustion in the main chamber. Depending on the local flow and mixing conditions close to the orifice exit, the hot jet can be broken into small kernels at low Tu or forms quickly a flame torch at hot conditions, strongly affecting ignition of the MC mixture and flame propagation in the main chamber. In the lean MC cases, the intense mixing with the stoichiometric PC mixture creates local compositions that are more favorable for ignition by the hot turbulent reactive jet that subsequently exits from the PC at a temperature that is significantly lower than the adiabatic flame temperature of the corresponding mixture. Despite the short residence time, the reactive state of the mixture is affected as it flows through the nozzle. The flame structures in the MC are described in terms of the progress variable and mixture fraction and compared to flamelet-type calculations. The local flame structure differs strongly from that of the 1-D unstrained premixed flame, particularly for the low Tu cases. The flamelet-type calculations show that ignition of the most reactive mixture is enhanced by the radicals in the hot reactive jet, while scalar dissipation rate accelerates the ignition of the whole mixture. The 2-D simulations show that ignition is significantly longer than what is predicted by the flamelet calculations.