Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data similarity”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

MSD CoP Webinar: Energy and AI

Context: This webinar was hosted by the MultiSector Dynamics Community of Practice (MSD CoP; https://multisectordynamics.org). Abstract: Projections of the need for new data centers to support Artificial Intelligence (AI) are large but highly uncertain. Recent projections indicate up to a 15% annual growth rate in data center electricity demand within the next 5-10 years. Given that most electric utilities are required to have a reserve margin of roughly the same magnitude as the projected growth in demand, these new data center loads could soon threaten resource adequacy and reliability unless data centers build their own generation, interruptible loads are negotiated, commensurate new capacity and/or transmission is built, or some combination of these options. Similarly, depending on the cooling technology and geographic location of new data centers, they could threaten water adequacy in water scarce regions. This webinar will provide an overview of the interactions between energy and AI and highlight two MSD projects exploring the grid and water implications of new data centers to support AI. Presenters : Dr. Casey Burleyson (Pacific Northwest National Laboratory); Dr. Stephanie Morris (Pacific Northwest National Laboratory); Kendall Mongird (Pacific Northwest National Laboratory) Moderator: Patrick M. Reed (MSD CoP Facilitation Team) This webinar was held on: June 16th, 2025 from 1-2 PM EST.

Artificial Intelligence↗
DOE OSTI2025

Energy and AI: Evaluating Future Grid and Water Stress Due to Data Centers

Projections of the need for new data centers to support Artificial Intelligence (AI) are large but highly uncertain. Recent projections indicate up to a 15% annual growth rate in data center electricity demand within the next 5-10 years. Given that most electric utilities are required to have a reserve margin of roughly the same magnitude as the projected growth in demand, these new data center loads could soon threaten resource adequacy and reliability unless data centers build their own generation, interruptible loads are negotiated, commensurate new capacity and/or transmission is built, or some combination of these options. Similarly, depending on the cooling technology and geographic location of new data centers, they could threaten water adequacy in water scarce regions. This presentation highlights the grid and water implications of new data centers to support AI.

Mongird, Kendall (ORCID:0000000328077088)↗
DOE OSTI2022

Long-run Marginal Emission Rates for Electricity - Workbooks for 2021 Cambium Data

These workbooks contain modeled estimates of long-run marginal emission rates (LRMER) for the contiguous United States. The LRMER is an estimate of the rate of emissions that would be either induced or avoided by a long-term (i.e., more than several years) change in electrical demand. It incorporates both the projected changes to the electric grid, as well as the potential for an incremental change in electrical demand to influence the structural evolution of the grid (i.e., the building and retiring of capital assets, such as generators and transmission lines). It is therefore distinct from the more-commonly-known short-run marginal, which treats grid assets as fixed. The Levelized LRMER worksheet within each workbook is set up to produce a levelized long-run marginal emission rate based on user-provided inputs. These levelized LRMER values are intended for analysts to use when estimating the emissions induced (or avoided) by a long-term change in end-use electricity demand. There are two workbooks that supply the data at two different geographic resolutions: states and GEA regions (20 regions that are similar to, but not exactly the same as, the US EPA's eGRID regions). For more data underlying these emissions factors, see the Cambium 2021 project at https://cambium.nrel.gov/. For more details on the inputs into the scenarios available in the workbooks, see the Standard Scenarios 2021 Report (https://www.nrel.gov/docs/fy22osti/80641.pdf). This data was produced as part of the Cambium project. For more details about the methodology, see the Cambium Documentation: Version 2021 (https://www.nrel.gov/docs/fy22osti/81611.pdf). This data is planned to be updated annually. Information on the latest versions can be found at https://www.nrel.gov/analysis/cambium.html.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Long-run Marginal Emission Rates for Electricity - Workbooks for 2022 Cambium Data

These workbooks contain modeled estimates of long-run marginal emission rates (LRMER) for the contiguous United States. A LRMER is an estimate of the rate of emissions that would be either induced or avoided by a change in electric demand, taking into account how the change could influence both the operation as well as the structure of the grid (i.e., the building and retiring of capital assets, such as generators and transmission lines). It is therefore distinct from the more-commonly-known short-run marginal, which treat grid assets as fixed. Long-run marginal emissions rates are generally appropriate to use when trying to comprehensively estimate the impact of a long-lived (i.e., more than several years) intervention. There are two workbooks that supply the data at two different geographic resolutions: states and GEA regions (20 regions that are similar to, but not exactly the same as, the US EPA's eGRID regions). For more data underlying these emissions factors, see the Cambium 2022 project at https://scenarioviewer.nrel.gov/. For more details on input assumptions and methodology see the associated report (Cambium 2022 Scenario Descriptions and Documentation, https://www.nrel.gov/docs/fy23osti/84916.pdf). This data is planned to be updated annually. Information on the latest versions can be found at https://www.nrel.gov/analysis/cambium.html.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Statistically-informed deep learning for gravitational wave parameter estimation

We introduce deep learning models to estimate the masses of the binary components of black hole mergers, $(m_1,m_2)$, and three astrophysical properties of the post-merger compact remnant, namely, the final spin, $a_\mathrm f$, and the frequency and damping time of the ringdown oscillations of the fundamental $\ell = m = 2$ bar mode, $(\omega_\mathrm R, \omega_\mathrm I)$. Our neural networks combine a modified WaveNet architecture with contrastive learning and normalizing flow. We validate these models against a Gaussian conjugate prior family whose posterior distribution is described by a closed analytical expression. Upon confirming that our models produce statistically consistent results, we used them to estimate the astrophysical parameters $(m_1,m_2, a_\mathrm f, \omega_\mathrm R, \omega_\mathrm I)$ of five binary black holes: GW150914, GW170104, GW170814, GW190521 and GW190630. We use PyCBC Inference to directly compare traditional Bayesian methodologies for parameter estimation with our deep learning based posterior distributions. Our results show that our neural network models predict posterior distributions that encode physical correlations, and that our data-driven median results and 90% confidence intervals are similar to those produced with gravitational wave Bayesian analyses. This methodology requires a single V100 NVIDIA GPU to produce median values and posterior distributions within two milliseconds for each event. Furthermore, this neural network, and a tutorial for its use, are available at the Data and Learning Hub for Science.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Angular clustering properties of the DESI QSO target selection using DR9 Legacy Imaging Surveys

ABSTRACT The quasar target selection for the upcoming survey of the Dark Energy Spectroscopic Instrument (DESI) will be fixed for the next 5 yr. The aim of this work is to validate the quasar selection by studying the impact of imaging systematics as well as stellar and galactic contaminants, and to develop a procedure to mitigate them. Density fluctuations of quasar targets are found to be related to photometric properties such as seeing and depth of the Data Release 9 of the DESI Legacy Imaging Surveys. To model this complex relation, we explore machine learning algorithms (random forest and multilayer perceptron) as an alternative to the standard linear regression. Splitting the footprint of the Legacy Imaging Surveys into three regions according to photometric properties, we perform an independent analysis in each region, validating our method using extended Baryon Oscillation Spectroscopic Survey (eBOSS) EZ-mocks. The mitigation procedure is tested by comparing the angular correlation of the corrected target selection on each photometric region to the angular correlation function obtained using quasars from the Sloan Digital Sky Survey (SDSS) Data Release 16. With our procedure, we recover a similar level of correlation between DESI quasar targets and SDSS quasars in two-thirds of the total footprint and we show that the excess of correlation in the remaining area is due to a stellar contamination that should be removed with DESI spectroscopic data. We derive the Limber parameters in our three imaging regions and compare them to previous measurements from SDSS and the 2dF QSO Redshift Survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

DOE CESER 6 GHz Interference Study

DOE CESER has sponsored Idaho National Laboratory (INL) to conduct an objective and independent study of potential 6 GHz interference from outdoor operation of unlicensed devices in the 6 GHz band on fixed service (FS) microwave communication links operated by electrical sector incumbents in that band. INL is collaborating with University of Notre Dame (UND), Electric Power Research Institute (EPRI), Lockard & White, Southern Company, and AT&T, to gather data with real-world 6 GHz interference experiments and identify (1) the potential for interference from unlicensed devices and (2) the interference necessary to cause harm to the incumbents. In addition to the functional assessment, a security assessment of the FCC mandated Automatic Frequency Coordination (AFC) System to regulate use of unlicensed 6 GHz standard power devices is also being conducted. A major objective is to create a science-based and defensible methodology used to produce the necessary data and to derive objective conclusions. This proven methodology can then be used to produce objective data and conclusions for other spectrum bands with similar incumbent uses including 4.4 – 4.9 GHz and 7.125 – 7.4 GHz, identified in the reconciliation bill that was adopted on July 4, 2025, as well as the National Spectrum Strategy discussions that are ongoing. This report contains 6 GHz field experiments and findings in the following real-world scenarios with commercial unlicensed standard power (SP) 6 GHz devices regulated by Automated Frequency Coordination (AFC): • University of Notre Dame (UND) Stadium with a capacity of 80,000 spectators, where Wi-Fi operating in 6 GHz has been deployed recently • Southern Company 6 GHz FS microwave link between Columbus and Fortson, Georgia Following are the following key findings from this study. 1. The AFC is under-protective of FS when line-of-sight exists along the path centerline. Data collected at Southern’s 6 GHz fixed link site shows significant erosion of as much as 21.4-24.4 dB of under-protection that can lead to potential service degradation under typical operating conditions. This first key finding is most likely the result of erroneous use of the RF propagation model. INL will collaborate with EPRI and the AFC Functional Requirements Working Group to submit a change request to the WinnForum TS-1014 standard towards correct use of the propagation model by the AFC. 2. There is additive interference effect of about 3 dB from nearly equal power interferers measured from simultaneous operation of two SP AP's operating co-channel with the FS receive from different locations along the path. This second key finding should be used to add the impact of additive interference of operation of multiple APs in the same geographical area, to the next generation of AFCs. INL will collaborate with FCC on the need for the AFC to consider additive interference. We also recommend that additional experiments are conducted on 6 GHz spectrum interference to further improve the AFC operation as the number of outdoor Wi-Fi devices continues to increase. These proposed steps and recommendations will make the co-existence of the incumbents and the 6 GHz outdoor Wi-Fi providers possible without any impact on the incumbents with a win-win outcome for all.

6 GHz↗
DOE OSTI2020

The good, the bad, and the ugly: Data-driven load profile discord identification in a large building portfolio

Reducing the overall energy consumption and associated greenhouse gas emissions in the building sector is essential for meeting our future sustainability goals. Recently, smart energy metering facilities have been deployed to enable monitoring of energy consumption data with hourly or subhourly temporal resolution. This unprecedented data collection has created various opportunities for advanced data analytics involving load profiles (e.g., building energy benchmarking programs, building-to-grid integration, and calibration of urban-scale energy models). These applications often need preprocessing steps to detect daily load profile discords, such as: 1) outliers due to system malfunctions (the bad) and 2) irregular energy consumption patterns, such as those resulting from holidays (the ugly) compared to normal consumption patterns (the good). However, current preprocessing methods predominantly focus on filtering using statistical threshold values, which fail to capture the contextual discords of daily profiles. In addition, discord detection algorithms in building research are often aimed at finding individual building-level discords, which are not suitable at a large scale. Thus, here, we develop a method for automated load profile discord identification (ALDI) in a large portfolio of buildings (more than 100 buildings). Specifically, ALDI 1) uses the matrix profile (MP) method to quantify the similarities of daily subsequences in time series meter data, 2) compares daily MP values with typical-day MP distributions using the Kolmogorov-Smirnov test, and 3) identifies daily load profile discords in a large building portfolio. We evaluate ALDI using the metering data of both an academic campus and a residential neighborhood. Our results demonstrate that ALDI efficiently discovers measurement errors by system malfunctions and low energy consumption days in the academic campus portfolio, and it detects unique load shape patterns likely driven by occupant behavior and extreme weather conditions in the residential neighborhood.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

NbZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-zirconium (Nb-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-zirconium (Ta-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ta and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2020

Technical Gaps in Hanford High-Level Waste Solids Settling Behavior and Settling Time Evaluation for Direct Feed High-Level Waste (DFHLW) Operations

Settling of high-level waste (HLW) solids in process vessels is a key conceptual process step in providing HLW feed directly to the Hanford Waste Treatment and Immobilization Plant (WTP) HLW Vitrification Facility. Direct Feed High-Level Waste (DFHLW) is a potential approach to initiating HLW vitrification prior to completing of the WTP Pretreatment Facility. Settling would be used with subsequent supernatant decant to concentrate HLW feed. To support planning for DFHLW, Washington River Protection Solutions (WRPS) requested support from the Pacific Northwest National Laboratory to evaluate the current data set available to predict the time needed for HLW solids to settle, to identify gaps in the understanding and predictive capability of HLW solids waste settling times, and to provide scoping estimates of the potential settling time. Eight technical gaps were identified for predicting settling times and characteristics of the formed sediment layers including: Gap 1: In-Tank Settling Rates Faster than Settling of Laboratory Samples, Gap 2: Effect of Sludge Leaching/Washing on Predicted Settling Times, Gap 3: Predicting Waste Settling from Waste Chemistry (Waste Type), Gap 4: Predicting Waste Settling from Particle Size and Density Distributions (PSDDs), Gap 5: Insufficient Laboratory and In-Tank Settling Data to Represent Hanford Waste, Gap 6: Methods for Real-Time, In-Tank Tracking of Settling, Gap 7: Prediction of Sediment Erosion Resistance as a Function of Settling Time, and Gap 8: Prediction of Sediment Solids Content as a Function of Settling Time. In addition to the data gaps, an overarching observation of the settling rate and settled layer data is the significant variation in behavior. At similar solids concentrations, settling rates can vary by as much as 3 orders of magnitude depending on the source waste tank, and significantly different settling rates are noted between laboratory and in situ tests for the same waste tank. The range of average solids concentration in existing HLW sediment, which may have been quiescent for decades, can vary from less than 7 wt% to greater than 74 wt% solids. The shear strengths (or yield stresses) measured on laboratory samples range from less than 27 Pa to greater than 6400 Pa. These variations can challenge process planning for the application of a settle/decant process for DFHLW. This report describes the significance of the gaps to the settle/decant process and presents uncertainties by way of examples. Potential technical approaches for resolving these gaps are described and the estimated difficulty in resolving these gaps is evaluated. Based on the significance of the gap and the difficulty of resolution, recommendations are made to address specific gaps. Scoping estimates of the potential settling times for DFHLW solids have been made based on the existing data set with its associated gaps. Depending on the process vessel depth and final sediment concentration, substantial fractions of the scoping estimate results for settling times for characterized HLW exceed the 2-week period that has been previously assumed for process planning. There is also significant disparity, potentially greater than a factor of 5000 difference, in the estimated settling times depending on process vessel depth and final sediment solid concentration. This variation in results underscores the significance of the identified gaps and uncertainties with respect to process planning for utilizing settle/decant operations for DFHLW.

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

Analysing the effects of heating and gas puffing in Proto-MPEX helicon and auxiliary heated plasmas

The Material Plasma Exposure Experiment (MPEX) is being constructed at Oak Ridge National Laboratory to investigate critical fusion reactor issues, such as plasma–material interactions (PMIs) under reactor-relevant conditions and time scales. The linear device Proto-MPEX was used as a test bed to address anticipated research and development issues associated with heating scenarios and establish the physics basis for MPEX. The SOLPS-ITER code suite has been applied to understand plasma and neutral transport in Proto-MPEX and to increase confidence in predictive simulations for MPEX. Coupling between COMSOL and SOLPS is performed to implement a 2D electron heating profile of the helicon source. The simulations show reasonable agreement with the experimental data for plasma with helicon and auxiliary electron cyclotron heating (ECH). Both Bohm and constant diffusion (D$_\perp$: 0.5 m 2 s –1 and $\chi_\perp$: 1 m 2 s –1 ) simulations show similar levels of agreement with respect to the sparse experimental data available, assuming a few per cent impurity concentration. ECH significantly increases the target electron temperature, however, the target electron density is reduced compared to helicon-only heated plasmas due to an increase in flow velocity and radial losses. The simulations show that further increasing the ECH power results in an increase in the target electron density due to increased recycling flux and ionization. The results indicate that ECH significantly enhances the target heat flux, with ECH power of 50 kW increasing the target heat flux from 0.4 to 17 MW m –2 . It is found that a small amount of gas puffing (GP) near the target plate can further increase the target heat fluxes at the higher ECH power cases, but the target heat flux is reduced at higher GP conditions due to a significant reduction in electron temperature via radiation. ECH and GP scenarios can generate a higher target flux, facilitating improved PMI studies with more reactor-relevant plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Multi-Species Complex and Standard Metabolomic Samples with Verified Truth Annotations Dataset

This dataset contains 4523251 (~6.35 GB) metabolite-spectra matches following identification with CoreMS. Data were manually curated as true positives, true negatives, or unknowns. Calculations for spectral similarity scores were carried out with two methods for a total of ~12.7 GB (6.35 * 2) of data. They are all .tsv files, though can easily be changed to .txt. The file types are: * human cerebrospinal fluid (CSF), human blood plasma human urine: already published here https://www.nature.com/articles/s41597-021-00894-y, • purchased FAMES standards • fungi species (A. niger, A. nidulans, T. reesei) • soil crust

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

An Efficient 1-D Thermal Stratification Model for Pool-Type Sodium-Cooled Fast Reactors

Investigating thermal stratification in the upper plenum of a sodium fast reactor (SFR) is presently a technology gap in SFR safety analysis. Understanding thermal stratification will promote safe operation of the SFR before its commercial deployment. Stratified layers of liquid sodium with a large vertical temperature gradient could be established in the upper plenum of an SFR during a down-power or a loss-of-flow transient. These stratified layers are unstable and could result in uncertainties for the core safety of an SFR. In order to predict the occurrence of the thermal stratification efficiently, we developed a one-dimensional (1-D) transport model to estimate the temperature profile of the ambient fluid in the upper plenum. This model demands much less computational effort than computational fluid dynamics (CFD) codes and provides calculations with higher fidelity than historical system-level codes. Two flow conditions were considered separately in the current study depending on if in-vessel components are presented in the upper plenum. For the condition where in-vessel components, specifically the upper internal structure, are presented, we assumed that the impinging sodium was evenly dispersed in the ambient fluid within the distance between the bottom of the in-vessel component and the jet inlet surface. For the condition where no in-vessel components are presented, we assumed that the impinging sodium was evenly dispersed in the ambient fluid within the jet length, which was determined through data-driven trainings. The newly developed 1-D model showed similar performance with the CFD model in both cases. However, due to the assumption of flat profiles of the impinging jet axial dispersion rate, nonnegligible discrepancies between the 1-D prediction and the measured data were observed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

ESS-DIVE reporting format for field measurements of soil respiration

This is a new data and metadata reporting format for soil respiration data, based on engagement with a wide range of researchers in the earth and ecological sciences as well as expert advisory panels. Reporting formats are intended to be ‘bottom-up’, community-developed documents and templates that can enable consistent, findable, accessible, interoperable, and reproducible (FAIR) data; our goal was a reporting format that would be relevant and useful for synthesis activities, optimizing data discoverability and usability while not placing an undue burden on data contributors. The format is split into two levels of information—metadata about chambers (used to measure soil respiration in the field), and the measured data themselves. The file data_reportingformat_template.xlsx contains information on how the data table should be structured, including field names, field descriptions, and units; the chammetadata_reportingformat_template.xlsx file contains similar information for the chamber-level metadata. Fundamentally, this format aims to enable FAIR data through an easy-to-use guide for data providers interested in formatting their data in a common way, and provide ‘future-proofing’ capabilities to support reanalyses of soil respiration data, a globally-important component of the carbon cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Evaluated 238 U(n,f) Average Prompt Fission Neutron Multiplicities Including the CGMF Model

This report documents an evaluation of the average prompt fission neutron multiplicity, $\overline{v}_p$, of 238 U from 800 keV to MeV. This evaluation had to be re-done from “scratch” as the input to previous $\overline{v}_p$ evaluations, specifically ENDF/B-VIII.0, was not found. That means that all available experimental data were re-analyzed and uncertainties were re-estimated. The new evaluated 238 U $\overline{v}_p$ based on only experimental data differs distinctly from ENDF/B-VIII.0 $\overline{v}_p$ from 2 to 4.5 MeV, and from 6 to 7 MeV, and is otherwise similar. The difference from 2 to 4.5 MeV stems from the fact that ENDF/B-VIII.0 was tweaked in this energy range to data of Frehaut, while two other, equally trustworthy, data sets would indicate an evaluated 238 U $\overline{v}_p$ that is up to 2% higher. Also, second chance fission in ENDF/B-VIII.0 was smoothed over from 6–7 MeV. Another major difference to ENDF/B-VIII.0 is that one of the evaluations presented here includes model information from the Hauser-Feshbach fission fragment decay code CGMF, while ENDF/B-VIII.0 is based purely on experimental data. CGMF links several fission quantities with each other; $\overline{v}_p$ is predicted by assumptions made on, e.g., pre-neutron emission yields as a function of mass, the total kinetic energy, or spin and parity of fission fragments. This allows to validate the new 238 U $\overline{v}_p$ by using CGMF parameters obtained from fitting to experimental 238 U $\overline{v}_p$ to predict yields as a function of mass, the average total kinetic energy, or the mean energy of the prompt fission neutron spectrum. These model-predicted values can then be compared to experimental and evaluated data. The model-predicted fission-observable values using evaluated parameters obtained here are reasonably close to experimental data for some observables, but are farther away from experimental data related to TKE observables. In addition to that, the evaluated 238 U(n,f) $\overline{v}_p$ shows similar deviations from ENDF/B-VIII.0 as for the evaluation with only experimental data. This difference is expected to lead to changes in simulated effective neutron multiplication factor, $k_{eff}$ of ICSBEP critical assemblies that are sensitive to 238 U in the fast range (BigTen, Flattop, Flattop-Pu). These changes in $k_{eff}$ need to be counter-balanced. Chi-Nu PFNS experimental data are expected to be released in the next few months that might lead to the needed changes in the PFNS. Until then, we hold off in benchmarking the new 238 U(n,f) $\overline{v}_p$ as well as submitting it to ENDF/B-VIII.1. Also, new high-precision 238 U $\overline{v}_p$ are expected to be measured by the CEA in the next two years that will shed further light on question on 238 U $\overline{v}_p$ from 2–4.5 and 6–7 MeV.

238U↗
DOE OSTI2025

The Importance of Being Adaptable: An Exploration of the Power and Limitations of Domain Adaptation for Simulation-Based Inference with Galaxy Clusters

The application of deep machine learning methods in astronomy has exploded in the last decade, with new models showing remarkably improved performance on benchmark tasks. Not nearly enough attention is given to understanding the models' robustness, especially when the test data are systematically different from the training data, or "out of domain." Domain shift poses a significant challenge for simulation-based inference, where models are trained on simulated data but applied to real observational data. In this paper, we explore domain shift and test domain adaptation methods for a specific scientific case: simulation-based inference for estimating galaxy cluster masses from X-ray profiles. We build datasets to mimic simulation-based inference: a training set from the Magneticum simulation, a scatter-augmented training set to capture uncertainties in scaling relations, and a test set derived from the IllustrisTNG simulation. We demonstrate that the Test Set is out of domain in subtle ways that would be difficult to detect without careful analysis. We apply three deep learning methods: a standard neural network (NN), a neural network trained on the scatter-augmented input catalogs, and a Deep Reconstruction-Regression Network (DRRN), a semi-supervised deep model engineered to address domain shift. Although the NN improves results by 17% in the Training Data, it performs 40% worse on the out-of-domain Test Set. Surprisingly, the Scatter-Augmented Neural Network (SANN) performs similarly. While the DRRN is successful in mapping the training and Test Data onto the same latent space, it consistently underperforms compared to a straightforward Yx scaling relation. These results serve as a warning that simulation-based inference must be handled with extreme care, as subtle differences between training simulations and observational data can lead to unforeseen biases creeping into the results.

Ntampaka, Michelle [Baltimore, Space Telescope Sci↗
DOE OSTI2026

Gearbox bearing crack growth prognostics and uncertainty quantification with physics-informed machine learning

This paper introduces the extreme theory of functional connections (X-TFC), a physics-informed machine learning algorithm, and tailors it to estimate the remaining useful life (RUL) of wind turbine gearbox bearings experiencing fatigue crack growth. Unlike purely data-driven methods, X-TFC embeds a physics model, based on Head's theory in this work, into its training objective. The core of X-TFC is a random-projection single-layer neural network trained via an extreme learning machine, which requires only limited damage progression data and solves for output weights with a least-squares optimization algorithm. A composite loss function balances the network's fit to observed degradation data against the residuals of the governing crack growth differential equation, ensuring the learned damage trajectory remains physically plausible. When applied to a vibration-based health-index (HI) dataset measured during the growth of a crack on the inner ring of a high-speed bearing in a wind turbine gearbox (Bechhoefer and Dubé, 2020), X-TFC achieves near-zero prediction bias. Even when trained on only the first 10 %–20 % of the damage progression data, with sufficient physics weighting its predictions remain monotonic and smooth, delivering high prognosability and trendability. To quantify the epistemic uncertainty, we employ a Monte Carlo ensemble of independently initialized X-TFC models trained on noise-perturbed data, which yields confidence intervals around each RUL estimate and captures both model-parameter and epistemic uncertainty. In addition to a vibration-based HI, we demonstrate that the proposed framework can be directly applied to a supervisory control and data acquisition (SCADA) data-based HI (Eftekhari Milani et al., 2026) measured during similar wind turbine gearbox bearing crack faults, preserving its accuracy and interpretability. This extension shows the versatility of our approach, which is applicable to bearings of multiple gearbox manufacturers, models, and ratings using only SCADA data. By integrating domain knowledge with machine learning, X-TFC offers a rapid, reliable tool for crack prognostics. Its adaptability to other bearing failure modes, such as pitch bearing ring cracks, positions X-TFC as a powerful enabler of data-driven, physics-informed asset management in the wind energy sector and beyond.

17 WIND ENERGY↗