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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Opposed Jet Burner Approach for Characterizing Flameholding Potentials of Hydrocarbon Scramjet Fuels

Opposed Jet Burner (OJB) tools have been used extensively by the authors to measure Flame Strength (FS) extinction limits of laminar H2/N2 air and (recently) hydrocarbon (HC) air Counterflow Diffusion Flames (CFDFs) at one atm. This paper details normalization of FSs of N2- diluted H2 and HC systems to account for effects of fuel composition, temperature, pressure, jet diameter, inflow Reynolds number, and inflow velocity profile (plug, contoured nozzle; and parabolic, straight tube). Normalized results exemplify a sensitive accurate means of validating, globally, reduced chemical kinetic models at approx. 1 atm and the relatively low temperatures approximating the loss of non-premixed idealized flameholding, e.g., in scramjet combustors. Laminar FS is defined locally as maximum air input velocity, U(sub air), that sustains combustion of a counter-jet of g-fuel at extinction. It uniquely characterizes a fuel. And global axial strain rate at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or (sub t)) can be compared directly with computed extinction limits, determined using either a 1-D Navier Stokes stream-function solution, using detailed transport and finite rate chemistry, or (better yet) a detailed 2-D Navier Stokes numerical simulation. The experimental results define an idealized flameholding reactivity scale that shows wide ranging (50 x) normalized FS s for various vaporized-liquid and gaseous HCs, including, in ascending order: JP-10, methane, JP-7, n-heptane, n-butane, propane, ethane, and ethylene. Results from H2 air produce a unique and exceptionally strong flame that agree within approx. 1% of a recent 2-D numerically simulated FS for a 3 mm tube-OJB. Thus we suggest that experimental FS s and/or FS ratios, for various neat and blended HCs w/ and w/o additives, offer accurate global tests of chemical kinetic models at the Ts and Ps of extinction. In conclusion, we argue the FS approach is more direct and fundamental, for assessing, e.g., idealized scramjet flameholding potentials, than measurements of laminar burning velocity or blowout in a Perfectly Stirred Reactor, because the latter characterize premixed combustion in the absence of aerodynamic strain. And FS directly measures a chemical kinetic characteristic of non-premixed combustion at typical flameholding temperatures. It mimics conditions where gfuels are typically injected into a subsonic flameholding recirculation zone that captures air, where the effects of aerodynamic strain and associated multi-component diffusion become important.

Pellett, Gerald L.↗
NASA NTRS2002

Evaluation of Transport in the Lower Tropical Stratosphere in a Global Chemistry and Transport Model

A general circulation model (GCM) relies on various physical parameterizations and provides a solution to the atmospheric equations of motion. A data assimilation system (DAS) combines information from observations with a GCM forecast and produces analyzed meteorological fields that represent the observed atmospheric state. An off-line chemistry and transport model (CTM) can use winds and temperatures from a either a GCM or a DAS. The latter application is in common usage for interpretation of observations from various platforms under the assumption that the DAS transport represents the actual atmospheric transport. Here we compare the transport produced by a DAS with that produced by the particular GCM that is combined with observations to produce the analyzed fields. We focus on transport in the tropics and middle latitudes by comparing the age-of-air inferred from observations of SF6 and CO2 with the age-of-air calculated using GCM fields and DAS fields. We also compare observations of ozone, total reactive nitrogen, and methane with results from the two simulations. These comparisons show that DAS fields produce rapid upward tropical transport and excessive mixing between the tropics and middle latitudes. The unrealistic transport produced by the DAS fields may be due to implicit forcing that is required by the assimilation process when there is bias between the GCM forecast and observations that are combined to produce the analyzed fields. For example, the GCM does not produce a quasi-biennial oscillation (QBO). The QBO is present in the analyzed fields because it is present in the observations, and systematic implicit forcing is required by the DAS. Any systematic bias between observations and the GCM forecast used to produce the DAS analysis is likely to corrupt the transport produced by the analyzed fields. Evaluation of transport in the lower tropical stratosphere in a global chemistry and transport model.

Douglass, Anne R.↗
DOE OSTI2024

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗
NASA NTRS1996

