Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ignition delay”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

Limitations of cetane number to predict transient combustion phenomena in high-pressure fuel sprays

Fundamental understanding of in-cylinder processes in diesel engines is important to screen emerging biofuels and advanced combustion modes that can reduce greenhouse gas emissions and regulated pollutants including soot. In this study, the role of fuel properties on spray development and combustion is investigated by systematically isolating chemical and thermophysical effects. Three different fuels are considered, two with similar chemical properties and two with similar thermophysical properties with one fuel common to both groups. Experiments are performed in a constant-pressure flow chamber de-signed to provide stable test conditions and facilitate acquisition of at least 150 injections in quick succession for each fuel under reacting conditions at high-pressure, high-temperature ambient conditions using a modified conventional diesel engine injector. Further, high speed optical diagnostics including rainbow schlieren deflectometry, OH* chemiluminescence, and a two-color pyrometry system are employed to simultaneously image the transient spray and reacting jets. Image analysis is performed to determine liquid length, vapor penetration length, timing and location of first stage and main ignition events, lift-off location, total soot mass, and more. Results show that fuels with similar chemical properties or cetane number (CN) exhibit similar delay times for first stage and main ignition events as may be expected, but very different liquid length, first stage and main ignition locations, lift-off length, apparent turbulent flame speed, and soot formation. As such, the ability to characterize candidate biofuels with CN or other parameters derived from simple flame configurations is called into question. In this study, thermo-physical properties controlling the liquid length are identified as the main contributing factor for the observed differences.

33 ADVANCED PROPULSION SYSTEMS↗

An experimental investigation of flame and autoignition behavior of propane

We report autoignition delay time data are one important means to develop, quantify, and validate fundamental understanding of combustion chemistry at low temperatures (T<1200 K). However, low-temperature chemistry often has higher uncertainties and scatter in the experimental data compared with high-temperature ignition data (T>1200 K). In this study, autoignition properties of propane and oxygen mixtures were investigated using the University of Michigan rapid compression facility in order to understand the effects of ignition regimes on low-temperature ignition data. For the first time for propane, autoignition delay times were determined from pressure histories, and autoignition characteristics were simultaneously recorded using highspeed imaging of the test section through a transparent end-wall. Propane mixtures with fuel-to-O 2 equivalence ratios of φ = 0.25 and φ = 0.5 and O 2 -to-inert gas molar ratios of 1:3.76 were studied over the pressure range of 8.9 to 11.3 atm and the temperature range of 930 – 1070 K. The results showed homogeneous or strong autoignition occurred for all φ = 0.25 experiments, and inhomogeneous or mixed autoignition occurred for all φ = 0.5 experiments. While a limited temperature range is covered in the study, importantly the data span predicted transitions in autoignition behavior, allowing validation of autoignition regime hypotheses. Specifically, the results agree well with strong-autoignition limits proposed based on the Sankaran Criterion. The autoignition delay time data at the strong-ignition conditions are in excellent agreement with predictions using a well-validated detailed reaction mechanism from the literature and a zerodimensional modeling assumption. However, the experimental data at the mixed autoignition conditions were systematically faster than the model predictions, particularly at lower temperatures (T< ~970 K). The results are an important addition to the growing body of data in the literature that show mixed autoignition phenomena are important sources of the higher scatter observed in the low-temperature autoignition data for propane and other fuels. The results are discussed in terms of different methods to capture the effects of pre-autoignition heat release associated with mixed autoignition conditions and thereby address some of the discrepancies between kinetic modeling and experimental measurements.

42 ENGINEERING↗

Modeling delayed thermal runaway in nitric acid-soaked cat litter mixed with radioactive waste

