SEARCH · Engineering Papers

Results for “Source term”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,585 records · Page 88

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Cell Free and Immobilization Technologies (CFIT)

Today, several key factors negatively impact the production of fuels and chemicals from renewable sources. Common hindrances in the biological production of biochemicals are: (1) end-product or intermediate toxicity to the microbial biocatalyst, (2) the diversion of carbon to biomass formation, and (3) co-production of undesired byproducts. A particularly attractive alternative is to eliminate the biocatalyst entirely and instead operate the desired metabolic pathways in isolation, thus circumventing the roadblocks of biological toxicity, lower yields, and lack of specificity. However, cell-free enzyme systems still suffer from low productivities owing in part to the effects of free diffusion of intermediates, lack of long term enzyme stability, cofactor cost or inefficient recycling rates, and finally, the cost of enzyme production/purification. This project represents a new effort to propose innovative and cost competitive routes to producing biochemicals from a variety of feedstocks using cell free approaches. These routes will help reduce the current risk and cost associated with classical cell free production. Cell free technologies show promise for application to the production of toxic/inhibitory products or products difficult to separate from microbial growth media and can help reduce the production barriers in multiple areas of biological conversion of feedstocks to biochemicals. More specifically, we are developing new technologies and routes that could be used to produce high value biochemicals such as 2,3 BDO and terpenes (among many others) from biomass derived C5/C6 sugars or lignin but also from waste byproducts such as glycerol. This project will lead to significant innovation and also lead to new concepts and rational design of pathways and enzymes. Within this project, we will develop new metabolic enzyme cascades that will represent natural or artificial combinations of enzymes to produce the desired biochemicals from a variety of feedstocks. We are also developing basic design principals for constructing synthetic metabolons, using fusion proteins and synthetic protein scaffolds, to promote substrate channeling and stability while conserving peak activity. Additionally, our efforts include a techno economic analysis (TEA) of cell free approaches to provide the sensitivities of the process to enzyme loading, activity, pH, reactor volumes, cofactor recycling rates. Finally, we are focusing on further increasing stability, operating lifetime and efficiency of the pathway enzymes by immobilization on support surfaces. We are also focus on immobilizing pathway enzymes or combinations of enzymes on several different conducting polymers and evaluate the effect on stability and operating lifetime. As more combinations become available we will conduct a more systematic study of the means of immobilizing these enzymes. This preliminary work will enable the in-depth study of cofactor recycling at these interfaces using mediators for electron transfer. Taken together, these approaches will enable process intensification, continuous operation, lower capital and separations costs and end-product flexibility, and thus has the potential to contribute significantly to BETO's goals of cost competitive biofuels and bioproducts. To date we have demonstrated 1) the conversion of pyruvate to 2,3 BDO (4 enzymes with cofactor recycling, no additional cofactors needed) at >100g/L (>3g/L/h) using our enzyme tethering approach, 2) Generated mutants of a key redox enzyme with >180- fold improvement in NADH utilization over WT to >85% of NADPH utilization. 3) Produced mevalonate at >10g/L from glucose, limonene at >10g/L from mevalonate and more than >5g/L from glucose with crude enzyme preparations and complete cofactor recycling. 4) successfully generated fully active cross-linked aggregates (CLEAs) of NOX, a cofactor-regenerating water-forming oxidase that is key to our process with increased stability 5) Engineered and identified enzymes able to use the biomimetic cofactor NMN and synthesized new synthetic cofactors.

biocatalysis↗

DOE OSTI2023

Geological Thermal Energy Storage Using Solar Thermal and Carnot Batteries: Techno-Economic Analysis: Preprint

