Search NASA⌕ Search

SEARCH · Search NASA

Results for “non-premixed”

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

Flame dynamics and kinetic coupling of ammonia and dimethyl-ether in non-premixed cool and warm flames at elevated pressure

Developing advanced low-temperature combustion engines with ammonia-biofuel blends requires a comprehensive understanding of low-temperature flame dynamics and kinetic interactions between ammonia and oxygenated fuels at elevated pressures. This work aims to study the dynamics and kinetics of non-premixed Dimethyl Ether (DME)/Ammonia (NH 3 ) cool and warm flames, and their reignition to hot flames. A counterflow burner is employed to establish DME/NH 3 cool/warm flames at pressures up to 5 atm. The extinction limits of cool flame and the reignition limits of warm flame to hot flame are measured by varying NH 3 concentrations and compared to simulations to quantitatively examine the effects on DME/NH 3 flames. It is found that NH 3 inhibits low-temperature DME oxidation and results in lower cool flame extinction limits. Warm flames in the presence of NH 3 are observed for the first time, revealing a non-monotonic effect of NH 3 addition: a small amount of NH 3 presence enhances warm flame chemistry and promotes reignition to hot flames, while a high NH 3 concentration weakens the warm flame. This trend is further explained by 0-D PSR kinetic simulations and 1-D S-curve flame dynamic calculations. Three flame transition regimes between cool flames (CF), warm flames (WF), and hot flames (HF) by different levels of NH 3 additions at a specific strain rate are identified, namely WFHF reignition, WF-CF transition, and WF extinction. Reaction sensitivity analyses of OH at low temperatures show that NH3 inhibits DME oxidation through OH consumption via H-abstraction and the kinetic couplings of RO 2 /NH 2 , RO 2 /NO x , R/NO x , and O 2 QOOH/NO x further suppress the low-temperature branching. At intermediate-temperatures, NH 2 /NO x /HO 2 coupling promotes warm flames via the pathway NH 2 → H 2 NO → HNO → NO by converting O 2 → HO 2 → OH. At even higher NH₃ concentrations, radical termination reactions of NH 2 + NO/NO 2 and excessive OH consumption via H-abstraction inhibit the flame. The insights into the kinetic coupling between NH 3 and low-temperature chemistry at elevated pressure and its impact on the dynamics of cool-warm-hot flame transitions will contribute to advancing combustion technologies with reduced emissions and improved energy-efficiency.

42 ENGINEERING↗
DOE OSTI2025

The effect on soot and its gas precursors of doping ethylene with 2,2,4,6,6-pentamethyl-heptane in the nitrogen-fuel stream of a laminar non-premixed Planar Mixing Layer Flame (PMLF)

Synthetic Aviation Turbine Fuels (SATFs) are promising for reducing soot emissions from the aviation sector and diversifying Jet Fuel (JF) sources. Accurately predicting the combustion and emissions behavior of SATFs (and other JFs) necessitates robust experimental databases to elucidate the chemistry of long-chain iso-paraffins, which can compose up to two-thirds of SATF blends and whose behavior is considered to be well-represented by that of iso-dodecane isomers. Here, this study characterizes two laminar non-premixed Planar Mixing Layer Flames (PMLFs) with mild soot loads fueled by nitrogen-diluted ethylene, pure and doped with 2,2,4,6,6-pentamethyl-heptane, respectively. The two PMLFs have the same stoichiometric mixture fraction and total hydrocarbon mole fraction in the fuel stream (X F,F =X C2H4,F +X C12H26,F = 0.260), resulting in nearly the same maximum temperature (T max ≈1800 K) and simple identification of the effects of doping. Importantly, any horizontal PMLF cross-section has a self-similar structure that can be modeled as an equivalent One-Dimensional Counterflow Flame (1D-CF) with vanishingly small strain rate (a). The cross-section at a Height Above the Burner (HAB) of 50 mm is characterized in terms of C 0 -C 18 gas species using capillary sampling followed by GC-MS analyses. Laser-Induced Emission Spectroscopy (LIES) quantifies the soot volume fraction (ƒ v ) profiles at HAB=25 and 50 mm where Elastic Laser Light Scattering (E-LLS) is performed to determine the a of the equivalent 1D-CFs and the profile of the E-LLS equivalent diameter ( d 6,3 ) of soot. The substitution of 1500 ppm of ethylene with 2,2,4,6,6-pentamethyl-heptane causes an increase of ≈1.5 in the concentrations of several polycyclic aromatic hydrocarbons and fv. Concurrently, the measured d 6,3 doubles in the oxidizer stream, yet remains the same in the fuel stream, at HAB=50 mm. Instead, at HAB= 25 mm, the iso-dodecane doping does not affect the d 6,3 profile in either stream. The experimental results partially validate the chemical reactions and soot formation kinetic model developed at Lawrence Livermore National Laboratory and provide directions to further improve its predictions.

