Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Programming”

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

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

At least 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of lead on an axisymmetric, single bead blown powder DED overhung geometry

Metal additive manufacturing can be utilized for the near net-shape manufacture of components in a layer wise technique. Traditionally, the buildup process is conducted in the direction of gravity or bottom to top in the vertical orientation; however, with the availability of commercial systems with additive heads or fixturing that can index from vertical, and the need to manufacture larger parts without modifying available systems, the limitation of manufacturing from bottom to top is removed enabling the buildup of larger components printed at an angle from vertical. The effect of gravity on the melt pool and as-deposited component quality when printing off vertical is unknown in literature, especially for axisymmetric components. This understanding is critical for the advancement of manufacturing components of increased size in 4 and 5-axis additive systems. This investigation utilizes blown powder directed energy deposition to evaluate the change in as-printed geometry when the start point is altered in relation to gravity while the part rotates to manufacture an axisymmetric component. The objectives of this work are to determine the impact of the deposition location on the geometric variability on axisymmetric components. This investigation tests the hypothesis that differences in layer height exist due to a change in catchment efficiency when the deposition location is moved to a tangent of the round geometry due to a change in the melt pool dynamics. It was found the ideal deposition location when printing at 26.6-degrees from vertical was at −90-degrees from the top center point of the component, along the tangent, where gravity was pushing the melt pool down, but the rotation of the part was pulling the deposited material towards the top center of the component. This work provides an understanding of layer height stability and catchment efficiency to guide print orientation strategy for high-aspect ratio components. It was found the −90-degree lead deposition had the best layer height stability at 0.6 % as compared to the top-center at 6.3 % and +90-degree deposition location at 9.2 % for the nominal programmed layer height. The change in layer height also effected the final diameter the greatest for the top center deposition location where the diameter diverged by 2.5 % during the overbuild condition and converged by 0.6 %. In conclusion, this finding will increase the manufacturing efficiency of axisymmetric components by increasing the passive stability of the printing process for the successful manufacture of parts without the need for perfectly calibrated manufacturing parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Complete wetting of Pt(111) by nanoscale liquid water films

The melting and wetting of nanoscale crystalline ice films on Pt(111) that are transiently heated above the melting point using nanosecond laser pulses are studied with infrared refection absorption spectroscopy (IRAS) and Kr temperature programmed desorption (TPD). The as-grown crystalline ice films consist of isolated nanoscale ice crystallites embedded in a hydrophobic water monolayer. Upon heating above the melting point, these ice crystallites rapidly melt to form nanoscale droplets of liquid water. Rapid cooling of the system to cryogenic temperatures after each laser pulse quenches the water films and allows them to be interrogated with IRAS, Kr TPD and other ultrahigh vacuum surface science techniques. With each successive heat pulse, these liquid drops spread across the surface until it is entirely covered with multilayer water films after several pulses. These results, which show that nanoscale water films completely wet Pt(111), are in contrast to molecular dynamics simulations predicting partial wetting of nanoscale water drops on a hydrophobic water monolayer. The results provide valuable new insights into the wetting characteristics of nanoscale water films on a clean, well-characterized single crystal surface.

Xu, Yuntao↗
DOE OSTI2017

Magnetic analysis of the Nb$_3$Sn low-beta quadrupole for the high luminosity LHC

As part of the Large Hadron Collider Luminosity upgrade (HiLumi-LHC) program, the US LARP collaboration and CERN are working together to design and build 150 mm aperture Nb 3 Sn quadrupoles for the LHC interaction regions. A first series of 1.5 m long coils were fabricated, assembled and tested in the first short model. This paper presents the magnetic analysis, comparing magnetic field measurements with the expectations and the field quality requirements. The analysis is focused on the geometrical harmonics, iron saturation effect and cold-warm correlation. Three dimensional effects such as the variability of the field harmonics along the magnet axis and the contribution of the coil ends are also discussed. Furthemore, we present the influence of the conductor magnetization and the dynamic effects.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Gaining Insights in Loading Events for Wind Turbine Drivetrain Prognostics

