A Conditional Moment Closure Study of Chemical Reaction Source Terms in SCCI Combustion

The objective of this study is to evaluate conditional moment closure (CMC) approaches to model chemical reaction rates in compositionally stratified, autoigniting mixtures, in thermochemical conditions relevant to stratified charge compression ignition (SCCI) engines. First-order closure, second-order closure and double conditioning are evaluated and contrasted as options in comparison to a series of direct numerical simulations (DNSs). The two-dimensional (2D) DNS cases simulate ignitions in SCCI-like thermochemical conditions with compositionally stratified n-heptane/air mixtures in a constant volume. The cases feature two different levels of stratification with three mean temperatures in the negative-temperature coefficient (NTC) regime of ignition delay times. The first-order closure approach for reaction rates is first assessed using hybrid DNS-CMC a posteriori tests when implemented in an open source computational fluid dynamics (CFD) package known as OpenFOAM. The hybrid DNS-CMC a posteriori tests are not a full CMC but a DNS-CMC hybrid in that they compute the scalar and velocity fields at the DNS resolution, thus isolating the first-order reaction rate closure model as the main source of modelling error (as opposed to turbulence model, scalar probability density function model, and scalar dissipation rate model). The hybrid DNS-CMC a posteriori test reveals an excellent agreement between the model and DNS for the cases with low levels of stratification, whereas deviations from the DNS are observed in cases which exhibit high level of stratifications. The a priori analysis reveals that the reason for disagreement is failure of the first-order closure hypothesis in the model due to the high level of conditional fluctuations. Second-order and double conditioning approaches are then evaluated in a priori tests to determine the most promising path forwards in addressing higher levels of stratification. The a priori tests use the DNS data to compute the model terms, thus directly evaluating the model assumptions. It is shown that in the cases with a high level of stratification, even the second-order estimation of the reaction rate source term cannot provide a reasonably accurate closure. Double conditioning using mixture-fraction and sensible enthalpy, however, provides an accurate first-order closure to the reaction rate source term.