Stochastic simulation of variations in the autoignition delay time of premixed methane and air

A mesoscopic stochastic particle model for homogeneous combustion is introduced. The model can be used to investigate the physical fluctuations in a system of coupled chemical reactions with energy (heat) release/consumption. In the mesoscopic model, the size of the homogeneous gas volume is an additional variable, which is eliminated in macroscopic continuum models by the thermodynamic limit N→∞. Thus, continuous homogeneous models are macroscopic models wherein fluctuations are excluded by definition. Fluctuations are known to be of particular importance for systems close to the autoignition limits. The new model is used to investigate the stochastic properties of the autoignition delay time in a homogeneous system with stoichiometric premixed methane and air. Temperature and species concentrations during autoignition of sub-macroscopic volumes, including physically meaningful fluctuations, are presented. It is found that different realizations mainly differ in the time when ignition occurs; besides this the development is similar. The mesoscopic range and the macroscopic limit are identified. Which range a specific system is assigned to is not only a question of the length scale or particle number, but also depends on the complete thermodynamic state. The stochastic algorithm yields the correct results for the macroscopic limit compared to the continuous balance equations. The sensitivity of the results to two different detailed reaction mechanisms (for the same system) is studied and found to be low. We show that when approaching the autoignition limit by decreasing the temperature, the fluctuations in the autoignition delay time increase and an increasing number of realizations will have exceedingly long ignition delay times, meaning they are in practice not autoignitable. With this result the mesoscopic simulations offer an explanation of the transition between autoignitable and non-autoignitable conditions. The calculated distributions were compared with ten repetitions of the same experiment. A mesoscopic distribution that matches the experimental results was found.