SEARCH · Engineering Papers

Results for “reactive burn model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Nonequilibrium effects in high-gain inertial confinement fusion

Recent experimental demonstrations of ignition and target gain in inertial confinement fusion (ICF) have stimulated interest in exploring the fundamental physics of violent deuterium-tritium (DT) burn in high-gain ICF targets. A significant DT-burn fraction is a necessary condition for high energy gain and large neutron yields (>100MJ). Using classical molecular-dynamics (MD) simulations and a hybrid fluid-kinetic model, we examine how a large fraction of low-energy 𝛼 particles can kick D and T ions out of equilibrium in high-gain ICF targets. The MD results suggest that (1) temperatures of 𝑇 𝐷 and 𝑇 𝑇 can differ by as much as ∼20% of their mean temperature and (2) the deviation of the DT energy distribution from the Maxwell-Boltzmann function can exceed ∼30%. Some of these MD observations, such as the preferential heating of D ions by low-energy 𝛼 particles and the temperature separation, can be explained by a proposed hybrid fluid-kinetic model. Furthermore, the implication of such nonequilibrium effects on the DT reactivity is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Aircraft measurements of aerosol and trace gas chemistry in the eastern North Atlantic

The Aerosol and Cloud Experiment in the Eastern North Atlantic (ACE-ENA) investigated properties of aerosols and subtropical marine boundary layer (MBL) clouds. Low subtropical marine clouds can have a large effect on Earth's radiative budget, but they are poorly represented in global climate models. In order to understand their radiative effects, it is imperative to understand the composition and sources of the MBL cloud condensation nuclei (CCN). The campaign consisted of two intensive operation periods (IOPs) (June–July 2017 and January–February 2018) during which an instrumented G-1 aircraft was deployed from Lajes Field on Terceira Island in the Azores, Portugal. The G-1 conducted research flights in the vicinity of the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) atmospheric observatory on Graciosa Island. An Aerodyne high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) and Ionicon proton-transfer-reaction mass spectrometer (PTR-MS) were deployed aboard the aircraft, characterizing chemistry of non-refractory aerosol and trace gases, respectively. The eastern North Atlantic region was found to be very clean, with an average non-refractory submicrometer aerosol mass loading of 0.6 µg m −3 in the summer and 0.1 µg m −3 in the winter, measured by the AMS. Average concentrations of the trace reactive gases methanol and acetone were 1–2 ppb; benzene, toluene and isoprene were even lower, <1 ppb. Mass fractions of sulfate, organics, ammonium and nitrate in the boundary layer were 69 %, 23 %, 7 % and 1 % and remained largely similar between seasons. The aerosol chemical composition was dominated by sulfate and highly processed organics. Particulate methanesulfonic acid (MSA), a well-known secondary biogenic marine species, was detected, with an average boundary layer concentration of 0.021 µg m −3 , along with its gas-phase precursor, dimethyl sulfide (DMS). MSA accounted for no more than 3 % of the submicron, non-refractory aerosol in the boundary layer. Examination of vertical profiles of aerosol and gas chemistry during ACE-ENA reveals an interplay of local marine emissions and long-range-transported aged aerosol. A case of transport of biomass burning emissions from North American fires has been identified using back-trajectory analysis. In the summer, the non-refractory portion of the background CCN budget was heavily influenced by aerosol associated with ocean productivity, in particular sulfate formed from DMS oxidation. Episodic transport from the continents, particularly of biomass burning aerosol, periodically increased CCN concentrations in the free troposphere. In the winter, with ocean productivity lower, CCN concentrations were overall much lower and dominated by remote transport. These results show that anthropogenic emissions perturb CCN concentrations in remote regions that are sensitive to changes in CCN number and illustrate that accurate predictions of both transport and regional aerosol formation from the oceans are critical to accurately modeling clouds in these regions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Chromium Nucleosynthesis and Silicon–Carbon Shell Mergers in Massive Stars

In this work, we analyze the production of the element Cr in galactic chemical evolution (GCE) models using the NuGrid nucleosynthesis yields set. We show that the unusually large [Cr/Fe] abundance at [Fe/H] ≈ 0 reported by previous studies using those yields and predicted by our Milky Way model originates from the merging of convective Si-burning and C-burning shells in a 20 ${M}_{\odot }$ model at metallicity Z = 0.01, about an hour before the star explodes. This merger mixes the incomplete burning material in the Si shell, including 51 V and 52 Cr, out to the edge of the carbon/oxygen (CO) core. The adopted supernova model ejects the outer 2 ${M}_{\odot }$ of the CO core, which includes a significant fraction of the Cr-rich material. When including this 20 M⊙ model at Z = 0.01 in the yields interpolation scheme of our GCE model for stars between 15 and 25 ${M}_{\odot }$, we overestimate [Cr/Fe] by an order of magnitude at [Fe/H] ≈ 0 relative to observations in the Galactic disk. This raises a number of questions regarding the occurrence of Si–C shell mergers in nature, the accuracy of different simulation approaches, and the impact of such mergers on the presupernova structure and explosion dynamics. According to the conditions in this 1D stellar model, the substantial penetration of C-shell material into the Si shell could launch a convective–reactive global oscillation if a merger does take place. In any case, GCE provides stringent constraints on the outcome of this stellar evolution phase.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2003

