Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional expectations”

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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Application of a Null-Space Monte Carlo Flow Model Set to the Composite Analysis Base Case Fate and Transport Modeling

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the Composite Analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3, documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application. The overall objective of the saturated zone modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within and downgradient of the Central Plateau of the Hanford Site. Specifically, the purpose of this environmental calculation is to describe the application of the hydraulic property fields and recharge parameters documented in ECF-HANFORD-20-0027, Null Space Monte Carlo Evaluation of the Plateau to River Model, to the CA flow and fate and transport simulation results to quantify the uncertainty in the simulated results due to input parameter selection. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, models lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. The result of NSCM analysis is a set of F&T simulations that provide an estimate of the range of possible outcomes that are used to quantify the uncertainty in simulated concentrations produced using the base case simulations. The simulated concentrations from all simulations will support calculation of the uncertainty of the total dose calculated in a separate calculation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗
DOE OSTI2026

A mechanistic, multiscale model for predicting Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles use silicon carbide (SiC) as the primary structural member and barrier against metallic fission product (FP) release. palladiums (PDs), produced by fission in the fuel kernel, can diffuse to and chemically interact with the SiC layer, degrading its structural integrity and ability to contain radioactive FPs. Existing temperature-dependent correlations for predicting Pd penetration rely on experimental data with significant scatter due to varying conditions, potentially complicating ongoing fuel qualification and licensing efforts for advanced reactors that would subject TRISO fuels to operating conditions outside of those examined in the experiments. A mechanistic model of Pd production, transport, and reaction is developed in this work to better understand and predict PD attack of SiC in TRISO particles. molecular dynamicss (MDs) simulations are utilized to calculate the Pd diffusivity in SiC grain bulk and grain boundaries. A mesoscale phase-field diffusion model, informed by the MD diffusivities, is used to develop a reduced order model (ROM) for the effect of SiC microstructure on the Pd penetration rate. The engineering scale BISON model calculates Pd production and transport, utilizing the ROM to predict penetration rates consistent with experimental data. This novel mechanistic ROM captures the effect of temperature, microstructure, and irradiation history on Pd penetration. In conclusion, these new capabilities are expected to support ongoing qualification and licensing efforts associated with near-term TRISO-fueled reactor applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Machine Learning Derived Dynamic Operating Reserve Requirements in High-Renewable Power Systems

Accurately forecasting wind and solar power output poses challenges for deeply decarbonized electricity systems. Grid operators must commit resources to provide reserves to ensure reliable operations in the face of forecast errors, a process which can increase fuel consumption and emissions. We apply neural network-based machine learning to expand the usefulness of median point forecast data by creating probabilistic distributions of short-term uncertainty in demand, wind, and solar forecasts that adapt to prevailing grid conditions. Machine learning derived estimates of forecast errors compare favorably to estimates based on incumbent methods. Reserves derived from machine learning are usually smaller than values derived using incumbent methods, which enables fuel savings during most hours. Machine learning reserves are generally larger than incumbent reserves during times of higher forecast error, potentially improving system reliability. Performance is tested using multi-stage production simulation modeling of the California Independent System Operator (CAISO) system. Machine learning reserves provide production cost and greenhouse gas (GHG) emission reductions of approximately 0.3% relative to historical 2019 requirements. Savings in the 2030 timeframe are highly dependent on battery storage capacity. At lower levels of battery capacity, savings of 0.4% from machine learning reserves are shown. Significant quantities of battery storage are expected to be added to meet California's resource adequacy needs and GHG reduction targets. Addition of these batteries saturate reserve needs and results in minimal within-hour balancing costs in 2030.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Opposing Hysteresis Patterns in Flow and Outflow Macroscopic Fundamental Diagrams and Their Implications

