Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Combustion simulation”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An experimental and computational analysis of combustion heat release transformation in dual fuel combustion

Dual fuel (DF) diesel-methane combustion, which employs a high-reactivity fuel (diesel) to ignite a low-reactivity fuel (methane), is a widely studied combustion strategy for internal combustion engines, with significant potential for engine-out emissions reductions without the need for major hardware modifications. A phenomenon, which has been reported in the DF literature, but not explained fully, is the transformation of the shape of the apparent heat release rate (AHRR) curve as the start of injection (SOI) of diesel is advanced beyond a certain threshold; coincidentally, this AHRR transformation is usually accompanied by a sharp decrease in engine-out emissions of oxides of nitrogen (NOx). The goal of the present work is to establish the underlying physical reason(s) that cause the AHRR transformation. The AHRR transformation was observed on a single cylinder research engine (SCRE) at an indicated mean effective pressure (IMEP) of 5 bar at a speed of 1500 rev/min. The transformation occurred over a range of SOIs from 330 to 320 crank angle degrees (CAD). While the 330 CAD SOI exhibited a typical two-stage AHRR curve, with a clearly definable first-stage peak followed by a second-stage AHRR with little-to-no low temperature heat release (LTHR) present and high engine-out NOx, the 320 CAD SOI exhibited a single-stage, Gaussian-like AHRR curve, with noticeable LTHR and at least one order-of-magnitude lower NOx emissions. Here, leveraging analysis of experimental data and three-dimensional computational fluid dynamic simulations, the authors show that the AHRR transformation is impacted mainly by differences in local equivalence ratio distributions within the cylinder at ignition onset for different diesel SOIs.

33 ADVANCED PROPULSION SYSTEMS↗
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 OSTI2021

Piston Bowl Geometry Effects on Gasoline Compression Ignition in a Heavy-Duty Diesel Engine

A design optimization campaign was conducted to search for improved combustion profiles that enhance gasoline compression ignition in a heavy-duty diesel engine with a geometric compression ratio of 17.3. A large-scale design of experiments approach was used for the optimization, employing three-dimensional computational fluid dynamics simulations. The main parameters explored include geometric features, injector specifications, and swirl motion. Both stepped-lip and re-entrant bowls were included in order to assess their respective performance implications. A total of 256 design candidates were prepared using the software package CAESES for automated and simultaneous geometry generation and combustion recipe perturbation. The design optimization was conducted for three engine loads representing light to medium load conditions. The design candidates were evaluated for fuel efficiency, emissions, fuel-air mixing, and global combustion behavior. Simulation results showed that the optimum designs were all stepped-lip bowls, due to improvements in fuel-air mixing, as well as reduced heat loss and emissions formation. Improvements in indicated specific fuel consumption of up to 3.2% were achieved while meeting engine-out NOx emission targets of 1-1.5 g/kW center dot h. Re-entrant bowls performed worse compared to the baseline design, and significant performance variations occurred across the load points. Specifically, the re-entrant bowls were on par with the stepped-lip bowls under light load conditions, but significant deteriorations occurred under higher load conditions. As a final task, selected optimized designs were then evaluated under full-load conditions.

CFD↗
DOE OSTI2023

A physics-driven Σ-Y atomization model for heavy-duty engine simulations

