Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Combustion modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Characterization of Flashback and Flame-Holding in a Jet-in-Crossflow Mixing Configuration With Methane–Hydrogen Fuel Blends

An approach that combines experimental and numerical analyses has been implemented to characterize the fundamentals of flashback events and flame-holding phenomena during high-hydrogen combustion in a jet-in-crossflow (JICF) configuration. Such flame dynamics are visualized experimentally using nanosecond (ns)-based hydroxyl planar laser-induced fluorescence (OH-PLIF) and chemiluminescence diagnostics techniques. The JICF burner has an optically accessible premixing tube allowing the optical diagnostics. The testing was conducted for varied premixer velocities (V) and equivalence ratio (ϕ) for 90–100% (H 2 , by mole) H 2 /CH 4 reactant mixtures at atmospheric temperature and pressure conditions. Two distinct flashback events were identified – conventional rich flashback and lean flashback, recorded while increasing ϕ and decreasing ϕ, respectively. The cause of lean flashback was attributed to lower momentum flux ratio which bends the jet sharply, closer toward the injection plane. The mean OH-PLIF images characterized the flame-holding behavior, where the flame was found to be stabilized on the leeward side only or on both windward and leeward sides as a lifted flame near the fuel port. A large eddy simulation (LES) with detailed chemistry combustion modeling approach was implemented along with an OH* submechanism, and it showed qualitative agreement with the integrated line-of-sight chemiluminescence results as well as the planar OH-PLIF measurements.

chemiluminescence diagnostic↗
DOE OSTI2020

Front propagation formulation for large eddy simulation of turbulent premixed flames

The front propagation formulation (FPF) for large eddy simulation (LES) of turbulent premixed flames is presented and discussed. The FPF method adopts a regularized Dirac delta function to minimize the spurious propagation of under-resolved filtered flame fronts in LES, while reproducing the reaction characteristics of premixed flame fronts. Here, the extension of the FPF method to LES and the relation to existing sub-filter combustion models are discussed. An algebraic model for the regularized Dirac delta function is presented and validated using direct numerical simulation (DNS) data for statistically planar premixed flames propagating in homogeneous isotropic turbulence. A posteriori validation for a filtered 1-D planar laminar flame is performed to demonstrate the characteristics of FPF in reproducing the propagation speed and the thickness of the under-resolved flame front. LES of a turbulent Bunsen flame is also performed to demonstrate the performance of the proposed method.

42 ENGINEERING↗
DOE OSTI2022

Quantification of Key Peroxy and Hydroperoxide Intermediates in the Low-Temperature Oxidation of Dimethyl Ether

Dimethyl ether (DME) oxidation is a model chemical system with a small number of prototypical reaction intermediates that also has practical importance for low-carbon transportation. Although it has been studied experimentally and theoretically, ambiguity remains in the relative importance of competing DME oxidation pathways in the low-temperature autoignition regime. To focus on the primary reactions in DME autoignition, we measured the time-resolved concentration of five intermediates, CH 3 OCH 2 OO (ROO), OOCH 2 OCH 2 OOH (OOQOOH), HOOCH 2 OCHO (hydroperoxymethyl formate, HPMF), CH 2 O, and CH 3 OCHO (methyl formate, MF), from photolytically initiated experiments. We performed these studies at P = 10 bar and T = 450-575 K, using a high-pressure photolysis reactor coupled to a time-of-flight mass spectrometer with tunable vacuum-ultraviolet synchrotron ionization at the Advanced Light Source. Our measurements reveal that the timescale of ROO decay and product formation is much shorter than predicted by current DME combustion models. The models also strongly underpredict the observed yields of CH 2 O and MF and do not capture the temperature dependence of OOQOOH and HPMF yields. Adding the ROO + OH $\rightarrow$ RO + HO 2 reaction to the chemical mechanism (with a rate coefficient approximated from similar reactions) improves the prediction of MF. Increasing the rate coefficients of ROO ↔ QOOH and QOOH + O 2 ↔ OOQOOH reactions brings the model predictions closer to experimental observations for OOQOOH and HPMF, while increasing the rate coefficient for the QOOH $\rightarrow$ 2 CH 2 O + OH reaction is needed to improve the predictions of formaldehyde. In conclusion, to aid future quantification of DME oxidation intermediates by photoionization mass spectrometry, we report experimentally determined ionization cross-sections for ROO, OOQOOH, and HPMF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

