Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing front”

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 163 records · Page 9

Modeling the microplastic distribution along the Delaware River Estuary: Accumulation patterns and hydrodynamic influences

Microplastic pollution is an escalating environmental concern, particularly in densely populated estuary regions, where it poses significant threats to aquatic life and human health. The dispersion patterns of microplastic particles along estuaries are influenced and complicated by multiple environmental factors such as river flow, tidal mixing, salt intrusion, and estuarine circulation. This study examines the accumulation and dispersion patterns by modeling three typical classes of microplastics in the Delaware River Estuary: synthetic fibers, sinking plastic films, and rising plastic pellets. Our findings reveal specific areas with high microplastic accumulation for each type. Notably, the upper estuary regions exhibit significant retention of rising microplastics, associated with a region with reduced along-thalweg velocities downstream of Trenton, NJ and upstream of Philadelphia, PA. Conversely, synthetic fibers and sinking plastic films accumulate in the flow convergence zone near the bottom salinity front, typically downstream of Philadelphia. All of the microplastic accumulation hot spot locations are controlled by the balance of river discharge and salinity intrusions. During the dry season, microplastic accumulation hot spots shift upstream in the estuary, whereas in the wet season, the strong river discharge pushes them downstream. Furthermore, on the other hand, tidal mixing, settling, and resuspension processes strongly impact the spreading of microplastics along the river.

Delaware River Estuary↗

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗

Water Vapor Column Measurements with Infrared Active Optical IPDA Lidar

Development of a novel triple-pulsed 2-μm direct detection Integrated Path Differential Absorption (IPDA) lidar to measure column average carbon dioxide (XCO2) and water vapor (XH2O) from an airborne platform provides a new tool to advance atmospheric science studies. Column water vapor measurements are used in global mapping of weather fronts, weather forecast, and in studies of atmospheric transport, radiation, clouds and climate. Water vapor interference is source of biases in XCO2 measurements due to overlap in absorption, line broadening, and retrieval of dry air mixing of carbon dioxide. Model calculations show a XH2O measurement capability of 0.5% accuracy over 1.0 km spatial scales from the NASA B-200 aircraft. Ground testing of the integrated IPDA was in progress during December 2017-January 2018, and airborne measurements and validation of XH2O measurements over land and ocean are planned during January-February 2018. Technology developments in the measurement of XH2O are presented along with performance projections, and results from ground and airborne testing and validations in this paper.

Singh, Upendra N.↗

Development and Application of Benchmark Examples for Mixed-Mode I/II Quasi-Static Delamination Propagation Predictions

The development of benchmark examples for quasi-static delamination propagation prediction is presented. The example is based on a finite element model of the Mixed-Mode Bending (MMB) specimen for 50% mode II. The benchmarking is demonstrated for Abaqus/Standard, however, the example is independent of the analysis software used and allows the assessment of the automated delamination propagation prediction capability in commercial finite element codes based on the virtual crack closure technique (VCCT). First, a quasi-static benchmark example was created for the specimen. Second, starting from an initially straight front, the delamination was allowed to propagate under quasi-static loading. Third, the load-displacement as well as delamination length versus applied load/displacement relationships from a propagation analysis and the benchmark results were compared, and good agreement could be achieved by selecting the appropriate input parameters. The benchmarking procedure proved valuable by highlighting the issues associated with choosing the input parameters of the particular implementation. Overall, the results are encouraging, but further assessment for mixed-mode delamination fatigue onset and growth is required.

Krueger, Ronald↗

Ionization wave propagation in a He plasma jet in a controlled gas environment

Characterizing ionization wave propagation in low temperature plasma jets is critical to predicting production of reactive species and plasma–surface interactions for biomedical applications and surface functionalization. In this work, results from optical emission and laser induced fluorescence measurements of the ionization wave in a He plasma jet operating in a controlled gas environment are discussed and used for comparison with numerical modeling. Additionally, the ionization wave was observed using ICCD (Intensified Charge Coupled Device) imaging and characterized by time and spatially resolved electron density measurements using laser-collision-induced fluorescence. The plasma jet was initially characterized using pure He (nominally at 200 Torr), while varying pressure and voltage. When operating in pure He, the ionization wave broadly expands exiting the plasma tube. Increasing the operating pressure reduces the speed and isotropic expansion of the ionization wave. Furthermore, the jet operated with a humid He shroud was also studied. The humid He shroud results in the electron density increasing and having an annular profile due to the lower ionization potential of H2O compared to He and localized photoionization in the mixing region. Numerical modeling highlighted the importance of resonance radiation emitted by excited states of He, photoelectron emission from the quartz tube, and the kinetic behavior of the electrons produced by photoionization ahead of the ionization front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-dimensional analysis of a postbuckled embedded delamination