Iso-dodecane (2,2,4,6,6-Pentamethyl-Heptane)↗
NASA NTRS1987

Analysis of non-premixed turbulent reacting flows

Studies of chemical reactions occurring in turbulent flows are important in the understanding of combustion and other applications. Current numerical methods are limited in their applications due to the numerical resolution required to completely capture all length scales, but, despite the fact that realistic combustion cannot be solved completely, numerical simulations can be used to give insight into the interaction between the processes of turbulence and chemical reaction. The objective was to investigate the effects of turbulent motion on the effects of chemical reaction to gain some insight on the interaction of turbulence, molecular diffusion, and chemical reaction to support modeling efforts. A direct turbulence simulation spectral code was modified to include the effects of chemical reaction and applied to an initial value problem of chemical reaction between non-premixed species. The influence of hydrodynamics on the instantaneous structure of the reaction was investigated. The local scalar dissipation rates and the local reaction rates were examined to determine the influence of vorticity or rate of strain on the reaction and the structure of the scalar field.

Leonard, Andy D.↗
NASA NTRS1992

Effect of finite-rate chemistry and unequal Schmidt numbers on turbulent non-premixed flames modeled with single-step chemistry

The interaction between a quasi-laminar flame and a turbulent flowfield is investigated through direct numerical simulations (DNS) of reacting flow in two- and three-dimensional domains. Effects due to finite-rate chemistry are studied using a single step global reaction A (fuel) + B (oxidizer) yields P (product), and by varying a global Damkoehler number, as a result of which the turbulence-chemistry interaction in the flame is found to generate a wide variety of conditions, ranging from near-equilibrium to near-extinction. Differential diffusion effects are studied by changing the Schmidt number of one reactive species to one-half. It is observed that laminar flamelet response is followed within the turbulent flowfield, except in regions where transient effects seem to dominate.

Chen, J. H.↗
NASA NTRS1992

Structure of turbulent non-premixed flames modeled with two-step chemistry

Direct numerical simulations of turbulent diffusion flames modeled with finite-rate, two-step chemistry, A + B yields I, A + I yields P, were carried out. A detailed analysis of the turbulent flame structure reveals the complex nature of the penetration of various reactive species across two reaction zones in mixture fraction space. Due to this two zone structure, these flames were found to be robust, resisting extinction over the parameter ranges investigated. As in single-step computations, mixture fraction dissipation rate and the mixture fraction were found to be statistically correlated. Simulations involving unequal molecular diffusivities suggest that the small scale mixing process and, hence, the turbulent flame structure is sensitive to the Schmidt number.

Chen, J. H.↗
NASA NTRS1987

On the modelling of non-reactive and reactive turbulent combustor flows

A study of non-reactive and reactive axisymmetric combustor flows with and without swirl is presented. Closure of the Reynolds equations is achieved by three models: kappa-epsilon, algebraic stress and Reynolds stress closure. Performance of two locally nonequilibrium and one equilibrium algebraic stress models is analyzed assuming four pressure strain models. A comparison is also made of the performance of a high and a low Reynolds number model for combustor flow calculations using Reynolds stress closures. Effects of diffusion and pressure-strain models on these closures are also investigated. Two models for the scalar transport are presented. One employs the second-moment closure which solves the transport equations for the scalar fluxes, while the other solves the algebraic equations for the scalar fluxes. In addition, two cases of non-premixed and one case of premixed combustion are considered. Fast- and finite-rate chemistry models are applied to non-premixed combustion. Both show promise for application in gas turbine combustors. However, finite rate chemistry models need to be examined to establish a suitable coupling of the heat release effects on turbulence field and rate constants.

Nikjooy, Mohammad↗
DOE OSTI2025

Progress in modeling hydrogen assisted ammonia oxidation with new experiments and a further reconciliation of the NH 3 + OH rate constant

Hydrogen-assisted oxidation of ammonia in a premixed, laminar flow tubular reactor under reducing conditions was investigated experimentally and through chemical kinetic modeling. Due to its impact on the competition for OH among ammonia and hydrogen, the rate constant for NH 3 + OH (R1) was determined through state-of-the-art theoretical kinetics calculations, employing composite energies that include the effects of higher order electronic excitations on the electronic energies along the variational reaction path and treating the limitations in the kinetics posed by the passage through a hydrogen-bonded complex. The resulting rate constant was in close agreement with the recent experimental value from Zaczek et al. (2025), settling a long-term dispute about the high-temperature value of k 1 and confirming within 20% the value previously used in modeling. The chemical kinetic model, with no other changes, captured well measured concentrations of NH 3 , H 2 , NO, and N 2 O from flow reactor oxidation of NH 3 /H 2 at slightly reducing conditions over a range of temperature (900-1350 K) and NH 3 /H 2 ratios (0.5-2.0). Comparison of the present results with reported data from a non-premixed setup indicates that for laminar flow tubular reactors, the reactor configuration may have implications for the observed H 2 consumption due to the possibility of preferential oxidation during mixing.