Correction to "Asian chemical Outflow to the Pacific in Spring: Origins, Pathways, and Budgets" by Isabelle Bey et al.

We analyze the Asian outflow of CO, ozone, and nitrogen oxides (NOx) to the Pacific in spring by using the GEOS-CHEM global three-dimensional model of tropospheric chemistry and simulating the Pacific Exploratory Mission-West (PEM-West B) aircraft mission in February-March 1994. The GEOS-CHEM model uses assimilated meteorological fields from the NASA Goddard Earth Observing System (GEOS). It reproduces relatively well the main features of tropospheric ozone, CO, and reactive nitrogen species observed in PEM-West B, including latitudinal and vertical gradients of the Asian pollution outflow over the western Pacific although simulated concentrations of CO tend to be too low (possibly because biogenic sources are underestimated). We use CO as a long-lived tracer to diagnose the processes contributing to the outflow. The highest concentrations in the outflow are in the boundary layer (0-2 km), but the strongest outflow fluxes are in the lower free troposphere (2-5 km) and reflect episodic lifting of pollution over central and eastern China ahead of eastward moving cold fronts. This frontal lifting, followed by westerly transport in the lower free troposphere, is the principal process responsible for export of both anthropogenic and biomass burning pollution from Asia. Anthropogenic emissions from Europe and biomass burning emissions from Africa make also major contributions to the Asian outflow over the western Pacific; European sources dominate in the lower troposphere north of 40 degrees N, while African sources are important in the upper troposphere at low latitudes. For the period of PEM-West B (February-March) we estimate that fossil fuel combustion and biomass burning make comparable contributions to the budgets of CO, ozone, and NO, in the Asian outflow. We find that 13% of NO, emitted in Asia is exported as NO, or PAN, a smaller fraction than for the United States because of higher aerosol concentrations that promote heterogeneous conversion of NOx to HNO3. Production and export of ozone from Asia in spring is much greater than from the United States because of the higher photochemical activity.

Bey, Isabelle↗

DOE OSTI2023

Modeling the Impact of the Organic Aerosol Phase State on Multiphase OH Reactive Uptake Kinetics and the Resultant Heterogeneous Oxidation Timescale of Organic Aerosol in the Amazon Rainforest

Accurate predictions of chemical lifetime, i.e., oxidation timescale and change in the mass of organic aerosols (OAs) in atmospheric models, are critical to quantify the impacts of OA on aerosol–cloud interactions, radiative forcing, and air quality. The heterogeneous oxidation of OA by hydroxyl (OH) radicals is a key process governing OA mass changes during their oxidation. Recently, laboratory observations of the OH uptake coefficient (γ OH ) and heterogeneous reaction rate (k het,OH ) for different OA systems with varying phase states have been combined to develop a new parameterization of γ OH and k het,OH as a function of the OA phase state, i.e., viscosity (η OA ). Here, in this work, we use a recently published viscosity prediction framework to analyze the new γ OH and k het,OH parameterization with a box model. Subsequently, we implement this box model to predict γ OH and k het,OH over the Amazon rainforest in the dry-to-wet transition season with both significant biomass burning and biogenic influences. Relevant parameters within our box model are specified based on detailed regional model simulations over the Amazon using the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem). η OA is predicted as a function of species volatility, OA composition (including water uptake) along with ambient relative humidity (RH), and temperature. Based on ambient conditions simulated by WRF-Chem over the entire atmospheric column, we use the box model to predict the upper bounds of the heterogeneous oxidation timescale of OA. Based on previous laboratory measurements, we assume that this heterogeneous OH oxidation causes fragmentation (carbon loss) of OA, followed by the evaporation of the fragmented molecules that results in the exponential decay of OA. We predict that the oxidation timescale of OA is ~1 month near the Earth’s surface in the pristine Amazonian background because of low OH concentrations. But OA is oxidized more rapidly within urban and wildfire plumes near the Earth’s surface with an order of magnitude higher OH concentrations compared to the pristine background, causing the simulated oxidation timescale of OA to be much shorter, ~3–4 days. At 3–5 km altitudes where biomass burning OA is predicted to be semisolid, the heterogeneous oxidation timescale is estimated as ~6 days to 3 weeks and decreases with increasing OH concentrations within plumes. We show that the simulated mass loss of OA is strongly size-dependent, where smaller particles oxidize more rapidly compared to larger particles due to their greater per-particle surface area-to-volume ratio. Increasing urbanization and deforestation in the Amazon in the future might increase OH concentrations in the background Amazon, causing faster oxidation of OA. At higher altitudes above liquid clouds, especially in the upper troposphere where temperatures approach ~250 K, future measurements are needed to reduce the uncertainties related to the mass loss of OA at colder temperatures and low relative humidity conditions due to its heterogeneous OH oxidation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