Delamination growth caused by local buckling of a delaminated group of plies was investigated. Delamination growth was assumed to be governed by the strain energy release rates, G(1), G(2) and G(3). The strain energy release rates were calculated using a geometrically nonlinear three-dimensional finite element analysis. The program is described and several checks of the analysis are discussed. Based on a limited parametric study, the following conclusions were reached: (1) the problem is definitely mixed mode (in some cases G(1) is larger than G(2), for other cases the opposite is true); (2) in general, there is a large gradient in the strain energy release rates along the delamination front; (3) the locations of maximum G(1) and G(2) depend on the delamination shape and the applied strain; (4) the mode 3 component was negligible for all cases considered; and (5) the analysis predicted that parts of the delamination would overlap. The results presented did not impose contact constraints to prevent overlapping. Further work is needed to determine the effects of allowing the overlapping.

Whitcomb, John D.↗

Three-dimensional analysis of a postbuckled embedded delamination

Delamination growth caused by local buckling of a delaminated group of plies was investigated. Delamination growth was assumed to be governed by the strain energy release rates, G(1), G(2) and G(3). The strain energy release rates were calculated using a geometrically nonlinear three-dimensional finite element analysis. The program is described and several checks of the analysis are discussed. Based on a limited parametric study, the following conclusions were reached: (1) the problem is definitely mixed mode (in some cases G(1) is larger than G(2), for other cases the opposite is true); (2) in general, there is a large gradient in the strain energy release rates along the delamination front; (3) the locations of maximum G(1) and G(2) depend on the delamination shape and the applied strain; (4) the mode 3 component was negligible for all cases considered; and (5) the analysis predicted that parts of the delamination would overlap. The results presented did not impose contact constraints to prevent overlapping. Further work is needed to determine the effects of allowing the overlapping.

Whitcomb, John D.↗

Transport Pathways for Asian Pollution Outflow Over the Pacific: Interannual and Seasonal Variations

The meteorological pathways contributing to Asian pollution outflow over the Pacific are examined with a global three-dimensional model analysis of CO observations from the Transport and Chemical Evolution over the Pacific (TRACE-P) aircraft mission (February-April 2001). The model is used also to place the TRACE-P observations in an interannual (1994-2001) and seasonal context. The major process driving Asian pollution outflow in spring is frontal lifting ahead of southeastward-moving cold fronts (the leading edge of cold surges) and transport in the boundary layer behind the cold fronts. Orographic lifting over central and eastern China combines with the cold fronts to promote the transport of Chinese pollution to the free troposphere. Outflow of seasonal biomass burning in Southeast Asia during spring takes place mostly by deep convection but also by northeastward transport and frontal lifting, mixing with the anthropogenic outflow. Boundary layer outflow over the western Pacific is largely devoid of biomass burning influence. European and African (biomass burning) plumes in Asian outflow during TRACE-P were weak (less than 60 ppbv and 20 ppbv CO, respectively) and were not detectable in the observations because of superposition of the much larger Asian pollution signal. Spring 2001 (La Nina) was characterized by unusually frequent cold surge events in the Asian Pacific rim and strong convection in Southeast Asia, leading to unusually strong boundary layer outflow of anthropogenic emissions and convective outflow of biomass burning emissions in the upper troposphere. The Asian outflow flux of CO to the Pacific is found to vary seasonally by a factor of 3-4 (maximum in March and minimum in summer). The March maximum results from frequent cold surge events and seasonal biomass burning emissions.

Liu, Hong-Yu↗

Observed Variability in Convective Cell Characteristics and Near-Storm Environments across the Sea- and Bay-Breeze Fronts in Southeast Texas