Ab initio theory↗
DOE OSTI2024

Direct numerical simulations of three-component Rayleigh–Taylor mixing and an improved model for multicomponent reacting mixtures

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

42 ENGINEERING↗
DOE OSTI2025

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗
DOE OSTI2025

CRN Modeling of Ammonia RQL Combustion using a Partially-Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However low flammability and a propensity for high nitrogen oxide (NOx) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NOx emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel, while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NOx concentrations below that of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing non-adiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and post-flame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors’ past work, this arrangement has demonstrated the significant sensitivity of exit NOx to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and non-premixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NOx emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have shown to create conditions which are conducive to increased N2O formation in the lean stage. The results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia..

advanced gas turbines↗
NASA NTRS1981

Reynolds number effects in combustion noise

Acoustic emission spectra have been obtained for non-premixed turbulent combustion from two small diameter laboratory gas burners, two commercial gas burners and a large gas burner in the firebox of a Babcock-Wilcox Boiler (50,000 lb steam/hr). The changes in burner size and firing rate represent changes in Reynolds number and changes in air/fuel ratio represent departure from stoichiometric proportions. The combustion efficiency was measured independently through gas analysis. The acoustic spectra obtained from the various burners exhibit a persistent shape over the Reynolds number range of 8200-82,000. The spectra were analyzed for identification of a predictable frequency domain that is most responsive to, and readily correlated with, combustion efficiency. A simple parameter (consisting of the ratio of the average acoustic power output in the most responsive frequency bandwidth to the acoustic power level of the loudest frequency) is proposed whose value increases significantly and unmistakably as combustion efficiency approaches 100%. The dependence of the most responsive frequency domain on the various Reynolds numbers associated with turbulent jets is discussed.

Seshan, P. K.↗
NASA NTRS1988

Analysis of homogeneous turbulent reacting flows

Full turbulence simulations at low Reynolds numbers were made for the single-step, irreversible, bimolecular reaction between non-premixed reactants in isochoric, decaying homogeneous turbulence. Various initial conditions for the scalar field were used in the simulations to control the initial scalar dissipation length scale, and simulations were also made for temperature-dependent reaction rates and for non-stoichiometric and unequal diffusivity conditions. Joint probability density functions (pdf's), conditional pdf's, and various statistical quantities appearing in the moment equations were computed. Preliminary analysis of the results indicates that compressive strain-rate correlates better than other dynamical quantities with local reaction rate, and the locations of peak reaction rates seem to be insensitive to the scalar field initial conditions.

Leonard, A. D.↗
NASA NTRS1990

Numerical analysis of bipropellant combustion in orbit maneuvering vehicle thrust chamber

The combustion and aerodynamic processes of liquid-liquid bipropellant; monomethylhydrazine and nitrogen tetroxide in the combustion chamber of a variable thrust engine (VTE) have been examined to assess the flow structure, combustion characteristics, and wall heat transfer under the simulated engine operations at two power levels. The Generalized Elliptic Multi-phase flow with Chemical reaction, Heat transfer and Interfacial Processes (GEM-CHIP) code and the Bipropellant Combustion (BICOMB) code, which is the GEMCHIP algorithm reformulated in the frame work of a body fitted coordinate system, have been used for the numerical study. Combustion flow structure with non-premixed type injection that simulates the impingement type injector is predicted and the results are compared with that of the premixed injector.

Chiu, H. H.↗
NASA NTRS1990

Full numerical simulation of coflowing, axisymmetric jet diffusion flames

The near field of a non-premixed flame in a low speed, coflowing axisymmetric jet is investigated numerically using full simulation. The time-dependent governing equations are solved by a second-order, explicit finite difference scheme and a single-step, finite rate model is used to represent the chemistry. Steady laminar flame results show the correct dependence of flame height on Peclet number and reaction zone thickness on Damkoehler number. Forced simulations reveal a large difference in the instantaneous structure of scalar dissipation fields between nonbuoyant and buoyant cases. In the former, the scalar dissipation marks intense reaction zones, supporting the flamelet concept; however, results suggest that flamelet modeling assumptions need to be reexamined. In the latter, this correspondence breaks down, suggesting that modifications to the flamelet modeling approach are needed in buoyant turbulent diffusion flames.

Mahalingam, S.↗
DOE OSTI2025

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