Two key aggregated traffic models are the relationship between average network flow and density (known as the network or flow macroscopic fundamental diagram [flow-MFD]) and the relationship between trip completion and density (known as network exit function or the outflow-MFD [o-FMD]). The flow- and o-MFDs have been shown to be related by average network length and average trip distance under steady-state conditions. However, recent studies have demonstrated that these two relationships might have different patterns when traffic conditions are allowed to vary: the flow-MFD exhibits a clockwise hysteresis loop, while the o-MFD exhibits a counter-clockwise loop. One recent study attributes this behavior to the presence of bottlenecks within the network. The present paper demonstrates that this phenomenon may arise even without bottlenecks present and offers an alternative, but more general, explanation for these findings: a vehicle’s entire trip contributes to a network’s average flow, while only its end contributes to the trip completion rate. This lag can also be exaggerated by trips with different lengths, and it can lead to other patterns in the o-MFD such as figure-eight patterns. A simple arterial example is used to demonstrate this explanation and reveal the expected patterns, and they are also identified in real networks using empirical data. Then, simulations of a congestible ring network are used to unveil features that might increase or diminish the differences between the flow- and o-MFDs. Finally, more realistic simulations are used to confirm that these behaviors arise in real networks.

42 ENGINEERING↗
DOE OSTI2021

Two-Phase Turbulence Statistics from High Fidelity Dispersed Droplet Flow Simulations in a Pressurized Water Reactor (PWR) Sub-Channel with Mixing Vanes

In the dispersed flow film boiling regime (DFFB), which exists under post-LOCA (loss-of-coolant accident) conditions in pressurized water reactors (PWRs), there is a complex interplay between droplet dynamics and turbulence in the surrounding steam. Experiments have accredited particular significance to droplet collision with the spacer-grids and mixing vane structures and their consequent positive feedback to the heat transfer recorded in the immediate downstream vicinity. Enabled by high-performance computing (HPC) systems and a massively parallel finite element-based flow solver—PHASTA (Parallel Hierarchic Adaptive Stabilized Transient Analysis)—this work presents high fidelity interface capturing, two-phase, adiabatic simulations in a PWR sub-channel with spacer grids and mixing vanes. Selected flow conditions for the simulations are informed by the experimental data found in the literature, including the steam Reynolds number and collision Weber number (Wec={40,80}), and are characteristic of the DFFB regime. Data were collected from the simulations at an unprecedented resolution, which provides detailed insights into the continuous phase turbulence statistics, highlighting the effects of the presence of droplets and the comparative effect of different Weber numbers on turbulence in the surrounding steam. Further, axial evolution of droplet dynamics was analyzed through cross-sectionally averaged quantities, including droplet volume, surface area and Sauter mean diameter (SMD). The downstream SMD values agree well with the existing empirical correlations for the selected range of Wec. The high-resolution data repository from the simulations herein is expected to be of significance to guide model development for system-level thermal hydraulic codes.

Saini, Nadish↗
DOE OSTI2023

Understanding of Low‐Porosity Sulfur Electrode for High‐Energy Lithium–Sulfur Batteries

Abstract The lithium–sulfur (Li–S) battery is a promising technology for large‐scale energy storage and vehicle electrification due to its high theoretical energy density and low cost. Reducing the sulfur cathode porosity has been identified recently as a viable strategy for improving the cell practical energy density and minimizing pore‐filling electrolytes to extend cell life at lean electrolyte conditions. Direct use of a low‐porosity cathode for Li–S battery results in poor electrode wetting, nonuniform electrode reactions, and thus early cell failure. To understand and mitigate the barriers associated with the use of low‐porosity electrodes, multiscale modeling is performed to predict electrode wetting, electrolyte diffusion, and their impacts on sulfur reactions in Li–S cells by explicitly considering the electrode wettability impacts and electrode morphologies. The study elucidates the critical impact of low tortuosity and large channel pore design for promoting electrode wetting and species diffusion. It is suggested that the secondary particle size should be comparable with the electrode thickness to effectively promote electrolyte wettability and sulfur reactivity. This study provides new insights into the low‐porosity electrode material and designs and is expected to accelerate the development of practical high‐energy Li–S batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Dynamic changes in evapotranspiration, canopy photosynthesis and expansive growth in open field under rapid fluctuating radiation

