Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulent combustion”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2026

Data-based filtered dissipation rate modelling for multi-modal turbulent combustion: evaluating a priori model generalizability

Manifold-based models offer a computationally efficient alternative to directly transporting the thermochemical state in computational simulations of turbulent reacting flows, projecting the high-dimensional thermochemical state-space onto a low-dimensional manifold. Recent efforts have yielded a manifold-based model applicable to multi-modal combustion, enabling reconstruction of the thermochemical state from solutions to two-dimensional manifold equations in mixture fraction and generalized progress variable that are parameterised by three scalar dissipation rates. In coarse-grained simulations such as Large Eddy Simulation (LES), closure of the multi-modal manifold equations and subfilter variances/covariance requires closure of three filtered scalar dissipation rates. Here, the present work adopts a data-based approach, providing closure for the three filtered scalar dissipation rates via deep neural networks (DNNs). High-fidelity datasets corresponding to an autoigniting n-dodecane jet flame and a bluff body swirl-stabilized confined lifted spray flame of two aviation fuels (Jet-A and C1) with different ignition propensities are leveraged to generate training data that spans a diverse range of thermodynamic conditions and combustion modes, including low- and high-temperature ignition regimes in addition to premixed and nonpremixed behaviour. A final DNN model is trained to enforce inherent physical constraints by learning nonlinear functional transformations of the three filtered scalar dissipation rates. The generalizability of this constrained DNN model is demonstrated a priori via conditional statistics evaluated on the lifted spray flame with C1–a configuration that had not been included in the training data. Excellent DNN agreement with conditional DNS statistics is observed, and integrated gradients are computed to identify the most sensitive input variables. The similarity of the marginal PDFs of the most informative input variables and outputs across configurations are quantified via the Wasserstein metric, demonstrating that data-based models may successfully generalize to unseen parametric conditions so long as the most informative input variables share similar distributions across training and testing datasets.

Data-based modelling↗
DOE OSTI2024

Hierarchically Informed Engineering Models for Predictive Modeling of Turbulent Premixed Flame Propagation in Pre- chamber Turbulent Jet Ignition

The goal of the project is to improve the predictive accuracy and efficiency of turbulent combustion sub-models for pre-chamber turbulent jet ignition (TJI). This goal is achieved through the development of a hierarchically informed engineering model for turbulent combustion in TJI. The model development starts with the highest level of model description of turbulent combustion with direct numerical simulation (DNS) from which fundamental characteristics and scaling properties of turbulent premixed flame propagation under TJI relevant conditions are obtained.

42 ENGINEERING↗
DOE OSTI2026

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↗
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↗
DOE OSTI2026

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗
DOE OSTI2025

In situ multi-tier auto-ignition detection applied to dual-fuel combustion simulations

Here we use an anomaly detection methodology that is centered on analyzing fourth-order joint moments (co-kurtosis), particularly focusing on its application in auto-ignition of combustion problems with large numbers of species. Unsupervised anomaly detection is challenging to generalize across problem types and domains. A recent technique, centered on analyzing information in the fourth-order joint moment co-kurtosis, has shown promise, especially for high-dimensional scientific data. In this work we present developments to the co-kurtosis based anomaly detection method needed to make it effective and scalable for large-scale distributed scientific data, such as those generated by massively parallel simulations. An in situ co-kurtosis algorithm is employed as the anomaly detection method for identifying ignition kernels in simulations of turbulent combustion. Here, we extend an existing methodology which identifies regions of the domain where anomalies are present, and add another tier of anomaly detection where the individual samples contributing to the anomaly are identified. We apply this algorithm on-the-fly to a variety of turbulent reacting flow problems and compare it to the widely used (but significantly more expensive) chemical explosive mode analysis (CEMA). We demonstrate the ability of the method to detect and identify the onset of low and high temperature ignition which can be used for computational steering, as chemical and combustion anomalies occur intermittently at spatio-temporal locations unknown a priori. Finally, we apply our lightweight in situ algorithm to an exascale high-fidelity simulation with a total of 2.4 Trillion degrees of freedom, performed using an adaptive mesh refinement solver. Furthermore, through a scalability analysis, we show that the relative computational cost of this in-situ anomaly detection algorithm compared to an iteration of the reacting flow solver is negligible.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Turbulent burning velocity of lean premixed hydrogen/air flames at engine conditions: Effects of turbulence intensity and length scale