Thermal ignition of radioactive waste within a 55-gallon drum was simulated by using a pressure-dependent waste decomposition model (Hobbs et al. in Process Saf Environ Prot https://doi.org/10.1016/j.psep.2022.09.047, 2022) calibrated with data from full-scale drum experiments (Parker et al. in The thermolytic response of a surrogate RNS waste mixture at the drum scale. Los Alamos National Laboratory Report LA-UR-16-21760, 2016) and validated with experiments from multiple laboratories (Hobbs et al. in Thermal analysis of aged nitric acid-soaked kitty litter in TRU waste drums-23370.WM2023 Conference, Phoenix, AZ, 2023). The acceleration of nitric acid chemistry reacting with an organic cat litter leading to thermal ignition was likely triggered by a restricted vent in the drum. Here, we address whether the form of the rate equation in (Hobbs et al. in Process Saf Environ Prot https://doi.org/10.1016/j.psep.2022.09.047, 2022) is sufficient to extrapolate thermal ignition within aged drums of similar content that have been stored in Texas for over nine years by investigating four different reaction rate forms for waste decomposition. A critical reaction rate reduction analysis is performed on each of these models to determine if delayed thermal runaway within vented aged waste is possible after nine years. We found that a pressure-dependent first-order rate expression not only predicted the accidental ignition of the waste drum, but the form also matches multiple experiments from different laboratories. Even though the waste composition decreases over time, the model predicts that acceleration leading to thermal runaway is possible if the waste is confined, even after 9 years. In conclusion, waste containing oxidizers such as nitric acid should not be mixed with organic adsorbents, especially if the waste is confined.

cookoff↗

Effects of non-thermal termolecular chemistry on detonation development in hydrogen (H2) / methane (CH4) - air mixtures

In combustion simulations, it is usually assumed that highly-energized collision complexes (either radical or stable species) formed in exothermic reactions are always under thermal equilibrium as they evolve through competing networks of reactions. However, in practical flames, non-trivial amounts of reactive radicals such as H, O and OH are present apart from O2. As a result, collisions of these reactive species with the energized collision complexes have been shown in recent studies to induce non-thermal reactivity. These studies have also demonstrated that such non-thermal reactions can be suitably represented in macroscopic kinetics models as chemically termolecular reactions. The present work was focused on identifying and quantifying the effects of including such reactions on the evolution of an initial deflagration front to a developing detonation in H2/CH4-air mixtures under boosted internal combustion (IC) engine conditions. Specifically, fully resolved simulations, with and without non-thermal reactivity, were performed for a constant volume reactor containing stoichiometric H2/CH4-air mixture. It was found that inclusion of termolecular chemistry resulted in a delayed onset of end-gas auto-ignition. Concurrently, the developing detonation intensity was observed to be significantly higher. Chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions responsible for the observed variation in the results.

Detonation↗

Advanced Intra-Cycle Detection of Pre-Ignition Events through Phase-Space Transforms of Cylinder Pressure Data

The widespread adoption of boosted, downsized SI engines has brought pre-ignition phenomena into greater focus, as the knock events resulting from pre-ignitions can cause significant hardware damage. Much attention has been given to understanding the causes of pre-ignition and identify lubricant or fuel properties and engine design and calibration considerations that impact its frequency. This helps to shift the pre-ignition limit to higher specific loads and allow further downsizing but does not fundamentally eliminate the problem. Real-time detection and mitigation of pre-ignition would thus be desirable to allow safe engine operation in pre-ignition-prone conditions. This study focuses on advancing the time of detection of pre-ignition in an engine cycle where it occurs. Furthermore, phase space transforms through time-delay embedding of cylinder pressure and principal component analysis were applied to same-cycle detection of pre-ignition and shown to enable detection on the order of a crank degree earlier than deviation in cylinder pressure can be identified through direct statistical observation of the pressure data. Additionally, it appears that the deviation of the trajectory in phase space may offer the opportunity to extend this method to further extend the detection window and allow more time for mitigation actions to occur.

42 ENGINEERING↗

Effects of non-thermal termolecular reactions on detonation development in hydrogen (H 2 )/methane (CH 4 ) - air mixtures

The binary fuel blend of H 2 /CH 4 is one of the most promising hydrogen-enriched hydrocarbon fuels in spark-ignition (SI) engines. Yet, the undesirable phenomenon of super-knock, which can severely and in-stantaneously damage an SI engine, limits its widespread adoption. Moreover, there is still a lack of con-sensus on the precise mechanism by which this phenomenon occurs i.e. via flame acceleration or spon-taneous ignition, despite numerous previous investigations. At the same time, recent studies [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078 -1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522] have demonstrated a high probability of occurrence of non-thermal reactions in premixed flames of such H 2 /CH 4 fuel blends with air due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 . The present study focuses on the evolution of an initial deflagration front to a detonation wave in H 2 /CH 4 - air mixtures under SI engine relevant conditions through fully resolved, constant volume 1D simulations with and without non-thermal reactivity. Non-thermal reactions were included in the macroscopic kinetics model as chemically termolecular reactions facilitated by the H + CH 3 and H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Further, the nonthermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions. Therefore, an additional set of simulations were performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522]. Compared to the baseline case, the onset of spontaneous ignition in the end-gas region was observed to be delayed in the presence of non-thermal termolecular reactions. Concurrently, the developing detonation was observed to be significantly stronger. In contrast, applying corrections to the recombi-nation/association rate constants resulted in a completely different behavior. Specifically, detonation was observed to occur due to self acceleration of the primary flame in the absence of spontaneous ignition in the end-gas region. Sensitivity analysis was performed to quantify the effects of non-thermal reactions on the duration of heat release rate and thereby the mechanism of detonation formation. In addition, chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions re-sponsible for the observed results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lightning Interferometry Uncertainty, Beam Steering Interferometry, and Evidence of Lightning Being Ignited by a Cosmic Ray Shower