Abstract In open fields, plants experience dynamic changes in environment, particularly radiation, temperature, wind, and humidity but their short-term responses have not been adequately characterized under natural conditions. In this study, we assessed causal effects of rapid radiation fluctuation on seven plant parameters in open fields of cotton and sweet corn canopies. The parameters are evapotranspiration ( ET ), canopy photosynthesis indicated by CO2 flux from air to canopy ( FCO2 ), sensible heat energy flux ( H ), canopy conductance for water vapor ( gcw ), canopy surface temperature ( Ts ), shoot and leaf elongation rate, and stem diameter change. The energy and CO2 fluxes were measured with Bowen ratio/energy balance/CO2 gradient (BREB+) technique, using averaging time of 5 min. Shoot + leaf elongation and stem diameter change were monitored with position transducers using averaging time of 1 min. All parameters were all found to respond to change in radiation and transpiration within minutes or sooner. While radiation effects on canopy gas exchanges are expected, illuminating are the indirect but immediate effects on shoot + leaf growth and stem diameter change through radiation effects on transpiration and plant water status. A novel finding is that gcw also responded within minutes or sooner to radiation fluctuations and that FCO2 was related almost linearly to gcw . Results are discussed in terms of soil-plant-atmosphere continuum, and interpreted in terms of dynamic interactions between transpiration and plant water status. The clear inverse relationship between ET and elongation rate or stem diameter changes provides additional evidence supporting the validity of 5-min averaging for the BREB + technique.

Hsiao, Theodore C.↗
DOE OSTI2022

Kinetic Model for the Reduction of Cu II Sites by NO + NH 3 and Reoxidation of NH 3 -Solvated Cu I Sites by O 2 and NO in Cu-SSZ-13

In this work, a kinetic model is developed for the reduction of Cu II sites by NO + NH 3 and the reoxidation of NH 3 -solvated Cu I sites by O 2 and NO in Cu-SSZ-13. Fourier transform infrared (FTIR) spectroscopy and spatially resolved capillary inlet mass spectrometry (SpaciMS) measurements during transient reactor experiments are utilized to identify the rate parameters associated with NO + NH 3 RHC (reduction half-cycle), proposed to occur via two distinct pathways involving adsorbed NH 3 and gas-phase NH 3 . The resulting NO + NH 3 RHC model is validated using spatiotemporal N 2 measurements covering a wide range of temperatures (200–450 °C) and space velocities (53,000–640,000 h –1 ). N 2 O formation is observed and modeled during NO + NH 3 RHC, with quantitative validation under standard selective catalytic reduction (SCR) conditions. Experimentally measured enthalpic and entropic changes associated with O 2 adsorption on NH 3 -solvated Cu I (ZCu(NH 3 ) 2 ) complexes [ Kamasamudram, K. Catal. Today 2010 , 151 (3–4), 212-222], along with activation energies estimated computationally for the intercage diffusion of ZCu(NH 3 ) 2 complexes [ Paolucci, C. Science 2017 , 357 (6 354), 898-903], are incorporated into a mean field kinetic model for the low-temperature oxidation half-cycle (OHC). Significant NH 3 release is observed during the isothermal oxidation of Cu I sites, attributed to desorption of NH 3 ligands from NH 3 -solvated Cu II dimers (Z 2 Cu 2 (NH 3 ) 4 O 2 ). Reduction of these dimeric complexes leads to the consumption of one NO/Cu II , contradicting the expected reduction stoichiometry. Inclusion of a global Arrhenius rate for the NO titration of Z 2 Cu 2 (NH 3 ) 4 O 2 complexes provides accurate representations of standard SCR on reduced and oxidized catalysts, predicting transient NO and NH 3 consumption between 150 and 250 °C as a function of hydrothermal aging. Deactivation of low-temperature standard SCR by NH 3 is observed at high NH 3 pressures, modeled via the formation of superoxo amino (ZCu(NH 3 ) 3 OO*) complexes during NH 3 titration of Z 2 Cu 2 (NH 3 ) 4 O 2 complexes [ Negri, C. J. Am. Chem. Soc. 2020 , 142 (37), 15884-15896]. The redox kinetic model presented here provides a foundational description of active site redox during low-temperature standard SCR, combining the recent kinetic, spectroscopic, and computational findings on the mechanism of standard SCR over Cu-SSZ-13.

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

