Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “work functions”

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 OSTI2022

Emergent Magnetic States and Tunable Exchange Bias at 3d Nitride Heterointerfaces

Abstract Interfacial magnetism stimulates the discovery of giant magnetoresistance (MR) and spin–orbital coupling across the heterointerfaces, facilitating the intimate correlation between spin transport and complex magnetic structures. Over decades, functional heterointerfaces composed of nitrides have seldom been explored due to the difficulty in synthesizing high‐quality nitride films with correct compositions. Here, the fabrication of single‐crystalline ferromagnetic Fe 3 N thin films with precisely controlled thicknesses is reported. As film thickness decreases, the magnetization dramatically deteriorates, and the electronic state changes from metallic to insulating. Strikingly, the high‐temperature ferromagnetism is maintained in a Fe 3 N layer with a thickness down to 2 u.c. (≈8 Å). The MR exhibits a strong in‐plane anisotropy; meanwhile, the anomalous Hall resistivity reverses its sign when the Fe 3 N layer thickness exceeds 5 u.c. Furthermore, a sizable exchange bias is observed at the interfaces between a ferromagnetic Fe 3 N and an antiferromagnetic CrN. The exchange bias field and saturation moment strongly depend on the controllable bending curvature using the cylinder diameter engineering technique, implying the tunable magnetic states under lattice deformation. This work provides a guideline for exploring functional nitride films and applying their interfacial phenomena for innovative perspectives toward practical applications.

Jin, Qiao↗
DOE OSTI2021

Historical Contingency in Microbial Resilience to Hydrologic Perturbations

Development of reliable biogeochemical models requires a mechanistic consideration of microbial interactions with hydrology. Microbial response to and its recovery after hydrologic perturbations (i.e., resilience) is a critical component to understand in this regard, but generally difficult to predict because the impacts of future events can be dependent on the history of perturbations (i.e., historical contingency). Fundamental issues underlying this phenomenon include how microbial resilience to hydrologic perturbations is influenced by historical contingency and how their relationships vary depending on the characteristics of microbial functions. To answer these questions, we considered a simple microbial community composed of two species that redundantly consume a common substrate but specialize in producing distinct products and developed a continuous flow reactor model where the two species grow with trade-offs along the flow rate. Simulations of this model revealed that (1) the history of hydrologic perturbations can lead to the shifts in microbial populations, which consequently affect the community's functional dynamics, and (2) while historical contingency in resilience was consistently predicted for all microbial functions, it was more pronounced for specialized functions, compared to the redundant function. As a signature of historical contingency, our model also predicted the emergence of hysteresis in the transitions across conditions, a critical aspect that can affect transient formation of intermediate compounds in biogeochemistry. This work presents microbial growth traits and their functional redundancy or specialization as fundamental factors that control historical contingencies in resilience.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2006

In-Space Radiator Shape Optimization using Genetic Algorithms

Future space exploration missions will require the development of more advanced in-space radiators. These radiators should be highly efficient and lightweight, deployable heat rejection systems. Typical radiators for in-space heat mitigation commonly comprise a substantial portion of the total vehicle mass. A small mass savings of even 5-10% can greatly improve vehicle performance. The objective of this paper is to present the development of detailed tools for the analysis and design of in-space radiators using evolutionary computation techniques. The optimality criterion is defined as a two-dimensional radiator with a shape demonstrating the smallest mass for the greatest overall heat transfer, thus the end result is a set of highly functional radiator designs. This cross-disciplinary work combines topology optimization and thermal analysis design by means of a genetic algorithm The proposed design tool consists of the following steps; design parameterization based on the exterior boundary of the radiator, objective function definition (mass minimization and heat loss maximization), objective function evaluation via finite element analysis (thermal radiation analysis) and optimization based on evolutionary algorithms. The radiator design problem is defined as follows: the input force is a driving temperature and the output reaction is heat loss. Appropriate modeling of the space environment is added to capture its effect on the radiator. The design parameters chosen for this radiator shape optimization problem fall into two classes, variable height along the width of the radiator and a spline curve defining the -material boundary of the radiator. The implementation of multiple design parameter schemes allows the user to have more confidence in the radiator optimization tool upon demonstration of convergence between the two design parameter schemes. This tool easily allows the user to manipulate the driving temperature regions thus permitting detailed design of in-space radiators for unique situations. Preliminary results indicate an optimized shape following that of the temperature distribution regions in the "cooler" portions of the radiator. The results closely follow the expected radiator shape.

Hull, Patrick V.↗
DOE OSTI2026

Soil physical and chemical measurements for topsoils collected during NEON campaign in East River, CO (06/14/2018-06/28/2018)