Energy storage is increasingly necessary as Variable Renewable Energy (VRE) technologies are deployed. Seasonal energy storage can shift energy generation from the summer to the winter, but these technologies must have extremely large energy capacities and low costs. Geological Thermal energy storage (GeoTES) is proposed as a solution for long-term energy storage. Excess thermal energy can be stored in permeable reservoirs such as aquifers and depleted hydrocarbon reservoirs for several months. In this article, we describe a techno-economic model that has been developed to evaluate GeoTES systems. The models are developed by combining the output of specialist models which enables the performance and cost of both the subsurface and surface systems to be captured. Off-design models are developed so that the performance can be evaluated at each hour of the year. GeoTES can be charged with two different energy sources: (1) concentrating solar thermal and (2) renewable electricity using heat pumps (henceforth known as a "Carnot Battery"). The stored thermal energy can be used to generate electricity and - uniquely - also directly produce heat that can be used by industrial processes. Furthermore, Carnot-Battery-GeoTES can also be used to form a cold storage reservoir. Preliminary results that quantify the technical and economic performance of these two GeoTES systems are presented.

carnot battery↗

DOE OSTI2023

Geological Thermal Energy Storage Using Solar Thermal and Carnot Batteries: Techno-Economic Analysis

Energy storage is increasingly necessary as variable renewable energy technologies are deployed. Seasonal energy storage can shift energy generation from the summer to the winter, but these technologies must have extremely large energy capacities and low costs. Geological thermal energy storage (GeoTES) is proposed as a solution for long-term energy storage. Excess thermal energy can be stored in permeable reservoirs such as aquifers and depleted hydrocarbon reservoirs for several months. In this article, we describe a techno-economic model that has been developed to evaluate GeoTES systems. The models are developed by combining the output of specialist models, which enables the performance and cost of both the subsurface and surface systems to be captured. Off-design models are developed so that the performance can be evaluated at each hour of the year. GeoTES can be charged with two different energy sources: (1) concentrating solar thermal and (2) renewable electricity using heat pumps (henceforth known as a "Carnot Battery"). The stored thermal energy can be used to generate electricity and, uniquely, also directly produce heat that can be used by industrial processes. Furthermore, Carnot Battery GeoTES can also be used to form a cold storage reservoir. Preliminary results that quantify the technical and economic performance of these two GeoTES systems are presented.

carnot battery↗

DOE OSTI2026

Measurement of $\nu_\mu/\bar\nu_\mu$ CC double-differential cross sections on MINERvA hydrocarbon target for the shallow inelastic scattering background region

Cross section measurements are essential for all neutrino oscillation experiments. In fact, uncertaintiesassociated to cross section model parameters constitute one of the dominant sources oferrors in current oscillation analyses. In particular, understanding neutrino-induced pion productionin the kinematic regime known as shallow inelastic scattering (SIS) is critical for improvingneutrino interaction modeling in event generators. In this study, 416,233 (237,468) muon neutrino(antineutrino) interactions are measured in a SIS background region, predominantly made ofbaryon resonances. The analyzed datasets were collected from 2013 to 2019, comprising neutrinosgenerated by the Fermilab NuMI facility, with mean energy of 6 GeV, and the interactions occurredon the MINERvA hydrocarbon target. The measurements are presented as double-differential crosssections in terms of the outgoing muon longitudinal and transverse momentum components, aswell as the Bjorken x and y variables. Comparisons between the extracted data and predictionsfrom several generators reveal significant discrepancies across most kinematic bins.

