Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive transport 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.

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

At least 433 records · Page 24

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗

Progress and Challenges in Quantifying Wildfire Smoke Emissions, Their Properties, Transport, and Atmospheric Impacts

Wildfire is a natural and integral ecosystem process that is necessary to maintain species composition, structure, and ecosystem function. Extreme fires have been increasing over the last decades, which have a substantial impact on air quality, human health, the environment, and climate systems. Smoke aerosols can be transported over large distances, acting as pollutants that affect adjacent and distant downwind communities and environments. Fire emissions are a complicated mixture of trace gases and aerosols, many of which are short‐lived and chemically reactive, and this mixture affects atmospheric composition in complex ways that are not completely understood. We present a review of the current state of knowledge of smoke aerosol emissions originating from wildfires. Satellite observations, from both passive and active instruments, are critical to providing the ability to view the large‐scale influence of fire, smoke, and their impacts. Progress in the development of fire emission estimates to regional and global chemical transport models has advanced, although significant challenges remain, such as connecting ecosystems and fuels burned with dependent atmospheric chemistry. Knowledge of the impact of smoke on radiation, clouds, and precipitation has progressed and is an essential topical research area. However, current measurements and parameterizations are not adequate to describe the impacts on clouds of smoke particles (e.g., CNN, INP) from fire emissions in the range of representative environmental conditions necessary to advance science or modeling. We conclude by providing recommendations to the community that we believe will advance the science and understanding of the impact of fire smoke emissions on human and environmental health, as well as feedback with climate system.

I N Sokolik↗

Development of a Computational Framework for Multiphysics Multiphase Species Tracking using NEAMS Tools

This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track isotopic species in Molten Salt Reactors (MSRs), with a specific focus on the 91-depletion chain within the Molten Salt Reactor Experiment (MSRE). The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the production, transport, and phase transitions of isotopes under steady-state and transient conditions. The main findings reveal that isotopes such as bromine-91 largely remain in the liquid phase, while others, including krypton-91and yttrium-91, transition to the gas phase, significantly influencing the reactor’s radiological source term. The study also shows that during transients, like a reactivity insertion transient, rapid void formation and the expansion of the liquid-gas interface led to substantial transfers of dissolved isotopes into the gas phase, altering isotope distribution and largely increasing the source term in the off-gas system. Additionally, the research highlights that short-lived isotopes dominate the initial off-gas response during transients, while longer-lived isotopes determine the equilibrium state, underscoring the necessity of dynamic simulations for accurate species tracking and reactor safety analysis. The developed methodology will be applied in the future to the tracking of a larger number of species and introduce other species tracking mechanisms, such as deposition and plating.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Computational design of microarchitected porous electrodes for redox flow batteries

Porous electrodes are used as the core reactive component across electrochemical technologies. In flowing systems, controlling the fluid distribution, species transport, and reactive environment is critical to attaining high performance. However, conventional electrode materials like felts and papers provide few opportunities for precise engineering of the electrode and its microstructure. To address these limitations, architected electrodes composed of unit cells with spatially varying geometry determined via computational optimization are proposed. Resolved simulation is employed to develop a homogenized description of the constituent unit cells. These effective properties serve as inputs to a continuum model for the electrode when used in the negative half-cell of a vanadium redox flow battery. Porosity distributions minimizing power loss are then determined via computational design optimization to generate architected porosity electrodes. The architected electrodes are compared to bulk, uniform porosity electrodes and found to lead to increased power efficiency across operating flow rates and currents. The design methodology is further used to generate a scaled-up electrode with comparable power efficiency to the bench-scale systems. Finally, the variable porosity architecture and computational design methodology presented here thus offers a novel pathway for automatically generating spatially engineered electrode structures with improved power performance.

25 ENERGY STORAGE↗

Following O and OH in He/O 2 and He/H 2 O gas mixtures—from the gas phase through the liquid phase to modifications on a biological sample

