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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2023

Multiple redox mechanisms for water-gas shift reaction on Fe 3 O 4 (1 1 1) surface: A density functional theory and mean-field microkinetic modeling study

In this work, we investigated redox mechanisms for the high temperature (HT) water–gas shift reaction (WGSR) on the Fe 3 O 4 (1 1 1) surface using density functional theory (DFT) calculations and mean-field microkinetic modeling (MKM). The redox pathways branched into three Langmuir-Hinshelwood processes (LH1, LH2 and LH3) and one Mars-van-Krevelen (MVK) process (in the presence of oxygen vacancy) from multidentate binding of CO and CO 2 over four available reactive sites (Fe oct2 , Fe tet1 , Fe bridge, and O1). We found that the LH1 and LH2 processes have CO adsorption at a single iron site (Fe oct2 or Fe tet1 ), while the LH3 and MVK processes have stronger chemisorption of CO or CO 2 by both the Fe oct2 site and the O 1 site. From the mechanistic study of these reaction path, we recognized that availability of O 1 sites was key to proceed to either CO oxidation by the single site (LH1, LH2) or by the dual site (LH3, MVK). We observed that the reaction energetics were significantly different in the CO oxidation steps, where the single or dual site results in markedly different apparent activation energy (73–281 kJ/mol) and reaction rates (10 -3 –10 -10 mol∙m -2 ∙s -1 ) among the four reaction mechanisms. The utilization of mean field MKM with DFT reaction energetics helps explain the experimental debates for the catalytic reaction details as well as providing a possible direction to engineer the catalyst with higher activity.

42 ENGINEERING↗
DOE OSTI2020

Near Neutral pH Redox Flow Battery with Low Permeability and Long-Lifetime Phosphonated Viologen Active Species

A highly stable phosphonate-functionalized viologen is introduced in this work as the redox-active material in a negative potential electrolyte for aqueous redox flow batteries (ARFBs) operating at nearly neutral pH. The solubility is 1.23 m and the reduction potential is the lowest of any substituted viologen utilized in a flow battery, reaching –0.462 V versus SHE at pH = 9. The negative charges in both the oxidized and the reduced states of 1,1'-bis(3-phosphonopropyl)-[4,4'-bipyridine]-1,1'-diium dibromide ( BPP–Vi ) effect low permeability in cation exchange membranes and suppress a bimolecular mechanism of viologen decomposition. A flow battery pairing BPP–Vi with a ferrocyanide-based positive potential electrolyte across an inexpensive, non-fluorinated cation exchange membrane at pH = 9 exhibits an open-circuit voltage of 0.9 V and a capacity fade rate of 0.016% per day or 0.00069% per cycle. Overcharging leads to viologen decomposition, causing irreversible capacity fade. Our research introduces extremely stable, extremely low-permeating and low reduction potential redox active materials into near neutral ARFBs.

25 ENERGY STORAGE↗
DOE OSTI2020

Instability and treatments of the coupled discrete element and lattice Boltzmann method by the immersed moving boundary scheme

The immersed moving boundary (IMB) scheme has been extensively used to couple the discrete element method (DEM) with the lattice Boltzmann method (LBM). In the literature, only the formulation of IMB for lattice nodal cells covered by a single-solid particle was given. The treatment of situations where a nodal cell is covered by two or more solid particles is seldom discussed. It is found that some numerical instability can occur for such situations due to an inappropriate computation of the weighting function in the IMB formulation. This work presents an enhanced treatment that can resolve the issue and validates it using some benchmark tests. Furthermore, to avoid the extra costs associated with the treatment and simplify the complicated procedure introduced, a simplified IMB scheme is proposed. The accuracy of both enhanced and simplified IMB schemes are validated by test cases including single-particle sedimentation, two-particle drafting-kissing-tumbling phenomenon, and multiple-particle sedimentation. Then, the robustness of both schemes is examined and discussed using a specially designed flow past cylinders test. Overall, the simplified IMB scheme is proved to be robust and sufficiently accurate and simpler and more effective than the enhanced scheme.

42 ENGINEERING↗
DOE OSTI2021

Electropolymerized-molecularly imprinted polymers (E-MIPS) as sensing elements for the detection of dengue infection

