Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Energy Transition Initiative”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dual functionality of over-lithiated NMC for high energy silicon-based lithium-ion batteries

Owing to their high specific capacity and suitably low operating potential, silicon-based anodes are an attractive alternative to graphite in next-generation lithium-ion batteries. However, silicon anodes suffer from low initial coulombic efficiency and fast capacity decay, limiting their widespread application. Pre-lithiation strategies are highly appealing to compensate for irreversible active lithium loss and to boost the cell energy density. In this work, we maximize the cell energy density by direct pre-lithiation of the NMC (LiNi 0.5 Mn 0.3 Co 0.2 O 2 ) cathode to Li 1+x NMCO 2 without introducing inactive deadweight to either electrode. First, we demonstrate that Li 1+x NMCO 2 can be synthesized chemically, via reaction between NMC and lithium napthalide, and electrochemically. The NMC cathode is tolerant of a one-time over-lithiation up to 60 mA h g NMC -1 , giving capacity retention on par with untreated NMC in half cell electrochemical cycling. Using synchrotron X-ray absorption spectroscopy (ex situ) and diffraction (in situ), we demonstrate that higher amounts of over-lithiation lead to local structure distortion – driven by transition metal reduction to Jahn–Teller active Mn 3+ and Co 2+ – as well as bulk structural hysteresis during over-lithiation and layer “buckling” that increases the amount of lithium extracted from the structure in the charged state. The Li 1+x NMCO 2 with low-to-moderate over-lithiation capacity (23, 46, and 70 mA h g NMC -1 ) is proven to be a highly effective dual-purpose lithium source and cathode material in full cell tests with a commercially relevant Si–graphite anode. These cells show higher capacity, superior cycle life, and improved coulombic efficiencies when compared to those with stoichiometric NMC cathodes. Finally, this study introduces a new and simple method to pre-lithiate layered transition metal oxide cathodes, opening up new possibilities for the development of high energy density lithium-ion batteries with next-generation anodes.

25 ENERGY STORAGE↗
DOE OSTI2020

Quench dynamics of two-component dipolar fermions subject to a quasiperiodic potential

Motivated by recent experiments in fermionic polar gases, we study the nonequilibrium dynamics of two-component dipolar fermions subject to a quasiperiodic potential. We investigate the localization of charge and spin degrees of freedom time evolving with a long-range spin-SU(2) symmetric fermionic Hamiltonian by calculating experimentally accessible dynamical observables. To study the nonequilibrium dynamics, we start the time evolution with two initial states at half-filling: (i) product state with doublons |↑↓0↑↓0↑↓0↑↓0↑↓$\rangle$ and (ii) product state with singlons |↑↓↑↓↑↓↑↓↑↓$\rangle$. We carried out the real-time evolution via the fermionic Hamiltonian using exact diagonalization (ED) and the time-dependent variational principle (TDVP) for finite matrix product states (MPS), within experimentally relevant timescales. For the product state with doublons, we observe a delocalized to localized phase transition varying disorder strengths by monitoring the decay of charge imbalance with time. For the long-range interacting Hamiltonian of our focus, and in the presence of strong enough disorder, starting the time evolution with singlons we find a strong reduction in the spin delocalization, contrary to results of previous studies using the disordered short-range (on-site) Hubbard model with SU(2) symmetry. In addition to the quench dynamics, we also demonstrate the localization of charge and spin in the full energy spectrum of the long-range spin-SU(2) symmetric Hamiltonian by monitoring spin and charge autocorrelation functions. Finally, our predictions for localization of both charge and spin should be observable in ultracold experiments with fermionic dipolar atoms subject to a quasiperiodic potential.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

NbTa_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-tantalum (Nb-Ta). 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. For Ta, p semi-core states as valence were chosen. 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 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 ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and 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. and 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

TaV_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-vanadium (Tanu-V). 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. 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 element V. For Ta, p semi-core states as valence were chosen. 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and Furthmuller, 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. and Furthmuller, 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