Wind energy is one of the largest sources of renewable energy in the world. To further reduce the operations and maintenance (O&M) costs of wind farms, it is essential to be able to accurately pinpoint the root causes of different failure modes of interest. An example of such a failure mode that is not yet fully understood is white etching cracks (WEC). This can cause the bearing lifetime to be reduced to 5–10% of its design value. Multiple hypotheses are available in literature concerning its cause. To be able to validate or disprove these hypotheses, it is essential to have historic high-frequency measurement data (e.g., load and vibration levels) available. In time, this will allow linking to the history of the turbine operating data with failure data. This paper discusses the dynamic loading on the turbine during certain events (e.g., emergency stops, run-ups, and during normal operating conditions). By combining the number of specific events that each turbine has seen with the severity of each event, it becomes possible to assess which turbines are most likely to show signs of damage.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2016

The Bee Microbiome: Impact on Bee Health and Model for Evolution and Ecology of Host-Microbe Interactions

ABSTRACT As pollinators, bees are cornerstones for terrestrial ecosystem stability and key components in agricultural productivity. All animals, including bees, are associated with a diverse community of microbes, commonly referred to as the microbiome. The bee microbiome is likely to be a crucial factor affecting host health. However, with the exception of a few pathogens, the impacts of most members of the bee microbiome on host health are poorly understood. Further, the evolutionary and ecological forces that shape and change the microbiome are unclear. Here, we discuss recent progress in our understanding of the bee microbiome, and we present challenges associated with its investigation. We conclude that global coordination of research efforts is needed to fully understand the complex and highly dynamic nature of the interplay between the bee microbiome, its host, and the environment. High-throughput sequencing technologies are ideal for exploring complex biological systems, including host-microbe interactions. To maximize their value and to improve assessment of the factors affecting bee health, sequence data should be archived, curated, and analyzed in ways that promote the synthesis of different studies. To this end, the BeeBiome consortium aims to develop an online database which would provide reference sequences, archive metadata, and host analytical resources. The goal would be to support applied and fundamental research on bees and their associated microbes and to provide a collaborative framework for sharing primary data from different research programs, thus furthering our understanding of the bee microbiome and its impact on pollinator health.

Microbiology↗
DOE OSTI2015

Uncertainty quantification of calculated temperatures for advanced gas reactor fuel irradiation experiments

This paper presents the quantification of uncertainty of the calculated temperature data for the Advanced Gas Reactor (AGR) fuel irradiation experiments conducted in the Advanced Test Reactor at Idaho National Laboratory in support of the Advanced Reactor Technology Research and Development program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR tests, the results of the numerical simulations are used in combination with statistical analysis methods to improve qualification of measured data. The temperature simulation data for AGR tests are also used for validation of the fission product transport and fuel performance simulation models. These crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. To quantify the uncertainty of AGR calculated temperatures, this study identifies and analyzes ABAQUS model parameters of potential importance to the AGR predicted fuel temperatures. The selection of input parameters for uncertainty quantification of the AGR calculated temperatures is based on the ranking of their influences on variation of temperature predictions. Thus, selected input parameters include those with high sensitivity and those with large uncertainty. Propagation of model parameter uncertainty and sensitivity is then used to quantify the overall uncertainty of AGR calculated temperatures. Expert judgment is used as the basis to specify the uncertainty range for selected input parameters. The input uncertainties are dynamic accounting for the effect of unplanned events and changes in thermal properties of capsule components over extended exposure to high temperature and fast neutron irradiation. The sensitivity analysis performed in this work went beyond the traditional local sensitivity. Using experimental design, analysis of pairwise interactions of model parameters was performed to establish sufficiency of the first-order (linear) expansion terms in constructing the response surface. To achieve completeness, uncertainty propagation made use of pairwise noise correlations of model parameters. Furthermore, time-dependent sensitivity over the test campaign duration was obtained using an interpolation scheme over the input parameter domain. This allows computation of uncertainty for the predicted peak fuel temperatures and the predicted graphite temperatures at TC locations over the whole AGR irradiation period.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Supplying materials needed for grain growth characterizations of nano-grained UO 2