Abstract During the DOE Atmospheric Radiation Measurement (ARM) Tracking Aerosol Convection Interactions Experiment (TRACER) IOP spanning June–September 2022, two fixed ARM sites and a mobile team concurrently sampled the airmass heterogeneity across sea- and bay-breeze fronts around the greater Houston metropolitan region. Here, we quantify the spatiotemporal variability between maritime (coastal/bay side of breeze fronts) and continental (inland side of breeze fronts) air masses over 15 IOP days characterized by strong sea-breeze forcing. We analyze environmental profile data from 177 radiosondes and use S- and C-band radar data to track and quantify the variability in attributes of more than 2300 shallow and transitioning cells across different air masses. The composite analysis of environmental profiles indicates that during the early afternoon, the sea-breeze maritime air mass exhibits lower convective available potential energy (CAPE) than the bay-breeze maritime air mass. As the sea breeze advances inland with time, CAPE within the maritime air mass exceeds that of the continental air mass to the north of the breeze fronts. In general, maritime cells have a larger mean composite reflectivity and cell widths than continental cells; however, the response varies between shallow and transitioning cells. Mean composite 20-dB Z echo-top heights, however, are similar across air masses for both shallow and transitioning cells. The continental and maritime inflow air mass for transitioning cells has significantly different mean values for mixed-layer entrainment CAPE, lifted condensation level, level of free condensation, boundary layer depth, and diluted equilibrium level. For shallow cells, only total precipitable water shows a significant difference. Significance Statement The greater Houston metropolitan area is a natural laboratory for understanding the individual impacts of background meteorology and aerosols on convective clouds. Due to its proximity to the Gulf Coast and Galveston Bay, the Houston region experiences a diurnal precipitation cycle in the summer, driven by convection triggered from sea- and bay-breeze fronts. These fronts act as a boundary between air masses with distinct thermodynamic and environmental characteristics. Convergence along these fronts and interactions between storm outflow and the fronts facilitate convection initiation in different mesoscale air masses. This study quantifies the heterogeneity among these air masses while investigating their influence on cloud microphysics. We find that the effect of airmass heterogeneity is more pronounced for the bulk microphysical properties in shallow clouds.

Sharma, Milind↗

Development and Application of Benchmark Examples for Mixed-Mode I/II Quasi-Static Delamination Propagation Predictions

The development of benchmark examples for quasi-static delamination propagation prediction is presented and demonstrated for a commercial code. The examples are based on finite element models of the Mixed-Mode Bending (MMB) specimen. The examples are independent of the analysis software used and allow the assessment of the automated delamination propagation prediction capability in commercial finite element codes based on the virtual crack closure technique (VCCT). First, quasi-static benchmark examples were created for the specimen. Second, starting from an initially straight front, the delamination was allowed to propagate under quasi-static loading. Third, the load-displacement relationship from a propagation analysis and the benchmark results were compared, and good agreement could be achieved by selecting the appropriate input parameters. Good agreement between the results obtained from the automated propagation analysis and the benchmark results could be achieved by selecting input parameters that had previously been determined during analyses of mode I Double Cantilever Beam and mode II End Notched Flexure specimens. The benchmarking procedure proved valuable by highlighting the issues associated with choosing the input parameters of the particular implementation. Overall the results are encouraging, but further assessment for mixed-mode delamination fatigue onset and growth is required.

Krueger, Ronald↗

FEED Study of CarbonCapture Inc DAC and CarbonCure Utilization Technologies Using United States Steel’s Gary Works Plant Waste Heat (Final Report)

The University of Illinois at Urbana-Champaign (UIUC) led this project to produce a front-end engineering design (FEED) study of an advanced Direct Air Capture and Utilization System (DACUS) system that can remove a minimum of 5,000 tonnes/yr net of carbon dioxide from air (based on cradle-to-gate LCA) and utilizing the CO 2 to produce low carbon intensity ready mix concrete. The designed system, if built, would be larger than any currently existing Direct Air Capture (DAC) collector in the U.S. Such carbon capture technologies are critical to meeting the goals of the DOE’s program to accelerate climate-critical technology. In addition to the power sector, industrial facilities for the manufacture of steel and cement/concrete are among the major sources of anthropogenic CO 2 . DAC is a promising new technology for reducing CO 2 , a potent greenhouse gas, in the atmosphere but is expensive, in part due to the energy required to adsorb and desorb captured CO 2 during cycles. By integrating CarbonCapture Inc. (CCI) DAC modules at United States Steel's Gary Works (USS) and utilizing the site's waste heat, energy, and location this project evaluates the feasibility of utilizing the captured CO 2 and the logistics of transportation. CarbonCapture Inc. has developed an innovative DAC system using novel adsorbents to cost-effectively capture CO 2 . The captured, liquified gas will be trucked to ready-mix concrete plants within the region, the closest of which is approximately 3.5 miles away, where CarbonCure will inject it into concrete during the mixing process at the facilities. The carbon dioxide reacts with concrete, mineralizing into calcium carbonate (CaCO 3 ), permanently locking the greenhouse gas into the matrix of the building material. This FEED study demonstrated a full CO 2 value chain for DACUS from industrial facilities. It also provided a means for Visage Energy Corp. (Visage) to assess the impact of this holistic approach on job creation, regional economic impact, and environmental justice. The project team also included Sargent & Lundy (S&L) to provide the constructability review and costing of the integration of the DAC with the steel plant. Ecotek Engineering USA, LLC designed the outside battery limit (OSBL) infrastructure to connect the DAC and the plant. Activities performed during the project included: (1) Project Management Plan; (2) Technology Maturation Plan (TMP); (3) Initial Workforce Readiness Plan; (4) Workforce Readiness Plan; (5) Front-End Engineering Design (FEED) Study; (6) Project Design Basis; (7) Hazards and Operability (HAZOP) Study; (8) Constructability Review; (9) Project Cost Assessment; (10) Logistics Analysis of CO 2 Transportation to the Utilization Site; (11) Business Case Analysis; (12) Life Cycle Analysis (LCA); (13) Environmental Health and Safety (EH&S) Analysis; (14) Environmental Justice Analysis; and (15) Economic Revitalization and Job Creation Outcomes Analysis. This report provides a summary of the outcomes and results of the project, which was performed between Oct. 1, 2022, through Sept. 30, 2024.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Plates and shells containing a surface crack under general loading conditions

