Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approach temperature”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Rapid attribution analysis of the extraordinary heat wave on the Pacific coast of the US and Canada in June 2021

Towards the end of June 2021, temperature records were broken by several degrees Celsius in several cities in the Pacific Northwest areas of the US and Canada, leading to spikes in sudden deaths and sharp increases in emergency calls and hospital visits for heat-related illnesses. Here we present a multi-model, multi-method attribution analysis to investigate the extent to which human-induced climate change has influenced the probability and intensity of extreme heat waves in this region. Based on observations, modelling and a classical statistical approach, the occurrence of a heat wave defined as the maximum daily temperature (TXx) observed in the area 45–52 °N, 119–123 °W, was found to be virtually impossible without human-caused climate change. The observed temperatures were so extreme that they lay far outside the range of historical temperature observations. This makes it hard to state with confidence how rare the event was. Using a statistical analysis that assumes that the heat wave is part of the same distribution as previous heat waves in this region led to a first-order estimation of the event frequency of the order of once in 1000 years under current climate conditions. Using this assumption and combining the results from the analysis of climate models and weather observations, we found that such a heat wave event would be at least 150 times less common without human-induced climate change. Also, this heat wave was about 2 °C hotter than a 1-in-1000-year heat wave would have been in 1850–1900, when global mean temperatures were 1.2 °C cooler than today. Looking into the future, in a world with 2 °C of global warming (0.8 °C warmer than today), a 1000-year event would be another degree hotter. Our results provide a strong warning: our rapidly warming climate is bringing us into uncharted territory with significant consequences for health, well-being and livelihoods. Adaptation and mitigation are urgently needed to prepare societies for a very different future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Direct Numerical Simulation of hydrogen combustion at auto-ignitive conditions: Ignition, stability and turbulent reaction-front velocity

Direct Numerical Simulations (DNS) are performed to investigate the process of spontaneous ignition of hydrogen flames at laminar, turbulent, adiabatic and non-adiabatic conditions. Mixtures of hydrogen and vitiated air at temperatures representing gas-turbine reheat combustion are considered. Adiabatic spontaneous ignition processes are investigated first, providing a quantitative characterization of stable and unstable flames. Results indicate that, in hydrogen reheat combustion, compressibility effects play a key role in flame stability and that unstable ignition and combustion are consistently encountered for reactant temperatures close to the mixture’s characteristic crossover temperature. Furthermore, it is also found that the characterization of the adiabatic processes is also valid in the presence of non-adiabaticity due to wall heat-loss. Finally, a quantitative characterization of the instantaneous fuel consumption rate within the reaction front is obtained and of its ability, at auto-ignitive conditions, to advance against the approaching turbulent flow of the reactants, for a range of different turbulence intensities, temperatures and pressure levels.

42 ENGINEERING↗
DOE OSTI2022

Evaluation of thermodynamic closure models for partially reacted two-phase mixture of condensed phase explosives

One of the key fundamental issues that is crucial in the continuum modeling of reactive flow phenomena is the thermodynamically consistent description of reaction mixture properties. To define the mixture properties, thermodynamic closure rules that relate the properties of the individual reaction components to the mixture properties are required. In the context of reactive two-phase modeling approaches, various strategies to define the thermodynamic closures have been adopted such as pressure temperature (PT) equilibrium between the individual reaction components, pressure (specific) volume (PV) equilibrium, etc. The choice of closure rules determines the relative distribution of specific volume and energy across the reaction components that comprise the mixture. Therefore, depending on the choice of the closure, the mixture thermodynamic behavior can vary. The present work examines the effect of different closure approaches on the thermodynamic properties of the reaction mixture. The analysis is performed for a condensed phase HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) based plastic bonded explosive (PBX) 9501 explosive using four different thermodynamic closures, viz., PT equilibrium, PV equilibrium, volume temperature (VT) equilibrium, and pressure (P) equilibrium with reactants on an isentrope. Furthermore, the relative variations in the thermodynamic properties of the mixture are analyzed and compared under both compression and expansion loading regimes. It is shown that out of the four closure models, only PT equilibrium and P equilibrium closures lead to a thermodynamically accurate description of the mixture under both compression and expansion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Doping engineering: Next step toward room temperature performance of terahertz quantum cascade lasers