NbV_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-vanadium (Nb-V). 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 V. 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and 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. and 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
NASA NTRS2024

Wave Rotor Enhanced Nuclear (WREN) Propulsion: NASA Innovative Advanced Concepts (NIAC) - Phase I Final Report

Nuclear Thermal Propulsion (NTP) is identified as one of the preferred propulsion technologies for manned missions throughout the solar system (NASA MSFC).[1, 2] The state-ofthe-art NTP cycle is based on a solid core Nuclear Engine for Rocket Vehicle Application (NERVA)[3] class technology (Fig. 1) that is envisioned to provide a specific impulse of 900 seconds doubling chemical rocket performance (450 seconds). Even with this impressive increase, the NTP NERVA designs still have issues providing adequate initial to final mass fractions for high ΔV missions.[4] Nuclear Electric Propulsion (NEP) can provide extremely high Isp (2,000 to over 10,000 seconds) but with only low thrust and limits on mass to power ratios. The need for an electric power source also adds the issue of heat rejection in space where thermal energy conversion is at best 30-40% under ideal conditions. NASA Space Technology Mission Directorate (STMD) has recently expressed interest in finding advanced nuclear propulsion technology through the NASA Go:Thrust RFI.[5, 6] A novel Wave Rotor (WR) topping cycle has been proposed for our NASA NIAC concept. It promises to deliver similar thrust as NERVA class NTP propulsion, but with Isp in the 1,200-2,000 second range. Coupled with an NEP cycle, the duty cycle Isp can further be increased (1,800-4,000 seconds) with minimal addition of dry mass. This bimodal design enables fast transit trajectories for manned missions to Mars and revolutionizes the deep space exploration of our solar system.

Nuclear Thermal Propulsion↗
NASA NTRS2018

Three Stage Cool Flame Droplet Burning Behavior of n-Alkane Droplets at Elevated Pressure Conditions: Hot, Warm and Cool Flame

Transient, isolated n-alkane droplet combustion is simulated at elevated pressure for helium-diluent substituted-air mixtures. We report the presence of unique quasi-steady, three-stage burning behavior of large sphero-symmetric n-alkane droplets at these elevated pressures and helium substituted ambient fractions. Upon initiation of reaction, hot-flame diffusive burning of large droplets is initiated that radiatively extinguishes to establish cool flame burning conditions in nitrogen/oxygen “air” at atmospheric and elevated pressures. However, at elevated pressure and moderate helium substitution for nitrogen ( X He > 20%), the initiated cool flame burning proceeds through two distinct, quasi-steady-state, cool flame burning conditions. The classical “Hot flame ” ( ∼1500 K) radiatively extinguishes into a “Warm flame ” burning mode at a moderate maximum reaction zone temperature ( ∼ 970 K), followed by a transition to a lower temperature ( ∼765 K), quasi-steady “Cool flame ” burning condition. The reaction zone (“flame”) temperatures are associated with distinctly different yields in intermediate reaction products within the reaction zones and surrounding near-field, and the flame-standoff ratios characterizing each burning mode progressively decrease. The presence of all three stages first appears with helium substitution near 20%, and the duration of each stage is observed to be strongly dependent on helium substitutions level between 20–60%. For helium substitution greater than 60%, the hot flame extinction is followed by only the lower temperature cool flame burning mode. In addition to the strong coupling between the diffusive loss of both energy and species and the slowly evolving degenerate branching in the low and negative temperature coefficient (NTC) kinetic regimes, the competition between the low-temperature chain branching and intermediate-temperature chain termination reactions control the “Warm” and “Cool ” flame quasi-steady conditions and transitioning dynamics.

Droplet Combustion↗
DOE OSTI2024

Development of a Computational Framework for Multiphysics Multiphase Species Tracking using NEAMS Tools