FLEX-FUEL MIXING CONTROLLED COMBUSTION ENABLED BY PRECHAMBER IGNITION

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

Zeman, Jared↗

DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through two methods, thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine the thermal quadrant powers. The N-16 system utilizes a beta chamber detector that outputs reactivity levels to calculate lobe power through an algorithm called lobe power calculation and indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. This is crucial for maintaining safety limits. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have assisted system engineers at INL with the objective to reach a more dynamic system. In return the system demonstrates an increase in the accuracy of power reading while maintaining safety margins. The new proposed coefficient and multiple method implements a cycle specific, and a physics-based model intended for changing coefficients during operation. This calibration experiment focuses on power splitting, outer shim and neck shim, as well as fuel burning into the RDAS weighting factors. Results show that the Physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. This proved true for both constrained and unconstrained testing. This is most likely due to the physics data fitting approach that resulted in favorable coefficients and multipliers for 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Towards an Agent-Based Blackboard System for Reactor Design Optimization

The field of reactor design is rich with opportunities for applications of computational optimization algorithms; these applications can range from preliminary core design to reactor shuffling patterns. Many of these schemes rely on sets of previously generated solutions (sometimes referred to as “generations”) to inform future decisions. While it is important to build upon prior knowledge, this process requires a full generation of solutions to be formed before future solutions can be examined. Rather than relying on a generational scheme to perform an optimization, we propose using an agent-based approach in conjunction with a blackboard framework for performing reactor design optimizations. Utilizing an agent-based approach allows agents to perform tasks independently, while retaining the ability to build off of previous solutions. We develop an agent-based blackboard system (ABBS) for determining the Pareto front (PF) in sodium fast reactor design optimization problems and compared this with the Non-Dominated Sorting Genetic Algorithm II (NSGA-II). Our goal is to evaluate the viability of the ABBS in producing a PF that is comparable with the NSGA-II algorithm. The design space consists of the fuel height, fuel smear, and plutonium fraction in the core, and we seek to minimize the reactivity swing and plutonium mass, while maximizing the burnup. The diversity, coverage, and spread of the PFs generated by the two methods are examined, and the ABBS is able to converge to the same PF as the NSGA-II algorithm. These results show that the ABBS is able to find optimal designs that are similar to those found by the NSGA-II algorithm. We conclude our study by applying the ABBS to the design of a sodium-cooled fast reactor to dispose of weapons-grade plutonium. The ABBS finds a core design that can burn upwards of 17.5 kg of weapons-grade plutonium per year and degrade an additional 195 kg of weapons-grade plutonium per year into non-weapons-grade material.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Investigating Secondary Aerosol Processes in the Amazon through Molecular-level Characterization of Semi-Volatile Organics