Abstract We present an uncertainty analysis on correlation‐based time delay estimate, the basis for broadband lightning interferometry. A normal interferometry would yield much higher uncertainties than the theoretically predicted lower bound due to finite correlation window lengths. However, if the signals are aligned in time, the uncertainties approach the lower bound and can be used to indicate the interferometric uncertainties. Based on this, we introduce a beam steering interferometry technique. It first estimates a direction centroid for the lightning sources with a normal interferometry and computes the time delays among the signals. It then shifts the raw data with the time delays to align the signals roughly in time and reprocess the aligned signals. It finally pushes the reprocessed results to their correct positions, based on the time delays estimated in the first step. We apply this technique on a fast positive breakdown (FPB) process that started a normal intracloud lightning. The FPB process is shown to have a much more complex structure and development than a normal interferometry would provide. More importantly, from both interferometric and polarization analyses, the FPB appears to be ignited by a cosmic ray shower (CRS). We estimate the radio frequency strength and frequency content related to a presumed CRS in a thunderstorm electric field and find that they are in agreement with the observations. We examine the electric field effect of the CRS front and find that it could raise the field above the threshold for positive breakdown and is capable of igniting the FPB (and the lightning).

Shao, Xuan‐Min↗

Controlling pollutant emissions in a high-pressure combustor with fuel-diluent blending

This paper investigates the formation of nitric oxide (NO) and carbon monoxide (CO) emissions as a function of nitrogen (N 2 ) or carbon dioxide (CO 2 ) diluent content, added to a premixed reacting jet in crossflow. Here, reaction characteristics of a rich methane-air jet injected into a lean vitiated crossflow were analyzed at an elevated pressure of 5 atm. The jet was pre-heated and enriched with 0%, 15%, 30%, and 50% mass fraction to quantify the effect on pollutant emissions. Simulated results of the full chemistry Star-CCM+CFD model were verified with data taken in the experimental high-pressure combustion facility, which provides pressure, temperature and velocity profiles, as well as line-of-sight chemiluminescence and exit emission measurements. The analysis revealed the significant influence of the diluent to delay axial combustion; this mixing delay increased the axial flame lift-off and controlled any thermal ignition limitations. Increasing the diluent content increased the timescale for the flame to stabilize, which allowed for greater entrainment of crossflow oxidizer into the axial jet stream, and led to decreased pollutant emissions. Hence, crossflow entrainment is a critical driving force at high diluent content. Local nitric oxide (NO) emission formation in the axial stage was predicted numerically, showing the correlation between diluent addition, axial heat release, and the formation of nitric oxide pollution. The high diluent levels leaned out the jet mixture and delayed axial combustion, while reducing prompt and thermal NOx by mitigating flame hotspots and minimizing the timescale that the products remain at high temperatures.