The package is part of the DOE Watershed Function Science Focus Area (SFA) project and includes soil physical and chemical measurements from topsoils collected at the East River, Colorado, in conjunction with the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) survey conducted in June 2018. The soil measurements include soil bulk density, soil volumetric water content, soil microbial biomass C (Carbon), N (Nitrogen) and C:N (C to N ratio), soil DNA yield, soil total extractable organic C, soil total extractable N, soil extractable nitrate, soil extractable ammonium, soil dissolved inorganic N, soil dissolved organic N, soil pH, soil TOC400 (total organic carbon at 400°C), soil ROC (residual oxidizable carbon), soil TIC (total inorganic carbon), soil TOC (total organic carbon), soil TC (total carbon), soil N, soil OM (organic matter) loss on ignition. Additional associated site metadata can be found in the ESS-DIVE package 10.15485/1618130. The dataset includes (1) 2018_NEON_soil_physical_chemical_measurements.csv: soil physical and chemical measurements indexed by soil sample IGSNs; (2) samples.csv: sample metadata file used to register International Generic Sample Numbers (IGSNs); (3) flmd.csv: file level metadata file; and (4) dd.csv: data dictionary file. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS (Catchment Hydrology and ↗
DOE OSTI2026

FAD-Toolset (Floating Array Design Toolset) [SWR-26-056]

The Floating Array Design (FAD) Toolset is a collection of tools for modeling and designing arrays of floating offshore structures. It was originally designed for floating wind systems but has applicability for many offshore applications. A core part of the FAD Toolset is the floating array model, which serves as a high-level library for efficiently modeling a floating array, such as a floating wind array. It combines site condition information and a description of the floating array design, and contains functions for evaluating the array's behavior considering the site conditions. For example, it combines information about site soil conditions, mooring line loads, and an array's anchor characteristics to estimate the holding capacity of each anchor. The library works in conjunction with the tools RAFT, MoorPy, and FLORIS to model floating platforms, wind turbines, mooring systems, power cables, and array wakes respectively. Layered on top of the floating array model is a set of design tools that can be used for algorithmically adjusting or optimizing parts of the a floating array. Specific tools existing for mooring lines, shared mooring systems, dynamic power cables, static power cable routing, and overall array layout. These capabilities work with the design representation and evaluation functions in the floating array model, and they can be applied by users in various combinations to suit different purposes. In addition to standalone uses of the FAD Toolset, a coupling has been made with Ard, (https://github.com/NLRWindSystems/Ard) a sophisticated and flexible wind farm optimization tool. This coupling allows Ard to use certain mooring system capabilities from FAD to perform layout optimization of floating wind farms with Ard's more advanced layout optimization capabilities. The FAD Toolset works with the IEA Wind Task 49 Ontology (https://github.com/IEAWindTask49/Ontology), which provides a standardized format for describing floating wind farm sites and designs. See example use cases in our examples folder (https://github.com/NLRWindSystems/FAD-Toolset/blob/main/examples/README.md) For working with the library, it is important to understand the floating array model structure, which is described more here: https://github.com/NLRWindSystems/FAD-Toolset/blob/main/fad/README.md.

Sirkis, Leah [National Laboratory of the Rockies (↗
DOE OSTI2020

Loschmidt-amplitude wave function spectroscopy and the physics of dynamically driven phase transitions

In this work, we introduce the Loschmidt amplitude as a powerful tool to perform spectroscopy of generic many-body wave functions and use it to interrogate the wave function obtained after ramping the transverse field quantum Ising model through its quantum critical point. Previous results are confirmed and a more complete understanding of the population of defects and of the effects of magnon-magnon interaction or finite-size corrections is obtained. The influence of quantum coherence is clarified.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1995

Analysis of surface cracks at hole by a 3-D weight function method with stresses from finite element method

Parallel with the work in Part-1, stress intensity factors for semi-elliptical surface cracks emanating from a circular hole are determined. The 3-D weight function method with the 3D finite element solutions for the uncracked stress distribution as in Part-1 is used for the analysis. Two different loading conditions, i.e. remote tension and wedge loading, are considered for a wide range in geometrical parameters. Both single and double surface cracks are studied and compared with other solutions available in the literature. Typical crack opening displacements are also provided.

Zhao, W.↗
DOE OSTI2021

Distributed Wind Controls: A Research Roadmap for Microgrids, Infrastructure Resilience, and Controls Launchpad (MIRACL)

The Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL) project’s controls research area aims to expand the benefits from distributed wind (DW) generation assets beyond solely providing low cost power directly to consumers. To make distributed wind turbines operate more effectively there is a need for more advanced ways to control them, allowing power companies, businesses, and energy consumers to take advantage of the unique characteristics of wind energy. DW can contribute to the safe and secure operation of grid by providing services such as voltage regulation, frequency regulation, spinning reserves, and even black start capabilities. In the larger DER context where much of the research has focused around solar technologies, the inertia available in wind technologies has generally not been considered. For wind turbines to provide these services in an optimal and expanded way, development and demonstration of control methods and communication interfaces within a grid and microgrid framework are required. In this document, NREL led the literature review in collaboration with SNL, of DER controls-focused integration research to establish a baseline for the controls research under the MIRACL project to identify specific control functions to be focused on throughout this project. This literature review focused primarily on the control functions of variable distributed generation, largely pulling from past solar PV and battery controls works, with a specific focus on applicability for distributed wind energy systems. The goal of this document is to identify a research roadmap based on the open literature and past national laboratory works to inform advanced wind turbine and power electronics control functions for four use cases: 1) distributed wind in isolated systems, 2) grid-connected microgrids (wind-hybrid systems and islanded operation), 3) behind-the-meter distributed wind applications in the power distribution systems, and 4) front-of-the-meter distributed wind applications in the power distribution systems.