In areas where biogenic emissions are oxidized in the presence of anthropogenic pollutants, it has become increasingly apparent that secondary organic aerosol (SOA) formation from biogenic volatile organic compounds (BVOCs) is substantially enhanced. The Amazon forest is the dominant source of BVOCs globally, and the forest is rapidly being converted to urban and agricultural uses. We participated in a collaborative field study located in the Amazon region around Manaus, Brazil, in 2014 (Green Ocean Amazon; GoAmazon). This study was designed to comprehensively examine how BVOC emissions (specifically, the most highly emitted non-methane hydrocarbon, isoprene, and lesser studied sesquiterpenes) and their interaction with anthropogenic pollution (i.e., nitrogen oxides, sulfur dioxide, and black carbon) alter the atmosphere’s oxidative capacity, influence secondary aerosol formation, atmospheric composition and, ultimately, affect the earth’s radiation balance and climate. We provide the first hourly, in-situ measurement of 30 sesquiterpenes and 4 diterpenes, roughly estimating that sesquiterpene oxidation contributes to 10-14% of ozone reactive loss by terpenes and at least 0.4–5% (median 1 %) of total submicron OA mass. However, this is likely a low-end estimate, as evidence for additional unaccounted sesquiterpenes and their oxidation products clearly exists. We also provide new perspectives on sulfate, NO x , and particle acidity influencing isoprene-derived SOA, comparing the central Amazon environment with Southeastern U.S.A. summer, representing clean to polluted conditions, respectively. We find that summed concentrations of isoprene-derived SOA tracers correlated with particulate sulfate spanning three orders of magnitude, suggesting that 1 μg m -3 reduction in sulfate corresponds with at least ~0.5 μg m -3 reduction in isoprene-derived SOA. We also find that SOA mass derived from isoprene oxidation in the presence of NO x is primarily comprised of aerosol sulfate bound with isoprene oxidation products (i.e. organosulfates), ~97% in the Amazon and ~55% in the Southeastern U.S. We infer under natural conditions in high isoprene emission regions, preindustrial aerosol sulfate was almost exclusively isoprene-derived organosulfates, which are traditionally thought as representative of anthropogenic influence. We further report the first field observations showing that particle acidity impacts the distribution of isoprene oxidation tracers in the gas and particle phases, providing physicochemical insight into isoprene SOA formation chemistry. Coupled with other co-located measurements by the DOE Atmospheric Radiation Measurement (ARM) team and collaborating PIs at our measurement site (called T3), on the G1 aircraft, and at additional field sites (T0, T1, T2) this data set has been used to understand emission, oxidation, and particle formation pathways at the molecular level, to elucidate how these pathways change when influenced by anthropogenic sources, and to evaluate their importance for the atmospheric budget of aerosols and the earth’s radiation balance on regional scales. Recommendations: Further constraint of the role of sesquiterpenes on ozone (O 3 ) reactivity and SOA formation is limited by the availability of authentic standards and synthesized compounds of sesquiterpenes and their oxidation products. Future research efforts should focus on making such standards available via custom synthesis to allow for accurate quantification and focused oxidation experiments that will fully elucidate the chemistry of these compounds in the atmosphere. Further, reductions in aerosol sulfate (via SO 2 emissions control) continues to be one of the most effective methods to reduce SOA formation from isoprene, but the concerted role of ammonia (NH 3 ) gas emissions from (agricultural) industry and other sources requires more investigation to better understand the role of aerosol acidity in SOA formation and potential controls. As GoAmazon particles can be more acidic compared to areas with greater anthropogenic NH 3 emissions, this promotes relatively greater organosulfate (OS) formation from isoprene even in the presence of NO x . Because OS and inorganic sulfate differ in water uptake properties, this implies preindustrial aerosol sulfate (albeit less abundant) may have been less reflective than current earth system models assume. Further research should investigate the radiative impacts of organic vs inorganic sulfate in aerosols to improve model representation. Finally, future work (currently underway) should include constraint of the contributions of monoterpene (C 10 H 16 ) BVOCs and biomass burning VOCs as precursors to SOA formation in the region as well as their impacts on radiative forcing in the climate system.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2008

Gaseous Surrogate Hydrocarbons for a Hifire Scramjet that Mimic Opposed Jet Extinction Limits for Cracked JP Fuels

This paper describes, first, the top-down methodology used to define simple gaseous surrogate hydrocarbon (HC) fuel mixtures for a hypersonic scramjet combustion subtask of the HiFIRE program. It then presents new and updated Opposed Jet Burner (OJB) extinction-limit Flame Strength (FS) data obtained from laminar non-premixed HC vs. air counterflow diffusion flames at 1-atm, which follow from earlier investigations. FS represents a strain-induced extinction limit based on cross-section-average air jet velocity, U(sub air), that sustains combustion of a counter jet of gaseous fuel just before extinction. FS uniquely characterizes a kinetically limited fuel combustion rate. More generally, Applied Stress Rates (ASRs) at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or t) can directly be compared with extinction limits determined numerically using either a 1-D or (preferably) a 2-D Navier Stokes simulation with detailed transport and finite rate chemistry. The FS results help to characterize and define three candidate surrogate HC fuel mixtures that exhibit a common FS 70% greater than for vaporized JP-7 fuel. These include a binary fuel mixture of 64% ethylene + 36% methane, which is our primary recommendation. It is intended to mimic the critical flameholding limit of a thermally- or catalytically-cracked JP-7 like fuel in HiFIRE scramjet combustion tests. Our supporting experimental results include: (1) An idealized kinetically-limited ASR reactivity scale, which represents maximum strength non-premixed flames for several gaseous and vaporized liquid HCs; (2) FS characterizations of Colket and Spadaccini s suggested ternary surrogate, of 60% ethylene + 30% methane + 10% n-heptane, which matches the ignition delay of a typical cracked JP fuel; (3) Data showing how our recommended binary surrogate, of 64% ethylene + 36% methane, has an identical FS; (4) Data that characterize an alternate surrogate of 44% ethylene + 56% ethane with identical FS and nearly equal molecular weights; this could be useful when systematically varying the fuel composition. However, the mixture liquefies at much lower pressure, which limits on-board storage of gaseous fuel; (5) Dynamic Flame Weakening results that show how oscillations in OJB input flow (and composition) can weaken (extinguish) surrogate flames up to 200 Hz, but the weakening is 2.5x smaller compared to pure methane; and finally, (6) FS limits at 1-atm that compare with three published 1-D numerical OJB extinction results using four chemical kinetic models. The methane kinetics generally agree closely at 1-atm, whereas, the various ethylene models predict extinction limits that average ~ 45% high, which represents a significant problem for numerical simulation of surrogate-based flameholding in a scramjet cavity. Finally, we continue advocating the FS approach as more direct and fundamental for assessing idealized scramjet flameholding potentials than measurements of "unstrained" premixed laminar burning velocity or blowout in a Perfectly Stirred Reactor.