Abstract Applied cold atmospheric plasma allows for the controlled delivery of reactive oxygen and nitrogen species tailored for specific applications. Through the manipulation of the plasma parameters, feed gases, and careful consideration of the environment surrounding the treatment target, selective chemistries that preferentially influence the target can be produced and delivered. To demonstrate this, the COST reference microscale atmospheric pressure plasma jet is used to study the generation and transport of O and ⋅ OH from the gas phase through the liquid to the biological model target cysteine. Relative and absolute species densities of ⋅ OH and O are measured in the gas phase through laser induced fluorescence (LIF) and two-photon absorption LIF respectively. The transport of these species is followed into the liquid phase by hydrogen peroxide quantification and visualized by a fluorescence assay. Modifications to the model biological sample cysteine exposed to ⋅ OH and H 2 O 2 dominated chemistry (He/H 2 O (0.25%)) and O dominated chemistry (He/O 2 (0.6%)) is measured by FTIR spectroscopy. The origin of these species that modify cysteine is considered through the use of heavy water (H 2 18 O) and mass spectrometry. It is found that the reaction pathways differ significantly for He/O 2 and He/H 2 O. Hydrogen peroxide is formed mainly in the liquid phase in the presence of a substrate for He/O 2 whereas for He/H 2 O it forms in the gas phase. The liquid chemistry resulting from the He/O 2 admixture mainly targets the sulfur moiety of cysteine for oxidation up to irreversible oxidation states, while He/H 2 O treatment leads preferentially to reversible oxidation products. The more O or OH/H 2 O 2 dominated chemistry produced by the two gas admixtures studied offers the possibility to select species for target modification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Automatic Volume Balancing for Online Refueling in Molten Salt Reactor Simulations with SCALE

This work introduces recently implemented capabilities in the SCALE code system’s TRITON reactor physics sequence that improve molten salt reactor (MSR) modeling: (1) a volume balancing option, which automatically balances the volumes of the fed material with a corresponding material removal to maintain fixed mixture volumes in the neutron transport model and accurate densities, and (2) continuous feed from mixtures, which enables users to define material feed streams directly from salt mixture definitions rather than individual nuclides. The new capabilities were demonstrated via SCALE/TRITON simulations of two representative MSR concepts. Simulations of the 180 MWth molten chloride fast reactor—a system with a fixed fuel salt volume in the reactor core—applied continuous refueling with fresh fuel salt. The results confirm approximately constant reactivity when the new volume balancing capability is used. Additionally, excellent agreement with a manual feed-and-drain method for volume balancing further verified the new implementation. Simulations of the 400 MWth EIRENE reactor, an integral MSR in which the salt volume may grow over time within the reactor core, demonstrated the capability to represent growing salt volume within the TRITON depletion calculation. Compared with reference solutions, the results show consistent trends in reactivity and isotopic evolution, confirming that these new user-friendly capabilities provide accurate, physically consistent approaches for modeling MSRs in SCALE.

Elzohery, Rabab [ORNL] (ORCID:0000000160043633)↗

2D-to-3D image translation of complex nanoporous volumes using generative networks

Image-based characterization offers a powerful approach to studying geological porous media at the nanoscale and images are critical to understanding reactive transport mechanisms in reservoirs relevant to energy and sustainability technologies such as carbon sequestration, subsurface hydrogen storage, and natural gas recovery. Nanoimaging presents a trade off, however, between higher-contrast sample-destructive and lower-contrast sample-preserving imaging modalities. Furthermore, high-contrast imaging modalities often acquire only 2D images, while 3D volumes are needed to characterize fully a source rock sample. In this work, we present deep learning image translation models to predict high-contrast focused ion beam-scanning electron microscopy (FIB-SEM) image volumes from transmission X-ray microscopy (TXM) images when only 2D paired training data is available. We introduce a regularization method for improving 3D volume generation from 2D-to-2D deep learning image models and apply this approach to translate 3D TXM volumes to FIB-SEM fidelity. We then segment a predicted FIB-SEM volume into a flow simulation domain and calculate the sample apparent permeability using a lattice Boltzmann method (LBM) technique. Results show that our image translation approach produces simulation domains suitable for flow visualization and allows for accurate characterization of petrophysical properties from non-destructive imaging data.

58 GEOSCIENCES↗

Quantifying Dissolution Dynamics in Porous Media Using a Spatial Flow Focusing Profile

Abstract The diverse range of patterns in porous media formed by dissolution processes depends on the relative magnitude of flow, transport, and chemical reactions at pore surfaces. However, distinguishing between regimes often relies solely on qualitative, visual comparisons of emergent structures. Here, we propose a quantitative measure capable of identifying different regimes using the concept of the spatial flow focusing profile, which segments the medium into cross sections along the flow direction to calculate the flow focusing index for each section. We employ this measure in numerical simulations of a dissolving porous medium using a pore network model. We obtain a morphological phase diagram of dissolution patterns, which we characterize using the flow focusing profile. In particular, we demonstrate that analyzing the temporal changes in the profile allows one to quantitatively distinguish between wormholing and channeling. The transition between them is shown to be affected by the heterogeneity of the system.