Souza Correia, Souza Correia, Daniel [Rio de Janei↗

DOE OSTI2021

Fatigue analysis of the Spallation Neutron Source 2 MW target design

Upgrades to the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory are under way with two major projects. The first project is the Proton Power Upgrade (PPU), which will double the power capacity of the SNS accelerator system to enable the operation of a future second neutron source target station—the second project. PPU also will increase power to the SNS first target station – which currently operates at 1.4 MW with 1.0 GeV protons – to 2 MW with 1.3 GeV protons. Final design of the PPU 2 MW target module is finished, including all necessary design and analysis calculations. The fatigue evaluation described herein is an assessment of the fitness of the PPU 2 MW target mercury vessel to resist failure from cyclic loading from the pulsed 60 Hz proton beam and thermal changes from disrupted operation. In addition to the reliability predictions, the method, detail, and depth of this fatigue evaluation are superior to those of past target design assessments. Furthermore, the thermal fatigue life is predicted to exceed 3125 thermal cycles per 1250 h of operation for all load cases. The design also meets the (relative) fatigue design goal for combined thermal and pulse loading of greater than 0.5 times the minimum fatigue life calculated for a jet-flow target design operating at 1.4 MW—the most robust design to date in terms of resistance to fatigue failure. From the onset of the PPU, it has been known that pulse loading requires effective application of helium gas injection into the target mercury to reduce beam pressure loading. For the PPU 2 MW target, none of the pulse load cases requires a maximum strain reduction from gas injection of more than 50% (for base material) to meet the fatigue design goal. This level of reduction is achievable from gas injection, according to historical strain measurements from 1.4 MW target operation. Overall, the PPU 2 MW target design is predicted to have superior resistance to fatigue failure compared with past target designs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 2. Quantitative Moisture Dynamics Estimation Model

Abstract The long‐term containment of high‐level radioactive waste in geological disposal repositories relies on Engineered Barrier Systems (EBS), with bentonite clay emerging as a candidate material due to its unique properties. Understanding moisture dynamics within bentonite buffers is crucial for EBS performance, as it directly influences the material's swelling capacity, thermal and hydraulic conductivity, mechanical properties, and long‐term evolution under complex thermal‐hydrological‐mechanical (THM) processes. This study develops an advanced Electrical Resistivity Tomography (ERT)‐based framework to quantitatively monitor moisture dynamics under THM conditions. Our framework extends the Waxman‐Smits model to incorporate the coupled effects of temperature, water content, fluid chemistry, and mechanical changes on bentonite's electrical properties. Utilizing HotBENT‐Lab data from our companion paper, which includes electrical conductivity, CT density, and thermocouple measurements, this study offers a novel methodological framework bridging different scales of the model. Our results show that the extended model can estimate water content from ERT data, capturing spatial and temporal variations in moisture distribution within bentonite columns. However, the model tends to overestimate water content compared to CT density‐derived measurements. We address this discrepancy by incorporating a simplified swelling effect model, which improves agreement between ERT and CT density‐based water content estimates. We also discuss model limitations, including simplified treatment of swelling and micropore effects, and propose a conceptual framework for transitioning from laboratory to field applications, addressing challenges such as parameter scalability, field validation methods, and integration of diverse data sources. This ERT‐based framework can potentially advance real‐world moisture monitoring of bentonite‐based EBS in nuclear waste repositories. Plain Language Summary Safely containing high‐level radioactive waste depends on barriers made from materials like bentonite clay, which is effective because it swells and seals in the waste. To ensure these barriers work well over time, it's important to understand how moisture moves through the clay. Our study developed a new method using ERT to monitor moisture levels in bentonite under conditions that mimic those in actual storage sites, including changes in temperature, water content, and mechanical stress. This study improved an existing model to better account for how these factors affect the clay, allowing us to create more accurate moisture maps. Initially, the proposed model overestimated the amount of water in the clay, but its accuracy was improved by factoring in how the clay swells when wet. This study also identified some limitations of the model and suggested ways to adapt it for use in real‐world waste storage sites. This new approach could lead to better monitoring and safety checks for nuclear waste storage systems, helping to ensure long‐term containment. Key Points This work develops an ERT‐based framework extending the Waxman‐Smits model to monitor bentonite moisture dynamics during coupled THM processes The extended model accurately estimates water content from Electrical Resistivity Tomography data, incorporating swelling effects to improve precision This work proposes a conceptual framework for transitioning from laboratory to field applications, advancing EBS monitoring in nuclear waste repositories

Chen, Hang↗

DOE OSTI2023

Versatile tape-drive target for high-repetition-rate laser-driven proton acceleration

We present the development and characterization of a high-stability, multi-material, multi-thickness tape-drive target for laser-driven acceleration at repetition rates of up to 100 Hz. The tape surface position was measured to be stable on the sub-micrometre scale, compatible with the high-numerical aperture focusing geometries required to achieve relativistic intensity interactions with the pulse energy available in current multi-Hz and near-future higher repetition-rate lasers ( > kHz). Long-term drift was characterized at 100 Hz demonstrating suitability for operation over extended periods. The target was continuously operated at up to 5 Hz in a recent experiment for 70,000 shots without intervention by the experimental team, with the exception of tape replacement, producing the largest data-set of relativistically intense laser–solid foil measurements to date. This tape drive provides robust targetry for the generation and study of high-repetition-rate ion beams using next-generation high-power laser systems, also enabling wider applications of laser-driven proton sources.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

GronOR: Scalable and Accelerated Nonorthogonal Configuration Interaction for Molecular Fragment Wave Functions

GronOR is a program package for nonorthogonal configuration interaction calculations. Electronic wave functions are constructed in terms of antisymmetrized products of multiconfiguration molecular fragment wave functions. The computational complexity of the nonorthogonal methodologies implemented in GronOR applied to large molecular assemblies requires a design that takes full advantage of massively parallel supercomputer architectures and accelerator technologies. This work describes the implementation strategy and resulting performance characteristics. In addition to parallelization and acceleration, the software development strategy includes aspects of fault resiliency and heterogeneous computing. The program was designed for large-scale supercomputers but also runs effectively on small clusters and workstations for small molecular systems. GronOR is available as open source to the scientific community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Deciphering the Role of Total Water Storage Anomalies in Mediating Regional Flooding

Regional floods result from various flood generation mechanisms. Traditional analyses mainly link flooding to extreme rainfall, with limited input from soil moisture. Total water storage (TWS) is a holistic measure of basin wetness, including additional storage components from surface water, snow, and groundwater. Utilizing a new 5-day Gravity Recovery and Climate Experiment and its Follow On (GRACE(-FO)) data set, we investigated the linkage between short-term TWS anomaly (TWSA) and regional flooding. The 5-day TWSA solutions revealed flood signals missed by monthly TWSA solutions. Global basins exhibit distinct storage-discharge co-evolution patterns, offering new insights into flood mechanisms and propensity. Our bivariate event analyses show the annual maximum river discharges co-occur more often with the TWSA maxima than with precipitation in many basins. Further analyses revealed TWSA's time-lagged effect on river discharge, particularly in basins susceptible to floods triggered by saturation-excess runoff. The 5-day TWSA provides a new source of information for enhancing global flood preparedness.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Techno-Economic Implications of Electrical Machine Scaling for Wave Energy Converters

The sizing of an electrical machine for a Wave Energy Converter (WEC) can have a substantial impact on the overall sizing, cost, and rating of the device. An electrical generator is typically part of the power take-off system, which is the mechanism by which the energy absorbed by the prime mover is transformed into usable electrical energy. For practically all WECs, the rate of change of actuation is predominantly determined by the wave resource (i.e., the wave height and frequency), and devices will see a sinusoidal varying velocity according to the wave conditions. The same can then be said for both directly and indirectly coupled power take-offs with electrical generators. This techno-economic study investigates electrical machine scaling and associated cost implications through core machine design theory, manufacturer data, supporting literature, and the Reference Model Project sponsored by the U.S. Department of Energy. The Reference Model Project was a partnered effort to develop open-source marine energy point designs as reference models to benchmark marine energy technology performance and costs, methods for design and analysis of marine energy technologies, estimations for capital costs, operational costs, and levelized cost of energy. The results from this study show torque is directly related to (1) the physical size of the machine required to increase the air-gap sheer stresses, (2)the amount of active material, (3) the support structure, (4) bearing size and rating, and (5) offshore cable rating, all of which have a significant effect on overall system costs in terms of both capital and operational expenditures. This paper aims to be a critical benchmark in helping determine an "optimal" nameplate rating for wave energy devices and their associated power take-offs. With an optimized rating and sizing process, WEC costs can be reduced and overall performance can be improved.

cables↗

DOE OSTI2023

Ice Cryo-Encapsulation Balloon (Project ICEBall) Field Campaign Report

The Ice Cryo-Encapsulation Balloon (ICEBall) field campaign was designed to sample the ice crystals that compose high-altitude cirrus with a passive device. The campaign made use of a new instrument, ICEBall, which is a balloon-borne ice crystal sampling system. The ice crystal sounding system is capable of measuring ice crystal concentration, temperature, atmospheric pressure, ice crystal habit, aerosol particle morphology, and residual composition. The 3-kg instrument is carried upwards at 5 m s-1 by a high-altitude balloon. The instrument can be cut down from the balloon at any altitude up to 20 km, and the apparatus returns to the surface by parachute. Basic measurements such as temperature and pressure are recorded onboard, and high-frequency Global Positioning System (GPS) records altitude and latitude/longitude. Ice crystal concentrations are measured through the use of a high-resolution video camera mounted on the device. Ice crystals are collected through an open aperture leading to insulated collection chambers cooled with dry ice. Upon exiting the top of the cloud system, the chamber aperture is closed, and the ~1 mm3 sample cell is magnetically sealed and isolated at -78 °C, ensuring that ice particles do not sublimate or grow after collection. Once the crystals are returned to the surface, they are double-sealed and immersed at liquid nitrogen temperature in “dry-cryo shippers” before being transported back to the laboratory. Dr. Magee’s laboratory at The College of New Jersey contains a cryo-stage scanning electron microscope (SEM), which was used to interrogate the crystals and aerosol particles. The main purpose of this pilot field campaign was to provide an unprecedented level of detail on the crystal habits and ice surface complexity in mid-latitude cirrus, which may help resolve issues associated with habit identification and classification in cirrus. The ICEBall campaign was originally scheduled to run from March 28 to April 18 of 2021 at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains observatory. The COVID-19 pandemic intervened and caused us to shift the dates of the experiment to October 16-November 6 of 2021. This period is also climatologically favorable for cirrus. The approximately six-month gap between our original field campaign dates and the actual dates afforded us the opportunity to build two new ICEBall payload instruments. These instruments were tested during an August 2021 trip to The College of New Jersey. During this field testing phase, we decided to launch the ICEBall payload upstream from the ARM SGP site with the goal of landing in the vicinity of the site. Our goal was to sample the ice crystals before the cirrus were advected over the remote-sensing instruments at the SGP site. The team assembled for the field campaign consisted of the Principle Investigator (PI) and Co-Principle Investigator (Co-PI) (Drs. Harrington and Magee), The Pennsylvania State University research scientist Dr. Alfred Moyle, and two graduate students (Ms. Marley Majetic and Gwenore Pokrifka). The team operated out of a house rented in Enid, Oklahoma. We successfully sampled seven cirrus cloud systems during the three-week field campaign (October 21, 23-26, 31, and November 1). This was a much higher success rate than either of the PIs anticipated (our goal was closer to sampling three or four cases). The balloon was typically launched from oil pads or farm fields northwest of Enid and the payload was typically retrieved somewhat north of the SGP site. We never landed directly at the SGP site, and so did not need regular access to the SGP facilities. Our greatest concern going into the field campaign was the longer-term storage of crystals in the -196°C cryo dry-shipper dewars and the subsequent transport across the country. We had tested storage and transport prior to the field campaign, but we had never stored crystals for a few weeks nor had we transported the dewars over long distances. To our great relief, the storage and transport worked flawlessly and we were able to image a large number of crystals from six of the seven cases. Working with the staff at the ARM SGP office was excellent. They not only helped us find the sources we needed for helium, liquid nitrogen, and other materials, but also helped with contacts within the Federal Aviation Administration (FAA) and Vance Air Force Base. One goal of our field project was to tie the in situ measurements of ice crystal habits to the radar signatures derived from the Ka-band ARM Zenith-pointing Radar (KAZR). Unfortunately KAZR was down for the duration of our experiment. However, the Ka-band Scanning ARM Cloud Radar (KASACR) was put into vertically pointing mode during the ICEBall campaign and those data, along with Doppler lidar measurements, have proved very useful.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Soil management legacy interacts with wheat genotype to determine access to organic N in a dryland system

Organic nutrient management through the application of compost and/or cover crops provides mineralizable sources of nutrients for plants while often building soil organic matter (SOM) and various aspects of soil health. Variability in nutrient acquisition strategies between crop genotypes may confer advantages under different soil health contexts and could be important for crop selection and breeding, but crop response under field conditions remains unexplored. We investigated the ability of different genotypes of winter wheat (Triticum aestivum L.) to access nitrogen (N) from newly added cover crop residues in two soils with contrasting levels of SOM and biological activity. We planted three previously characterized wheat genotypes in a long-term dryland compost amendment field trial: 1) Byrd (modern, deep roots, low exudation), 2) Cheyenne (historic, drought susceptible, intermediate exudation), and 3) Snowmass (modern, drought-susceptible, high exudation). 15 N-labelled cover crop residue was added to each plot and traced into wheat tissue. In the low SOM soil, the high exudate genotype Snowmass and historic genotype Cheyenne took up the most residue-derived N (6.4–8.1 kg N ha −1 ) compared to the low-exudate genotype Byrd (4.4 kg N ha −1 ), suggesting a strong exudate effect in the more carbon-limited soil. However, in the high SOM soil, the low-exudate, deep rooted genotype, Byrd, took up the most residue N (4.6 kg N ha −1 vs. 2.8 and 3.3 hg N ha −1 for Cheyenne and Snowmass, respectively), which indicated higher native N cycling activities and greater importance of drought resistance. Enzyme activity, inorganic N, and microbial communities were not influenced by genotype, though did show strong effects of compost application legacy. Furthermore, our results show that belowground allocation strategies that favor microbial stimulation may be less successful under water limitation, especially when high native SOM and biological activity can support mineralization of residue N without added investment in root inputs. Increased soil health through SOM-building management likely enhances nutrient cycling, and may better support root strategies that invest less in microbial stimulation in favor of other limiting resources.