On the subject of large-scale pool fires and turbulent boundary layer interactions

The role to which a realistic inflow turbulent boundary layer (TBL) influences transient and mean large-scale pool fire quantities of interest (QoIs) is numerically investigated. High-fidelity, low-Mach large-eddy simulations that activate low-dissipation, unstructured numerics are conducted using an unsteady flamelet combustion modeling approach with mutiphysics coupling to soot and participating media radiation transport. Three inlet profile configurations are exercised for a large-scale, high-aspect rectangular pool that is oriented perpendicular to the flow direction: a time-varying, TBL inflow profile obtained from a periodic precursor simulation, the time-mean of the transient TBL, and a steady power-law inflow profile that replicates the mean TBL crosswind velocity of 10.0 m/s at a vertical height of 10 m. Results include both qualitative transient flame evolution and quantitative flame shape with ground-level temperature and convective/radiative heat flux profiles. While transient fire events, which are driven by burst-sweep TBL coupling, such as blow-off and reattachment are vastly different in the TBL case (contributing to increased root mean square QoI fluctuation prediction and disparate flame lengths), mean surface QoI magnitudes are similar. Quadrant analysis demonstrates that the TBL configuration modifies burst-sweep phenomena at windward pool locations, while leeward recovery is found. In conclusion, positive fluctuations of convective heat flux correlate with fast moving fluid away from the pool surface due to intermittent combustion events.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

On the modes of nanosecond pulsed plasmas for combustion ignition of quiescent CH 4 -air mixtures

The effect of transient plasma modes on ignition kernel development are discussed here for a quiescent CH 4 -air combustion model system. A 10 ns high-voltage pulse was applied to a pin-to-pin electrode in lean fuel-air mixtures at room temperature and atmospheric pressure. High-impedance streamer, transient spark and low-impedance spark discharges were identified based on pulse waveforms of voltage and current. A sustained ignition kernel expansion was observed when the plasma discharge transitioned into a transient spark or spark discharge. The minimum ignition energy was obtained at the transient spark mode, which has less than a third of the energy or Coulomb transfer compared to the low-impedance spark. Employing repetitive 10-pulse sequence at 10 kHz, the lean-fuel limit was extended from an equivalence ratio of 0.6 for the single pulse ignition to 0.5. The use of repetitive pulses also allowed streamer breakdown or spark initiation to occur at a lower voltage.

CH 4↗
DOE OSTI2024

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Conjugate heat transfer↗
DOE OSTI2021

Transformation of an Embedded Five-Membered Ring in Polycyclic Aromatic Hydrocarbons via the Hydrogen-Abstraction–Acetylene-Addition Mechanism: A Theoretical Study