Here, a straightforward in situ detection method for dengue infection was demonstrated through the molecular imprinting of a dengue nonstructural protein 1 (NS1) epitope into an electropolymerized molecularly imprinted polyterthiophene (E-MIP) film sensor. The key enabling step in the sensor fabrication is based on an epitope imprinting strategy, in which short peptide sequences derived from the original target molecules were employed as the main template for detection and analysis. The formation of the E-MIP sensor films was facilitated using cyclic voltammetry (CV) and monitored in situ by electrochemical quartz crystal microbalance (EC-QCM). Surface properties were analyzed using different techniques including atomic force microscopy (AFM), X-ray photoelectron spectroscopy (XPS), and polarization modulation-infrared reflection-adsorption (PM-IRRAS). The standard calibration curve (R = 0.9830) was generated for the detection of the epitope, Ac-VHTWTEQYKFQ-NH 2 , with a linear range of 0.2 to 30 μg/mL and detection limit of 0.073 μg/mL. A separate calibration curve (R = 0.9786) was obtained using spiked buffered solutions of dengue NS1 protein, which resulted in a linear range of 0.2 to 10 μg/mL and a detection limit of 0.056 μg/mL. The fabricated E-MIP sensor exhibited long-term stability, high sensitivity, and good selectivity towards the targeted molecules. These results indicated that the formation of the exact and stable cavity imprints in terms of size, shape, and functionalities was successful. In our future work, we aim to use our E-MIP sensors for NS1 detection in real-life samples such as serum and blood.

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

Pretreatment of lignocellulosic feedstocks for cellulose nanofibril production

Cellulose nanofibrils (CNFs) have attracted a great deal of research interest in recent years attributable to the low cost and abundance of lignocellulosic biomass from which they can be extracted. These materials have potential applications in a wide array of areas because of their unique properties such as ultra-high aspect ratios and specific strengths. However, the high energy required to extract CNFs from biomass through fibrillation often makes them prohibitively expensive or negates their inherent sustainability. As such, a variety of biomass treatments prior to fibrillation have been investigated by researchers to reduce the energy requirements of CNF extraction, improve the efficacy of biomass fibrillation and subsequent processes, and/or impart functionality in resulting nanofibrils. In this work, both widely used and emerging mechanical, chemical, and enzymatic pretreatments used prior to fibrillation of lignocellulosic biomass for CNF extraction are reviewed. Attention is given to the effect of these various pretreatments on the properties of the resulting CNFs. Finally, the energy consumption in fibrillation processes with and without the use of pretreatments is compared, and future perspectives on challenges and opportunities in lignocellulosic feedstock pretreatments are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

SuperCam Calibration Targets: Design and Development

Abstract SuperCam is a highly integrated remote-sensing instrumental suite for NASA’s Mars 2020 mission. It consists of a co-aligned combination of Laser-Induced Breakdown Spectroscopy (LIBS), Time-Resolved Raman and Luminescence (TRR/L), Visible and Infrared Spectroscopy (VISIR), together with sound recording (MIC) and high-magnification imaging techniques (RMI). They provide information on the mineralogy, geochemistry and mineral context around the Perseverance Rover. The calibration of this complex suite is a major challenge. Not only does each technique require its own standards or references, their combination also introduces new requirements to obtain optimal scientific output. Elemental composition, molecular vibrational features, fluorescence, morphology and texture provide a full picture of the sample with spectral information that needs to be co-aligned, correlated, and individually calibrated. The resulting hardware includes different kinds of targets, each one covering different needs of the instrument. Standards for imaging calibration, geological samples for mineral identification and chemometric calculations or spectral references to calibrate and evaluate the health of the instrument, are all included in the SuperCam Calibration Target (SCCT). The system also includes a specifically designed assembly in which the samples are mounted. This hardware allows the targets to survive the harsh environmental conditions of the launch, cruise, landing and operation on Mars during the whole mission. Here we summarize the design, development, integration, verification and functional testing of the SCCT. This work includes some key results obtained to verify the scientific outcome of the SuperCam system.

Manrique, J. A. (ORCID:0000000220532819)↗
DOE OSTI2020

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