The atomization of a liquid jet is a multi-scale and multi-physics problem of interest for many engineering applications. Particularly, it drives the fuel-oxidizer mixing and dictates the efficiency of combustion engines. High-fidelity multi-phase simulations remain challenging due to the excessive computational cost required to capture all the atomization scales. Thus, atomization models are necessary to represent sub-grid liquid structures. Here, in this work, a modified Σ-Y model in the context of Eulerian-Lagrangian Spray Atomization (ELSA) is used to transport the surface area density of the spray. The model's predictive performance is assessed under various operating conditions relevant to heavy-duty engines using the Engine Combustion Network (ECN) Spray C and Spray D research-grade injectors. The classic droplet collision formulation of the Σ-Y model alone does not replicate the response of the measured spray surface area to changes in injector, ambient pressure, and injection pressure, requiring individual tuning of the model's parameters. Instead, a transition between dense and dilute spray breakup mechanisms is proposed in terms of the average droplet spacing. The collision breakup mechanism represents the dilute spray, whereas the droplet size in the dense spray is driven by a competition between the integral scale of turbulence and the balance between the turbulent kinetic energy and the surface energy of the droplets. Such an approach minimizes the requirements for model tuning. Moreover, the role of the constants of a compressible standard k-ε RANS model is assessed in the injector's internal flow and external spray simulation framework, and an updated set is proposed. The results are validated against x-ray radiography and Ultra-Small Angle X-ray Scattering (USAXS) data and highlight the predictive capabilities of the proposed physics-driven Σ-Y model, which is compatible with engine simulation turnaround times.

42 ENGINEERING↗
DOE OSTI2025

Assessment of flamelet/progress variable methods for supersonic combustion

Tabulated chemistry models, including the flamelet/progress variable approach, have been successfully used for a variety of turbulent flame simulations. The progress variable describes the progress of reactions in a system and parameterizes a lookup table of thermochemical variables. This approach reduces the cost of simulations, transporting only one scalar (progress variable) instead of the many species mass fractions required for detailed chemistry. Originally developed for low Mach number flame simulations, recent works have focused on extensions of this approach to compressible flames, supersonic combustion, and detonations, with applications such as scramjet combustors and rotating detonation engines. Unlike low Mach simulations, compressible flow simulations require solving the energy transport equation, which is coupled to the equation of state. This leads to additional modeling challenges regarding the thermodynamics and its impact on the chemistry. The validity of modeling assumptions, for example the relationship between energy and temperature, also varies with the combustion regime. The present work provides a detailed assessment of the existing strategies for chemistry tabulation for compressible/supersonic combustion, including detonations. A priori analysis indicates that approximations which are reasonable for weakly compressible flames may break down for shock-induced combustion. Furthermore, the analysis identifies specific assumptions and approximations that do not hold for detonations, emphasizing that care must be taken when applying tabulated chemistry models outside their intended combustion regimes.

Detonations↗
DOE OSTI2022

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗
DOE OSTI2022

Detailed chemistry modeling of rotating detonations with dilute n -heptane sprays and preheated air