03 NATURAL GAS↗

Autoignition Delay Time Measurements of Dimethyl Ether at 110 Bar Inside a Shock Tube

Fuel mixtures of propane and dimethyl ether have been proposed as potential diesel substitutes for use in compression-ignition engines. There is an interest in phasing out diesel fuel because its combustion byproducts are harmful to the environment. Propane and dimethyl ether burn cleaner and produce significantly less carbon dioxide than diesel. The combustion of these fuels at engine-relevant conditions must be characterized to create an engine capable of running on propane and dimethyl ether. In this work, IDT measurements of stoichiometric neat DME diluted in 75% argon were completed for temperatures spanning 800-1200K at 110 bar. A shock tube will be used for all experiments in this work. Experimental IDTs were then compared to simulated values using several well-established mechanisms in the literature.

Guenther, Aaron C. [University of Central Florida]↗

Thermite and intermetallic projectiles examined experimentally in air and inert gas environments

Intermetallic (aluminum and zirconium) and thermite (aluminum and molybdenum trioxide) projectiles were launched using a high velocity impact ignition testing system. The experiments were designed to simulate reactivity in high (argon) and low (air) altitude environments. The projectiles were launched into a chamber that included a steel target plate for projectile penetration before impacting a rear witness plate. The chamber was semi-sealed and instrumented for quasi-static pressure data. The results provide an understanding of energy release from the projectile materials and of the environmental influence on performance. The transient pressure traces provide insight into reaction kinetics. A bifurcation in transient pressure rise was an indication of a shift in reaction kinetics from the inherent reactive material to metal oxidation with the environment. The bifurcation was delayed by about 0.15 ms for the intermetallic relative to the thermite, evidence that the thermite reaction proceeded faster upon impact than the intermetallic. The two-step process (impact ignition of the reactive material followed by metal oxidation) was shown to produce higher energy conversion efficiencies than projectiles composed of pure fuel (i.e., aluminum) reported previously. Both reactive materials showed energy conversion efficiencies greater than 30% (for air) and 50% (for argon), and an explanation of underestimated efficiency and energy losses is provided. These results have implications for advancing formulations for ballistic applications. Structural reactive materials can be used to modify the effective reactivity of metal-containing formulations in varied atmospheric environments.

36 MATERIALS SCIENCE↗

Forward modeling of Doppler velocity interferometer system for improved shockwave measurements

We report that when near-instantaneous shocks are recorded by using a Doppler velocity interferometer velocity interferometer system for any reflector (VISAR), they typically exceed the detector’s ability to react and “skipped fringes” result, where its visibility briefly reduces. Traditionally, replacing skipped fringes required guesswork in analysis, which increased arrival time errors. Second, the use of long but velocity-sensitive interferometer delays with fast detectors which can resolve the delay has traditionally been avoided because of the fear of confusing the arrival time signal. However, shorter delays produce smaller fringe phase shifts per velocity and, thus, can decrease velocity precision. We realize that while some loss of fringe information occurs at shock events, this is often just a partial loss and the residual fringe information can still hold valuable information. We describe a forward model (FM) of the interferometer action and detector blurring that assists with VISAR fringe analysis at skipping events: (1) more precise shock arrival times, (2) arrival time precision not limited by long delays, and (3) improved ghost subtraction, which improves accuracy over a broad time region. Lastly, we demonstrate the utility of a FM on National Ignition Facility and Laboratory for Laser Energetics Omega shots and discuss data process methods to increase the high time resolution of VISAR systems despite their limitations and complexities.

47 OTHER INSTRUMENTATION↗

Actinide Crystallization and Fission Reactions in Cooling White Dwarf Stars