Compost↗

DOE OSTI2024

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Modeling impurity concentrations in liquid argon detectors

Impurities in noble liquid detectors used for neutrino and dark matter experiments can significantly impact the quality of data. We present an experimentally verified model for describing the dynamics of impurity distributions in liquid argon (LAr) detectors. The model considers sources, sinks, and transport of impurities within and between the gas and liquid argon phases. Measurements of oxygen concentrations in a 20-L LAr multi-purpose test stand are compared to calculations made with this model to show that an accurate description of the concentrations under various operational conditions can be obtained. A result of this analysis is a determination of Henry’s coefficient for oxygen in LAr. Additionally, these calculations also show that some processes have small effects on the impurity dynamics and excluding them yields a solution as a sum of two exponential terms. This solution provides a simple way to extract Henry’s coefficient with negligible approximation error. It is applied to the data and the Henry’s coefficient for oxygen in LAr is obtained as 0.84 +0.09 -0.05 , consistent with literature results. Finally, based on the analysis of the data with the model, we further suggest that, for a large liquid argon detector, barriers to flow (“baffles”) installed in the gas phase to restrict flow can help reduce the ultimate impurity concentration in the LAr.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

A Quasi-Ordered Mn-Rich Cathode with Highly Reversible Oxygen Anion Redox Chemistry

Anionic oxygen redox chemistry in Li-rich Mn-based layer oxide cathodes represents a transformative approach for boosting the energy density of next-generation lithium-ion batteries. However, conventional oxygen redox reactions often induce oxygen dimerization at high voltages, leading to irreversible lattice oxygen loss and a rapid voltage fade. Herein, we achieve highly reversible oxygen redox chemistry through a new quasi-ordered structural design that incorporates both intra- and interlayer cation disorder configurations. This unique structure significantly enhances lattice oxygen stability, effectively stabilizes oxidized oxygen, and inhibits the formation of peroxo- or superoxol-like species, thereby enabling anionic redox reactions to proceed reversibly even at deep delithiation states. The quasi-ordered design mitigates irreversible phase transitions and preserves the structural integrity throughout extended cycling. Consequently, the proposed cathode demonstrates exceptional cyclability with negligible capacity and voltage fade, retaining 99% capacity and 98% average voltage after long-term cycling. Finally this work provides fresh insights into addressing issues related to lattice oxygen instabilities and reforming strategies for developing long-life, high-energy-density anionic redox cathode materials for advanced batteries.

chemical structure↗