Five-membered rings are constituents of many polycyclic aromatic hydrocarbons (PAHs) and their presence on the edges of large PAHs has been repeatedly observed experimentally. However, modern kinetic combustion models often do not consider the growth of PAHs through the transformation of the five-membered rings. Here, in connection with the above, we carried out a theoretical study of the mechanism of Hydrogen-Abstraction-aCetylene-Addition (HACA) transformation of an embedded five-membered ring on the armchair PAH edge to a six-membered ring, considering cyclopenta[d,e,f]phenanthrene (4,5-methylenephenanthrene) as a prototype system for this process. The potential energy surface for the reactions of cyclopenta[d,e,f]phenanthrenyl radicals produced by direct H abstractions from cyclopenta[d,e,f]phenanthrene with acetylene has been compiled at the G3(MP2,CC)//B3LYP/6-311G(d,p) level of theory including zero-point vibrational energy corrections. The computed energies and molecular parameters were then used to solve the Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) in order to calculate reaction rate at various pressures and temperatures, which were fitted to modified Arrhenius equation for further kinetic modeling. The results show that the HACA transformation of the embedded five-membered ring to a six-membered ring is possible, albeit slow. The most viable reaction mechanism involves the R2 + C 2 H 2 reaction, where the acetylene molecules adds to a σ-radical in the six-membered ring adjacent to the five-membered ring via a low entrance barrier. The predominant product of R2 + C 2 H 2 is predicted to be 3 ethynyl 4H cyclopenta[def]phenanthrene Pr5 via immediate H elimination from the initial addition complex. Next, Pr5 undergoes H-assisted isomerization to 4aH-pentaleno[4,3,2,1-cdef]phenanthrene Pr4, and the latter adds an H atom eventually forming the 1-pyrenylmethyl radical Pr3: R2 + C 2 H 2 ⇌ 3 ethynyl 4H cyclopenta[def]phenanthrene (Pr5) + H ⇌ 4aH-pentaleno[4,3,2,1-cdef]phenanthrene (Pr4) + H Pr5 + H ⇌ Pr4 + H Pr4 + H → 1-pyrenylmethyl (Pr3) This HACA sequence may be competitive with the methyl radical addition to the R1 radical formed by H abstraction from the CH 2 group in the five-membered ring of cyclopenta[d,e,f]phenanthrene, which provides a pathway to pyrene following two H atom losses. Relative contributions of the two mechanisms of the five- to six-membered ring transformation would strongly depend on the branching ratios of the R1 and R2 radicals produced by the H abstractions and the available concentration of C 2 H 2 vs. CH 3 and hence differ in different flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

The Reaction of Atomic Hydrogen with an Imine: Direct Kinetic Measurements and Computational Studies

Imines are of interest in the combustion modeling of nitrogen-based fuels such as ammonia, as intermediates in atmospheric oxidation of amines, and may offer a new class of environmentally benign substitutes for perhalogenated reagents. Elementary gas-phase rate constants for the consumption of H atoms by 1,1,1,3,3,3-hexafluoro-2-propanimine were measured over 294–736 K using the laser flash photolysis/atomic resonance fluorescence technique. The results are summarized as 3.0 × 10 –11 exp(−14.5 kJ mol –1 /RT) cm 3 molecule –1 s –1 with a confidence interval of ±13%. These are the first temperature-dependent kinetic experiments on an imine of the form R 2 C═NH and thereby provide the only validation of transition-state modeling for such systems. Several pathways are feasible; comparison of theory and experiment suggests the dominant channel is addition of H atoms to the nitrogen atom. On this basis, the observations are rationalized quantitatively. Computations also yield the thermochemistry of (CF 3 ) 2 CNH, with Δ f H 298 = −1232 kJ mol –1 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Active Thermochemical Tables: the thermophysical and thermochemical properties of methyl, CH 3 , and methylene, CH 2 , corrected for nonrigid rotor and anharmonic oscillator effects

The thermophysical properties (isobaric heat capacity, entropy, enthalpy increment) of two prominent radicals, methyl, CH 3 , and methylene, CH 2 , were computed using the Nonrigid Rotor Anharmonic Oscillator (NRRAO) approach and compared to their RRHO counterparts, demonstrating significant differences between the results from the two approaches. Methylene presents a typical case in which the NRRAO thermophysical properties have significantly higher values than their RRHO counterparts at higher temperatures. In the case of methyl, the positive anharmonicity of the umbrella motion causes an opposite effect, and the NRRAO corrected thermophysical properties have lower values than their RRHO counterparts. The NRRAO corrected thermophysical properties, in turn, affect the resulting thermochemical properties. Here, two reactions important in combustion modelling were tested: the recombination of methyl radical with hydrogen atoms to form methane, and the recombination of two methyl radicals to form ethane. The related NRRAO equilibrium constants differ significantly from their RRHO analogs, and the consequences for chemical modelling are discussed. Also reported are the most current ATcT enthalpies of formation for CH n (n = 4-0) species and for C 2 H 6 , together with the tightly related sequential bond dissociation enthalpies along with the CHn series.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Simulations of Multi-Mode Combustion Regimes Realizable In a Gasoline Direct Injection Engine