We hereby offer a comprehensive analysis of various factors that could potentially enable terahertz quantum cascade lasers (THz QCLs) to achieve room temperature performance. We thoroughly examine and integrate the latest findings from recent studies in the field. Our work goes beyond a mere analysis; it represents a nuanced and comprehensive exploration of the intricate factors influencing the performance of THz QCLs. Through a comprehensive and holistic approach, we propose novel insights that significantly contribute to advancing strategies for improving the temperature performance of THz QCLs. This all-encompassing perspective allows us not only to present a synthesis of existing knowledge but also to offer a fresh and nuanced strategy to improve the temperature performance of THz QCLs. We draw new conclusions from prior works, demonstrating that the key to enhancing THz QCL temperature performance involves not only optimizing interface quality but also strategically managing doping density, its spatial distribution, and profile. This is based on our results from different structures, such as two experimentally demonstrated devices: the spit-well resonant-phonon and the two-well injector direct-phonon schemes for THz QCLs, which allow efficient isolation of the laser levels from excited and continuum states. In these schemes, the doping profile has a setback that lessens the overlap of the doped region with the active laser states. Our work stands as a valuable resource for researchers seeking to gain a deeper understanding of the evolving landscape of THz technology. Furthermore, we present a novel strategy for future endeavors, providing an enhanced framework for continued exploration in this dynamic field. This strategy should pave the way to potentially reach higher temperatures than the latest records reached for T max of THz QCLs.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Future Climate Projections for South Florida: Improving the Accuracy of Air Temperature and Precipitation Extremes With a Hybrid Statistical Bias Correction Technique

Projecting future climate variables is essential for comprehending the potential impacts on hydroclimatic hazards like floods and droughts. Evaluating these impacts is challenging due to the coarse spatial resolution of global climate models (GCMs); therefore, bias correction is widely used. Here, we applied two statistical methods—standard empirical quantile mapping (EQM) and a hybrid approach, EQM with linear correction (EQM-LIN)—to bias correct precipitation and air temperature simulated by nine GCMs. We used historical observations from 20 weather stations across South Florida to project future climate under three shared socioeconomic pathways (SSPs). Compared to the EQM, the hybrid EQM-LIN method improved R 2 of daily quantiles by up to 30% over the historical period and improved MAE up to 70% in months that contain most extreme values. Projected extreme precipitation at the weather stations showed that, compared to the EQM-LIN, the EQM method underestimates the high quantiles by up to 26% in SSP585. The projected changes in annual maximum precipitation from historical period (1985–2014) to near future (2040–2069) and far future (2070–2100) were between 2% and 16% across the study area. Projected future precipitation suggested a slight decrease during summer but an increase in fall. This, along with rising summer temperatures, suggested that South Florida can experience rapid oscillations from warmer summers and increased flooding in fall under future climate. Additionally, our comparative analyses with globally and nationally downscaled studies showed that such coarse scale studies do not represent the climatic extremes well, particularly for high quantile precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Cryogenic-Refined MOSFET Modeling for Oscillator, Frequency Divider, and Amplifier Designs Below 4 K

Capturing device characteristic changes at cryogenic temperatures is crucial for cryo-CMOS circuit designs. In this work, we present an isothermal cryogenic-refined modeling approach for CMOS transistors that is simple, low overhead, and easy to implement while offering the required accuracy for predicting circuit performance at the designated temperatures. Guided by die-level measurement data and circuit design principles, the model introduces corrections to only five critical parameters: threshold voltage, carrier mobility, elevated low-frequency flicker noise, dominant high-frequency shot noise, and subthreshold swing (SS). These refinements are implemented around the foundry-provided SPICE model, which is typically validated only down to about 200 K. With these adjustments, the proposed cryogenic-refined model achieves less than 5% error in both large-signal metrics (I–V characteristics) and small-signal parameters (e.g., transconductance) when compared with device measurements at deep-cryogenic temperatures. The methodology is validated in two advanced technologies: TSMC 40-nm CMOS and GlobalFoundries (GF) 45-nm RF-SOI. We further demonstrate its applicability in three representative RF circuits: a 30-GHz LC oscillator, a high-speed current-mode-logic (CML) frequency divider (FD), and a subthreshold Gb/s amplifier, all showing close agreement between simulated predictions and measurements performed at 4 and 2.5 K. Finally, we believe that the proposed approach is implementation-friendly and can significantly accelerate the development of cryo-CMOS integrated circuits.

circuit modeling↗
DOE OSTI2024