This activity is supported by the US Nuclear Energy Advanced Modeling and Simulation (NEAMS) Fuels Product Line (FPL) and aims at providing experimental data for the validation of the mesoscale simulation code MARMOT. MARMOT is a mesoscale multiphysics code that predicts the coevolution of microstructure and properties within reactor fuel during its lifetime in the reactor. It is an important component of the Moose-Bison-Marmot (MBM) code suite that has been developed by Idaho National Laboratory (INL) to enable next generation fuel performance modeling capability as part of the NEAMS Program FPL. In order to ensure the accuracy of the microstructure based materials models being developed within the MARMOT code, extensive validation efforts must be carried out. In this report, we summarize our preliminary synchrotron radiation experiments at APS to determine the grain size of nanograin UO 2 . The methodology and experimental setup developed in this experiment can directly apply to the proposed in-situ grain growth measurements. The investigation of the grain growth kinetics was conducted based on isothermal annealing and grain growth characterization as functions of duration and temperature. The kinetic parameters such as activation energy for grain growth for UO 2 with different stoichiometry are obtained and compared with molecular dynamics (MD) simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Synchrotron characterization of nanograined UO 2 grain growth

This activity is supported by the US Nuclear Energy Advanced Modeling and Simulation (NEAMS) Fuels Product Line (FPL) and aims at providing experimental data for the validation of the mesoscale simulation code MARMOT. MARMOT is a mesoscale multiphysics code that predicts the coevolution of microstructure and properties within reactor fuel during its lifetime in the reactor. It is an important component of the Moose-Bison-Marmot (MBM) code suite that has been developed by Idaho National Laboratory (INL) to enable next generation fuel performance modeling capability as part of the NEAMS Program FPL. In order to ensure the accuracy of the microstructure based materials models being developed within the MARMOT code, extensive validation efforts must be carried out. In this report, we summarize our preliminary synchrotron radiation experiments at APS to determine the grain size of nanograin UO 2 . The methodology and experimental setup developed in this experiment can directly apply to the proposed in-situ grain growth measurements. The investigation of the grain growth kinetics was conducted based on isothermal annealing and grain growth characterization as functions of duration and temperature. The kinetic parameters such as activation energy for grain growth for UO 2 with different stoichiometry are obtained and compared with molecular dynamics (MD) simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Lower-length-scale modeling of chemical additions, corrosion and fission products in select MSR base salts

Molten salts play a crucial role in advancing Generation IV nuclear energy technology, with chloride-based salts like NaCl-UCl 3 garnering significant attention due to their distinctive properties. However, the corrosive nature of molten salts causes the dissolution of chromium (Cr), leading to the formation of CrCl 2 and CrCl 3 species in molten chlorides. Moreover, the radioactive decay of nuclear fuel gives rise to fission products, including Cs, Sr, and I chlorides. The first part of this report presents a comprehensive study utilizing ab initio molecular dynamics (AIMD) simulations to investigate the properties of eutectic NaCl-UCl 3 molten salt in the presence of corrosion products (CrCl 2 and CrCl 3 ) and fission product (CsCl). The study focuses on essential structural and thermophysical properties such as density, mixing energy, coordination numbers (CN), and Radial Distribution Functions (RDF) of the salts with varying compositions of corrosion products and fission product, covering a range from 0% to 13.5%. The results offer valuable insights into the behavior of corrosion and fission products in uranium-based molten salts, providing essential data that can be used as input to the MTDB-TC and MSTDB-TP property databases being developed by the NEAMS program. Due to their favorable characteristics such as low melting points, high boiling points, and low costs, MgCl 2 +NaCl+KCl (MNK) eutectic salts have recently attracted attention as high-temperature heat transfer fluids. The incorporation of LiCl into MNK salts can further reduce their melting points and increase their specific heat capacities, which is de- sirable for high-temperature heat transfer applications. The second part of this report presents the development and validation of a new shell-model potential for the MgCl 2 +NaCl+KCl+LiCl system, which captures the polarization of Cl anions. The extensive comparison with experimental data and AIMD simulations demonstrates the accuracy and reliability of the potential. Furthermore, using this potential, we elucidate the intricate network structure in MNK eutectic salts that is formed through polyhedron sharing. This research contributes to a fundamental understanding of the atomic structures and thermophysical properties of multi-component molten salts, which is crucial for future development of heat transfer fluids for applications in molten salt reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanical properties, failure mechanisms, and scaling laws of bicontinuous nanoporous metallic glasses