A Burst Mode, Ultrahigh Temperature UF4 Vapor Core Reactor Rankine Cycle Space Power System Concept

Static and dynamic neutronic analyses have been performed on an innovative burst mode (100's of MW output for a few thousand seconds) Ulvahigh Temperature Vapor Core Reactor (UTVR) space nuclear power system. The NVTR employs multiple, neutronically-coupled fissioning cores and operates on a direct, closed Rankine cycle using a disk Magnetohydrodynamic (MHD) generater for energy conversion. The UTVR includes two types of fissioning core regions: (1) the central Ultrahigh Temperature Vapor Core (UTVC) which contains a vapor mixture of highly enriched UF4 fuel and a metal fluoride working fluid and (2) the UF4 boiler column cores located in the BeO moderator/reflector region. The gaseous nature of the fuel the fact that the fuel is circulating, the multiple coupled fissioning cores, and the use of a two phase fissioning fuel lead to unique static and dynamic neutronic characteristics. Static neutronic analysis was conducted using two-dimensional S sub n, transport theory calculations and three-dimensional Monte Carlo transport theory calculations. Circulating-fuel, coupled-core point reactor kinetics equations were used for analyzing the dynamic behavior of the UTVR. In addition to including reactivity feedback phenomena associated with the individual fissioning cores, the effects of core-to-core neutronic and mass flow coupling between the UTVC and the surrounding boiler cores were also included in the dynamic model The dynamic analysis of the UTVR reveals the existence of some very effectlve inherent reactivity feedback effects that are capable of quickly stabilizing this system, within a few seconds, even when large positive reactivity insertions are imposed. If the UTVC vapor fuel density feedback is suppressed, the UTVR is still inherently stable because of the boiler core liquid-fuel volume feedback; in contrast, suppression of the vapor fuel density feedback in 'conventional" gas core cavity reactors causes them to become inherently unstable. Due to the strength of the negative reactivity feedback in the UTVR, it is found that external reactivity insertions alone are inadequate for bringing about significant power level changes during normal reactor operations. Additional methods of reactivity control such as variations in the gaseous fuel mass flow rate, are needed to achieve the desired power level oontrol.

Dugan, E. T.↗
DOE OSTI2022

Optimization of Canister Loading Patterns in Dual Purpose Canisters for Criticality Suppression

This report documents work performed in support of the US Department of Energy Office of Nuclear Energy (NE) Spent Fuel and Waste Disposition, Spent Fuel and Waste Science and Technology, under work breakdown structure element 1.08.01.03.05, “Direct Disposal of Dual Purpose Canisters.” In particular, this report fulfills milestone M3SF-22PN010305094, “Application of AI/ML techniques to optimize DPC loading,” within work package SF-22PN01030509, “Direct Disposal of Dual Purpose Canisters - PNNL.” This report continues the process of examining the potential for using loading optimization as a disposal criticality suppression technique. This latest update to the report added: 1. A validation of the artificial neural network (ANN) reactivity prediction tool against as-loaded dual-purpose canisters from the UNF-ST&DARDS database 2. Refinement of the algorithm to incorporate an additional DPC design to improve the predictions 3. Investigation of a storage and transportation focused loading optimization algorithm. A key result from this year’s work is that the ANN performance was significantly improved by including the additional model for the NUHOMS canisters and it is apparent that fuel type specific modeling considerations should be incorporated into the ANN in future work. The GRASP-enabled adaptive multi-objective memetic algorithm with partial clustering (GAMMA-PC) algorithm was evaluated as a tool for storage and transportation oriented optimization. The GAMMA-PC optimization routine considers decay heat at loading, the time between loading and when the canister is eligible for transportation and minimizes the number of casks loaded. Future work will develop a dose minimization optimization routine.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Deep Learning Potential For Investigation of Chemical Conversion of CO2 at Air-Reactive Deep Eutectic Solvent Interface

Capturing CO2 from air is urgent to battle against the climate change. Traditional CO2 capture by alkaline solution or aqueous amine has high regeneration cost, equipment corrosion, and amine leakage. Deep eutectic solvents (DESs) are promising alternatives with low corrosion, non-toxicity, and biodegradable nature. Functionalized DESs can capture and convert CO2 at low partial pressures. Developing reactive DESs for efficient CO2 capture needs a foundational understanding of the reaction-transport coupling of CO2 at the air-DES interface. In this study, we develop a deep learning potential (DLP) to investigate the absorption and reaction of CO2 at the air-DES interface (DES: 1-ethyl-3-methylimidazolium 2-cyanopyrrolide ([Emim][2-CNpyr]) and ethylene glycol (EG)). The DLP model allows us to simulate chemical reactions at a lower computational cost with ab initio accuracy. By analyzing reaction free energy surfaces, molecular interactions among CO2, DES components and reaction products, and transport of CO2 through the interface, we identify three main reaction pathways of CO2 at the interface: forming carboxylate with [Emim], carbamate with [2-CNpyr], and carbonate with EG. The mechanistic understanding of CO2 chemisorption at the interface will facilitate the development of novel direct air capture technology using DESs.

Nguyen, Manh Tien↗
DOE OSTI2022

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Formulation, Implementation and Validation of a 1D Boundary Layer Inflow Scheme for the QUIC Modeling System

Recent studies have highlighted the importance of accurate meteorological conditions for urban transport and dispersion calculations. In this work, we present a novel scheme to compute the meteorological input in the Quick Urban & Industrial Complex () diagnostic urban wind solver to improve the characterization of upstream wind veer and shear in the Atmospheric Boundary Layer (ABL). The new formulation is based on a coupled set of Ordinary Differential Equations (ODEs) derived from the Reynolds Averaged Navier–Stokes (RANS) equations, and is fast to compute. Building upon recent progress in modeling the idealized ABL, we include effects from surface roughness, turbulent stress, Coriolis force, buoyancy and baroclinicity. We verify the performance of the new scheme with canonical Large Eddy Simulation (LES) tests with the GPU-accelerated FastEddy"Equation missing" solver in neutral, stable, unstable and baroclinic conditions with different surface roughness. Furthermore, we evaluate QUIC calculations with and without the new inflow scheme with real data from the Urban Threat Dispersion (UTD) field experiment, which includes Lidar-based wind measurements as well as concentration observations from multiple outdoor releases of a non-reactive tracer in downtown New York City. Compared to previous inflow capabilities that were limited to a constant wind direction with height, we show that the new scheme can model wind veer in the ABL and enhance the prediction of the surface cross-isobaric angle, improving evaluation statistics of simulated concentrations paired in time and space with UTD measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Strain-Driven Surface Reactivity in Magnesium-Ion Battery Cathodes

Successful deployment of a Mg-ion battery requires cathodes that can achieve reversible Mg intercalation and high energy density. Recent theoretical and experimental studies indicated that the overall transport is likely limited by sluggish Mg transport at the cathode-electrolyte interface and not Mg diffusion through bulk. Here, in this work, we investigated the surface electrochemical activity of Mg ions by using a spinel-structured manganese oxide thin-film model system and in situ X-ray scattering. In combination with post-mortem microscopy analysis, we found that magnesium insertion was more favorable than subsequent extraction near the surface of the Mg x Mn 2 O 4 film, resulting in overmagnesiation, and eventually amorphization of the surface. This structural irreversibility and high overpotential required for Mg extraction could explain significant voltage hysteresis and Mg surface enrichment previously observed in bulk cathodes. Density functional theory calculations suggested that the tendency for the Mg surface enrichment could be associated with Mg diffusion kinetics, which varies with the strain state evolved due to constrained film volume change during Mg insertion and extraction. Particularly, out-of-plane Mg migration was predicted to be favorable in the tensile strain rather than in the compressive case.

25 ENERGY STORAGE↗
DOE OSTI2020

Review of Tools for Modeling Core Radial Expansion in Liquid Metal-Cooled Fast Reactors

Core radial expansion in liquid-metal cooled fast reactor systems is a well-known phenomenon that produces strong reactivity feedback effects. It is crucial from a safety standpoint to design a reactor with the proper core restraint system to guide expansion of the fuel assemblies into a formation that produces negative reactivity feedback in accident conditions. However, the modeling of core radial expansion and its associated reactivity feedback is extremely challenging. Liquid metal-cooled fast reactor designers (both industry and commercial) have pointed to improved modeling of core radial expansion as a key modeling challenge for these types of reactors. The Department of Energy Nuclear Energy Advanced Modeling and Simulation (DOE-NEAMS) program has an important role to play in helping to develop a robust, predictive tool for this phenomenon. A literature review has been performed to identify current simulation capabilities applicable to the prediction of core radial expansion in liquid metal-cooled fast reactors and the resultant reactivity feedback effects. An overview as well as a detailed description of the multi-physics phenomena underlying core radial expansion is given in the first two chapters. The objective of this report is to identify existing software that has been or could be applied to the core radial expansion problem, and then to identify possible paths to develop a robust, coupled Multiphysics tool for advanced analysis, in which all three physics (radiation transport, structural mechanics, thermal fluids) are taken into account considering feedback effects from each physics to the other. Several computer codes and workflows have been developed in the United States and abroad to perform individual aspects of radial expansion calculations or to perform loosely coupled analysis. No code system currently exists which tightly and robustly couples the neutronics, thermal hydraulics, and thermal mechanical physics with enough detail to fully resolve the complex core radial expansion reactivity feedback effects. The future availability of such a code system is of vital importance to fully understanding the reactivity feedback effects that occur due to radial expansion, and consequently to optimizing the design of the core restraint system. A potential code development path forward using MOOSE-based tools is suggested in order to leverage the natural tight coupling and robustness of MOOSE-based applications. Additionally, potential gaps for modeling and common assumptions are listed which require further investigation to ensure accuracy for prediction of deformations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Hydration-Controlled Proton Transport in Respiratory Complex I

Proton pumping by respiratory complex I is one essential element for generating the proton motive force that drives ATP synthesis in mitochondria. Although it is understood that electrons from NADH reduce ubiquinone at the peripheral arm and that four protons are transferred in the membrane domain, the mechanism by which this redox reaction initiates proton translocation remains unclear. A lateral pathway linking the quinone binding site to the membrane domain via ND1, ND3, and ND4L subunits has been proposed as a possible initial path of an excess proton. However, experimental structures indicate that the hydration connectivity between D66 ND3 and E34 ND4L is comparatively weaker than in neighboring segments, suggesting a potential regulatory point for proton transfer. Using multiscale reactive molecular dynamics (MS-RMD) and a water wire connectivity metric, we directly simulate proton transport through this region as coupled to the hydration by water molecules. Our results reveal that proton transfer is thermodynamically feasible when transient hydration aligns with the presence of an excess proton, revealing the strong coupling between hydration and proton transfer (PT) in this region of Complex I. These findings support a model where proton injection enhances local hydration, dynamically opening the pathway for proton transfer and regulating the onset of proton pumping in Complex I.

Chemical structure↗
NASA NTRS2006

On the Quality of the Nimbus 7 LIMS Version 6 Ozone for Studies of the Middle Atmosphere

The Nimbus 7 Limb Infrared Monitor of the Stratosphere (LIMS) radiance profile dataset of 1978/79 was reconditioned and reprocessed to Version 6 (V6) profiles of temperature and species that are improved significantly over those from Version 5 (V5). The LIMS V6 dataset was archived for public use in 2002. Improvements for its ozone include: (1) a more accurate accounting for instrument and spacecraft motion effects in the radiances, (2) the use of better spectroscopic line parameters for its ozone forward model, (3) retrievals of all its scans, (4) more accurate and compatible temperature versus pressure profiles (or T(p)), which are needed for the registration of the ozone radiances and for the removal of temperature effects from them, and (5) a better accounting for interfering species in the lower stratosphere. The retrieved V6 ozone profiles extend from near cloud top altitudes to about 80 km and from 64S to 84N latitude with better sampling along the orbit than for the V5 dataset. Calculated estimates of the single-profile precision and accuracy are provided for the V6 ozone from this study. Precision estimates based on the data themselves are of order 3% or better from 1 to 30 hPa. Estimates of total systematic error for a single profile are hard to generalize because the separate sources of error may not all be of the same sign and they depend somewhat on the atmospheric state. It is estimated that the V6 zonal mean ozone distributions are accurate to within 9% to 7% from 50 hPa to 3 hPa, respectively. Effects of a temperature bias can be significant and may be present at 1 to 2 hPa though. There may be ozone biases of order 10% at those levels due to possible biases of up to +2 K, but there is no indication of a similar problem elsewhere in the stratosphere. Simulation studies show that the LIMS retrievals are also underestimating slightly the small amplitudes of the atmospheric temperature tides, which affect its retrieved day/night ozone differences. There are small biases in the middle to lower stratosphere for the ascending versus descending node LIMS ozone, due principally to not accounting for the asymmetric weighting of its radiances across the tangent layer. The estimates of total accuracy were assessed by comparing the daily zonal mean LIMS ozone distributions against those from the Nimbus 7 SBUV Version 8 (V8) dataset for the same period. Generally, the LIMS V6 ozone agrees well with SBUV, except perhaps in the tropical lower stratosphere where the LIMS ozone is less. Still, the accuracy for LIMS V6 ozone in the lower stratosphere is improved over that found for LIMS V5, as indicated by several LIMS comparisons with ECC ozonesonde profiles. The LIMS V6 ozone is considered especially suitable for detailed studies of large-scale stratospheric processes above the 100-hPa level. Comparison of diurnal, photochemical model calculations with the monthly-averaged, upper stratospheric ozone observed from LIMS V6 indicates only a slight ozone deficit for the model at about 2 hPa. However, that deficit exhibits little to no seasonal variation and is in good agreement with similar model comparisons for a seasonal time series of ozone obtained with ground-based microwave instruments. Because the LIMS V6 ozone in the lower stratosphere has improved accuracy and sampling versus that of V5, it should now be possible to conduct quantitative studies of ozone transport and chemistry for the northern hemisphere, polar stratospheric winter of 1978/79 a time period when the catalytic loss of ozone due to reactive chlorine should not have been a major factor for the Arctic region.

Remsberg, Ellis↗
DOE OSTI2022

Development of the TREAT M2 experiment as a transient benchmark

A transient benchmark based on the Transient Reactor Test Facility (TREAT) M2 Calibration experiment (M2-CAL) is under development. TREAT, at Idaho National Laboratory, is a graphite moderated air-cooled research reactor which has been used extensively for fuel material testing under extreme and accident conditions. Accurate benchmark models are a beneficial component in the operations, experimental planning, and development of TREAT. In this work, we present a transient benchmark model for the M2-CAL experiment core loading. The benchmark model incorporates the coupling between the Monte Carlo Code SERPENT and the computational fluid dynamics code OpenFOAM to capture the temperature feedback mechanism. The M2-CAL transient 2580 was simulated in this work. A pre-transient analysis was performed to determine the optimum core composition and conditions before the beginning of the transient. The analysis was validated against historic TREAT kinetic measurements and the worth of the transient rod T-2. The axial power distribution in the flux wire for the M2-CAL experiment was determined and contrasted with the experiment. The model was then used to simulate the M2-CAL transient 2580 experiment based on the reported pre-transient and transient conditions. Several transient observables were calculated and compared to the experiment. The model shows good agreement with the experimental power traces, period, and average increase of power at the power ramp (time > 7.8 s). Inverse point kinetics analysis was performed during the period of the power ramp for the model and the experiment based on that, the temperature feedback component was isolated and contrasted with the experimental feedback. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1997

On the Relation between Stratospheric Chlorine/Bromine Loading and Short-Lived Tropospheric Source Gases

Current methods for estimating the concentrations of inorganic chlorine/bromine species Cl(y)/Br(y) in the stratosphere due to decomposition of tropospheric source gases assume that the Cl(y)/Br(y) concentration in the stratosphere is determined mainly by the balance between production from in situ oxidation of the source gases in the stratosphere and removal by transport of Cl(y)/Br(y) out of the stratosphere. The rationale being that for source gases whose lifetimes are of the order of several months or longer the concentration of Cl(y)/Br(y) in the troposphere is small because they are produced at a relatively slow rate and also removed efficiently by washout processes. As a result of the small concentration, the rate at which Cl(y)/Br(y) is transported to the stratosphere is expected to be small compared to the in situ stratospheric production. Thus the transport of Cl(y)/Br(y) from the troposphere contributes little to the stratospheric concentration. In contrast, the origin of stratospheric Cl(y)/Br(y) from reactive source gases with tropospheric lifetimes comparable to the washout lifetime of Cl(y)/Br(y) (of the order of 10-30 days) in the troposphere is distinctly different. The in situ source in the stratosphere is expected to be significantly smaller because only a small portion of the source gas is expected to survive the troposphere to be transported into this region. At the same time these short-lived source gases produce appreciable amounts of Cl(y)/Br(y) in the troposphere such that transport to the stratosphere offers a larger source for stratospheric Cl(y)/Br(y) than in situ production. Thus, for reactive source species, simple methods of estimating the concentration of stratospheric Cl(y)/Br(y) that ignore the tropospheric contribution will seriously underestimate the loading. Therefore estimation of the stratospheric Cl(y)/Br(y) loading requires not only measurements of tropospheric source gases but also measurements of Cl(y)/Br(y) at the tropopause. This paper illustrates the mechanism by using results from a two-dimensional chemistry-transport model. However, in view of the importance of tropospheric transport on stratospheric loading the detailed values should be further evaluated using a three-dimensional model with appropriate treatment of convective transport.

Ko, Malcolm K. W.↗
DOE OSTI2024

Pebble Tanker Model for Nuclear Criticality Safety Needs

This report documents a study performed to investigate the requirements for criticality safety benchmark experiments for high-assay, low-enriched uranium (HALEU) fuel in transportation applications. In this work, an exploratory application model, the “Pebble Tanker,” was developed to represent TRISO fuel in a transportation scenario for an analysis of the validation basis in industrial quantities. An aspect of the criticality validation process involves assessing the “similarity” between application and experimental benchmark systems through an integral index parameter evaluation. Here, this includes propagating nuclear data uncertainties and calculating a correlation coefficient (hereinafter referred to as “c k ”) to evaluate the similarity of benchmark experiments compared with the application Pebble Tanker model. Finding sufficient critical benchmark experiments allows for the evaluation of bias and bias uncertainty, thus determining the upper subcritical limit (USL) of the transportation package. A target k eff of ~0.94 was used in this work to establish appropriate modeling conditions, reflecting a reasonable estimate for a USL. Two container models were investigated: one with the Hermes-type pebble and one with the Pebble Bed Modular Reactor (PBMR)–type pebble. The models were simplified, considering only fuel, containment structure, and either water or air. This allows a focus on the underlying physics of applications involving TRISO fuel pebbles using the Pebble Tanker model. A crucial consideration is the transport package's ability to safely hold pebbles while flooded, maintaining subcritical conditions. Tools available in the SCALE 6.3.1 suite—the CSAS6-Shift, TSUNAMI-3D-Shift, and TSUNAMI-IP sequences—were employed for neutronics and sensitivity and uncertainty (S/U) analysis of the Pebble Tanker. Findings demonstrated sufficient available critical experiment benchmarks to perform a validation of the Pebble Tanker in the most reactive state, i.e., when the Tanker is flooded.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2015

Quantifying Sources and Sinks of Reactive Gases in the Lower Atmosphere Using Airborne Flux Observations

Atmospheric composition is governed by the interplay of emissions, chemistry, deposition, and transport. Substantial questions surround each of these processes, especially in forested environments with strong biogenic emissions. Utilizing aircraft observations acquired over a forest in the southeast U.S., we calculate eddy covariance fluxes for a suite of reactive gases and apply the synergistic information derived from this analysis to quantify emission and deposition fluxes, oxidant concentrations, aerosol uptake coefficients, and other key parameters. Evaluation of results against state-of-the-science models and parameterizations provides insight into our current understanding of this system and frames future observational priorities. As a near-direct measurement of fundamental process rates, airborne fluxes offer a new tool to improve biogenic and anthropogenic emissions inventories, photochemical mechanisms, and deposition parameterizations.

emission↗
DOE OSTI2023

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↗