17 WIND ENERGY↗
NASA NTRS2020

Electrical Cable Design for Urban Air Mobility Aircraft

Urban Air Mobility (UAM) describes a new type of aviation focused on efficient flight within urban areas for moving people and goods. There are many different configurations of UAM vehicles, but they generally use an electric motor driving a propeller or ducted fan powered by batteries or a hybrid electric power generation system. Transmission cables are used to move energy from the storage or generation system to the electric motors. Though terrestrial power transmission cables are well established technology, aviation applications bring a whole host of new design challenges that are not typical considerations in terrestrial applications. Aircraft power transmission cable designs must compromise between resistance-per-length, weight-per-length, volume constraints, and other essential qualities. In this paper we use a multidisciplinary design optimization to explore the sensitivity of these qualities to a representative tiltwing turboelectric UAM aircraft concept. This is performed by coupling propulsion and thermal models for a given mission criteria. Results presented indicate that decreasing cable weight at the expense of increasing cable volume or cooling demand is effective at minimizing maximum takeoff weight (MTO). These findings indicate that subsystem designers should update their modeling approach in order to contribute to system-level optimality for highly-coupled novel aircraft. Mobility (UAM) vehicles have the potential to change urban and intra-urban transport in new and interesting ways. In a series of two papers Johnson et al.1 and Silva et al.2 presented four reference vehicle configurations that could service different niches in the UAM aviation category. Of those, this paper focuses on the Vertical Take-off and Landing (VTOL) tiltwing configuration shown in Figure 1. This configuration uses a turboelectric power system, feeding power from a turbo-generator through a system of transmission cables to four motors spinning large propellers on the wings. Previous work on electric cable subsystems leaves much yet to be explored, especially in the realm of subsystem coupling. Several aircraft optimization studies1, 3, 4 only considered aircraft electrical cable weight and ignored thermal effects. Electric and hybrid-electric aircraft studies by Mueller et al.5 and Hoelzen et al.6 selected a cable material but did not investigate alternative materials. Advanced cable materials have been examined by a number of authors: Alvarenga7 examined carbon nanotube (CNT) conductors for low-power applications. De Groh8, 9 examined CNT conductors for motor winding applications. Behabtu et al.,10 and Zhao et al.11 examined CNT conductors for a general applications. There were some studies that examined the thermal effects of cables but they did not allow the cable material to change; El-Kady12 optimized ground-cable insulation and cooling subject constraints. Vratny13 selected cable material based on vehicle power demand, and required resulting cable heat to be dissipated by the Thermal Management System (TMS). None of these previous studies allowed for the selection of the cable material based on a system level optimization goal. Instead, they focused on sub-system optimality such as minimum weight, which comes at the expense of incurring additional costs for other subsystems. Dama14 selected overhead transmission line materials using a weighting function and thermal constraints. However, that work was not coupled with any aircraft subsystems like a TMS. The traditional aircraft design approach, which relies on assembling groups of optimal subsystems, breaks down when considering novel aircraft concepts like the tiltwing vehicle. In a large part, this is because novel concepts have a much higher degree of interaction or coupling between subsystems. For example, when a cable creates heat, this heat needs to be dissipated by the TMS, which needs power supplied by the turbine, and delivering the power creates more heat. The cable, the TMS, and the turbine are all coupled. A change to one subsystem will affect all the other subsystems, much to the consternation of subsystem design experts. Multidisciplinary optimization is the design approach that can address these challenges. However, to fully take advantage of this, we must change the way we think about subsystem design. Specifically, we must move away from point design, and focus on creating solution spaces. The work presented in this paper uses the multidisciplinary optimization approach with aircraft level models to study the system-level sensitivity of cable traits: weight-per-length and resistance-per-length. Additionally, we examined the effects of vehicle imposed volume constraints on these traits. This is useful for three purposes: (1) to demonstrate a framework that can perform a coupled analysis between the aircraft thermal and propulsion systems, (2) to provide a method by which future cable designs can be evaluated against each other given a system-level design goal, (3) to provide insight into what cable properties may be promising for future research. This last element is explored given the caveat that the models contained in this analysis do not represent high-fidelity systems. Thus, while we can demonstrate coupling in between systems, the exact system-level sensitivity to a given parameter may change if a subsystem model or the assumptions governing that model change. The organization of this paper is as follows, in Sec II we outline a method to combine the VTOL vehicle design and cable information in order to produce cables sensitivity studies. Results analysis and discussion are contained in Sec III. Conclusions are presented in Sec IV.