This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track isotopic species in Molten Salt Reactors (MSRs), with a specific focus on the 91-depletion chain within the Molten Salt Reactor Experiment (MSRE). The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the production, transport, and phase transitions of isotopes under steady-state and transient conditions. The main findings reveal that isotopes such as bromine-91 largely remain in the liquid phase, while others, including krypton-91and yttrium-91, transition to the gas phase, significantly influencing the reactor’s radiological source term. The study also shows that during transients, like a reactivity insertion transient, rapid void formation and the expansion of the liquid-gas interface led to substantial transfers of dissolved isotopes into the gas phase, altering isotope distribution and largely increasing the source term in the off-gas system. Additionally, the research highlights that short-lived isotopes dominate the initial off-gas response during transients, while longer-lived isotopes determine the equilibrium state, underscoring the necessity of dynamic simulations for accurate species tracking and reactor safety analysis. The developed methodology will be applied in the future to the tracking of a larger number of species and introduce other species tracking mechanisms, such as deposition and plating.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Density functional theory investigations into the magnetic ordering of U 3 O 8

Density functional theory (DFT) has been highly successful in supporting experimental materials science; however, a correct electronic ground state is required to realize the full theoretical capacity of DFT. The uranium oxides, α–U 3 O 8 in particular, are simultaneously technologically important materials and theoretically challenging for DFT because the uranium magnetic ground state is not obvious. This is true for both experiment and theory—magnetic susceptibility measurements indicate an antiferromagnetic (AFM) ground state with transitions near 4.2 and 8.0 K, but the ordering itself is not known. Theoretical literature reports are in contradiction, with independent studies finding paramagnetic, ferromagnetic (FM), and AFM states as the lowest energy configuration. However, recent inelastic neutron scattering experiments suggested an uninvestigated magnetic structure with ordering along the [0.5 1 1] plane, motivating a theoretical reinvestigation. Using this insight, we calculated the relative energy of FM and AFM orderings along [0.5 1 1], [0.5 0 0], [0 1 0], and [0 0 1] using noncollinear DFT calculations with spin-orbital coupling. In this work, we found that the [0.5 1 1] AFM structure is lower in energy than FM or AFM orderings along the low Miller index directions. We also investigated polarization of the magnetic moment along each lattice vector and found that polarization along the out-of-plane direction is the energetically preferred orientation for the AFM structures. Additionally, we found in all calculations that moments initially pointing along the in-plane lattice vectors significantly relax until they point along the coordinate between the two crystallographically distinct uranium sites with complex noncollinear magnetic configurations. The new [0.5 1 1] AFM magnetic structure provides an additional path forward toward understanding the electronic structure of α-U 3 O 8 and lends theoretical credibility to recent neutron scattering results.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Cybersecurity Threat Profile for a Connected Lighting System

