Non-equilibrium effects on thermal ignition using hard sphere molecular dynamics simulations


Abstract in English

The present study addresses the role of molecular non-equilibrium effects in thermal ignition problems. We consider a single binary reaction of the form A+B -> C+C. Molecular dynamics calculations were performed for activation energies ranging between RT and 7.5RT and heat release of 2.5RT and 10RT. The evolution of up to 10,000 particles was calculated as the system undergoes a thermal ignition at constant volume. Ensemble averages of 100 calculations for each parameter set permitted to determine the ignition delay, along with a measure of the stochasticity of the process. A well behaved convergence to large system sizes is also demonstrated. The ignition delay calculations were compared with those obtained at the continuum level using rates derived from kinetic theory: the standard rate assuming that the distribution of the speed of the particles is the Maxwell-Boltzmann distribution, and the perturbed rates by Prigogine and Xhrouet [1] for an isothermal system, and Prigogine and Mahieu [2] for an energy releasing reaction, obtained by the Chapman-Enskog perturbation procedure. The molecular results were found in very good agreement with the latter at low temperatures, confirming that non-equilibrium effects promote the formation of energetic particles, that serve as seeds for subsequent reaction events: i.e., hot spots. This effect was found to lower the ignition delay by up to 30%. At high temperatures, the ignition delay obtained from the standard equilibrium rate was found to be up to 60% longer than the molecular calculations. This effect is due to the rapidity of the reactive collisions that do not allow the system to equilibrate. For this regime, none of the perturbation solutions obtained by the Chapman-Enskog procedure were valid. This study thus shows the importance of non-equilibrium effects in thermal ignition problems, for most temperatures of practical interest.

Download