Numerical modeling of lean hydrogen spark-ignition engines: On the role of intrinsic instabilities
Hydrogen-fueled internal combustion engines (H 2 -ICEs) hold strong potential as a pathway toward CO 2 -neutral propulsion. To reduce emissions, H 2 -ICEs are usually operated under fuel-lean conditions, where the flames are prone to thermo-diffusive instabilities (TDIs). These TDIs govern both local and global flame propagation, but their impact on full-scale engine combustion remains an open question. In this study, high-fidelity three-dimensional large-eddy simulations (LES) are performed at multiple mesh resolutions, with the finest grid sufficiently resolved to directly characterize flame front dynamics relevant to engine-scale combustion. The simulations reveal cellular and finger-like flame structures characteristic of TDIs throughout the entire combustion process. Analysis of the local thermo-chemical state demonstrates that differential diffusion induces pronounced mixture stratification and elevates reaction rates, resulting in super-adiabatic temperatures that strongly correlate with flame curvature. Building on these findings, the performance of the baseline artificially thickened flame (ATF) model and a recently developed thermo-diffusive (TD)-aware extension is assessed. Unlike the state-of-the-art ATF model, which suffers from grid dependence and underestimates the experimental pressure trace, the TD-aware formulation captures experimental trends more accurately and provides consistent, grid-independent integrated heat-release (IHR) traces. For the operating condition considered here, the results show that TD effects represent sub-grid-scale contributions that need to be accounted for to obtain consistent predictions of global combustion behavior under the investigated lean H 2 -ICE conditions. • LES captured cellular and finger-like thermo-diffusive flame structures. • Coarser grids suppressed fine-scale instabilities resolved at high resolution. • Local mixture stratification enhanced reactivity and caused super-adiabatic states. • ATF model showed grid bias from missing thermo-diffusive instability treatment. • Thermo-diffusive-aware ATF model reduced grid bias and improved predictive accuracy.
- Research Article
21
- 10.1080/13647830.2015.1132010
- Feb 23, 2016
- Combustion Theory and Modelling
The detailed flame structure of laminar premixed cellular flames in the tubular domain is simulated in 2D using a fully-implicit primitive variable finite difference formulation that includes multicomponent transport and detailed chemical kinetics. Numerical results for H2/air flames are presented and compared against spatially resolved experimental measurements of temperature and chemical species including atomic H and OH. The experimental results compare well for flame structure and cell number, despite the numerical model under-predicting the peak temperature by 200 K. Numerical experiments were performed to assess the ability for cellular tubular flames to impact experimental and numerical investigations of practical flames. The cellular flame structure is found to provide a highly sensitive geometry that is useful for validating diffusive transport modelling approximations. This capability is exemplified through the development of a simple and accurate approximation for thermal diffusion (i.e. the Soret effect) that is suitable for practical combustion codes.
- Research Article
56
- 10.1016/0360-1285(89)90009-9
- Jan 1, 1989
- Progress in Energy and Combustion Science
Selective diffusional demixing: Occurrence and size of cellular flames
- Research Article
9
- 10.1016/j.combustflame.2009.12.011
- Jan 12, 2010
- Combustion and Flame
Three-dimensional simulations of cellular non-premixed jet flames
- Research Article
165
- 10.1016/j.combustflame.2022.112254
- Jun 24, 2022
- Combustion and Flame
Synergistic interactions of thermodiffusive instabilities and turbulence in lean hydrogen flames
- Research Article
3
- 10.1016/j.fuel.2023.130620
- Dec 17, 2023
- Fuel
The effect of mixture inhomogeneity and turbulence on the flame front curvature and flame surface density of turbulent planar flames of natural gas
- Research Article
57
- 10.1016/j.ijhydene.2016.07.086
- Jul 30, 2016
- International Journal of Hydrogen Energy
High hydrogen content syngas fuel burning in lean premixed spherical flames at elevated pressures: Effects of preferential diffusion
- Abstract
- 10.1016/s0959-6526(97)88902-0
- Jan 1, 1996
- Journal of Cleaner Production
5538594 Method for producing a blade coated paper from recycled, high lignin content, waste paper: Hank Mark A; Mulcahy Leo T; Peterson Ralph S; Streisel Robert C, Covington, VA, United States assigned to Westvaco Corporation
- Research Article
55
- 10.1063/1.870201
- Nov 1, 1999
- Physics of Fluids
The influence of hydrodynamic instability on the structure of two-dimensional (2D) and three-dimensional (3D) cellular flames is numerically investigated. The equation used is the compressible Navier–Stokes equation including a one-step irreversible chemical reaction. We superimpose an infinitesimal disturbance on the stationary plane flame and calculate the evolution of the disturbed flame front to obtain the relation between the growth rate and the wave number, i.e., the dispersion relation. With an increase in flame temperature, the growth rate increases since hydrodynamic instability becomes stronger. The unstable range normalized by the preheat zone thickness hardly changes, even though the flame temperature increases. The critical wave number, which corresponds to the maximum growth rate, is almost constant. Therefore, the normalized spacing between cells of the cellular flame is independent of the flame temperature. Moreover, we superimpose the disturbance with the critical wave number to investigate the structure of cellular flames. The stationary cellular flame is obtained when the inlet-flow velocity is set to the flame velocity of the cellular flame. The higher the flame temperature, the deeper the cell and the broader the flame surface. In addition, the cell depth and the surface area of the 3D flame are larger than those of the 2D flame. This is caused by the difference in the disposition of cells.
- Research Article
1
- 10.1299/kikaib.64.3491
- Jan 1, 1998
- TRANSACTIONS OF THE JAPAN SOCIETY OF MECHANICAL ENGINEERS Series B
The influence of hydrodynamic instability on the structure of two-dimensional (2D) and three-dimensional (3D) cellular flames is numerically investigated. The equation used is the compressible Navier–Stokes equation including a one-step irreversible chemical reaction. We superimpose an infinitesimal disturbance on the stationary plane flame and calculate the evolution of the disturbed flame front to obtain the relation between the growth rate and the wave number, i.e., the dispersion relation. With an increase in flame temperature, the growth rate increases since hydrodynamic instability becomes stronger. The unstable range normalized by the preheat zone thickness hardly changes, even though the flame temperature increases. The critical wave number, which corresponds to the maximum growth rate, is almost constant. Therefore, the normalized spacing between cells of the cellular flame is independent of the flame temperature. Moreover, we superimpose the disturbance with the critical wave number to investigate the structure of cellular flames. The stationary cellular flame is obtained when the inlet-flow velocity is set to the flame velocity of the cellular flame. The higher the flame temperature, the deeper the cell and the broader the flame surface. In addition, the cell depth and the surface area of the 3D flame are larger than those of the 2D flame. This is caused by the difference in the disposition of cells.
- Research Article
8
- 10.1016/j.combustflame.2021.111442
- Apr 15, 2021
- Combustion and Flame
Cellular flame structures and thermal characteristics of axi-symmetric ceiling fires: An experimental study and scaling analysis
- Research Article
- 10.1299/kikaib.65.3475
- Jan 1, 1999
- TRANSACTIONS OF THE JAPAN SOCIETY OF MECHANICAL ENGINEERS Series B
The two-dimensional (2-D) and three-dimensional (3-D) unsteady reactive flows are calculated to study the structure of downwards/upwards propagating cellular flames at the Lewis number unity. We superimpose an infinitesimal disturbance on the stationary plane flame to obtain the growth rate depending on the wave number, i.e., the dispersion relation. When premixed flames are propagated downwards/upwards, the growth rate decreases/increases and the unstable range becomes narrower/broader with an increase in acceleration. Since we have positive growth rates, the cellular shape of flame fronts appears owing to intrinsic instability. To study the formation of a cellular flame front, we superimpose the disturbance with the peculiar wave number which corresponds to the maximum growth rate. The disturbance on a flame is evolved, and eventually the cellular flame front is formed. When the flame is propagated downwards, the spacing between cells is almost constant and the depth of cells becomes smaller as the acceleration increases. When the flame is propagated upwards, on the other hand, the former becomes smaller and the latter becomes larger.
- Research Article
9
- 10.1016/0010-2180(94)00223-f
- Jun 1, 1995
- Combustion and Flame
Hydrodynamic aspects of premixed flame stripes in two-dimensional stagnation-point flows
- Research Article
36
- 10.1016/j.combustflame.2017.05.031
- Jun 22, 2017
- Combustion and Flame
A consistent Artificially Thickened Flame approach for spray combustion using LES and the FGM chemistry reduction method: Validation in Lean Partially Pre-vaporized flames
- Research Article
61
- 10.1002/aic.11180
- Apr 23, 2007
- AIChE Journal
Non‐premixed methane‐oxygen flame dynamics and structures confined within an alumina combustor are described. Non‐stabilized, transient flame dynamics and phenomena leading to the formation of a stable edge‐like flame and distinct cellular structures that take place within combustor channels of dimensions 35 mm long, 5 mm wide, and 0.75 mm high (combustor volume ∼ 130 mm3) are discussed. These confined flames are surveyed by measuring the external wall temperatures, high‐speed and still‐frame visual flame imaging, recordings of emitted acoustics from the combustor, and capturing visible, CH*, and OH* chemiluminescence through a sapphire window. The observed dynamic flame structure is an oscillating edge‐like flame accompanied by ignition‐extinction events that precede the formation of a stable edge‐like flame and flame cells in the reaction channel. The cellular flame structures in all cases exhibit a confined tribrachial structure with a folded or extinguished rich branch while the lean branch survives. © 2007 American Institute of Chemical Engineers AIChE J, 2007
- Research Article
6
- 10.1016/s0167-2789(96)00136-4
- Dec 1, 1996
- Physica D: Nonlinear Phenomena
Structure and dynamics of kink and cellular flames stabilized on a rotating burner