58 GEOSCIENCES↗

Corrosion Response of Ni-19Cr-5Fe in Molten NaCl-MgCl 2 Salt

Ni-based alloys are strong candidates for use in high-temperature molten salt reactors due to their superior corrosion resistance and mechanical stability compared to stainless steels. In this study, the static corrosion behavior of a Ni-19Cr-5Fe model alloy was systematically evaluated in a purified molten NaCl-MgCl 2 salt at 700°C for 30, 240, and 500 h. Post-exposure analyses were conducted to assess microstructural evolution, corrosion depth, and elemental depletion profiles. Corrosion rates, quantified by chromium depletion depth, followed an inverse power-law trend with increasing exposure time, indicating diffusion-limited kinetics. This trend is attributed to limitations in mass transport in the static salt and to the progressive local depletion of reactive chromium species. Coupled electron backscatter diffraction and energy dispersive X-ray spectroscopy analysis further revealed that the grain boundary character significantly influences corrosion susceptibility: high-angle grain boundaries exhibited pronounced Cr depletion and pitting, while low-angle and Σ3 boundaries remained comparatively resistant. These results offer valuable insight into the role of microstructure in corrosion processes and reinforce the importance of time-dependent material evaluation in molten salt environments relevant to advanced reactor designs.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling a generic TRISO-fueled heat pipe microreactor using SCALE: Depletion, transportation criticality, and shielding

This paper demonstrates the applicability of the SCALE code system to tristructural-isotropic (TRISO)-fueled heat pipe microreactors through depletion, transportation criticality, and shielding analyses of a generic reference design. The study conducted supports US Nuclear Regulatory Commission code readiness efforts for advanced non–light-water reactor technologies and is intended as a code capability demonstration rather than as an optimization of a specific microreactor design. The modeled reactor employs high-assay low-enriched uranium (HALEU) uranium oxycarbide (UCO) TRISO fuel and beryllium oxide (BeO) reflectors and operates at 7.5 MWth with a nominal lifetime of about 3 effective full power years. Representative cases for fresh and irradiated cores were selected to exercise SCALE methods relevant to reactor operation and post-irradiation transport. The discharged-core decay heat is approximately 6% of operating power immediately after shutdown. Transportation criticality calculations show that internal water ingress is the dominant reactivity effect, with fully flooded fresh core and irradiated core configurations remain above the subcriticality criterion, even with the available control mechanisms. Shielding calculations for a simplified transportation package indicate that normal-condition dose rates are governed mainly by shielding thickness and cooling time, whereas the breached hypothetical accident case is governed primarily by cooling time. Overall, the study shows that SCALE supports depletion, transportation criticality, and shielding evaluations efficiently for TRISO-fueled heat pipe microreactors within a single code system.

Criticality↗

The effect of water vapor on surface oxygen exchange kinetics of thin film (La,Sr)(Co,Fe)O 3-δ

Oxygen reduction and transport in electrochemical applications are strongly affected by the reactant gas composition, including impurities such as water. Understanding the relationship between oxygen reactivity and humidity is key for attaining stability and high efficiency of solid oxide fuel cells, which are especially attractive for transportation applications as low emission power generation devices. Both short- and long-term effects of moisture have previously been associated with use of ambient air. To better understand these effects, we study oxygen exchange kinetics at the cathode surface and the cathode/electrolyte interface using La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3-δ (LSCF) epitaxial thin films as a model-surface system. In situ synchrotron X-ray techniques evaluate the oxygen reduction reaction (ORR) kinetics in response to environmental variables. The results suggest a clear enhancement of the ORR rate upon short term exposure of the perovskite-structured films to water. Finally, on the basis of these measurements, along with in situ X-ray characterization of Sr segregation and computational first-principles studies, we suggest a model for this increased ORR activity that may lead to further improvements by stabilizing active cathode surface configurations.

25 ENERGY STORAGE↗

Florida Thunderstorms: A Faucet of Reactive Nitrogen to the Upper Troposphere