Molecular dynamics simulations are employed to study the mechanical properties of nanoporous Cu x Zr 1-x metallic glasses (MGs) with five different compositions, x = 0.28, 0.36, 0.50, 0.64, and 0.72, and porosity in the range 0.1 < Φ < 0.7. Results from tensile loading simulations indicate a strong dependence of Young's modulus, E, and Ultimate Tensile Strength (UTS) on porosity and composition. By increasing the porosity from Φ = 0.1 to Φ = 0.7, the topology of the nanoporous MG shifts from closed cell to open-cell bicontinuous. The change in nanoporous topology enables a brittle-to-ductile transition in deformation and failure mechanisms from a single critical shear band to necking and rupture of ligaments. Genetic Programming (GP) is employed to find scaling laws for E and UTS as a function of porosity and composition. Here, a comparison of the GP-derived scaling laws against existing relationships shows that the GP method is able to uncover expressions that can predict accurately both the values of E and UTS in the whole range of porosity and compositions considered.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Automated Work Packages Prototype: Initial Design, Development, and Evaluation. Light Water Reactor Sustainability Program

The goal of the Automated Work Packages (AWP) project is to demonstrate how to enhance work quality, cost management, and nuclear safety through the use of advanced technology. The work described in this report is part of the digital architecture for a highly automated plant project of the technical program plan for advanced instrumentation, information, and control (II&C) systems technologies. This report addresses the DOE Milestone M2LW-15IN0603112: Describe the outcomes of field evaluations/demonstrations of the AWP prototype system and plant surveillance and communication framework requirements at host utilities. A brief background to the need for AWP research is provided, then two human factors field evaluation studies are described. These studies focus on the user experience of conducting a task (in this case a preventive maintenance and a surveillance test) while using an AWP system. The remaining part of the report describes an II&C effort to provide real time status updates to the technician by wireless transfer of equipment indications and a dynamic user interface.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Progress on Grizzly Development for Reactor Pressure Vessels and Reinforced Concrete Structures

The Grizzly code has been under development for the U.S. Department of Energy’s Light Water Reactor Sustainability program to provide predictive tools for the evolution of material degradation in critical light water reactor structural components due to long-term exposure to the environmental conditions of normal reactor operation, and for the effects of this degradation on the safety of these components. This development has primarily targeted reactor pressure vessels and reinforced concrete structures. Work performed during Fiscal Year 2019 added important features for both of these types of structures. For probabilistic fracture mechanics analysis of reactor pressure vessels, capabilities were added to account for the effects of warm prestressing and residual stresses, as well as to model crack initiation, growth, and arrest. For modeling microstructure evolution in reactor pressure vessels, a cluster dynamics model for evolution of Mn-Ni-Si phases has been implemented. For modeling degraded reinforced concrete structures, capabilities to consider the combined effects of damage and creep in the concrete constitutive model and model nonlinear behavior of reinforcing bars has been added. These newly-developed features fill important gaps in Grizzly’s feature set for engineering-scale analysis of degraded reactor pressure vessels and reinforced concrete structures. Grizzly development in both of these areas has reached a point where it is ready for extensive testing on relevant engineering problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

A research program to measure the lifetime of spin polarized fuel

The use of spin polarized fuel could increase the deuterium-tritium (D-T) fusion cross section by a factor of 1.5 and, owing to alpha heating, increase the fusion power by an even larger factor. Issues associated with the use of polarized fuel in a reactor are identified. Theoretically, nuclei remain polarized in a hot fusion plasma. The similarity between the Lorentz force law and the Bloch equations suggests polarization can be preserved despite the rich electromagnetic spectrum present in a magnetic fusion device. The most important depolarization mechanisms can be tested in existing devices. The use of polarized deuterium and 3 He in an experiment avoids the complexities of handling tritium, while encompassing the same nuclear reaction spin-physics, making it a useful proxy to study issues associated with full D-T implementation. 3 He fuel with 65% polarization can be prepared by permeating optically-pumped 3 He into a shell pellet. Dynamically polarized 7 Li-D pellets can achieve 70% vector polarization for the deuterium. Cryogenically-frozen pellets can be injected into fusion facilities by special injectors that minimize depolarizing field gradients. Alternatively, polarized nuclei could be injected as a neutral beam. Once injected, the lifetime of the polarized fuel is monitored through measurements of escaping charged fusion products. Multiple experimental scenarios to measure the polarization lifetime in the DIII-D tokamak and other magnetic-confinement facilities are discussed, followed by outstanding issues that warrant further study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