Stability of the fcc phase in shocked nickel up to 332 GPa

Despite making up 5-20 wt.% of Earth’s predominantly iron core, the melting properties of elemental nickel at core conditions remain poorly understood, due largely to a dearth of experimental data. We present here an in situ X-ray diffraction study performed on laser shock-compressed samples of bulk nickel, reaching pressures up to ~ 500 GPa. Hugoniot states of nickel were targeted using a flat-top laser drive, with in situ X-ray diffraction data collected using the Linac Coherent Light Source. Rietveld methods were used to determine the densities of the shocked states from the measured diffraction data, while peak pressures were determined using a combination of measured particle velocities, shock transit times, hydrodynamic simulations, and laser intensity calibrations. We observed solid compressed face-centered cubic (fcc) Ni up to at least 332 ± 30 GPa along the Hugoniot—significantly higher than expected from the majority of melt lines that have been proposed for nickel. We also bracket the partial melting onset to between 377 ± 38 GPa and 486 ± 35 GPa.

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

A practicable estimation of opening angle of dust torus in Type-1.9 AGN with double-peaked broad Hα

In this manuscript, an independent method is proposed to estimate opening angle of dust torus in AGN, through unique properties of Type-1.9 AGN with double-peaked broad Hα (Type-1.9 DPAGN) coming from central accretion disc. Type-1.9 AGN without broad Hβ can be expected by the commonly accepted unified model of AGN, considering central BLRs seriously obscured by dust torus with its upper boundary in the line of sight. For the unique Type-1.9 DPAGN, accretion disc originations of double-peaked broad Hα can be applied to determine the inclination angle of the central accretion disc, which is well accepted as substitute of the half opening angle of the central dust torus. Then, among low redshift Type-1.9 DPAGN in SDSS, SDSS J1607+3319 at redshift 0.063 is collected, and the half opening angle of the central dust torus is determined to be around 46 ± 4°, after considering the disfavoured BBH system to explain the double-peaked broad Hα through long-term none variabilities and disfavoured local physical conditions to explain the disappearance of broad Hβ through virial BH mass properties. The results indicate that more detailed studying on dust torus of AGN can be appropriately done through Type-1.9 DPAGN in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

20 years of the CEA/DAM NNSA/DP Agreement

For the past twenty years, there has been a very active, productive International Agreement between France and the United States of America for Cooperation on Fundamental Science supporting Stockpile Stewardship. Under this Agreement the scientists at the nuclear weapons laboratories in both countries have collaborated on many unclassified research projects in areas such as Materials in Extreme Conditions, Nuclear Physics, and Atomic and Plasma Physics. The results of their efforts have not only been published in the open literature but have enhanced the physics base of the computer codes essential to the mission of ensuring that the nuclear stockpiles are safe, secure and effective. Work on these collaborations is extremely important. Not only does such joint work bring more brilliant minds to work on pressing research problems for both countries but also the collaborative effort sharpens the scientific skills of and presents scientific challenges to the scientific and technical staff of the laboratories. As we reach the 20th Anniversary of the formal signing of the Agreement, we feel it is important to thank all of the individuals who have contributed to its continuing success. Our expectations for brilliant, exciting, challenging joint research projects under this Agreement are at an all-time high. We toast the achievement of this milestone and look forward to the research results to come.

36 MATERIALS SCIENCE↗
DOE OSTI2025

ALARA Review for Venting of Flanged Tritium Waste Containers (FTWCs) at TA-54

