Multi-scale modeling of dynamics and ignition to flame transitions of high pressure stratified n -heptane/toluene mixtures

Abstract The transitions from ignition to flames as well as the combustion dynamics in stratified n-heptane and toluene mixtures are numerically modeled by a correlated dynamic adaptive chemistry method coupled with a hybrid multi-timescale method (HMTS/CO-DAC) in a one-dimensional constant volume chamber. The study attempts to answer how the kinetic difference between n-alkanes and aromatics leads to different ignition to flame transitions and knocking-like acoustic wave formation at low temperature and engine pressure conditions with fuel stratification. It is found that the low temperature chemistry (LTC) and fuel stratification of n-heptane leads to the formation of multiple ignition fronts. Four different combustion wave fronts, a low temperature ignition (LTI) front followed by a high temperature ignition (HTI) front, a premixed flame (PF) front, and a diffusion flame (DF) front, are demonstrated. The fast LTI and HTI wave front propagation leads to a shock-like strong acoustic wave propagation, thus strongly modifying the dynamics of the subsequent diffusion and premixed flame fronts. On the other hand, for the toluene mixture, due to the lack of LTC, only two combustion wave fronts are formed, a HTI front and a premixed flame front, exhibiting stable flow field and no formation of shock-like acoustic wave. The dynamics of transition from combustion to shock waves is further analyzed by using a modified Burgers’ equation. The analysis for n-heptane/air mixture indicates that both the onset of LTI and the strong dependency of HTI on the equivalence ratio can either promote or attenuate the transition from strong acoustic wave to shock wave. However, the toluene/air mixture exhibits no coupling with acoustic wave, suggesting that the rich LTC reactivity with fuel stratification, specific to the n-alkane chemistry, can lead to knocking and acoustic formation.

[1]  J. Burgers The Nonlinear Diffusion Equation , 2013 .

[2]  H. Im,et al.  Direct numerical simulation of ignition front propagation in a constant volume with temperature inhomogeneities: I. Fundamental analysis and diagnostics , 2006 .

[3]  Philippe Pierre Pebay,et al.  Direct numerical simulation of ignition front propagation in a constant volume with temperature inhomogeneities. II. Parametric study , 2006 .

[4]  Y. Zeldovich,et al.  Regime classification of an exothermic reaction with nonuniform initial conditions , 1980 .

[5]  John B. Heywood,et al.  Internal combustion engine fundamentals , 1988 .

[6]  Amanda Ramcharan,et al.  Relative Reactivities of the Isomeric Butanols and Ethanol in an Ignition Quality Tester , 2011 .

[7]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[8]  Y. Ju,et al.  Multi-timescale modeling of ignition and flame propagation of diesel surrogate fuel mixtures , 2013 .

[9]  Chi-Wang Shu,et al.  High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems , 2009, SIAM Rev..

[10]  Y. Ju,et al.  Direct numerical simulations of NOx effect on multistage autoignition of DME/air mixture in the negative temperature coefficient regime for stratified HCCI engine conditions , 2014 .

[11]  Chi-Wang Shu,et al.  High Order ENO and WENO Schemes for Computational Fluid Dynamics , 1999 .

[12]  Xiaolong Gou,et al.  A path flux analysis method for the reduction of detailed chemical kinetic mechanisms , 2010 .

[13]  M. P. Burke,et al.  Effect of cylindrical confinement on the determination of laminar flame speeds using outwardly propagating flames , 2009 .

[14]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

[15]  Thomas A. Litzinger,et al.  The experimental evaluation of a methodology for surrogate fuel formulation to emulate gas phase combustion kinetic phenomena , 2012 .

[16]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[17]  Johannes Martinus Burgers,et al.  The Nonlinear Diffusion Equation: Asymptotic Solutions and Statistical Problems , 1974 .

[18]  John E. Dec,et al.  Advanced compression-ignition engines—understanding the in-cylinder processes , 2009 .

[19]  John Dinwoodie,et al.  SAE Technical Paper Series , 1987 .

[20]  R. Reitz Directions in internal combustion engine research , 2013 .

[21]  M. P. Burke,et al.  Multi-timescale modeling of ignition and flame regimes of n-heptane-air mixtures near spark assisted homogeneous charge compression ignition conditions , 2011 .

[22]  Jacqueline H. Chen,et al.  A numerical study of the autoignition of dimethyl ether with temperature inhomogeneities , 2013 .

[23]  Xiaolong Gou,et al.  A dynamic multi-timescale method for combustion modeling with detailed and reduced chemical kinetic mechanisms , 2010 .

[24]  C. Westbrook,et al.  A Comprehensive Modeling Study of n-Heptane Oxidation , 1998 .

[25]  Y. Ju,et al.  Kinetic effects of toluene blending on the extinction limit of n-decane diffusion flames , 2010 .

[26]  Zheng Chen,et al.  On the extraction of laminar flame speed and Markstein length from outwardly propagating spherical flames , 2011 .

[27]  Y. Ju,et al.  Direct numerical simulations of exhaust gas recirculation effect on multistage autoignition in the negative temperature combustion regime for stratified HCCI flow conditions by using H2O2 addition , 2013 .

[28]  Thomas A. Litzinger,et al.  A jet fuel surrogate formulated by real fuel properties , 2010 .

[29]  Yiguang Ju,et al.  Effects of Lewis number and ignition energy on the determination of laminar flame speed using propagating spherical flames , 2009 .

[30]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[31]  D. Splitter,et al.  Fuel reactivity controlled compression ignition (RCCI): a pathway to controlled high-efficiency clean combustion , 2011 .