FAST.Farm Development and Validation of Structural Load Prediction Against Large Eddy Simulations

FAST.Farm is a mid-fidelity engineering tool developed by the National Renewable Energy Laboratory targeted at accurately and efficiently predicting wind turbine power production and structural loading in wind farm settings, including wake interactions between turbines. FAST.Farm is based on several principles of the dynamic wake meandering (DWM) model, but also addresses limitations of previous DWM implementations. Previous FAST.Farm studies have shown the similarities and differences between FAST.Farm and large eddy simulations for rigid turbine cases. The objective of this work is to quantify the ability of FAST.Farm to accurately predict turbine structural response in a small wind farm. This is done by comparing FAST.Farm structural response to SOWFA-OpenFAST results for three laterally-aligned turbines. The purpose of this study is to characterize the similarities and differences between FAST.Farm and a higher fidelity model for predicting turbine structural response in a wind farm for differing atmospheric inflows. Strong statistical agreement was found between FAST.Farm and SOWFA-OpenFAST structural response for the non-waked upstream turbine, and good agreement was found for the downstream turbines for most structural quantities. Higher differences were seen for downstream turbines with low ambient turbulence intensity or yawed turbines, suggesting areas for FAST.Farm wake dynamics modeling improvements. For all cases and turbines, small statistical differences were seen between blade deflections and bending moments, with larger differences for tower-top and tower-base bending moments. Overall, the results establish confidence for applying FAST.Farm to wind farm power and loads analyses and identify areas where further model validations and model improvements should be targeted.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
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 OSTI2016

CDAC Student Report: Summary of LLNL Internship

Multiple objectives motivated me to apply for an internship at LLNL: I wanted to experience the work environment at a national lab, to learn about research and job opportunities at LLNL in particular, and to gain greater experience with code development, particularly within the realm of high performance computing (HPC). This summer I was selected to participate in LLNL's Computational Chemistry and Material Science Summer Institute (CCMS). CCMS is a 10 week program hosted by the Quantum Simulations group leader, Dr. Eric Schwegler. CCMS connects graduate students to mentors at LLNL involved in similar re- search and provides weekly seminars on a broad array of topics from within chemistry and materials science. Dr. Xavier Andrade and Dr. Erik Draeger served as my co-mentors over the summer, and Dr. Andrade continues to mentor me now that CCMS has concluded. Dr. Andrade is a member of the Quantum Simulations group within the Physical and Life Sciences at LLNL, and Dr. Draeger leads the HPC group within the Center for Applied Scientific Computing (CASC). The two have worked together to develop Qb@ll, an open-source first principles molecular dynamics code that was the platform for my summer research project.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Accelerator Real-time Edge AI for Distributed Systems (READS) (Proposal)

Over the last decade, Machine Learning (ML) technologies have slowly made their way into the accelerator community. Rapid advances in recent years in deep learning, particularly reinforcement learning for control system applications and the accessibility of deep learning in embedded hardware, have generated renewed interest and spawned a number of applications. The Fermilab Accelerator Complex, shown in Fig. 1, has provided High Energy Physics (HEP) experiments with proton beams for nearly fifty years. The current focus of the laboratory is its world-class experimental program at the intensity frontier. While increasing beam intensity certainly presents its own challenges, preserving beam size while minimizing beam losses – particles lost through interactions with the beam vacuum pipe – turns out to be, in many ways, the main challenge. The accelerator is controlled via a complex system of hundreds of thousands of devices. Enabling fine tuning and real-time optimization of their parameters using ML methods and stepping beyond experience-based reasoning of human operators are key to the success of future intensity upgrades. Our objective will be to integrate ML into accelerator operations and furthermore, provide an accessible framework, which can also be used by a broad range of other accelerator systems with dynamic tuning needs.

43 PARTICLE ACCELERATORS↗