Pellett, Gerald L.↗

DOE OSTI2023

Quantifying the Contributions of Aerosol- and Snow-Produced ClNO 2 through Observations and 1D Modeling

Nitryl chloride (ClNO 2 ) is a radical reservoir that forms and accumulates in the nocturnal atmospheric boundary layer influenced by combustion emissions and chloride (e.g., sea salt and road salt). Upon sunrise, ClNO 2 rapidly photolyzes to generate highly reactive chlorine radicals (Cl • ) that affect the air quality by generating secondary air pollutants. Recent studies have shown road salt aerosols and saline snowpack to be sources of ClNO 2 in the wintertime urban environment; however, the quantitative contributions of each chloride source are not known. In this study, we examine the vertically resolved contributions of aerosol particles and saline snowpack as sources of ClNO 2 by using an observationally constrained snow–atmosphere coupled one-dimensional model applied to wintertime Kalamazoo, Michigan, U.S. Model simulations show that ClNO 2 emitted from urban snowpack can be vertically transported throughout the entire atmospheric boundary layer and can be a significant source of ClNO 2 , contributing up to ~60% of the ClNO 2 budget near the surface. Here, modeled snowpack ClNO 2 emission rates were 6 (±7) times higher than the observationally derived emission rates, suggesting that not all snow chloride is available for reaction. ClNO 2 production from both aerosol particles and snow emissions are required to best simulate the observed surface-level ClNO 2 . Using the traditional bulk parameterization for ClNO 2 produced from particles significantly overestimated ClNO 2 due to the assumption of having equivalent dinitrogen pentoxide (N 2 O 5 ) uptake and chloride availability for the entire particle population. In comparison, the chemically resolved surface area-based parameterization slightly underestimated the observations and had uncertainties deriving from ClNO 2 production from residential wood burning particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

The reactive flow evolution of the polymer-bonded explosive PBX 9502: Experiments and model validation in extreme pressure regimes

The shock-to-detonation transition properties of the triaminotrinitrobenzene based PBX 9502 high explosive (HE) are experimentally and computationally explored in extremely high input pressure conditions. These include both slightly sub-Chapman–Jouguet and overdriven input pressure conditions, namely, ∼25 and ∼31 GPa, respectively. Our experiments capture the transient buildup of a shock-induced reaction via measurement of HE and polymethyl methacrylate window interface particle velocity profiles for a variety of sample thicknesses for this insensitive HE. These observations necessitate extremely thin explosive samples, and the high rates of reaction provide a considerable challenge to optical diagnostics. Samples at these thicknesses also provide an opportunity for evaluation of potential micro-structure effects on the resulting shock-to-detonation-transition measurements. To address this, the thin samples are also characterized via x-ray micro-computed tomography. Finally, a pair of previously established continuum-level detonation performance modeling approaches for PBX 9502 were used to analyze the experiments. The employed model variants crucially differ in their definition of each model’s empirical reaction rate functional form, utilization of shock state quantities, and local flow variable dependencies. As a result, the present experiments provide a novel platform to evaluate the quantitative and qualitative consequences stemming from these modeling choices in a challenging initiation scenario, largely beyond the chosen calibration range of either model. This new experimental information will provide a platform for both improved physics and model parameterizations for this well-studied explosive.

36 MATERIALS SCIENCE↗

NASA NTRS2011

Reactive Nitrogen, Ozone and Ozone Production in the Arctic Troposphere and the Impact of Stratosphere-Troposphere Exchange