Various through and part-through crack problems in plates and shells are considered. The line-spring model of Rice and Levy is generalized to the skew-symmetric case to solve surface crack problems involving mixed-mode, coplanar crack growth. Compliance functions are introduced which are valid for crack depth to thickness ratios at least up to .95. This includes expressions for tension and bending as well as expressions for in-plane shear, out-of-plane shear, and twisting. Transverse shear deformation is taken into account in the plate and shell theories and this effect is shown to be important in comparing stress intensity factors obtained from the plate theory with three-dimensional solutions. Stress intensity factors for cylinders obtained by the line-spring model also compare well with three-dimensional solution. By using the line-spring approach, stress intensity factors can be obtained for the through crack and for part-through crack of any crack front shape, without recalculation integrals that take up the bulk of the computer time. Therefore, parameter studies involving crack length, crack depth, shell type, and shell curvature are made in some detail. The results will be useful in brittle fracture and in fatigue crack propagation studies. All problems considered are of the mixed boundary value type and are reducted to strongly singular integral equations which make use of the finite-part integrals of Hadamard. The equations are solved numerically in a manner that is very efficient.

Joseph, Paul F.↗

Seamless transitions from early prototypes to mature operational software - A technology that enables the process for planning and scheduling applications

New and innovative software technology is presented that provides a cost effective bridge for smoothly transitioning prototype software, in the field of planning and scheduling, into an operational environment. Specifically, this technology mixes the flexibility and human design efficiency of dynamic data typing with the rigor and run-time efficiencies of static data typing. This new technology provides a very valuable tool for conducting the extensive, up-front system prototyping that leads to specifying the correct system and producing a reliable, efficient version that will be operationally effective and will be accepted by the intended users.

Hornstein, Rhoda S.↗

Modeling Autonomous Vehicle-Targeted Aggressive Merging Behaviors in Mixed Traffic Environment

Promising advances in Autonomous Vehicle (AV) technology have fueled industry and research fields to dedicate significant effort to the study of the integration of AVs into the traffic network. This study focuses on the transition phase between all Human Driven Vehicles (HDVs) in the network to all AVs, where these different vehicle types coexist in a mixed traffic environment. This paper investigates the potential impacts of aggressive merging behaviors by human drivers on traffic performance in a mixed environment. For this, three vehicle types – AVs, HDVs, and Aggressive HDVs (AHDVs) – are modeled in an open-source microscopic traffic simulation model, SUMO. In the developed simulation, the AHDVs are modeled to emulate aggressive merging behaviors in front of AVs at a merge section of a freeway exit ramp. Several experiments are used to study the impact of such behavior. Results show travel time gains by AHDVs at the expense of AVs and HDVs.

Bae, JongIn↗

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗

Three-Dimensional Distribution of a Major Desert Dust Outbreak over East Asia in March 2008 Derived from IASI Satellite Observations