Lean and dilute gasoline compression ignition (GCI) operation in spark ignition (SI) engines are an attractive strategy to attain high fuel efficiency and low NOx levels. However, this combustion mode is often limited to low-load engine conditions due to the challenges associated with autoignition controllability. In order to overcome this constrain, multi-mode (MM) operating strategies, consisting of advanced compression ignition (ACI) at low load and conventional SI at high load, have been proposed. In this three-dimensional computational fluid dynamics study, the concept of multi-mode combustion using two RON98 gasoline fuel blends (Co-Optima Alkylate and E30) in a gasoline direct injection (GDI) engine were explored. To this end, a new reduced mechanism for simulating the kinetics of E30 fuel blend is introduced in this study. To cover the varying engine load demands for multi-mode engines, primary combustion dynamics observed in ACI and SI combustion modes was characterized and validated against experimental measurements. In order to implement part-load conditions, a strategy of mode transition between SI and ACI combustion (i.e., mixed-mode combustion) was then explored numerically by creating a virtual test condition. The results obtained from the mixed-mode simulations highlight an important feature that deflagrative flame propagation regime coexists with ignition-assisted end-gas autoignition. This study also identifies a role of turbulent flow property adjacent to premixed flame front in characterizing the mixed-mode combustion. The employed hybrid combustion model was verified to perform simulations aiming at suitable range of multi-mode engine operations.

Kim, Sayop↗
DOE OSTI2023

An Integrated Approach to Predicting Ash Deposition and Heat Transfer in Coal-Fired Boilers

The overall goal of this project is to develop via measurements and simulations an advanced online technology to predict, monitor and manage fireside ash deposition in a coal-fired boiler allowing for more efficient operations under a range of load conditions and fuel property variability. With this in place fuel sorting and blending can be done upstream and operations can be optimized to compensate for load and fuel properties. In support of this objective, three experimental campaigns were undertaken during the course of the project to measure ash deposition rates within the boiler at different fuel flow rates and its ash composition variability. Simulations of the experimental conditions representing actual geometry, operational scenarios in terms of air flow rates, coal flow rates as well as coal compositions, heating values, and particle size distributions were also carried out. Deposition rates were predicted using a unique particle kinetic energy and viscosity based ash deposition methodology whose validity was ascertained by comparing against deposition rate measurements for widely varying operating conditions and ash compositions in a lab-scale furnace. With a unique end-to-end combustion modeling methodology established and different simulation scenarios carried out, the results from our computational fluid dynamic (CFD) simulations in conjunction with the plant data summarized in this report were used to refine Microbeam Technology Incorporated’s MTI CSPI-CT Tool to predict and monitor fire-side ash deposition under a range of load conditions and fuel property variability in real time.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Independent Fuel Property Effects of Fuel Volatility on Low Temperature Heat Release and Fuel Autoignition (Final Report)