We analyze the aircraft observations obtained during the Arctic Research of the Composition of the Troposphere from Aircraft and Satellite (ARCTAS) mission together with the GEOS-5 CO simulation to examine O3 and NOy in the Arctic and sub-Arctic region and their source attribution. Using a number of marker tracers and their probability density distributions, we distinguish various air masses from the background troposphere and examine their contribution to NOx, O3, and O3 production in the Arctic troposphere. The background Arctic troposphere has mean O3 of approximately 60 ppbv and NOx of approximately 25 pptv throughout spring and summer with CO decreases from approximately 145 ppbv in spring to approximately 100 ppbv in summer. These observed CO, NOx and O3 mixing ratios are not notably different from the values measured during the 1988 ABLE-3A and the 2002 TOPSE field campaigns despite the significant changes in the past two decades in processes that could have changed the Arctic tropospheric composition. Air masses associated with stratosphere-troposphere exchange are present throughout the mid and upper troposphere during spring and summer. These air masses with mean O3 concentration of 140-160 ppbv are the most important direct sources of O3 in the Arctic troposphere. In addition, air of stratospheric origin is the only notable driver of net O3 formation in the Arctic due to its sustainable high NOx (75 pptv in spring and 110 pptv in summer) and NOy (approximately 800 pptv in spring and approximately 1100 pptv in summer) levels. The ARCTAS measurements present observational evidence suggesting significant conversion of nitrogen from HNO3 to NOx and then to PAN (a net formation of approximately 120 pptv PAN) in summer when air of stratospheric origin is mixed with tropospheric background during stratosphere-to-troposphere transport. These findings imply that an adequate representation of stratospheric O3 and NOy input are essential in accurately simulating O3 and NOx photochemistry as well as the atmospheric budget of PAN in tropospheric chemistry transport models of the Arctic. Anthropogenic and biomass burning pollution plumes observed during ARCTAS show highly elevated hydrocarbons and NOy (mostly in the form of NOx and PAN), but do not contribute significantly to O3 in the Arctic troposphere except in some of the aged biomass burning plumes sampled during spring. Convection and/or lightning influences are negligible sources of O3 in the Arctic troposphere but can have significant impacts in the upper troposphere in the continental sub-Arctic during summer.

Liang, Q.↗

DOE OSTI2020

Combined Experimental/Numerical Study of the Soot Formation Process in a Gasoline Direct-Injection Spray in the Presence of Laser-Induced Plasma Ignition

Here, combustion issued from an eight-hole, direct-injection spray was experimentally studied in a constant-volume pre-burn combustion vessel using simultaneous high-speed diffused back-illumination extinction imaging (DBIEI) and OH* chemiluminescence. DBIEI has been employed to observe the liquid-phase of the spray and to quantitatively investigate the soot formation and oxidation taking place during combustion. The fuel-air mixture was ignited with a plasma induced by a single-shot Nd:YAG laser, permitting precise control of the ignition location in space and time. OH* chemiluminescence was used to track the high-temperature ignition and flame. The study showed that increasing the delay between the end of injection and ignition drastically reduces soot formation without necessarily compromising combustion efficiency. For long delays between the end of injection and ignition (1.9 ms) soot formation was eliminated in the main downstream charge of the fuel spray. However, poorly atomized and large droplets formed at the end of injection (dribble) eventually do form soot near the injector even when none is formed in the main charge. The quantitative soot measurements for these spray and ignition scenarios, resolved in time and space, represents a significant new achievement. Reynolds-averaged Navier-Stokes (RANS) simulations were performed to assess spray mixing and combustion. An analysis of the predicted fuel-air mixture in key regions, defined based upon experimental observations, was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing. The mixture analysis indicates that soot production can be avoided if the flame propagates into regions where the equivalence ratio (Φ) is already below 2. Reactive RANS simulations have also been performed, but with a poor match against the experiment, as the flame speed and heat-release rate are largely over estimated. This modeling weakness appears related to a very high level of turbulent viscosity predicted for the high-momentum spray in the RANS simulations, which is an important consideration for modeling ignition and flame propagation in mixtures immediately created by the spray.

42 ENGINEERING↗

DOE OSTI2025

Enabling Low-Temperature (LTP) Ignition Technologies for Multi-Mode Engines through the Development of a Validated High-Fidelity LTP Model for Predicative Simulations Tools