We describe the daily evolution of the three-dimensional (3D) structure of a major dust outbreak initiated by an extratropical cyclone over East Asia in early March 2008, using new aerosol retrievals derived from satellite observations of IASI (Infrared Atmospheric Sounding Interferometer). A novel auto-adaptive Tikhonov-Phillips-type approach called AEROIASI is used to retrieve vertical profiles of dust extinction coefficient at 10 microns for most cloud-free IASI pixels, both over land and ocean. The dust vertical distribution derived from AEROIASI is shown to agree remarkably well with along-track transects of Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) spaceborne lidar vertical profiles (mean biases less than 110 meters, correlation of 0.95, and precision of 260 meters for mean altitudes of the dust layers). AEROIASI allows the daily characterization of the 3D transport pathways across East Asia of two dust plumes originating from the Gobi and North Chinese deserts. From AEROIASI retrievals, we provide evidence that (i) both dust plumes are transported over the Beijing region and the Yellow Sea as elevated layers above a shallow boundary layer, (ii) as they progress eastward, the dust layers are lifted up by the ascending motions near the core of the extratropical cyclone, and (iii) when being transported over the warm waters of the Japan Sea, turbulent mixing in the deep marine boundary layer leads to high dust concentrations down to the surface. AEROIASI observations and model simulations also show that the progression of the dust plumes across East Asia is tightly related to the advancing cold front of the extratropical cyclone.

Cuesta, Juan↗

DNS of ignition and flame stabilization in a simplified gas turbine premixer

With the increasing need for fuel flexibility, mitigation of auto-ignition (AI) inside gas turbine (GT) premixers becomes crucial. They must be designed to yield a sufficiently homogeneous fuel-air mixture to achieve low emissions while at the same time avoiding the occurrence of AI and subsequent flame stabilization. This challenge requires a detailed understanding of turbulent mixing and chemistry interactions. In the present work, a direct numerical simulation (DNS) of an array of jets in crossflow (JICF), representative of an industrial GT premixer, is reported to shed light on these complex phenomena. It is found that AI kernels form in the aft part of the premixer and coalesce into a flame front that then propagates upstream, mainly through the boundary layer, and successively engulfs the jets. This, therefore, suggests a significant role of the jet array pattern on the flame stabilization. It is noted that AI kernels continue to form independently during the whole time of the simulation. To clarify the contribution of AI and diffusion in the ignition kernels and the main flame, chemical explosive mode analysis (CEMA) is employed jointly with a kernel tracking algorithm. It is found that during the initial formation of the flame, many ignition kernels form in mixtures with low scalar dissipation rate and large contribution from AI mode. As they quickly grow, they merge into a single flame front that becomes increasingly more diffusion-assisted over time, balancing the AI mode. Turbulence is shown to have a significant enhancing effect in lean premixed flames, but further analysis is required to fully characterize it. These findings are relevant for the industrial premixer studied, and also for novel micromix concepts that may be used in the next generation of GT combustion systems.

ADVANCED PROPULSION SYSTEMS↗

Evaluation of Minimum NOx Emission From Ammonia Combustion

Abstract Ammonia (NH3) is being explored as a hydrogen carrier with no carbon emissions. However, if burned directly as NH3, rather than being completely decomposed back to N2/H2, the fuel-bound nitrogen comes with a potentially significant NOx emissions penalty. Indeed, several existing studies are showing ammonia combustion NOx emissions that exceed current natural gas fueled, DLN technologies by one to two orders of magnitude. Therefore, it is important to establish the theoretical minimum NOx emissions for an ammonia combustor, to determine how much NOx levels can be reduced via further technology development. In other words, the purpose of this work is not to analyze the performance of a specific combustor but, rather, the fundamental limits of what is achievable. This study quantifies this minimum NOx level for a two-stage combustor system for a given combustor exit temperature and residence time, with a constraint on unburned fuel levels. As expected, the optimum configuration is a rich front end combustor to burn and crack ammonia with significant H2 production, followed by an NO relaxation reactor, followed by a lean stage that consumes the remaining H2. The optimum residence time and stoichiometry of each zone are determined in the fast mixing limit, which essentially balances between NOx production in the primary and secondary zones. These results show minimum NOx levels are in 200–400 ppm range at 1 bar, but drop to levels of ∼25 ppm at 20 bar. These NOx emissions are dominated by NOx production in the primary stage which relaxes to equilibrium levels quite slowly. As processes controlling NOx relaxation to equilibrium in the primary stage dominate overall NO emission levels, combustor NOx sensitivities are essentially opposite that of natural gas fired, DLN systems. Specifically, NOx values drop with increased combustor residence time, increased pressure, and increased combustor exit temperature. These results also suggest that the most important strategy for NOx minimization is to provide sufficient relaxation time after the primary zone for NOx to approach equilibrium—this can be done via kinetic means to accelerate this relaxation rate, such as enhancing pressure or temperature, or increasing residence times. Indeed, this work shows that low pressure combustors specifically optimized for ammonia will have residence times that are one to two orders of magnitude larger than current natural gas systems. By doing so, NOx levels below 10 ppm may be achievable. Finally, we discuss the sensitivity of these values to uncertainties in ammonia kinetics.

Engineering↗