In anticipation of improved energy performance and cost savings, cities and building owners are increasingly considering “smart lighting initiatives” that aim to convert their collection of simple luminaires (i.e., lighting fixtures) into an intelligent connected lighting system (CLS) capable of remotely monitoring energy consumption and fault conditions, and possibly implementing adaptive lighting schemes. The U.S. Department of Energy (DOE) has set an national goal of tripling the energy efficiency and demand flexibility of the buildings sector by 2030, relative to 2020 levels 1. It is forecast that connected lighting systems can contribute to that goal by delivering 125 TWh of annual energy savings by 2035 2, equivalent to the annual output of 50 typical (500 MW) power plants. However, these energy savings and the DOE goal are put at significant risk if connected technologies are not adopted due to real or perceived cybersecurity concerns. Connected IoT devices such as these have historically been rife with vulnerabilities which sometimes put security considerations secondary to functionality and operability. What are the cybersecurity threats that will impact these systems, as formerly banal luminaires transition into intelligent connected devices that collect information about themselves, their surrounding environment, and possibly us? In this paper we analyze a threat profile performed on a fault-detection use case for streetlights. A threat profile establishes security requirements, justifies security measures, yields actionable controls, and effectively communicates risk to stakeholders. This effort provides critical information for making threat-based decisions to increase security at a reasonable cost, and can effectively be used by development teams, software architects, and managers to make cybersecurity a part of their ongoing culture of awareness, training, and prevention. This leads to more secure systems and better-understood security. On-premise, cloud, and hybrid architectures with different authentication mechanisms were modeled and later categorized using the Microsoft STRIDE framework. An analysis of the recommended controls for each threat was performed to determine which controls could and should be put in place by manufacturers or third-party suppliers, and which controls need to be left up the end-user to implement. Fifty-seven threats were identified. Among our key findings: (1) 65% (37/57) of the threats did not involve the luminaires, but rather the other components needed to communicate with and manage them; (2) 63% (36/57) of the threats could have been mitigated through manufacturer-implemented defensive techniques or “controls”; and (3) 23% (13/57) of the threats were dependent on the network configuration. Recommendations based on the results of this work are made to key stakeholder groups. Notably, lighting technology developers are advised to address all threats that can be reasonably controlled with baked-in technology solutions (e.g., encryption or authentication controls), and employ some form of secure supply chain management and tracking where other parts (e.g., sensors, microprocessors) of a luminaire must also be built and manufactured with the proper security controls in place. Developers should also review threats involving assets not developed in-house to understand how connectivity with other devices will affect their product during system operation and determine if a compensating control for a defense-in-depth strategy will be needed. Finally, those interested in deploying CLS should compare the differences between cloud and on-premise models to determine which is more suitable for their needs and the abilities of their security team.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1987

Low power arcjet thruster pulse ignition

An investigation of the pulse ignition characteristics of a 1 kW class arcjet using an inductive energy storage pulse generator with a pulse width modulated power converter identified several thruster and pulse generator parameters that influence breakdown voltage including pulse generator rate of voltage rise. This work was conducted with an arcjet tested on hydrogen-nitrogen gas mixtures to simulate fully decomposed hydrazine. Over all ranges of thruster and pulser parameters investigated, the mean breakdown voltages varied from 1.4 to 2.7 kV. Ignition tests at elevated thruster temperatures under certain conditions revealed occasional breakdowns to thruster voltages higher than the power converter output voltage. These post breakdown discharges sometimes failed to transition to the lower voltage arc discharge mode and the thruster would not ignite. Under the same conditions, a transition to the arc mode would occur for a subsequent pulse and the thruster would ignite. An automated 11 600 cycle starting and transition to steady state test demonstrated ignition on the first pulse and required application of a second pulse only two times to initiate breakdown.

Sarmiento, Charles J.↗
NASA NTRS1987

Low power arcjet thruster pulse ignition

An investigation of the pulse ignition characteristics of a 1 kW class arcjet using an inductive energy storage pulse generator with a pulse width modulated power converter identified several thruster and pulse generator parameters that influence breakdown voltage including pulse generator rate of voltage rise. This work was conducted with an arcjet tested on hydrogen-nitrogen gas mixtures to simulate fully decomposed hydrazine. Over all ranges of thruster and pulser parameters investigated, the mean breakdown voltages varied from 1.4 to 2.7 kV. Ignition tests at elevated thruster temperatures under certain conditions revealed occasional breakdowns to thruster voltages higher than the power converter output voltage. These post breakdown discharges sometimes failed to transition to the lower voltage arc discharge mode and the thruster would not ignite. Under the same conditions, a transition to the arc mode would occur for a subsequent pulse and the thruster would ignite. An automated 11 600 cycle starting and transition to steady state test demonstrated ignition on the first pulse and required application of a second pulse only two times to initiate breakdown.

Sarmiento, Charles J.↗
NASA NTRS2001

Baseline Testing of the EV Global E-Bike with Ultracapacitors

