Large Eddy simulation of the turbulent spark ignition and of the flame propagation in spark ignition engines
Large Eddy simulation of the turbulent spark ignition and of the flame propagation in spark ignition engines
复制标题
火花点火发动机中湍流火花点火和火焰传播的大涡模拟
DOI:
--
复制
发表时间:
2016
期刊:
影响因子:
--
通讯作者:
S. Mouriaux
中科院分区:
文献类型:
--
作者:
S. Mouriaux
The use of lean equivalence ratios or high EGR rates in spark ignition engines (SIE) enables to optimize CO2 and NOx emissions; however too important dilution rates leads to increased cycle-to-cycle variability. These latter are mostly due to the ignition phase, which becomes critical when dilution rates are important and requires high ignition energy. The ECFM-LES model currently used in IFPEN, which is based on the flame surface density concept, is not sufficient to describe ignition in these critical conditions. The TF-LES approach was chosen in this study, principally because it directly resolved chemistry and can thus model ignition via a local raise of the temperature. The present work defines and evaluates a simulation strategy for TF-LES in SIE configurations, that enables a fine prediction of critical ignitions and of the turbulent flame propagation.In the first part, DNS of turbulent ignition were performed. The ignition phase was modeled using a thermal energy deposit (ED model, Lacaze et al.). Simulations reproduced the ignition experiments of Cardin et al. who determined the minimum ignition energy (MIE) of lean premixed methane/air mixtures, for different turbulence characteristics. The main purpose of the study was to determine the numerical and physical model parameters, which enable to reproduce Cardin et al. experiments. Two types of kinetic schemes were evaluated: a simplified kinetic scheme and an analytical kinetic scheme (ARC), that can predict both the auto-ignition delays and the laminar flame speed, while keeping affordable CPU times. Results analysis enabled to define ignition criteria and to highlight the differences in terms of ignition prediction using the two kinetic schemes. Results also demonstrated that the chosen approach could recover correct levels of ignition energy for laminar and low Karlovitz number cases (Ka<10). For higher Karlovitz number cases, the ED model was found to be insufficient to predict the ignition and a finer description of the energy deposit is required.In the second part, a dynamic wrinkling model (Wang et al., 2012) was studied to describe the out-of-equilibrium behavior of the flame during the propagation phase. Studies on laminar spherical flames were first performed, to assess the laminar degeneration of the model. Then, as first tests in an engine configuration have revealed incompatibilities of the model, modifications were proposed. The modified dynamic model was finally tested in the ICAMDAC engine configuration. Results of the simulations were compared against previous results of Robert et al. obtained with the ECFM-LES model using a transport equation for the flame surface density that can describe the out-of-equilibrium wrinkling of the flame. Results obtained with the dynamic model are in very good agreement with the ones of Robert et al., thus demonstrating the ability of the dynamic model to predict out-of-equilibrium values in the engine configuration. Besides, the dynamic model self-adapts to the turbulence conditions, hence does not require any model parameter adjustment, as is it the case for models based on the flame surface density transport equation.