Aretskin-Hariton, Eliot D.↗
DOE OSTI2024

Phase-field modeling for restructuring in the dark zone of high burnup UO 2

This report summarizes the mesoscale modeling work performed in fiscal year 2024 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstructural evolution and restructuring observed in the dark regions of high burnup UO 2 nuclear fuel. This is the first attempt to realistically simulate the restructuring behavior observed in different region of a high burnup fuel. We employ a grand-potential based phase-field model to concurrently evaluate the formation of subgrains and growth of fission bubbles within the fuel. A energy-based subgrain formation criteria is introduced to simulate the restructuring process. Effect of different initial conditions and different modeling parameters are studies systematically to capture how each of these parameters influence the characteristics of the restructured fuel. It is observed that the subgrain formation begins around existing fission gas bubbles and then proceeds towards triple junctions, grain boundaries and grain interiors. It is demonstrated that restructuring is influenced by a combination of initial dislocation densities, subgrain formation rate, and temperature. Rate of restructuring increases with increase in fuel temperature. A restructuring bias is observed within the microstructure due to variation in defect accumulation among different grains. Furthermore, bubble sizes and distribution does not have a significant effect on rate of restructuring. The predicted microstructures resembles the characteristics of the restructured regions as observed in experiments. Finally, a correlation is presented that demonstrates the evolution of the restructuring volume fraction as a function of local effective burnup. This work provides a first of its kind restructuring model for darkzone that can be used by BISON for performance prediction of high burnup UO 2 fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Modeling Pb(II) Adsorption on Mineral Surfaces: Bridging Density Functional Theory and Experiment with Thermodynamic Insights

Despite decades of work on aqueous lead (Pb) adsorption on a-Fe2O3 (hematite) and a-Al2O3 (alumina), gaps between measurements and modeling obscure molecular-level understanding. Achieving well-matched geometries between theory and experimental for mineral-water interfaces is a hurdle, as surface functional group type and distribution must be accounted for in determining mechanisms. Additionally, computational methods that can describe the substrate are often not appropriate to capture aqueous effects. Progress requires focusing on well-studied and relevant systems, such as key facets (001),(012), and (110) of hematite and alumina, and ubiquitous contaminants such as aqueous Pb. In the past, bulk-parametrized bond-valence principles were used to rationalize Pb(II) adsorption trends. These approaches can break down at surfaces, where flexible bonding environments and adsorption-induced surface relaxations play a critical role. Here, we adapt and apply a density functional theory (DFT) and thermodynamics framework, integrating DFT-calculated energies with experimental data and electrochemical principles, to predict Pb(II) adsorption. Our model results capture trends across the full set of surfaces and predict that inner-sphere Pb(II) sorption on (001) alumina varies from unfavorable to weakly favorable across a range of pH conditions. This aligns with experiment insights that Pb(II) interacts at that surface through outer-sphere interactions. Extending to Fe(II) adsorption, we demonstrate a coverage-dependent site preference, potentially explaining disorder in overlayers grown by the oxidative adsorption of Fe(II) on hematite (001).

lead contamination↗
DOE OSTI2025

Montane Conifer, Aspen, Meadow, and Sagebrush Metagenome Resolved Genomes and Traits in East River Watershed, Colorado, USA