Localized Plasmonic Heating for Single-Molecule DNA Rupture Measurements in Optical Tweezers

To date, studies on the thermodynamic and kinetic processes that underlie biological function and nanomachine actuation in biological- and biology-inspired molecular constructs have primarily focused on photothermal heating of ensemble systems, highlighting the need for probes that are localized within the molecular construct and capable of resolving single-molecule response. Here we present an experimental demonstration of wavelength-selective, localized heating at the single-molecule level using the surface plasmon resonance of a 15 nm gold nanoparticle (AuNP). Our approach is compatible with force-spectroscopy measurements and can be applied to studies of the single-molecule thermodynamic properties of DNA origami nanomachines as well as biomolecular complexes. We further demonstrate wavelength selectivity and establish the temperature dependence of the reaction coordinate for base-pair disruption in the shear-rupture geometry, demonstrating the utility and flexibility of this approach for both fundamental studies of local (nanometer-scale) temperature gradients and rapid and multiplexed nanomachine actuation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Rapid and efficient room-temperature serial synchrotron crystallography using the CFEL TapeDrive

Serial crystallography at conventional synchrotron light sources (SSX) offers the possibility to routinely collect data at room temperature using micrometre-sized crystals of biological macromolecules. However, SSX data collection is not yet as routine and currently takes significantly longer than the standard rotation series cryo-crystallography. Thus, its use for high-throughput approaches, such as fragment-based drug screening, where the possibility to measure at physiological temperatures would be a great benefit, is impaired. On the way to highthroughput SSX using a conveyor belt based sample delivery system – the CFEL TapeDrive – with three different proteins of biological relevance (Klebsiella pneumoniae CTX-M-14 β-lactamase, Nectria haematococca xylanase GH11 and Aspergillus flavus urate oxidase), it is shown here that complete datasets can be collected in less than a minute and only minimal amounts of sample are required.

36 MATERIALS SCIENCE↗
DOE OSTI2020

New heat flux model for non-axisymmetric divertor infrared structures

A convective heat flux model for perturbed plasmas, based on guiding center ion drift in vacuum elds (A. Wingen, et al., Phys. Plasmas 21 (2014) 012509), has been updated. The old model only considered ion heat flux, while here also electron heat flux is included. The updated model predicts divertor heat flux distributions in non-axisymmetric (3D) plasmas with applied Resonant Magnetic Perturbation (RMP) fields, and includes electric scalar potentials. It is found that a radial electric field in the near Scrape-o Layer (SOL) can considerably shift the footprints toroidally, leading to a smearing out eff ect of the incident heat flux, while a simple model for sheath potential has little impact on footprints. Various approaches to model electron heat flux are studied. A convective electron model, based on collisionless free streaming, is found to yield the best agreement with measurements, while a conductive model requires a at temperature gradient inside lobes to yield acceptable peak heat ux values. A heuristic heat flux layer approach, based on a fixed layer width also requires a limited heat flux inside the last closed flux surface (LCFS); by selecting various locations of the LCFS, the results of the conductive or convective model can be recovered respectively. The sum of ion and electron heat fluxes, both obtained by the convective model, is compared to experimental data for multiple time slices in DIII-D. Strike point splitting is observed with peak heat fluxes and layer widths that compare well to infrared camera (IR) measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Thermodynamics and kinetics of solution combustion synthesis: Ni(NO 3 ) 2 + fuels systems

Solution combustion synthesis (SCS) utilizes exothermic self-propagating reactions to prepare nanoscale materials that can be used widely in energy, electronics, and biomedical technologies and other applications. SCS is a specific variety of a more general combustion synthesis (CS) method. Investigations of the thermodynamics, kinetics, and the mechanisms of SCS reactions, are not as well studied as the other CS processes. This work reports on a systematic study of the thermodynamics and kinetics of SCS reactions involving Ni(NO 3 ) 2 , an oxidizer, and either glycine (C 2 H 5 NO 2 ) or hexamethylenetetramine (HMT, C 6 H 12 N 4 ) as fuels. A thermodynamic modeling approach, based on the Gibbs free energy minimization principle, is applied to the simultaneous calculations of the adiabatic temperatures and compositions of the equilibrium products. Our calculations reveal the influence of fuel-to-oxidizer ratio, amount of water, and the oxygen in air on the combustion temperature under adiabatic conditions and the composition of the resulting products. We have, in turn, measured the combustion temperature and phase composition of products and compared them with the calculations. Variations of drying times for the solutions yield precursor gels with varying water contents. This approach enables the manipulation of combustion parameters and confirms the use of calculated activation energies for reactions using the Merzhanov-Khaikin method. Here, the results show that SCS reactions in fuel-lean solutions producing NiO have higher activation energy in contrast to reactions with fuel-rich solutions that form Ni. Reduction of activation energies due to the increase in the fuel-to-oxidizer ratio could be related to the observed change of the rate-limiting stages of the endothermic decomposition of the individual reactants to the exothermic decomposition of coordinate compounds formed between the reactants.

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