This Cooperative Research and Development Agreement (CRADA) project between Argonne National Laboratory (ANL), Oak Ridge National Laboratory (ORNL), and Shell Global Solutions (Shell) was initiated as part of a Directed Funding Opportunity (DFO) call for proposals from the Co-Optimization of Fuels and Engines (Co-Optima) initiative. Shell had observed that volatile fuels suppress low temperature heat release (LTHR) more than expected based on conventional gasoline autoignition metrics: research octane number (RON) and motor octane number (MON). The role of LTHR contributes to autoignition phenomena for both boosted spark ignition (BSI) and advanced compression ignition (ACI) combustion modes. ACI combustion modes are applicable to large engines in the hard-to-electrify applications such as off-road, rail, and marine. Thus, having a reliable understanding of autoignition phenomena, including being able to accurately account for the effects of fuel volatility, is particularly important as new synthetic and bio-fuel compositions are considered.This CRADA project aimed to test the hypothesis that the decreased LTHR is due to preferential evaporation of multicomponent fuels when using direct injection (DI) fueling technology, creating composition and reactivity stratification. A custom set of fuels was designed and blended to test this hypothesis by Shell, with experimental engine studies at ORNL and engine combustion modeling by ANL. However, the initial experimental findings did not show the expected effect of fuel volatility suppressing LTHR. Instead, the LTHR propensity observed was independent of the fuel volatility. Due to the unexpected experimental result, the remainder of the experimental effort was redirected to study the effect of fuel volatility on emissions under spark-ignited cold-start conditions. However, as with the LTHR experiments, the cold start effort did not show a meaningful effect of fuel volatility on cold start emissions. Meanwhile, improved engine CFD models have been developed for both LTHR and cold start operations for the Shell fuels with different volatilities. While the simulation efforts were not pursued further due to the insignificant effects of fuel volatility as shown in experiments, the models developed can be easily retooled for off-road, rail, and marine applications.

09 BIOMASS FUELS↗
DOE OSTI2024

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Characterization of Flashback and Flame-Holding in a Jet-in-Crossflow Mixing Configuration with Methane-Hydrogen Fuel Blends

An approach that combines experimental and numerical analyses has been implemented to characterize the fundamentals of flashback events and flame-holding phenomena during high-hydrogen combustion in a jet-in-crossflow (JICF) configuration. Such flame dynamics are visualized experimentally using nanosecond (ns)-based hydroxyl planar laser-induced fluorescence (OH-PLIF) and chemiluminescence diagnostics techniques. The JICF burner has an optically accessible pre-mixing tube allowing the optical diagnostics. The testing was conducted for varied pre-mixer velocities (V) and equivalence ratio (ϕ) for 90%-100% (H2, by mole) H2/CH4 reactant mixtures at atmospheric temperature and pressure conditions. Two distinct flashback events were identified – conventional rich flashback and lean flashback, recorded while increasing ϕ and decreasing ϕ, respectively. The cause of lean flashback was attributed to the lower momentum flux ratio which bends the jet sharply, closer toward the injection plane. The mean OH-PLIF images characterized the flame-holding behavior where the flame was found to be stabilized on the leeward side only or on both windward and leeward sides as a lifted flame near the fuel port. A Large Eddy Simulation (LES) with detailed chemistry combustion modeling approach was implemented along an OH* sub-mechanism, and it showed qualitative agreement with the integrated line-of-sight chemiluminescence results as well as the planar OH-PLIF measurements.

flame holding↗
DOE OSTI2025

Accelerating computational fluid dynamics simulation of post-combustion carbon capture modeling with MeshGraphNets

Packed columns are commonly used in post-combustion processes to capture CO 2 emissions by providing enhanced contact area between a CO 2 -laden gas and CO 2 -absorbing solvent. To study and optimize solvent-based post-combustion carbon capture systems (CCSs), computational fluid dynamics (CFD) can be used to model the liquid–gas countercurrent flow hydrodynamics in these columns and derive key determinants of CO 2 -capture efficiency. However, the large design space of these systems hinders the application of CFD for design optimization due to its high computational cost. In contrast, data-driven modeling approaches can produce fast surrogates to study large-scale physics problems. We build our surrogates using MeshGraphNets (MGN), a graph neural network framework that efficiently learns and produces mesh-based simulations. We apply MGN to a random packed column modeled with over 160K graph nodes and a design space consisting of three key input parameters: solvent surface tension, inlet velocity, and contact angle. Our models can adapt to a wide range of these parameters and accurately predict the complex interactions within the system at rates over 1700 times faster than CFD, affirming its practicality in downstream design optimization tasks. This underscores the robustness and versatility of MGN in modeling complex fluid dynamics for large-scale CCS analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Evaluation of residual gas fraction estimation methods for cycle-to-cycle combustion variability analysis and modeling