This project will vent headspace hydrogen and oxygen from specialized high-pressure storage vessels called flanged tritium waste containers (FTWCs). There are four of these containers located in Los Alamos National Laboratory’s (LANL’s) Technical Area 54 (TA-54), Building 1028. This building is in the southwest corner of LANL’s Material Disposal Area G. The vented headspace gas is expected to contain tritium in the form of water vapor, elemental hydrogen or a combination of both. The venting operation is expected to take place in 2025. There will be a series of operational readiness reviews prior to venting activities commencing. The FTWCs at TA-54 contain tritium-contaminated metal parts and molecular sieve media, which is a pebble-like material used to absorb water vapor from the air. This molecular sieve media inside the FTWCs is contained in metal canisters, along with some loose media material in bags. Over time, tritiated water vapor that had been adsorbed onto the media can become liberated into the FTWC headspace gas. Radiolysis can cause separation of the water vapor into its hydrogen and oxygen components, resulting in the potentially hazardous gas mixture within the FTWC. LANL has determined that continued tritium storage in these containers can pose an unsafe condition due to possible hydrogen and oxygen gas buildup within the FTWCs with a potential for explosion if sparks are generated within the FTWCs. To mitigate this hazard, the FTWCs will be stabilized by venting them in-place to remove hazardous gases.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Development of Multiresolution Capabilities for the Holistic Energy Resource Optimization Network (HERON) tool A progress update

INL researchers work on technoeconomic analyses for integrated energy systems (IES) using the Framework for Optimization of ResourCes and Economics (FORCE). Within FORCE, researchers use the Holistic Energy Resource Optimization Network (HERON) tool to conduct optimization of grid portfolios under uncertain market conditions. These optimizations determine optimal capacities for all IES components and strategies for resource dispatch which maximize some economic metric (e.g., net present value). Resource dispatch occurs on finer timescales (typically hours) and thus are asked to respond to a given time series (e.g. hourly load demand profiles for a grid, or pre-determined electricity prices). Volatile and complex bidding dynamics as well as poorly forecasted weather events within deregulated markets add uncertainty to the time series; FORCE can address this uncertainty by training a reduced order model on historical time series and generate unique synthetic time series which represent individual scenarios or realizations of the market. The IES configuration can be simulated under these different sampled realizations and a stochastic optimization is conducted which optimizes the expected value of the desired economic metric. The training of a synthetic time series generator is limited by the chosen time resolution; dynamics can occur on different time scales. Seasonal demand trends can dominate faster dynamical events (such as power outages from certain sectors or severe weather events) which might not get captured correctly by the trained model. In this report, we investigate different ways of addressing the training and generation of time series on multiple time scales using three main algorithms: wavelet decomposition, dynamic mode decomposition, and generative adversarial networks for time series. We demonstrate a time series analysis that yields information on not just the frequency space but also temporal space: where a fast Fourier transform can provide what frequencies dominate, the new algorithms can provide when the frequencies dominate as well. These analyses can help improve IES optimization by allowing researchers to couple simulations at different timescales when it is most needed - seasonal, day-ahead, and real time optimization - with greater computational efficiency. Future work will include implementation of a subset of the proposed algorithms into the FORCE toolset and application of these analyses into multiple timescale optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Assessing mass balance-based inverse modeling methods via a pseudo-observation test to constrain NO x emissions over South Korea

This study constrains NO x emissions by utilizing mass balance–based inverse modeling methods and examines the optimization of inverse modeling conditions through a pseudo-observation test for South Korea, which has complex topography. Here, we applied mass balance–based inverse modeling methods such as Basic Mass Balance (BMB), Finite Difference Mass Balance (FDMB), and Iterative Finite Difference Mass Balance (IFDMB). We performed a numerical simulation of air quality using the Community Multi-scale Air Quality (CMAQ) model to calculate the NO 2 column density required for the inverse modeling. The pseudo-observation test was performed according to season, modeling resolution, and regridding methodology of satellite observation data to identify various conditions while applying the inverse modeling for South Korea. Comparing the inverse modeling results from the BMB, FDMB, and IFDMB methods, IFDMB was the most effective method in constraining NO x emissions in the South Korean region since it minimized smearing effects (i.e., transport-induced errors) through iterative calculations. The accuracy of the constrained NO x emissions using mass balance–based inversions in South Korea was the highest in the summer due to the minimized smearing effects, and the 9 km resolution modeling was the most efficient for inverse modeling of the region. In addition, the results of inverse modeling varied depending on the regridding methods, implying the importance of using a suitable regridding method and modeling resolution. This study used pseudo-observations, but the inversions are expected to be applied based on actual satellite data in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The impact of metastability on the high-pressure behavior of cerium