Correct Interpretations of ENDF-102 Definitions for Resonance Effects

My Uncle Willie circa 1600 wrote “What’s in a name; a rose by any other name would smell as sweet.” I fear in this case we have a somewhat similar problem in that we may be using the same word but are not using the same definition; specifically, the word Unresolved. The simplest physics definition as it applies to neutron resonances, is the energy point where we can no longer see/measure ALL – let me repeat that – ALL - of the individual resonances. That seems simple and clear, but the question is: how to represent resonances beyond this point in order to accurately reproduce the effects we have seen in measurements and expect/need to reproduce in our applications. We know there are more, unseen resonances, otherwise we wouldn’t say Unresolved. The ENDF approach is well defined in ENDF-102 and simple: for ENDF data the only way to represent Unresolved data is by using a theoretical model to define the distribution of resonances, including those that are too narrow to measure (i.e., are unresolved). It is important to note that in ENDF this is the one and only Unresolved model, e.g., there is no provision in ENDF to accurately define individually ALL resonances above the Resolved energy range – by ALL here I mean both those that we can measure and those that we cannot individually measure, but that theory and integral measurements tells us are present. An alternative approach, which would appear to be equally valid, would be to include the latest measured data as tabulated energy expendent data extending upwards in energy above the Resolved energy range. In this approach the evaluation would not include an ENDF style Unresolved energy range; it would only include a Resolved resonance region, followed by tabulated higher energy points, representing the resonances that could be measured beyond the Resolved range. But an important point to note: By listing these resonances above the resolved energy one admits that at least some resonances in this energy range are missing as Unresolved; i.e., they are too narrow or overlapping to measure. The purpose of this paper is to illustrate that the later approach, while done with good intentions, and appearing to be valid/adequate in plots, does not meet the need of our engineering applications. Why? As we will see below, of these two possible approaches, only the ENDF use of a model to statistically include the missing, i.e., unresolved, resonances, can meet our engineering needs to reproduce the integral effects we have measured and understand. Only with this statistical model can we predict and include in our calculated results the important effects of temperature (Doppler broadening), and energy integrals (self-shielding). Below I will first present results using two ENDF/B-VIII.1 evaluations, U235 and U238, that use the correct ENDF-102 definition of an Unresolved resonance region, using a statistical model to include the effects of resonances that theory predicts are present, but are too narrow to measure. These two evaluations reproduce the expected temperature (Doppler) and energy integral (self-shielding) effects that we expect. Next I will present results using one ENDF/B-VIII.1 evaluation, 26-Fe-56, that does not use an ENDF-102 Unresolved resonance region; instead above its Resolved energy range it lists many tabulated energy points, that look like measured data, but by definition, since they are included above the ENDF Resolved energy range there are missing Unresolved resonances, i.e., there are missing the resonances that are too narrow to resolve, i.e., are unresolved. My conclusion, and I hope yours, is that the below figures illustrate that this approach does not reproduce the temperature and energy integrals that we expect and need to accurately calculate results for our fission reactor calculations. As such this approach should not be used in ENDF formatted evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Creating high-precision reference gas standards of 85Kr for groundwater age-dating