Climate change is driving vegetation shifts in mountain watersheds, with unknown impacts on biogeochemical cycles. We hypothesize that these shifts will reshape soil microbiomes and associated biogeochemical processes. As a part of Lawrence Berkeley National Laboratory (LBNL) Watershed Science Focus Area (SFA), we assessed microbiome and microbial functional trait differences between soils under conifer, aspen, forby meadows, and sagebrush across the East River Watershed, CO, controlling for elevation and aspect.Here we present metagenome assembled genomes (MAGs) for the bacterial and archaeal communities from soils 0-20cm in depth across three locations in the watershed—Headwaters, Upper Reaches, and Lower Reaches from August 3-11th 2016. Each location was further subdivided into two blocks, with one block on a west facing aspect, and two on the east aspect of the valley. Within blocks, two samples per vegetation type were taken (one at each depth). This resulted in 66 samples, which were sequenced at JGI and can be found under the Joint Genome Institute (JGI) Genomes Online Database (GOLD) sequencing project Gs0118068. Metagenomes were assembled through an inhouse pipeline (see methods), binned using four autobinners (concoct, maxbin2, metabat2, and vamb) and consolidated using dastool. The consolidated bins from all metagenomes were pooled, filtered by completeness (>75%) and contamination (<25%), and dereplicated at 95% ANI using drep. The dataset includes a zip file of 687 genomes (Vegtype_MAGS.zip), the accession numbers for the underlying metagenomes, a csv file with MAG quality metrics and taxonomy from Genome Taxonomy Database (GTDB) and National Center for Biotechnology Information (NCBI) taxonomic representative genome proteins (EastRiver_Vegtype_drep_genome_info.csv), and a file containing MAG quality metrics and taxonomy (gtdb_drep_bin_taxonomy.csv). The dataset additionally includes a sample metadata file (EastRiver_Vegtype_sample_metadata.csv), a metadata file used to register associated samples with IGSNs (International Generic Sample Numbers) (samples.csv), a Google KML file for the sampled locations (sample_collection_sites.kml), a location metadata file (locations.csv), a file-level metadata file (flmd.csv), and a data dictionary (dd.csv) file.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Comparative study of Minnesota functionals performance on ferroelectric BaTiO 3 and PbTiO 3

Density functional theory based simulations are the leading tool for the computational investigation of ferroelectrics and the parametrization of their classical potentials. However, the predictions often depend strongly on the exchange correlation functional. The most popular choices, LDA and GGA, tend to underestimate or overestimate some structural, electric, and energy properties. These impede development of highly accurate classical potentials that extend the reach of first-principles simulations to finite temperatures and realistic sizes. In this work, we investigate the performance of recently developed Minnesota exchange correlation functionals on prototypical ferroelectrics BaTiO 3 , and PbTiO 3 , in comparison with some popular ones. We find that there exists strong correlation between predictions for some properties (tetragonality, phase energy difference, and polarization) by different functionals. Along the correlation line, we find a range of functionals (including some from the Minnesota suite) whose predictions fall in between those of LDA and GGA and, therefore, offer an improvement. A way to relatively rank functional performance with respect to chosen benchmarks is proposed and applied to identify the top performers for BaTiO 3 and PbTiO 3 . Here, the performance is found to be material dependent. Therefore, we propose that the performance assessment carried out in this work is employed to "screen" functionals prior to their use on ferroelectrics. Quick and computationally inexpensive, it is likely to lead to improved descriptions, especially for classical potential parametrization. The Minnesota functional suite is found to be suitable for this task.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Metallic Fuel Performance Code Requirements for the Versatile Test Reactor Project

Metallic nuclear fuels have been proposed for use in several emerging nuclear reactor designs, and a number of codes have been developed to model these fuels and assess their performance. Qualification of these metallic nuclear fuels will ultimately require monitored irradiation of lead test assemblies, but the use of fuel performance codes can reduce the uncertainty associated with these efforts by quantifying uncertainties and estimating margins to failure ahead of time. In this work, metallic fuel performance code requirements are defined for the Versatile Test Reactor (VTR) project using input from the Experimental Breeder Reactor II functional requirements, operational requirements, and design criteria. This work focuses on the thermomechanical responses and irradiation behaviors of metallic fuel elements that are representative of the proposed VTR driver fuel design concepts. As such, many of the code requirements and physics/modeling discussions in this work are expected to be broadly applicable to metallic fuel applications outside of the VTR project.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1997

Characterization of Ionospheric Flow Patterns During Times of Unusual Auroral Activity Using DE and DMSP Data

This final report describes the work done by Dr. Marc Hairston and Dr. Rod Heelis on NASA SR&T grant NAGW-4411 studying the theta aurora using DE-1 ultraviolet imager data and DMSP particle data. This report covers the period from summer 1995 through summer 1996 along with a review of the previous work. Previous work on this grant looked at the time period from the launch of DMSP-F8 in June 1987 through the end of mission of DE-1 in summer 1991. Despite the sporadic and decreasing frequency of observations from DE-1 over this time period, we were able to identify six events for study where the DE-1 imager observed a theta aurora occurring during a period in which the DMSP-F8 satellite flew through the region of the aurora at an altitude of 800 km. We focused on the best two events where the theta aurora persisted for an extended period so that we could observe the DMSP particle signatures in both the hemisphere observed by DE-1 and in the other polar hemisphere immediately before or afterwards. These results were presented at the Fall 1994 meeting of the AGU. Initially we had hoped to expand on this work. However, further work showed that none of these events gave a clear enough signature in the DMSP data for us to identify a theta aurora in the hemisphere opposite to the hemisphere imaged by DE-1. Without that, there were no results from this work that were new enough to warrant publication. So instead we used the final year of the grant to work with our colleague, J. A. Cumnock on a similar project using DE data to study the evolution of theta auroras as a function of the IMF A paper from that work was published which acknowledged this grant and a copy of that paper is included with this final report.