The NASA John H. Glenn Research Center initiated baseline testing of the EV Global E-Bike SX with ultracapacitors as a way to reduce pollution in urban areas, reduce fossil fuel consumption, and reduce operating costs for transportation systems. The E-Bike provides an inexpensive approach to advance the state of art in hybrid technology in a practical application. The project transfers space technology to terrestrial use via nontraditional partners, and provides power system data valuable for future space applications. The work was done under the Hybrid Power Management (HPM) Program, which includes the Hybrid Electric Transit Bus (HETB). The E-Bike is a state of the art, ground up, hybrid electrical bicycle. Unique features of the vehicle's power system include the use of an efficient, 400 W electric hub motor, and a seven-speed derailleur system that permits operation as fully electric, fully pedal, or a combination of the two. Other innovative features, such as regenerative braking through ultracapacitor energy storage, are planned. Regenerative braking recovers much of the kinetic energy of the vehicle during deceleration. A description of the E-bike, the results of performance testing, and future vehicle development plans are given in this report. The report concludes that the E-Bike provides excellent performance, and that the implementation of ultracapacitors in the power system can provide significant performance improvements.

Eichenberg, Dennis J.↗
NASA NTRS2004

Effects of Kinetic Roughening and Liquid-Liquid Phase Transition on Lysozyme Crystal Growth Velocities

We measured the growth velocities of the (110) face of tetragonal lysozyme, V (centimeters per second), at four different concentrations, c (milligrams per milliliter), as the solution temperature, T (Centigrade), was reduced. For a broad range of T dependent on c, we find that the growth velocities increased as the solution temperature was reduced. The initial increase in V is well characterized by the 2D nucleation model for crystal growth, yielding the magnitude of an effective barrier for growth, gamma(sub s) = 1.2 plus or minus 0.1 x 10(exp -13) erg/molecule. Below certain temperatures, T(sub cr), dependent on c, however, a kinetic roughening hypothesis that considers the continuous addition of molecules anywhere on the crystal surface better describes the observed growth velocities. The application of the continuous growth model, up to the solution cloud-point temperatures, T(sub cl), enabled the determinations of the crossover concentration, c(sub r), from estimated values of T(sub cr). For all conditions presented, we find that the crossover from growth by 2D nucleation to continuous addition occurs at a supersaturation, sigma (sub c), = 2.0 plus or minus 0.1. Moreover, we find the energy barrier for the continuous addition, E(sub c), within the temperature range T(sub cl) less than T less than T less than T (sub cr), to be 6 plus or minus 1 x 10(exp -13) erg/molecule. Further reduction of T below approximately 2-3 C of T(sub cl), also revealed a rapid slowing of crystal growth velocities. From quasi-elastic light scattering investigations, we find that the rapid diminishment of crystal growth velocities can be accounted for by the phase behavior of lysozyme solutions. Namely, we find the reversible formation of dense fluid proto-droplets comprised of lysozyme molecules to occur below approximately 0.3 C of T(sub cl). Hence, the rapid slowing of growth velocities may occur as a result of the sudden depletion of "mobile" molecules within crystal growth solutions as dense fluid proto-droplets form.

Gorti, Sridhar↗
DOE OSTI2024

Early photometric and spectroscopic observations of the extraordinarily bright INTEGRAL-detected GRB 221009A