Cycle-to-cycle combustion variability in spark-ignition engines during normal operation is mainly caused by random perturbations of the in-cylinder conditions such as the flow velocity field, homogeneity of the air-fuel distribution, spark energy discharge, and turbulence intensity of the flame front. Such perturbations translate into the variability of the energy released observed at the end of the combustion process. During normal operating conditions, the cycle-to-cycle variability (CCV) of the energy release behaves as random uncorrelated noise. However, during diluted combustion, in either the form of exhaust gas recirculation (EGR) or excess air (lean operation), the CCV tends to increase as dilution increases. Moreover, when the ignition limit is reached at high dilution levels, the combustion CCV is exacerbated by sporadic occurrences of incomplete combustion events, and the uncorrelation assumption no longer holds. The low or null energy released by partial burns and misfires has an impact on the following combustion event due to the residual gas that carries burned and unburned gases, which contributes to the deterministic coupling between engine cycles. Many residual gas fraction estimation methods, however, only address the nominal case where complete combustion occurs and combustion events are uncorrelated. Here we evaluate the efficacy of such methods on capturing the effects of partial burns and misfires on the residual gas estimate for high-EGR operation. The advantages and disadvantages of each method are discussed based on their ability to generate cycle-to-cycle estimates. Finally, a comparison between the different estimation techniques is presented based on their usefulness for control-oriented modeling.

42 ENGINEERING↗
DOE OSTI2023

SootLib: A soot model library for combustion simulation

Soot formation in combustion is an important process that affects radiative heat transfer, flame temperatures, and emissions with health and environmental impacts. Soot formation involves complex chemistry for nucleation, growth, oxidation, and coagulation processes. The soot particles vary widely in size and accurate modeling requires representation of the particle size distribution (PSD). Modeling soot is not trivial, and is only one of several physical processes active in combustion systems. This paper presents a software package called SootLib, which is an open-source library for modeling soot formation and other aerosol systems. SootLib is written in C++, is documented with Doxygen, and is available on GitHub. The library includes several models for soot chemistry and coagulation, and it represents the PSD using either a sectional model or the method of moments (MOM). Four closure approaches for the MOM are implemented allowing up to eight moments: monodispersed, an assumed-shape lognormal distribution, the quadrature method of moments, and the method of moments with interpolative closure. SootLib provides an interface for inclusion in other combustion packages including CFD or reacting flow solvers. The range of models allows comparisons and sensitivity studies, and the modularity facilitates extension to other soot models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

CFD modeling of non-catalytic, partial-oxidation engine reformer for flare mitigation

Flaring associated natural gas is commonly employed in the oil and gas industry to reduce methane (CH 4 ) emissions but generates carbon dioxide (CO 2 ) and harmful pollutants, significantly contributing to air pollution and posing risks to public health. To mitigate this impact, M2X Energy Inc. has developed a small-scale, modular gas-to-methanol system. This system features an engine reformer that performs fuel-rich partial oxidation of wellhead gas to produce syngas—a mixture of carbon monoxide (CO) and hydrogen (H 2 )—followed by a downstream reactor for methanol synthesis. This study focused on computational fluid dynamics (CFD) modeling of the engine reformer to simulate partial oxidation chemistry, predict the rich-burn operating limit, and assess syngas quality, ultimately aiding in design and operational optimization. The CFD model, developed within a Reynolds-Averaged Navier-Stokes (RANS) turbulence framework, incorporated sub-models for turbulent combustion, a chemical mechanism with polycyclic aromatic hydrocarbon (PAH) pathways, and soot emissions to accurately capture the fuel-rich, turbulent jet ignition and combustion processes. Model validation against experimental data showed good agreement across pre- and main-chamber pressures, apparent heat release rates, and exhaust gas concentrations of key species (H 2 , CO, CO 2 , CH 4 ) for varying intake equivalence ratios. Here, the model identified a rich-burn operating limit near a fuel-air equivalence ratio of 2.35, consistent with experimental observations. Furthermore, syngas quality analysis revealed that extending the rich-burn limit through engine reformer optimization could enhance syngas production, contributing to higher methanol synthesis efficiency.

Computational Fluid Dynamics↗