During the NASA Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) enhanced mixing ratios of nitric oxide were measured in the anvils of thunderstorms and in clear air downwind of storm systems on flights of a Wl3-57F high-altitude aircraft. Mixing ratios greater than l0 - 20 times background were readily observed over distances of 25-120 km due to lightning activity. In many of the Florida storms deposition of NO occurred up to near the tropopause but major deposition usually occurred 1 - 2 km below the tropopause, or mostly within the visible anvil volume formed prior to storm decay. Observations from two storms of very different anvil size and electrical activity allowed estimates of the total mass of NO, vented to the middle and upper troposphere. Using the cloud-to ground (CG) flash accumulations from the National Lightning Detection Network, climatological intra-cloud (IC) to CG ratios, and assuming that CG and IC flashes were of equivalent efficiency for NO production, the ranges of production per flash for a moderate-sized and a large storm were (0.51 - 1.0) x l0(exp 26) and (2.3 - 3.1) x 10(exp 26) molecules NO/flash, respectively. Using the recently determined average global flash rate of 44 8, a gross extrapolation of these two storms to represent possible global annual production rates yield 1.6 - 3.2 and 7.3 - 9.9 Tg(N)/yr, respectively. If the more usual assumption is made that IC efficiency is l/l0th that of CG activity, the ranges of production for the moderate-sized and large storm were (1.3 - 2.7) x l0(exp 26) and (6.0 - 8.1) x l0(exp 26) molecules NO/CG flash, respectively. The estimates from the large storm may be high because there is indirect evidence that the IC/CG ratio was larger than would be derived from climatology. These two storms and others studied did not have flash rates that scaled as approx. H(sup 5) where H is the cloud top altitude. The observed CG flash accumulations and NO(x) mass production estimate for the month of July over the Florida area were compared with a representative 3D global Chemistry-Transport Model (CTMJ that uses the Price et al. lightning parameterization. For two land grid points representing the Florida peninsula the model compared well with the observations: CG flash rates were low by only a factor of approx. 2. When the model grid points included the coastal regions of Florida the flash accumulations were lower than observed by a factor of 3.4 - 4.6. It is recommended that models using the Price et al. parameterization allow any global coastal grid point to maintain the land rather than the marine flash rate parameterization. The convection in this CTM underestimated the actual cloud top heights over Florida by 1 - 2 km and thus the total lightning flash rates and the altitude range of reactive nitrogen deposition. Broad scale (20 - 120 km) median mixing ratios of NO within anvils over Florida were significantly larger than in storms previously investigated over Colorado and New Mexico.

Ridley, B.↗

Boost-phase discrimination research activities

Theoretical research in two areas was performed. The aerothermodynamics research focused on the hard-body and rocket plume flows. Analytical real gas models to describe finite rate chemistry were developed and incorporated into the three-dimensional flow codes. New numerical algorithms capable of treating multi-species reacting gas equations and treating flows with large gradients were also developed. The computational chemistry research focused on the determination of spectral radiative intensity factors, transport properties and reaction rates. Ab initio solutions to the Schrodinger equation provided potential energy curves transition moments (radiative probabilities and strengths) and potential energy surfaces. These surfaces were then coupled with classical particle reactive trajectories to compute reaction cross-sections and rates.

Cooper, David M.↗

An overview of the Evaluation of Oxygen Interactions with Materials 3 experiment: Space Shuttle Mission 46, July-August 1992

The Evaluation of Oxygen Interactions with Materials 3 (EOIM-3) flight experiment was developed to obtain benchmark atomic oxygen reactivity data and was conducted during Space Transportation System Mission 46 (STS-46), July 31 to August 7, 1992. In this paper, we present an overview of EOIM-3 and the results of the Lyndon B. Johnson Space Center (JSC) materials reactivity and mass spectrometer/carousel experiments. Mass spectrometer calibration methods are discussed briefly, as a prelude to a detailed discussion of the mass spectrometric results produced during STS-46. Mass spectrometric measurements of ambient O-atom flux and fluence are in good agreement with the values calculated using the MSIS-86 model of the thermosphere as well as estimates based on the extent of O-atom reaction with Kapton polyimide. Mass spectrometric measurements of gaseous products formed by O-atom reaction with C(13) labeled Kapton revealed CO, CO2, H2O, NO, and NO2. Finally, by operating the mass spectrometer so as to detect naturally occurring ionospheric species, we characterized the ambient ionosphere at various times during EOIM-3 and detected the gaseous reaction products formed when ambient ions interacted with the C(13) Kapton carousel sector. By direct comparison of the results of on-orbit O-atom exposures with those conducted in ground-based laboratory systems, which provide known O-atom fluences and translational energies, we have demonstrated the strong translational energy dependence of O-atom reactions with a variety of polymers. A 'line-of-centers' reactive scattering model was shown to provide a reasonably accurate description of the translational energy dependence of polymer reactions with O atoms at high atom kinetic energies while a Beckerle-Ceyer model provided an accurate description of O-atom reactivity over a three order-of-magnitude range in translational energy and a four order-of-magnitude range in reaction efficiency. Postflight studies of the polymer samples by x-ray photoelectron spectroscopy and infrared spectroscopy demonstrate that O-atom attack is confined to the near-surface region of the sample, i.e. within 50 to 100 A of the surface.

Koontz, Steven L.↗

Using Radiogenic Noble Gas Nuclides to Identify and Characterize Rock Fracturing

Abstract Fracture‐released radiogenic noble gas nuclides are used to identify locations and constrain the volume of new fracture creation during subsurface detonations. Real‐time, in situ noble gases and reactive gases were monitored using a field‐deployed mass spectrometer and automated sampling system in a multilevel borehole array. Released gases were measured after two different detonations having distinct energy, pressure, and gas volume characteristics. Explosive‐derived gases (N 2 O, CO 2 ) and excess radiogenic 4 He and 40 Ar above atmospheric background are used to identify locations of gas transport and new fracture creation after each detonation. Fracture‐released radiogenic 4 He is used to constrain the volume of newly created fractures with a model of helium release from fracturing. Explosive by‐product gas was observed in multiple locations both near and distal to the shot locations for both detonations. Radiogenic 4 He and 40 Ar release from rock damage was observed in locations near the detonation after the second, more powerful detonation. Observed 4 He response is consistent with a model of diffusive release from newly created fractures. Volume of new fractures estimated from the 4 He release ranges from 1 to 5 m 2 with apertures ranging from 0.1 to 1 m. Our results provide evidence that radiogenic noble gases released during fracture creation can be identified at the field scale in real time and used to identify timing and location of fracture creation during deformation events. This technique could be useful in subsurface science and engineering problems where the location and amount of newly created rock fracturing is of interest including fault rupture, mine safety, subsurface detonation monitoring and reservoir stimulation.

58 GEOSCIENCES↗

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

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

Nikjooy, Mohammad↗

Non-equilibrium chemistry in the solar nebula and early solar system: Implications for the chemistry of comets

Theoretical models of solar nebula and early solar system chemistry which take into account the interplay between chemical, physical, and dynamical processes have great utility for deciphering the origin and evolution of the abundant chemically reactive volatiles (H, O, C, N, S) observed in comets. In particular, such models are essential for attempting to distinguish between presolar and solar nebula products and for quantifying the nature and duration of nebular and early solar system processing to which the volatile constituents of comets have been subjected. The diverse processes and energy sources responsible for chemical processing in the solar nebula and early solar system are discussed. The processes considered include homogeneous and heterogeneous thermochemical and photochemical reactions, and disequilibration resulting from fluid transport, condensation, and cooling whenever they occur on timescales shorter than those for chemical reactions.

Fegley, Bruce, Jr.↗

Radicals and Reservoirs in the GMI Chemistry and Transport Model: Comparison to Measurements

We have used a three-dimensional chemistry and transport model (CTM), developed under the Global Modeling Initiative (GMI), to carry out two simulations of the composition of the stratosphere under changing halogen loading for 1995 through 2030. The two simulations differ only in that one uses meteorological fields from a general circulation model while the other uses meteorological fields from a data assimilation system. A single year's winds and temperatures are repeated for each 36-year simulation. We compare results from these two simulations with an extensive collection of data from satellite and ground-based measurements for 1993-2000. Comparisons of simulated fields with observations of radical and reservoir species for some of the major ozone-destroying compounds are of similar quality for both simulations. Differences in the upper stratosphere, caused by transport of total reactive nitrogen and methane, impact the balance among the ozone loss processes and the sensitivity of the two simulations to the change in composition.

Douglass, Anne R.↗