The first solids that form as a cooling white dwarf (WD) starts to crystallize are expected to be greatly enriched in actinides. This is because the melting points of WD matter scale as Z5/3 and actinides have the largest charge Z. We estimate that the solids may be so enriched in actinides that they could support a fission chain reaction. This reaction could ignite carbon burning and lead to the explosion of an isolated WD in a thermonuclear supernova (SN Ia). Our mechanism could potentially explain SN Ia with sub-Chandrasekhar ejecta masses and short delay times.

79 ASTRONOMY AND ASTROPHYSICS↗

Optical Investigation of Catalyst Heating Operation in a Light-Duty Compression-Ignition Engine

It is well known that cold start engine operation contributes a majority of tail-pipe emissions during a vehicle drive cycle. Aftertreatment catalysts have poor conversion efficiency during this period due to low conversion efficiency at temperatures less than their light-off temperature. A late post injection strategy, with a delayed injection to increase exhaust enthalpy, can help decrease the time required to heat-up the catalyst. However, late injections into conditions with decreasing temperature and density can result in high emissions and combustion instability. In this study a three-injection strategy with a late-post injection was investigated in an optically-accessible engine. For the first two injections (pilot and main) timing and duration were held constant as the post injection timing was swept and duration varied to hold load constant, at a relevant light-load catalyst heating operating condition. Simultaneous double-pass schlieren and OH chemiluminescence data were acquired in an optically-accessible 0.4 L single-cylinder engine to understand spray characteristics and combustion behavior for #2 diesel fuel. Optical images analysis revealed higher liquid penetration and lower vaporization of the liquid fuel jet as post injection is retarded. It is also observed that retarding post injection timing to +25 CAD resulted in lower observed soot concentration in the piston bowl.

99 GENERAL AND MISCELLANEOUS↗

Influence of Operating Conditions on Ethanol Passive Prechamber Cold-start Operation

Cold-start operation in spark ignition engines is characterized by many drawbacks such as poor fuel vaporization, substantial wall film formation, and weak ignition which can result in unstable combustion and high concentrations of carbon monoxide (CO), nitrogen oxides (NOx), and unburnt hydrocarbons (UHC) prior to catalyst light-off. These challenges are particularly important for ethanol as the fuel’s high latent heat of vaporization and low vapor pressure can hinder reliable ignition during cold-start. In this context, passive prechamber ignition systems offer a viable pathway to improve ignitability by producing hot turbulent jets that can promote combustion in the main chamber. In this study, the influence of operating conditions on jet formation and combustion characteristics in an ethanol-fueled passive prechamber spark ignition engine is investigated. Two intake pressures (40 kPa and 60 kPa) and two engine speeds (450 and 600 rpm) are studied to represent various stages of the cold-start ramp-up. Results indicate that lower intake pressures lead to delayed and asymmetric turbulent jet formation, and lower peak heat release rates. Lower engine speeds, on the other hand, produced higher peak heat release and faster combustion due to greater residence time in the engine. This greater residence time also increases the likelihood of autoignition in the main chamber, potentially promoting quicker heat release for lower engine speeds. These results provide insight into the role of operating conditions on cold-start ramp-up of passive prechamber engines.

Banagiri, Shrikar [ORNL] (ORCID:0000000239745099)↗

Simulation of Premixed and Partially Premixed Jet-in-Crossflow Flames at High Pressure

In this paper, we explore the operational map of a lean axial-staged combustor of premixed and partially premixed reacting jet-in-crossflow flames at high -pressure (5 atm). This study attempts to expand the data to relatively high pressure and could significantly aid scaling to real gas turbine engine conditions at 20–30 atm. High-speed camera, particle image velocimetry (PIV), CH* chemiluminescence, temperature, and pressure measurements were taken and processed to allow accurate reconstruction of six operating points relative to computational fluid dynamics (CFD) simulations under minimal adjustments. Variation of lean main stage (φ = 0.575 and 0.73) and rich jet (φ = 1.1, 4, and 8) equivalence ratio has been investigated for a four mm axial jet. The fully premixed flames were found to be controlled by the crossflow temperature before ignition and the crossflow oxygen content during combustion. Analysis of flame shape and position for the partially premixed operating points describes a lee stabilized as well as a more unsteady windward flame branch. Adjustment of added jet fuel and crossflow temperature along with its corresponding oxygen level is required to attain a compact flame body. Here, the risk of delaying combustion progress is significantly increased at a richer jet φ = 8 and an overshooting, spatially divided flame was attained with a main stage φ = 0.73. Control toward a compact flame body is critical to allow combustion at reasonable reaction rate.