Hairston, Marc↗
DOE OSTI2023

FUN-BioCROP model with litter decomposition parameters derived from the LIDET dataset

This repository contains the code and data necessary to run the FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms, and Protection) model with litter decomposition parameters derived from a modified Monte Carlo simulation that used the Long-term Intersite Decomposition Experiment Team dataset. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of Files: FUNBioCROP_LIDET Study.Rmd R code with FUN-BioCROP model that can be run with 10 different sets of parameters for litter decomposition (Baseline, LIDET, or eight other best parameter sets identified in the modified Monte Carlo simulation. CORPSE Functions_Bioenergy_V2.R Code with CORPSE model functions, called by FUNBioCROP_LIDET Study.Rmd Model Input Data: bulk.csv, bulk_till.csv, rhizo.csv, rhizo_till.csv, litter.csv Initial C and N (kg C or N/m2) pool values for each soil compartment, final values from spin up. All five files have the same columns: (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 Model Input Data: FluxTower_AvgSoilT.csv: Average daily soil temperature (oC) at 10 cm depth at University of Illinois Urbana-Champaign (UIUC) Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: FluxTower_AvgSoilVWC.csv: Average daily soil volumetric water content (VWC) at 10 cm depth at UIUC Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: input_CCS_LIDET Study.csv: This file has daily data to run FUN-BioCROP (Column - Description - Units): yr - calendar year - year doy - day of year (1 to 365) (no leap year) - day anpp - aboveground NPP (DayCent) - kg C/m2/day bnpp - belowground NPP (DayCent) - kg C/m2/day aglivc - live aboveground biomass carbon (DayCent) - kg C/m2 bglivcj - live juvenile fine root biomass carbon (DayCent) - kg C/m2 bglivcm - live mature fine root biomass carbon (DayCent) - kg C/m2 aglivn - live aboveground biomass nitrogen (DayCent) - kg N/m2 bglivnj - live juvenile fine root biomass nitrogen (DayCent) - kg N/m2 bglivnm - live mature fine root biomass nitrogen (DayCent) - kg N/m2 nyr - simulation year - year cult - indicates a cultivation event (0 or 1) crop - indicates a new crop (0 or 1) fert - indicates a fertilizer event (0 or 1) frst - indicates the first day of the growing season (0 or 1) harv - indicates a harvest event (0 or 1) last - indicates the end of the growing season (0 or 1) croptype - crop type (0=none; 1=alfalfa; 2=corn; 3=grass clover pasture; 4=soybean; 5=wheat) cropsrl - crop specific root length - mm/g root cultrhizmix - fraction of rhizosphere mixed with bulk soil during cultivation (0.0-1.0) - fraction cultlitmix - fraction of litter mixed with bulk soil during cultivation (0.0-1.0) - fraction harvremov - fraction of above ground biomass removed during harvest (0.0-1.0) - fraction fertamt - fertilization amount - g N/m2 lifehist - plant life history (0 = annual, 1 = perennial) froot_turnover_c - amount of C in fine root turnover - kg C/m2 froot_turnover_n - amount of N in fine root turnover - kg N/m2 agrd_turnover_c - amount of C in aboveground biomass turnover - kg C/m2 agrd_turnover_n - amount of N in aboveground biomass turnover - kg N/m2 leaf_litter_fastfrac - Fast decomposing fraction of leaf litter (0.0-1.0) - fraction root_litter_fastfrac - Fast decomposing fraction of root litter (0.0-1.0) - fraction root_diameter - root diameter - mm root_length - root length - mm root/m2 rhizo_frac - fraction of total soil volume that is rhizosphere (0.0 - 1.0) - fraction date - date in format YYYY-MM-DD Instructions: Save the model code ("FUN-BioCROP_LIDET Study.Rmd") and accompanything files (data streams and CORPSE function code) in the same folder. In model code "Chunk 3: Load CORPSE Data Streams" set the working directory (setwd) to the folder with the files saved in step #1. In "Chunk 5: Define LIDET parameter sets" select the litter decomposition parameter set to be used in the run, and comment out all other sets. If changing any parameter values, edit them in "Chunk 6: Load parameters." Run all chunks up to and including "Chunk 10: Prepare Data for Export." In "Chunk 11: Export Output Data" edit data frames for export and filenames, as necessary. "Chunk 12: Graph Total Soil C" makes a figure of C remaining over the model run period. Description of each model chunk (in file FUN-BioCROP_LIDET Study.Rmd): Chunk 1: Remove all functions, clear memory. Removes all functions from R environment, clears the memory. Chunk 2: Load Packages. Loads packages necessary to run the code. Chunk 3: Load CORPSE Data Streams. Sets the working directory and loads the data files necessary to run CORPSE. Chunk 4: Load CORPSE Functions. Loads the R script with CORPSE functions from the working directory, "CORPSE Functions_Bioenergy_V2.R". Chunk 5: Define LIDET parameter sets. Has ten different parameter sets for litter decomposition tested in this study: Baseline parameters, LIDET parameters, and the other 8 best performing parameter sets identified in the modified Monte Carlo. To run the model, all but one parameter set must be commented out. Chunk 6: Load Parameters. Loads all fixed parameters to run the model. Data frame with definitions of parameters is in the CORPSE function script "CORPSE Functions_Bioenergy_V2.R" Chunk 7: Prepare Data Streams. Takes data streams loaded in Chunk 3 and puts them in the format necessary to run the model. The model is coded to run at least two sites at a time, so if only one site is being run it must be run in duplicate. Individual data tables of daily values are created in this chunk from the input data file. Chunk 8: Set Initial Conditions. Creates data tables of soil C and N pools for each soil compartment (rhizo_till, rhizo, bulk_till, bulk, litter) and loads initial values into the data tables. Creates lists for each soil compartment to hold model output. Chunk 9: Load FUN Data and Set Up Matrices. Uses DayCent data to calculate FUN input data: root and leaf N demand, total N demand, plant CN, leaf N available for retranslocation, and litter production. Creates matrices for FUN model outputs. Chunk 10: Run Model. Runs the model. Chunk 11: Prepare Data for Export. Combines data from each day saved as lists into data frames for each soil compartment. Adds values from all soil compartments together to calculate total soil values, creates separate data frames for each soil C and N pool (e.g., protected slow C) for the total soil value. Adds different C and N pools together to calculate total soil C and N for all layers. Creates data frame of ratio of protected to unprotected SOC. Organizes FUN data for export. Chunk 12: Export Results. Exports CSV files of model results to the working directory. Chunk 13: Graph Total Soil C. Makes figure of C remaining over time. Related Links: Original FUN-BioCROP model: https://github.com/BrzostekEcologyLab/FUN-BioCROP LIDET dataset: https://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023

Juice, Stephanie↗
DOE OSTI2024

Multiple Functional Bonds Integrated Interphases for Long Cycle Sodium-Ion Batteries

Sodium-ion batteries (SIBs) have garnered significant interest as one of the most promising energy suppliers for power grid energy storage. However, the poor electrode/electrolyte interfacial stability leads to continual electrolyte decomposition and transition metal dissolution, resulting in rapid performance degradation of SIBs. In this work, we propose a strategy integrating multiple functional bonds to regulate electrode/electrolyte interphase by triple-coupling of succinonitrile (SN), sodium hexafluorophosphate (NaPF 6 ) and fluorinated ethylene carbonate (FEC). Theoretical calculation and experiment results show that the solvation structure of Na + and ClO 4 – is effectively reconfigured by the solvated FEC, SN and PF 6 – in PC-based carbonate electrolyte. The newly developed electrolyte demonstrates increased Na + -FEC coordination, weakened interaction of Na + -PC and participation of SN and PF 6 – anions in solvation, resulting in the formation of a conformal interfacial layer comprising of sodium oxynitrides (NaN x O y ), sodium fluoride (NaF) and phosphorus oxide compounds (NaP x O y ). Consequently, a 3 Ah pouch full cell of hard carbon//NaNi 1/3 Fe 1/3 Mn 1/3 O 2 exhibits an excellent capacity retention of 90.4 % after 1000 cycles. Detailed postmortem analysis of interface chemistry is further illustrated by multiple characterization methods. Finally, this study provides a new avenue for developing electrolyte formulations with multiple functional bonds integrated interphases to significantly improve the long-term cycling stability of SIBs.

25 ENERGY STORAGE↗
DOE OSTI2026

Meteorological and Soil Data from Ecohydrology Sensor Towers at Pump House and Snodgrass Mountain in East River Watershed, Colorado, 2019-2025

This data package includes hourly meteorological and soil sensor data at eight ecohydrology monitoring sites in East River Watershed, Colorado as part of the Watershed Function Scientific Focus Area (WFSFA) research led by Lawrence Berkeley National Lab (LBNL). Four field sites were located on the hillslope of East River (ER) near Pump House (PH) at Mount Crested Butte (ER-PHS1 to 4), and the other four are in the Snodgrass Mountain (SG) area (SG-EHS5 to 8). In terms of vegetation cover, three sites are in montane grasslands (ER-PHS1, ER-PHS2, and SG-EHS5), three are below evergreen conifer canopy (ER-PHS3, SG-EHS6, and SG-EHS7), and two are below deciduous aspen canopy (ER-PHS4 and SG-EHS8). The monitoring period began in October 2019 at the East River sites, in October 2020 at SG-EHS5 and SG-EHS6, and in October 2021 at SG-EHS7 and SG-EHS8. In September 2024, all four East River sites were fully retired. The four Snodgrass Mountain sites remain active. Each site is equipped with a comprehensive suite of meteorological sensors on a tripod and soil sensors that measure weather, energy fluxes, and soil variables. This data package includes measurements from ten different types of sensors and up to thirteen individual sensors per site, including (1) a weather station (measurement height ranges from 2.8~3.8 meters (m) above ground), (2) a quantum sensor for photosynthetic active radiation (PAR) (2.4~3.3m), (3) a net radiometer (1.7~2.1m), (4) an infrared radiometer (1.6~2.2m), (5) a sonic distance sensor (1.5~1.9m), (6) a soil carbon dioxide (CO2) flux chamber (0m), (7) a soil heat flux plate (-0.05m below ground), (8) a soil oxygen sensor (-0.3m), (9) a soil water potential sensor (-0.3m), and (10) soil water content sensors at 3~4 depths (-1.15 ~ -0.1m). A total of twenty-three variables is reported in this data package, including (1) atmospheric variables: air temperature (TA), atmospheric pressure (PA), vapor pressure (VP), and vapor pressure deficit (VPD), (2) precipitation variables: rain precipitation (P) and snow depth (D_SNOW), (3) energy fluxes variables: four-component net radiation (NETRAD) (shortwave/longwave incoming/outgoing radiation, SW_IN, SW_OUT, LW_IN, LW_OUT), photosynthetic photon flux density (PPFD), and soil heat flux (G), (4) soil variables: soil water content (SWC), soil water potential (SWP), soil temperature (TS), soil bulk electrical conductivity (COND_SOIL), and soil gaseous oxygen concentration (O2_SOIL), (5) wind variables: two-dimensional wind speed (WS), gust speed (WS_MAX), and wind direction (WD), and (6) surface variables: surface infrared temperature (T_CANOPY) and soil CO2 flux (CO2_SOIL). Please see the Methods section for data processing and QA/QC steps taken to generate the hourly datasets. The following files are included in this data package (notes on version: v{x}-{y}, where x is the metadata version, and y is the data version, when applicable): (1) “metadata_site_v{x}-{y}.csv” - a site metadata file that summarizes location information of all sites, including site ID, description, coordinates, timeframe, elevation, and vegetation cover, (2) “metadata_instrument_v{x}-{y}.csv” - an instrument metadata file that summarizes sensor information of all sites, including sensor manufacturer and model, measurement height, and sampling and averaging interval of all variables, (3) "data_{SITE_ID}_v{x}-{y}.csv" - eight data files that contain hourly data of each site indicated by {SITE_ID} in the filename, (4) “/figure/data_{SITE_ID}_v{x}-{y}.png" - eight figures that help visualize data of each site indicated by {SITE_ID} in the filename, (5) “/photo/*” - photos of each site indicated by {SITE_ID} in the filename, and (6) four file level metadata (flmd.csv) and data dictionary (*_dd.csv) files that summarize file, header, column, and variable information of all files. Notes: (1) Measurement height: Each variable name is followed by conventional positional qualifiers “H_V_R”, where H indicates the relative horizontal positions of that specific variable, V the vertical positions, and R the replicates. In this data package, only the vertical qualifier V varies, and V increases from the highest vertical position (V=1) to the lowest. Variables with the same qualifier are not necessarily measured by the same sensor, and the same variable with the same qualifier across different sites are not necessarily measured at the same height. Please refer to “metadata_instrument.csv” for the sensor information and measurement heights, and whether a variable is measured below the canopy. (2) Variable availability: Snow depth is not available at ER-PHS3 and SG-EHS7. SWC, soil temperature, and soil bulk EC at the deepest depth (<-1m) are not available at SG-EHS6 and SG-EHS7. The missing value code for numeric variables is -9999, except for SWP. For SWP, the missing value code is +9999, because SWP values are negative. (3) Sampling frequency: Please refer to “metadata_instrument.csv” for the increase of sampling frequency of some variables from 30-min to 1-min at ER-PHS1 to 4 in July 2020. (4) Sensors: While the methods of each sensor are not detailed, all sensors are commercially available, and their methods can be found in their manuals. Please refer to “metadata_instrument.csv” for the sensor manufacturer and model information. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