Context. GRB 221009A, initially detected as an X-ray transient by Swift, was later revealed to have triggered the Fermi satellite about an hour earlier, marking it as a post-peak observation of the event’s emission. This GRB distinguished itself as the brightest ever recorded, presenting an unparalleled opportunity to probe the complexities of GRB physics. The unprecedented brightness, however, challenged observation efforts, as it led to the saturation of several high-energy instruments.Aims. Our study seeks to investigate the nature of the INTEGRAL-detected GRB 221009A and elucidate the environmental conditions conducive to these exceptionally powerful bursts. Moreover, we aim to understand the fundamental physics illuminated by the detection of teraelectronvolt (TeV) photons emitted by GRB 221009A.Methods. We conducted detailed analyses of early photometric and spectroscopic observations that span from the Fermi trigger through to the initial days following the prompt emission phase in order to characterize GRB 221009A’s afterglow, and we complemented these analyses with a comparative study.Results. Our findings from analyzing INTEGRAL data confirm GRB 221009A as the most energetic event observed to date. Early optical observations during the prompt phase negate the presence of bright optical emissions with internal or external shock origins. Spectroscopic analyses enabled us to measure GRB 221009A’s distance and line-of-sight properties. The afterglow’s temporal and spectral analysis suggests prolonged activity of the central engine and a transition in the circumburst medium’s density. Finally, we discuss the implications for fundamental physics of detecting photons as energetic as 18 TeV from GRB 221009A.Conclusions. Early optical observations have proven invaluable for distinguishing between the potential origins of optical emissions in GRB 221009A, underscoring their utility in GRB physics studies. However, the rarity of such data underscores the need for dedicated telescopes capable of synchronous multiwavelength observations. Additionally, our analysis suggests that the host galaxies of TeV GRBs share commonalities with those of long and short GRBs. Expanding the sample of TeV GRBs could further solidify these findings.Key words: techniques: photometric / techniques: spectroscopic / gamma-ray burst: general / gamma-ray burst: individual: GRB 221009A

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Molecular View of Methane Activation on Ni(111) through Enhanced Sampling and Machine Learning

A combination of machine learned interatomic potentials (MLIPs) and enhanced sampling simulations is used to investigate the activation of methane on a Ni(111) surface. The work entails the development and iterative refinement of MLIPs, initially trained on a dataset constructed via ab initio molecular dynamics (AIMD) simulations, supplemented by adaptive biasing forces, to enrich the sampling of catalytically relevant configurations. Our results reveal that by incorporating collective variables that capture the behavior of the reactant molecule, as well as additional frames that describe the dynamic response of the catalytic surface, it is possible to enhance considerably the accuracy of predicted energies and forces. By employing enhanced sampling schemes in the refinement of the MLIP, we systematically explore the potential energy surface, leading to a refined MLIP capable of predicting DFT-level energies and forces and replicating key geometric characteristics of the catalytic system. The resulting free energy landscapes at several temperatures provide a detailed view of the thermodynamics and dynamics of methane activation. Specifically, as methane approaches and dissociates on the catalytic surface, the process involves the dynamic interplay of CH 4 and the Ni catalyst that includes both enthalpic and entropic contributions. The progression towards the transition state involves an CH 4 moiety that is increasingly restrained in its ability to rotate or translate, while the stage following the transition state is characterized by a notable rise of the Ni atom that interacts with the cleaved C–H bond. Furthermore, this leads to an increase in the mobility of the adsorbed species, a feature that becomes more pronounced at higher temperatures.

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

Translational diffusion in supercooled water at and near the glass transition temperature—136 K

The properties of amorphous solid water at and near the calorimetric glass transition temperature, T g , of 136 K have been debated for years. One hypothesis is that water turns into a “true” liquid at T g (i.e., it becomes ergodic) and exhibits all the characteristics of an ergodic liquid, including translational diffusion. A competing hypothesis is that only rotational motion becomes active at T g , while the “real” glass transition in water is at a considerably higher temperature. To address this dispute, we have investigated the diffusive mixing in nanoscale water films, with thicknesses up to ∼100 nm, using infrared (IR) spectroscopy. The experiments used films that were composed of at least 90% H 2 O with D 2 O making up the balance and were conducted under conditions where H/D exchange was essentially eliminated. Because the IR spectra of multilayer D 2 O films (e.g., thicknesses of ∼3–6 nm) embedded within thick H 2 O films are distinct from the spectrum of isolated D 2 O molecules within H 2 O, the diffusive mixing of (initially) isotopically layered water films could be followed as a function of annealing time and temperature. The results show that water films with total thicknesses ranging from ∼20 to 100 nm diffusively mixed prior to crystallization for temperatures between 120 and 144 K. The translational diffusion had an Arrhenius temperature dependence with an activation energy of 40.8 ± 3.5 kJ/mol, which indicates that water at and near T g is a strong liquid. The measured diffusion coefficient at 136 K is 6.25 ± 1.4 × 10 −21 m 2 /s.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