The goal of multi-mode engine architectures is to extend current lean-burn dilution limits with renewable fuels, which requires spark plugs to deposit high energies (hundreds of mJ) in order to initiate ignition and complete combustion. At elevated energy deposition rates, spark plugs experience increased electrode erosion and thermal losses, which ultimately shortens the spark-plug lifetime and lowers ignition efficiency. As such, in order to safeguard the efficiency gains of multi-mode concepts, new and improved ignition technologies are required. Recently, non-equilibrium low-temperature plasmas (LTP) have been shown to promote energy-efficient ignition via quenching and transport of electronically excited atoms and molecules, selective radical production and fast heating of hydrocarbon/air mixtures [1-2]. Thus, LTP is seen as a technology that can potentially improve the energy extraction efficiency of fuels, while enabling kinetically controlled combustion modes towards fuel leaner conditions to realize current DOE VTO goals of improving the sustainability of future mobility [3]. Although many previous studies have demonstrated the efficacy of plasma-assisted ignition to enhance combustion, the detailed enhancement mechanisms remain largely unknown, especially for oxygenated fuels and at elevated pressures that are most relevant to practical engine conditions. These barriers hinder the development of accurate and comprehensive numerical models that seek to describe LTP-based ignition in existing engine design software tools and methods. Current state-of-the-art simulation capabilities for LTP ignition systems are in need of improvements since they deliver qualitative results only due to important limitations of existing approaches. Firstly, validated kinetic models with elementary steps for plasma discharges in oxygenated fuel/air mixtures of relevance to the transportation sector are required. Such kinetic models do not exist at present and will be developed and validated within this project. Secondly, plasma discharges and reactive mixture ignition are multi-scale, unsteady processes requiring high-performance numerical methods and software that execute efficiently on DOE supercomputers. Such software does not exist at present and will be developed and applied to practical LTP ignition scenarios as part of this project. Thirdly, experimental databases that are tailored to serve as benchmark in support of the development of predictive computational models of LTP ignition do not exist and will be part of this project.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Plasma assisted NH 3 /H 2 /air ignition in nanosecond discharges with non-equilibrium energy transfer

Ammonia (NH 3 ), with its high energy density and easiness to store and transport as a hydrogen carrier, has become a promising alternative green fuel. However, its adoption in power generation is hindered by challenges such as low burning velocity, slow low-temperature oxidation, high NO x emissions, and ignition difficulty. Here, this work computationally investigates the effects of non-equilibrium energy transfer by nanosecond discharges on NH 3 ignition and flame propagation in an NH 3 /H 2 /air flow at 700 K and 1 atm. The simulation results demonstrate that NH 3 /air mixtures require a large ignition energy due to their large critical ignition radius. It is shown that adding 30 % hydrogen significantly reduces the critical ignition radius and minimum ignition energy. Two-dimensional modeling further shows a non-monotonic dependence of ignition kernel volume on the applied voltage and reduced electric field. The optimum ignition enhancement occurs at 200 Td where the generation of electronically excited species and radicals including N 2 (B), O( 1 D) and OH becomes most efficient. Higher voltages divert electron energy toward ionization, which makes it less effective for NH 3 ignition. The study also identifies an optimal electrode gap size for a given pulse energy. Smaller gap sizes increase deposited energy density, raising temperature and radical concentrations. However, excessive reduction of the gap distance reduces flame propagation speed due to the flame stretch effect in rich mixtures with the effective Lewis number greater than unity. A nonlinear relationship between pulse repetition frequency and ignition kernel volume is observed in a nanosecond pulsed high frequency discharge (NPHFD). An optimal frequency range of 200 kHz to 2 MHz is found when two pulses are used. In addition, an optimal number of pulses exists for each pulse repetition frequency, with higher frequencies requiring more pulses to maximize the overlap region. These findings provide critical insights on developing controlled plasma discharge techniques for efficient NH 3 ignition in reactive flows within internal combustion engines and gas turbines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1999

Planar Strain-Rate-Free Diffusion Flames: Initiation, Properties, and Extinction