The structures adopted by solids under pressure are often assumed to reflect thermodynamic equilibrium, yet in many materials phase selection is strongly influenced by kinetic pathways and microstructural inheritance. Elemental cerium (Ce) exemplifies this challenge, with decades of conflicting reports describing different high-pressure crystal structures emerging under nominally identical conditions. Here we use neutron diffraction from large ( ~ 60 mm 3 ) sample volumes to follow the structural evolution in ultra-high-purity Ce during controlled pressure-temperature cycling between 85 and 295 K and up to 8 GPa. We find that the crystal structure formed at high pressure depends on the compression pathway: slow compression ( ~ 0.25 GPa hr −1 ) at room temperature favors an orthorhombic phase (α'), whereas slow ( ~ 0.25 GPa hr −1 ) and also moderately faster ( ~ 0.5 GPa hr −1 ) compression at low temperature stabilizes a pure monoclinic phase (α"). The low-temperature phase persists metastably over a wide temperature range but transforms irreversibly upon heating above ~ 280 K, or modest pressure cycling. Remarkably, the lower-pressure γ phase remains trapped far beyond the expected stability field, persisting to the highest pressures studied. These observations show that phase selection in Ce is governed by kinetics and microstructural memory rather than equilibrium thermodynamics or sample morphology alone, establishing path dependence as a defining feature of its high-pressure behavior.

Ridley, Christopher J. [Oak Ridge National Laborat↗
DOE OSTI2021

Hydrogenous Content in Residues from Aqueous Chloride Pu and Am Recovery Operations

Byproduct residues generated from Aqueous Chloride Pu and Am recovery operations at Los Alamos National Laboratory have been investigated and determined to contain high weight % chloride salts and minimal hydrogen content. Hydrogen content is the result of limited surface area adsorption of ambient moisture. The assumption used in technical analysis NCS-TECH-19-029 to model potentially hydrogenous items as mixtures of polyethylene is shown to be overly conservative. Our results and literature precedent indicate that water adsorption of 1-3 % by weight is expected for hydroxide cakes and Pu oxides, and a conservative bounding assumption for these materials is 5% water by weight. LOH results for a dissolution heel stored long term showed 6.25% weight loss, which was conservatively assumed to be water. A conservative bounding assumption for dissolution heels is 8 wt.% water content. This equates to 0.55 % hydrogen by weight for hydroxide cakes, and 0.88% hydrogen by weight for dissolution heels. The bulk composition of dissolution heels reflects the composition of the original byproduct feed material, chloride salts of potassium, sodium, magnesium and calcium. The bulk composition of calcined hydroxide precipitation residues is composed of potassium chloride as a result of potassium hydroxide neutralization of hydrochloric acid solutions. Metal hydroxides generated during the hydroxide precipitation process are decomposed under the conditions of calcination to metal oxide forms.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Removing numerical pathologies in a turbulence parameterization through convergence testing

Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e. clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, this noise can be amplified by an advective term in CLUBB's background diffusion. Smoothing the limiters and removing the advective component of the background diffusion reduces the noise and restores the expected first-order convergence in CLUBB's solutions. These model reformulations improve the results at coarser, near-operational grid spacing and time step in cumulus cloud and dry turbulence tests. In addition, convergence testing is proved to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

Mathematics and Computing, Geosciences↗