Utilization of liquid fuels is crucial to enabling commercialization of rotating detonation engines (RDE) in the near future. In this study, Eulerian-Lagrangian simulations are conducted for rotating detonative combustion with dilute n-heptane sprays and preheated air. Two-dimensional flattened configuration is used and a skeletal chemical mechanism with 44 species and 112 elementary reactions for n-heptane combustion is adopted. The flow structure, droplet distribution, and thermochemical parameters in the refill zone are first analyzed. It is shown that the mixture in the refill zone is heterogeneous, including evaporating droplets, vapor, and air. When the total temperature is below 950 K, the average equivalence ratio increases with the total temperature. When it is higher than 950 K, the average equivalence ratio is almost constant. Subsequently, the chemical explosive mode analysis is applied to identify the controlling reactions and dominant combustion modes in the fuel refill zone and reaction fronts. Results demonstrate that the initiation reaction (R104: n-C 7 H 16 + O 2 2-C 7 H 15 + HO 2 ) and low-temperature reaction (R107: RO 2 R'O 2 H) are dominant in the upstream and downstream of the refill zone, respectively. The intermediate species from low-temperature chemistry, R'O 2 H, is found to be important for the chemical explosive mode in the undetonated mixture. The influence of species diffusion and dispersed droplets is further analyzed. Results show that vapor autoignition facilitated by droplet evaporation occurs in the refill zone. Finally, the effects of the air total temperature on the detonation propagation speed and RDE propulsion performance are investigated. Further, it is found that the detonation propagation speed and specific impulse increase with air total temperature. The total pressure ratio first increases and then decreases as the air total temperature increases. Moreover, when the total temperature of the preheated air is above 1,300 K, the effects of low-temperature chemistry are negligible.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Enhancing Molecular Isotope Detection with a Chirp Tuned Optical Centrifuge

This paper gives a brief experimental overview on creating a chirp-tunable optical centrifuge with the goal of amplifying the isotope shift of a given molecule. Results for this experiment were recorded with a CO 2 plume as the probed material. Beyond amplification of the isotope shift, this technology could be modified to alter magnetic field characteristics of common molecules as well as simulating high temperature gases found in combustion environments. Our specific design uses double mirror and grating stretcher and will include a unique tuning feature for customization of the beam’s red and blue components' temporal chirp, allowing for quick reactivity to varying task requirements. The design and parameters will be discussed, as well as the methods used to model the optical centrifuge. Although this experiment is ongoing, results have been gathered previously that display the efficacy of the optical centrifuge. For the probed CO 2 sample, Raman Shifts centered around a wavenumber of 60 and 275 were recorded with durations of 1000 and 800 picoseconds, respectively.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Foreword: Message From the MTS Journal Editor

Here, it is my pleasure to present this general issue to the MTS community. This general issue features manuscripts on a variety of topics of interest to the MTS community including articles, reviews, and commentaries: technology gaps for monitoring birds and marine mammals at offshore wind facilities; adaptive tracking of the Barents Sea polar front using an autonomous underwater vehicle; advanced simulation and environmental impact assessment of combustion in maritime energy systems; performance of moored real time ocean observations during Cyclones in the Bay of Bengal; recent advances in assessing environmental effects of marine renewable energy; and the new blue economy and future for marine renewable energy.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

High Fidelity Simulations of a Methane-Oxygen Rotating Detonation Rocket Engine

The rotating detonation combustor (RDC) is an important realization of pressure gain combustion for practical applications. The RDC system is characterized by a highly unsteady flow field, with multiple reflected pressure waves following the strong detonations and the existence of partially-burnt gases entrained within the detonation annulus during the injector recovery process. While experimental efforts have provided macroscopic properties of an RDC with a given fuel-oxidizer mixture and geometry, limited optical access and reduced lifespan of flow-field diagnostic equipment constrain understanding of the mechanisms that lend to wave stability, controllability, and sustainability. To this end, high-fidelity numerical simulations of the Air Force Research Laboratory (AFRL) methane-oxygen rotating detonation rocket engine (RDRE) are performed to provide detailed insight into the detonation physics and anomalous behavior within the combustor. The results highlight the merit of high-fidelity numerical studies in analyzing the combustion processes and injector dynamics of a complex combustion system. The simulations capture the general wave behavior observed within the experimental studies. Three different operating conditions, with varying global equivalence ratio and system mass flow rate, are considered: (1) baseline configuration, (2) high equivalence ratio, and (3) high mass flow rate. The high mass flow rate case results in increased detonation wave strength whereas the high equivalence ratio results in increased deflagrative heat release. The high equivalence ratio case results in chaotic behavior with two pairs of counter-rotating detonation waves, whereas a single pair of co-rotating waves is observed in the baseline and high mass flow rate cases. Further, parasitic combustion of the reactant mixture ahead of the detonation wave diminishes its strength and reduces propagation speed to approximately 50-60\% of the Chapman-Jouguet (CJ) condition. In this work, a detailed analysis of the detonation wave behavior, injector dynamics, and the statistics of combustion parameters are presented.

Prakash, Supraj↗
DOE OSTI2021

A Transported Livengood–Wu Integral Model for Knock Prediction in Computational Fluid Dynamics Simulation

This work describes the development of a transported Livengood–Wu (L–W) integral model for computational fluid dynamics (CFD) simulation to predict autoignition and engine knock tendency. The currently employed L–W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics, and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition (HCCI) combustion for three different fuels, showing good accuracy in prediction of autoignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L–W integral model is coupled with G-equation model to indicate end-gas autoignition and knock tendency in CFD simulations of a direct-injection spark-ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely, Co-Optima alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multivariable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity, and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency—indicated by the fuel energy contained in the autoignited region—of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

CMLM (Co-Optimized Machine-Learned Manifolds) [SWR-23-41]

Co-optimized Machine-Learned Manifolds (CMLM) is a data-driven approach for developing reduced-order manifold models for high-dimensional chemically reacting systems. It involves a specially designed neural network, the training of which simultaneously optimizes linear combinations of species that define the manifold, nonlinear mapping to outputs of interest such as reaction rates, and (optionally) subfilter closure for large eddy simulation. This software package provides an implementation of the CMLM approach in Python using the PyTorch machine learning library. A few example cases are included, showing how the tool can be applied to different types of data from 0D and 1D reacting simulations performed using Cantera. The neural networks can be saved in a format that is readable by the Pele suite of combustion solvers for use in reacting computational fluid dynamics simulations. This software repository contains several python scripts to perform various tasks associated with the Co-optimized Machine Learned Manifolds (CMLM) model, which is described in Perry, Henry de Frahan, and Yellapantula, CNF, 2022 (https://doi.org/10.1016/j.combustflame.2022.112286). This includes not only the code that defines the CMLM model, but also scripts to generate suitable training data, scripts to pre-process the data, scripts to train the CMLM model, and scripts to plot the output, as well as various other helper files. The scripts depend on several commonly used python libraries for data analysis and chemical reaction computations. The trained models that result from this tool are designed to work with the an interface being implemented in the Pele suite of reacting flow solvers (https://github.com/AMReX-Combustion).

Perry, Bruce↗
DOE OSTI2023

Validating Connected, Automated, and Electric Vehicle Models and Simulation - Research Performance Progress Report

The objective of this project is to test connected and automated vehicles with both electrified and internal combustion engine powertrains to support updates and validation of modeling and simulation tools. This includes the development of the components and network architecture to execute and collect empirical data for multiple scenarios and traffic interactions. Specific program objectives include: • Translate Lab algorithms into vehicle and infrastructure controls • Conduct physical testing at realistic scale • Evaluate system performance • Improve models using empirical data • Improve control algorithms from lessons learned • Identify system and algorithm assumptions which need refinement It is important to note that the objective of this project was not to demonstrate the efficacy of the selected algorithms to improve energy efficiency but rather to validate and improve modeling and simulation tools using empirical data. While it is a desirable outcome to concurrently demonstrate improved energy efficiency through use of these algorithms, and in most cases that was the outcome, the success of this project was not predicated on the performance of the algorithm towards improving energy efficiency across all scenarios and test matrices.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Contrasting Stratospheric Smoke Mass and Lifetime From 2017 Canadian and 2019/2020 Australian Megafires: Global Simulations and Satellite Observations

Stratospheric injections of carbonaceous aerosols and combustion gases by extreme wildfires have become increasingly common. Recent “megafires,” particularly large and intense fires, delivered particulate burdens to the lower stratosphere comparable to those of moderate volcanic eruptions. The 2017 Canadian megafire generated four large Pyrocumulonimbi (pyroCbs), injecting up to ≈0.3 Tg of smoke in the lower stratosphere. Even more extreme, the 2019/2020 Australian event produced a pyroCb activity resulting in stratospheric smoke intrusions of ≈1 Tg. To understand their contrasting behavior, we present global climate simulations of the atmospheric response to these events, applying smoke burdens informed by remote observations. Model outcomes, compared to satellite data of smoke transport, reproduce reasonably well the initial plume rise, at 0.2–0.3 km/day, attaining heights of ≈20 km in Canada and above 30 km in Australia. Global dispersal of the plume occurs within about 3 weeks in both cases, consistent with observations. Smoke removal timescales, ≈5 months for the Canadian megafire, agree with remote measurements. During the Australian megafire, observations indicate stratospheric injections three times as large, and models predict comparatively longer smoke lifetimes, ≈16 months. After the latter event, atmospheric optical depths and radiative cooling achieved values close to those measured following the Pinatubo eruption. Sensitivity tests of model assumptions indicate, in accord with prior studies, that smoke burden, injection heights, and black carbon content can determine plume evolution and possible climate impacts. An empirical relation between peak heights of stratospheric plumes and lifetimes is derived that can help assess megafire impacts on the stratosphere, climate and the Earth system

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Effects of Turbulence and Temperature Fluctuations on Knock Development in an Ethanol/Air Mixture

The effects of turbulence on knock development and intensity for a thermally inhomogeneous stoichiometric ethanol/air mixture at a representative end-gas autoignition condition in internal combustion engines are investigated using direct numerical simulations with a skeletal reaction mechanism. Two- and three-dimensional simulations are performed by varying the most energetic length scale of temperature, \(l_T\) , and its relative ratio with the most energetic length scale of turbulence, \(l_T/l_e\) , together with two different levels of the turbulent velocity fluctuation, \(u'\) . It is found that \(l_T\) / \(l_e\) and the ratio of ignition delay time to eddy-turnover time, \(\tau _{ig}/\tau _t\) , are the key parameters that control the detonation development. An increase in either \(l_T\) or \(l_e\) enhances the detonation propensity by allowing a longer run-up distance for the detonation development. The characteristic length scale of the temperature field, \(l_T\) , is significantly modified by high turbulence intensity achieved by a large \(l_e\) and \(u'\) . The intense turbulence mixing effectively distributes the initial temperature field to broader scales to support the developing detonation waves, thereby increasing the likelihood of the detonation formation. On the contrary, high turbulence intensity with a short mixing time scale, achieved by a small \(l_e\) and a large \(u'\) , reduces the super-knock intensity attributed to the finer broken-up structures of detonation waves. Either \(\tau _{ig}/\tau _t\) less than unity or \(l_e = l_T\) even with a large \(u'\) is found to have no significant effect on super-knock mitigation. Finally, high turbulent intensity may induce high-pressure spikes comparable to the von Neumann spike. Overall, increased temperature and pressure by combustion heating, noticeably after the peak of heat release rate, significantly enhance the collision and interaction of multiple emerging autoignition fronts near the ending combustion process, resulting in localized high-pressure spikes.

42 ENGINEERING↗
DOE OSTI2021

Predicting large-scale pool fire dynamics using an unsteady flamelet- and large-eddy simulation-based model suite

A low-Mach, unstructured, large-eddy-simulation-based, unsteady flamelet approach with a generalized heat loss combustion methodology (including soot generation and consumption mechanisms) is deployed to support a large-scale, quiescent, 5-m JP-8 pool fire validation study. The quiescent pool fire validation study deploys solution sensitivity procedures, i.e., the effect of mesh and time step refinement on capturing key fire dynamics such as fingering and puffing, as mesh resolutions approach O(1) cm. A novel design-order, discrete-ordinate-method discretization methodology is established by use of an analytical thermal/participating media radiation solution on both low-order hexahedral and tetrahedral mesh topologies in addition to quadratic hexahedral elements. Here, the coupling between heat losses and the flamelet thermochemical state is achieved by augmenting the unsteady flamelet equation set with a heat loss source term. Soot and radiation source terms are determined using flamelet approaches for the full range of heat losses experienced in fire applications including radiative extinction. The proposed modeling and simulation paradigm are validated using pool surface radiative heat flux, maximum centerline temperature location, and puffing frequency data, all of which are predicted within 10% accuracy. Simulations demonstrate that under-resolved meshes predict an overly conservative radiative heat flux magnitude with improved comparisons as compared to a previously deployed hybrid Reynolds-averaged Navier–Stokes/eddy dissipation concept-based methodology.

42 ENGINEERING↗
DOE OSTI2025

Chemical Reactor Network Modeling of Ammonia Rich-Quench-Lean 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 (NO x ) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NO x 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 NO x concentrations below those of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing nonadiabatic 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 postflame 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 NO x 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 nonpremixed 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 NO x emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have been shown to create conditions that are conducive to increased N 2 O formation in the lean stage. In conclusion, the results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia.

Combustion↗