Absolute gas counting (AGC) was applied to two gas blends of 85Kr in argon-methane (P10) counting gas to establish a high-precision specific activity (Bq/cm3) reference value for characterizing 85Kr detection efficiency for groundwater age dating measurements. The AGC or length-compensated technique has been utilized by the metrology community for decades and is an accepted method for developing radioactive gas standards. The AGC capability at Pacific Northwest National Laboratory (PNNL) uses a set of nine unequal-length proportional counters with precisely-measured internal volumes, and a gas loading system with high-precision pressure and temperature sensors. A series of AGC measurements were collected at multiple pressures to determine the inverse pressure relationship (1/P) for 85Kr and define a wall-effect correction that accounts for events decaying into the detector wall and not depositing sufficient energy in the gas to be detected. In addition to the wall-effect, two additional corrections were evaluated and are discussed in detail. Specifically, the threshold effect which accounts for events deposited below the analysis threshold and a detection efficiency as a function of detector volume effect that was observed during analysis. A robust uncertainty model was developed using the Guide to the expression of Uncertainty in Measurements (GUM) approach. The combination of carefully scrutinized correction factors, precise measurements of pressure, temperature and detector volume, and robust counting statistics resulted in the determination of high-precision specific activity values with 0.50% or less total combined uncertainty for two Kr-in-P10 reference gas standards (KP10) that will enable new groundwater age-dating measurements at PNNL.

Kr-85↗
DOE OSTI2024

Heat transfer optimization of uo 2 -mo fuel using genetic algorithms

Two genetic algorithm (GA) methods were applied to thermal finite element models to optimize the heat transfer efficacy of a UO 2 -Mo composite fuel pellet with typical pressurized water reactor fuel geometry. Mo additions to UO 2 have been shown to increase the thermal conductivity, thus reducing centerline temperatures and temperature gradients. Previous studies evaluated uniformly dispersed Mo or continuous Mo internal geometries (e.g., fins, plates, discs) that were selected using engineering intuition. The current study uses two different implementations of the same GA to optimize Mo placement and minimize the fuel temperature with the only constraint being a maximum 10% Mo volume fraction. One approach superimposed Mo line elements onto the monolithic UO 2 pellet model, and the other converted entire UO 2 volume elements to Mo. The former method generated 1D heat transfer connections between nodes, whereas the latter method allowed for the formation of 3D structures. Features of the optimal fuel design produced by the GAs included dispersed Mo near the centerline that shifted the peak fuel temperature outward by 0.6 mm, Mo chains in the high-heat-flux region in the mid-to-outer radial zone, and a large continuous structure that spanned the full radius and height of the pellet and accounted for 87.7 % of the total Mo in the pellet. Analysis of this design indicates that the optimal Mo configuration is a balance between creating continuous heat transfer pathways and optimally dispersing Mo to minimize the heat transfer distance through UO 2 . This architecture ultimately produced an effective thermal conductivity of 11.3 W/m·K under the assumed boundary conditions. This result is higher than any previous values from the literature. In conclusion, potential fabrication methods and challenges are discussed in addition to the implications on fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles

In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.

Adsorption-desorption process↗
DOE OSTI2025

Low-loss superconducting resonators fabricated from tantalum films grown at room temperature

The use of α-tantalum in superconducting circuits has enabled a considerable improvement in the coherence time of transmon qubits. The standard approach to grow α-tantalum thin films on silicon involves heating the substrate, which takes several hours per deposition and prevents the integration of this material with wafers containing temperature-sensitive components. We report a detailed experimental study of an alternative growth method of α-tantalum on silicon, which is achieved at room temperature through the use of a niobium seed layer. Despite a substantially higher density of oxygen-rich grain boundaries in the films sputtered at room temperature, resonators made from these films are found to have state-of-the-art quality factors, comparable to resonators fabricated from tantalum grown at high temperature. This finding challenges previous assumptions about correlations between material properties and microwave loss of superconducting thin films, and opens a new avenue for the integration of tantalum into fabrication flows with limited thermal budget.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantifying and simulating the weather forecast uncertainty for advanced building control

Weather forecast uncertainty is unavoidable despite technological advancements. Accurately quantifying and modelling this uncertainty is essential for developing and comparing advanced building controllers. In this study, we present a structured approach using a first-order autoregressive model (AR(1)) to model uncertainty in ambient temperature and global solar irradiation (GHI) forecasts. We analyzed weather data from four cities and employed Jensen–Shannon divergence (JSD) to evaluate the similarity between synthetic and actual forecast errors. The average JSD values for temperature are 0.027 (Berkeley), 0.021 (Leuven), 0.018 (Berlin), and 0.008 (Oslo), and for GHI, the average JSD values are 0.016 (Berkeley), 0.058 (Leuven), and 0.013 (Berlin). The low JSD values indicate a high similarity between the synthetic and real forecast error distributions. Further, our approach successfully generates synthetic weather forecasts that mirror the statistical properties of actual forecasts. The implementation of our method for uncertain forecast generation is being added to the BOPTEST framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

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↗