For turbulent lean premixed hydrogen flames with strong thermodiffusively instabilities, most previous studies have focused on the influence of turbulence intensity, whereas the role of turbulence length scale is less well understood. Here, this study addresses this gap by conducting direct numerical simulations (DNS) of statistically planar turbulent premixed flames for a lean (ϕ=0.35) hydrogen/air mixture under independently varied turbulence intensity (u') and length scale (l T ) at engine-relevant thermodynamics conditions. Results show that as u' increases, the flame front becomes increasingly wrinkled, forming smaller cellular structures. In contrast, l T variations do not significantly alter the size of these structures. For the turbulent burning velocity (S T ), the normalized S T (i.e., S T /S L , where S L is the laminar flame speed) increases linearly with u', driven by both enhanced flame surface wrinkling (i.e., increased A T /A L ) and enhanced local burning rate (i.e., increased I 0 ). However, increasing l T reduces I 0 , despite a continued increase in A T /A L , resulting in only a marginal increase in S T /S L . To reveal the underlying mechanisms, especially the decreasing trend of I 0 with l T , local flame dynamics analyses are performed. It is found that as l T increases, the interaction between thermodiffusive effects and turbulence weakens due to the reduced tangential strain rate, while the flame curvature remains largely unchanged. This suppresses local reactivity enhancement and thus decreases I 0 , In contrast, an increase in u' enhances the interaction by amplifying both curvature fluctuation and tangential strain rate, leading to increased local reactivity (increased I 0 ). Finally, based on the DNS data, several new scaling models are proposed for the three global properties, S T /S L , A T /A L , and I 0 , and showed improvements compared to existing models. These findings provide new insights into the flame-turbulence interactions in thermodiffusively unstable hydrogen flames. The DNS dataset is also useful for the development of turbulent combustion models applicable to practical engine simulations.

Engine-relevant condition↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Enclosed Sydney Swirl Burner Experimental Data Set - Release 1.0

This submission includes a comprehensive data set for the NETL Enclosed Sydney Swirl Burner, generated between 2012 and 2014. The data set includes global parametric data and advanced diagnostic data for (3) flames of interest to CFD model validation. The submission includes a companion summary document describing the data and associated references.

Diagnostics↗
DOE OSTI2025

Towards stabilization of pressurized oxy-coal diffusion flames

Designing an effective burner is vital for the development of pressurized oxy-fuel combustion technologies. Turbulent jet diffusion burners are adopted for the pressurized oxy-combustor in this research and a bluff-body is employed to create a recirculation zone to stabilize the flame in pressurized oxy-combustor. The objective of this numerical work is to perform a systematic analysis of the characteristics of this pressurized non-premixed flame. Specifically, a 15-bar POFC combustor of power 100 kWth is modeled employing the Ansys FLUENT commercial platform, using Reynolds-averaged Navier-Stokes (RANS) modeling. The present computational work focuses on identifying the aerothermodynamic features of the isothermal and reacting flow with a disk-shape bluff-body. Results show that fuel-oxidizer stream momentum ratio has great impacts on the temperature profile of the down-fired, co-axial, pressurized oxy-coal, diffusion flame. Parametric study of the width of the burner tip identifies an optimal range for this pressurized burner.

Li, Lei↗
DOE OSTI2026

Pixel-Resolved Long-Context Learning for Turbulence at Exascale: Resolving Small-scale Eddies Toward the Viscous Limit

Turbulence plays a crucial role in multiphysics applications, including aerodynamics, fusion, and combustion. Accurately capturing turbulence's multiscale characteristics is essential for reliable predictions of multiphysics interactions, but remains a grand challenge even for exascale supercomputers and advanced deep learning models. The extreme-resolution data required to represent turbulence, ranging from billions to trillions of grid points, pose prohibitive computational costs for models based on architectures like vision transformers. To address this challenge, we introduce a multiscale hierarchical Turbulence Transformer that reduces sequence length from billions to a few millions and a novel RingX sequence parallelism approach that enables scalable long-context learning. We perform scaling and science runs on the Frontier supercomputer. Our approach demonstrates excellent performance up to 1.1 EFLOPS on 32,768 AMD GPUs, with a scaling efficiency of 94\%. To our knowledge, this is the first AI model for turbulence that can capture small-scale eddies down to the dissipative range in three-dimensional turbulence at high Reynolds numbers.

Yin, Junqi [ORNL] (ORCID:0000000338435520)↗
DOE OSTI2025

Exploring the Effects of Varying Pre-Chamber Geometry in a Heavy-Duty Natural Gas Optical Engine under Dilution Conditions

Pre-chamber combustion is an advanced ignition strategy that has been shown to enhance spark ignition (SI) combustion stability in natural gas (NG) engines by providing distributed ignition sites from turbulent jets and enhancing main-chamber turbulence. Pre-chamber combustion has been proven especially advantageous compared to SI in ultra-lean and dilute operating conditions. This work involves experimental investigation of the effects of varying passive pre-chamber nozzle configuration on pre-chamber and main chamber combustion under simulated exhaust gas recirculation (EGR) dilution (0 and 20%) conditions in a heavy-duty, single-cylinder, optically accessible NG engine at stoichiometric fuel-air ratio. Pre-chamber nozzle configurations include four pre-chambers with constant nozzle area to pre-chamber volume ratio (A/V) with different nozzle sizes and orientations and one configuration with larger nozzles. The optical engine is operated in a skip-fire sequence consisting of 18 motored cycles followed by two consecutive fired cycles to elucidate the effect of combustion residuals (internal EGR) on combustion evolution. Pressure-based diagnostics are used to monitor pre-chamber and subsequent main chamber combustion, and optical diagnostics include high-speed OH* chemiluminescence to visualize the development of pre-chamber jets and the resulting ignition of the main chamber charge. Heat release analysis of the in-cylinder pressure data indicates that a faster pre-chamber pressure rise does not always translate into faster main-chamber combustion. The pre-chamber with the smallest nozzle diameter produced the highest pre-chamber pressure rise and fastest combustion under non-diluted conditions. However, dilution delays the main chamber ignition for smaller nozzles despite a comparable rise in pre-chamber pressure compared to configurations with larger nozzles. This effect is more pronounced for cycles with in-cylinder combustion residuals in addition to external dilution. Additionally, it was observed that pre-chambers with swirling nozzles have a faster pressure rise in the pre-chamber and main chamber under dilute conditions. Optical diagnostics suggest that the main reason for the delay between the pre-chamber pressure rise and main-chamber combustion lies in jet quenching and delayed re-ignition, which can even lead to misfire if jets emitted from small nozzles combined with dilution fail to re-ignite.

Dhotre, Akash [University Of Minnesota-Twin Cities↗
DOE OSTI2025

FLEX-FUEL MIXING CONTROLLED COMBUSTION ENABLED BY PRECHAMBER IGNITION

There is an imminent need to displace fossil diesel fuel with cleaner burning, domestically produced, renewable fuels for use in heavy-duty engines. Bioethanol is a prime candidate as it widely adopted in the U.S. as a gasoline additive ranging in volume percentage from 10% (E10) up to 85% (E85). Direct substitution of market available ethanol-gasoline blends for diesel fuel is not plausible as the stark reactivity differences would not constitute the same ignition quality nor achieve auto-ignition at all. This work focuses on the development of prechamber enabled mixing-controlled combustion (PC-MCC) as an advanced combustion strategy to facilitate reliable ignition and diffusion style combustion ethanol-gasoline fuel blends. PC-MCC involves integration of an actively fueled prechamber (PC) into a conventional compression ignition combustion system. When ignited, the PC ejects hot turbulent jets into the main combustion chamber that then interact with the direct injected fuel, prompting immediate ignition. The PC jet flames provide a robust thermal ignition source that allows the engine to operate agnostic of fuel composition, or flex-fuel. Computational fluid dynamics (CFD) modeling was used to assess critical design features of the PC while garnering insights into the ignition strategies that facilitate robust performance. A key finding was the ignition performance benefits of fuel-rich PC operation which yield exothermic jets. Based on the numerical findings, a prototype igniter was tested experimentally on both single and multi-cylinder engine platforms at a variety of operating conditions. The experimental results indicate flex-fuel PC-MCC is well capable of diesel-like combustion processes by demonstrating matched or improved gross thermal efficiencies and load variability within 2%. Fuel grade ethanol (E98) exhibited consistently lower NOx and immeasurable soot across the load space. E98 also demonstrated a significant improvement in thermal efficiency at light loads.

Zeman, Jared↗
DOE OSTI2025

Multi-modal effects on indirect noise induced by turbulent entropy fields

Planar entropy waves are commonly assumed for predicting indirect combustion noise. However, the non-planar and turbulent nature of flows found in most practical combustors challenges this assumption. In the present paper, we examine the indirect noise generated by non-planar and turbulent entropy fields through subsonic nozzles. Firstly, we introduce a new transfer function framework that accounts for the contribution of non-planar Fourier modes of the entropy field to the indirect noise spectra. When applied to a turbulent flow field, this method demonstrates a significant improvement in spectral predictions compared with a conventional approach that only considers the planar mode. Secondly, simulations show that non-planar Fourier modes become significant above a threshold frequency ƒ thresh , found in the mid- to high-frequency range. This contribution of non-planar modes is explained by two-dimensional shear effects that distort the entropy waves. A scaling relation that uses residence times along streamlines is developed for ƒ thresh , showing good agreement with simulation results. Finally, we show that the indirect noise from non-planar entropy modes found in aviation combustors can be significant at frequencies below 1 kHz, which might be relevant in situations of thermo-acoustic instabilities coupled to indirect noise.

42 ENGINEERING↗
DOE OSTI2026

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)↗