03 NATURAL GAS↗

Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR

Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.

computational fluid dynamics simulation↗

Silicon alloying enhances fast heating rate combustion of aluminum particles

Alloys and composites of aluminum (Al) have shown promise in regulating and enhancing particle combustion for energy generation applications. Recent work on aluminum-silicon (Al-Si) spherical alloy particles has shown improved combustion at low heating rates through enhanced diffusion accompanying a lower melting temperature. This study extends reactivity analysis to higher heating rates, comparing oxidation of Al-Si with Al. Flame speeds of Al-Si powder mixed with molybdenum trioxide (MoO 3 ) powder (Al-Si+MoO 3 ) exhibited a faster transition to steady propagation relative to Al+MoO 3 . Bomb calorimetry experiments revealed up to 5.8% greater early temperature rise for Al-Si powder. Rapid steady propagation for Al-Si mixture and faster temperature rise for Al-Si particles were attributed to accelerated kinetics evidenced in thermal equilibrium analysis of the mixtures using a differential scanning calorimeter (DSC). Larger DSC exotherms in the early stages of oxidation (i.e., 480–720 °C) correlated with early heat release and promoted steady flame propagation for Al-Si+MoO 3 compared with Al+MoO 3 . Furthermore, in multiple heating rate DSC studies, Al-Si+MoO 3 ignited at a lower heating rate (i.e., 15 ⁰C/min) than Al+MoO 3 (i.e., 20 ⁰C/min). Both mixtures ignited at temperatures less than 625 °C which is above the melting temperature of Al-Si (574 °C) and below Al (659 °C). Thus, alloy particle fuels with enhanced diffusion-controlled kinetics promote steady flame propagation and show promise for energy generation applications. This is especially promising for technologies driven by enhanced condensed phase combustion such as some propellant additives or applications dependent on reliable burn rates such as pyrotechnic time delay formulations and primers.

36 MATERIALS SCIENCE↗

CFD modeling of near-wall combustion and unburned methane prediction in natural gas spark ignition engines

Natural gas-powered engines play a critical role in gas drilling, compression, and transmission sectors, but methane (CH 4 ) from engine combustion slip can be significant over their lifespan, contributing to atmospheric pollution and signaling reduced engine efficiency. Here, to address this challenge, computational fluid dynamics (CFD) simulations offer valuable insights into the in-cylinder combustion process, enabling the optimization of combustion strategies and engine designs to minimize unburned CH 4 slip. This study aims to evaluate and improve combustion models for simulating the combustion process and predicting unburned CH 4 concentrations in natural gas spark-ignition (SI) engines, including engines that are part of combined reformer-engine systems. Specifically, the performance of two flamelet-based combustion models—the Extended Coherent Flame Model (ECFM) and the G-equation model—was assessed using experimental engine data collected under varying excess-air ratio (λ) conditions and fuel compositions, including natural gas and syngas blends. In addition, to enhance the predictive capabilities of the G-equation model, a flame-wall interaction (FWI) sub-model was integrated into its framework. The effects of its model parameters, such as quenching and influence distance, on combustion behavior and unburned methane predictions were analyzed in detail. The ECFM tended to predict delayed combustion phasing under diluted mixture conditions, resulting in overprediction of unburned CH 4 concentrations. In contrast, the G-equation model provided reasonable predictions of combustion pressure, while representing higher the CH 4 reduction rate across the operating condition compared to experimental data. Incorporating the FWI sub-model—with the quenching distance calculated based on a pressure-dependent relation (P -0.48 ) and a fixed influence distance of 1.5 mm—further improved the G-equation model’s accuracy in predicting CH 4 reduction rates without compromising its ability to simulate the combustion process.

Combustion model↗