An effectively strain-rate-free diffusion flame constitutes the most vigorous laminar combustion of initially unmixed reactive gases. Such a diffusion flame is characterized by a relatively long residence time and by a relatively large characteristic length scale. If such a flame were also planar, providing high symmetry, it would be particularly suitable for experimental and theoretical investigations of key combustion phenomena, such as multicomponent diffusion, chemical kinetics, and soot inception, growth, and oxidation. Unfortunately, a planar strain-rate-free diffusion flame is highly disrupted in earth-gravity (e.g., in a counterflow-diffusion-flame apparatus) because of the very rapid onset (approx. 100 ms) of gravity-induced instability. Accordingly, a specially dedicated apparatus was designed, fabricated, and initially checked out for the examination of a planar strain-rate-free diffusion flame in microgravity. Such a diffusion flame may be formed within a hollowed-out squat container (initially configured as 25 cm x 25 cm x 9 cm), with isothermal, noncatalytic, impervious walls. At test initiation, a thin metallic sheet (approx. 1 mm in thickness) that separates the internal volume into two equal portions, each of dimensions 25 cm x 25 cm x 4.5 cm, is withdrawn, by uniform translation (approx. 50 cm/s) in its own plane, through a tightly fitting slit in one side wall. Thereupon, diluted fuel vapor (initially confined to one half-volume of the container) gains access to diluted oxygen (initially with the same pressure, density, and temperature as the fuel, but initially confined to the other half-volume). After a brief delay (approx. 10 ms), to permit limited but sufficient-for-flammability diffusional interpenetration of fuel vapor and oxidizer, burning is initiated by discharge of a line igniter, located along that side wall from which the trailing edge of the separator withdraws. The ignition spawns a triple-flame propagation across the 25 cm x 25 cm centerplane. When a diffusion flame is emplaced in the centerplane, any subsequent travel, and change in temperature, of that planar diffusion flame may be tracked, along with the effectively spatially uniform but temporally evolving pressure within the container. Eventually, nearly complete depletion of the stoichiometrically deficient reactant, along with heat loss to the container surfaces, effects extinction. These data afford an opportunity to check theoretical models of diffusion and chemical kinetics under conditions ranging from intense burning to flame out, or, alternatively, to evolve simple empirical representations of these phenomena. Thus, the project sought to utilize microgravity testing to elucidate commonly encountered phenomenology, arising in the commonly-encountered mode of combustion (whether related to heating, manufacturing, boiling, and propulsion, or to uncontrolled, free-burning fire in structures and wildland vegetation), of those commonly utilized fuels usually categorized as gaseous fuels (such as hydrogen, natural gas, and propane, which are gaseous under atmospheric conditions).

Fendell, Francis↗

DOE OSTI2023

Coupled Decay Heat and Thermal Hydraulic Capability for Loss-of-Coolant Accident Simulations

As the nuclear energy industry considers ways to achieve improved economics in the current fleet of light-water reactors (LWRs), one possible approach is to operate each cycle for longer durations. This causes a greater portion of the fuel to be burned and reduces the frequency of outages, which ultimately reduces the cost to operate the reactor. However, this also leads to higher burnup fuels than has traditionally been allowed in these reactors. Thus, there are concerns about integrity of high-burnup (HBu) fuel, especially during accident conditions such as loss-of-coolant accidents (LOCAs), as shown by Capps et al.. To investigate these concerns, advanced modeling and simulation capabilities are being leveraged to determine the susceptibility of HBu fuel to fuel fragmentation, relocation, and dispersion (FFRD). Improvements have previously been made to fuel performance capabilities to more accurately model these phenomena; multiphysics simulations have also been conducted to determine the power and burnup histories of the HBu fuel, which are needed as inputs for the fuel performance calculations. Most recently, new statistical approaches have been developed to identify a subset of fuel rods that have greater FFRD susceptibility, reducing the total number of fuel performance simulations required. Prior LOCA simulations have relied on the TRACE systems code, which can model the core and primary loop during accident conditions. TRACE includes many models for various aspects of the primary loop, but two sets of models are important for this report. First, TRACE uses a lumped-fuel approach for modeling the core. This approximates the ~50,000 fuel rods in the core with a much smaller number of rods. The rods can be lumped in various ways as determined by the user. For example, one lumped rod may be used to represent all rods in an assembly, sometimes with an additional rod representing the hottest fuel rod. However, due to runtime constraints and complexity of modeling, a more common approach is to group several assemblies or larger regions of the core into single lumped rods. These lumping schemes apply not only to fuel rods but to flow channels as well. Second, TRACE has several different models for treating decay heat, ranging from pregenerated decay heat curves based on an ANSI/ANS-5.1 standard (hereinafter abbreviated simply as ANSI) to explicit time-dependent heat inputs from the user. None of these models account for differences in isotopics between different rods, which is an approximation the work in this report seeks to eliminate. This report focuses on the implementation of coupled decay heat capabilities in the Virtual Environment for Reactor Applications (VERA) code suite to address a gap identified in previous LOCA simulations. This constitutes an improvement for both the lumped-fuel and decay heat models in TRACE. VERA has been developed to perform high-fidelity, whole-core multiphysics simulations for LWRs. Previously, during the Consortium for Advanced Simulation of LWRs (CASL) program, the emphasis was on providing accurate steady-state analysis—with a secondary focus on reactivity insertion accident (RIA) analysis—to address operational challenges in the nuclear energy industry. Under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, these capabilities are being extended to a broader range of transient analyses with the goal of quantifying the risk of fuel failures such as FFRD. To properly model such conditions with VERA, decay heat calculations have been integrated with the multiphysics to enable rod-by-rod thermal hydraulic (TH) conditions to be driven by the decay heat in long-running accidents such as LOCAs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