Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model 3”

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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Feature Extraction for Pipeline Defects Inspection Based Upon Distributed Acoustic Fiber Optic Sensing Data

Fiber-optic distributed acoustic sensing (DAS) is becoming an increasingly important tool for real-time monitoring of energy and civil infrastructure structural health such as pipelines. We present a systematic theoretical study of the potential for DAS to be directly coupled with guided ultrasonic waves typically used in conventional acoustic non-destructive evaluation (NDE) methods for real-time pipeline health monitoring. We are referring to this innovative new NDE technique as ultrasonic guided wave and optical fiber sensor fusion. In the practical application of DAS coupled with guided ultrasonic waves, the structural design of (1) the specific guided waves excited, (2) the physical installation of the acoustic transducers and the fiber optic sensors, and (3) the functional performance specifications (gauge length, sensitivity, Etc.) of fiber optic DAS have an important influence on overall capabilities of the monitoring system. Meanwhile, physics-based analysis of acoustic waves is still a challenge due to the complex nature of the Lamb wave when it propagates, scatters, and disperses in the presence of structural defects. In this work, we simulate carbon steel pipes relevant for oil and gas pipeline applications with diameters of approximately 6-12” and wall thickness of 0.5” as the objects to be monitored. By establishing and implementing these capabilities, we seek to pursue an in-depth study on structural parameter optimization of DAS network, measurement range, and signal processing with an ultimate goal of increasing the sensitivity and efficacy of DAS to defect identification for various modes of corrosion expected in practice. To study the characteristics of scattered acoustic waves and performance of DAS for defect identification, we simulated the response of DAS for multiple pipe structures, defect types, and DAS sensor network configuration using finite element software Ansys, then the properties of signal response are extracted to construct defect-sensitive features. The raw data simulated, and the associated features extracted can ultimately be utilized as annotated training data to benchmark various designs for DAS applications, guided acoustic excitation sources, and learning model parameters to enhance early detection of potentially problematic defects.

Pipeline Defects Inspection, Fiber-optic sensors, ↗
DOE OSTI2020

Analysis of Respirator Cartridge Performance Testing on Hanford Tank A-101

Washington River Protection Solutions (WRPS) conducted tests on two types of chemical cartridges for use in air purifying respirators (APR) to determine the period of time that the cartridges would provide adequate performance1 to protect workers when exposed to a mixture of Chemicals of Potential Concern (COPCs) from vapors emanating from the headspace of tank A-101 on the Hanford Site. The Occupational Safety and Health Administration (OSHA) identifies cartridge testing as a valid approach for establishing a cartridge change schedules.[3] Testing is commonly applied in situations where mixtures of COPCs exist, and where other approaches, such as manufacturer recommendations and modeling, are less reliable. The tests were designed and conducted to ensure measurement and/or control of the key variables OSHA identified as important to estimate cartridge service life, including temperature, humidity, COPC concentration, breathing rate, and cartridge adsorption capacity. Testing was conducted over a period from July 22-24, 2016, using headspace vapors from Hanford tank A-101 under static conditions2 fed to a respirator cartridge test stand developed by WRPS in collaboration with HiLine Engineering (Richland, Washington). Multipurpose respirator cartridges, SCOTT 7422-SD1 and 7422-SC1 (SCOTT Safety, Monroe, North Carolina) were assessed on separate days with A-101 headspace vapors. Sample media (sorbent tubes) were used to collect samples of the vapor stream entering and exiting the respirator cartridge, and were subsequently analyzed for COPC concentrations. Pacific Northwest National Laboratory was tasked with conducting an independent analysis of the analytical results and making recommendations based on the results for respiratory cartridge performance and service life. Key conclusions from the assessment of the 59 COPCs in this study are described below: • Based on d measured cartridge inlet vapor concentrations from tank A-101 headspace, both ammonia and N-Nitrosodimethylamine (NDMA) exceeded their Occupational Exposure Limits (OEL).3 These measurements were generally consistent with maximum A-101 headspace measurements previously obtained for these two compounds. • Ammonia had respirator cartridge outlet concentrations that exceeded 10% of the OEL for both cartridges tested, indicating breakthrough for each. For the SCOTT 7422-SD1 cartridge, ammonia appeared to breakthrough above 10% of its OEL after 2 hours. For the SCOTT 7422-SC1 cartridge, ammonia breakthrough appeared to occur after 8 hours. 1 “Adequate performance” refers to the breakthrough criteria used in this analysis—that is, sustained cartridge outlet concentrations above 10% of the OEL. For some COPCs, an alternate threshold has been applied when necessary due to higher detection limit (DL)/reporting limit (RL) values for specific compounds. Ultimately, Industrial Hygiene professionals will use these results along with specific hazard assessments to determine service life and cartridge use that provides the necessary performance. 2 These tests were conducted under static conditions absent waste-disturbing activities in the subject tank or tank farm. 3 Occupational Exposure Limits accepted for Hanford Tank Farm use are based on OELs established by a U.S. governmental agency or national professional organization (e.g., OSHA, National Institute for Occupational Safety and Health, American Conference of Governmental Industrial Hygienists), or if no U.S. OEL exists, standard toxicological practices are applied to develop OELs using non-U.S. exposure limits, other established OELs for chemical surrogates when available, or other standard procedures. The OEL for NDMA was established in 2005 based on the MAK (Maximale Arbeitsplatzkonzentration) Commission standard adopted in Europe. iv • Despite respirator inlet measurements for NDMA that were in excess of its OEL, all corresponding outlet measurements, from both respirator cartridges, were below analytical DL 1 for the

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Synthetic Design of Microorganisms for Lignin Fuel

Despite the proof-of-the-concept for lignin bioconversion at the time when the project was first funded, the conversion efficiency and titer were extremely low for bioconversion of lignin into lipid. The low conversion efficiency is due the inherent recalcitrance of lignin and limited understanding of lignin chemistry. In nature, lignin is a complex aromatic heteropolymer composed of phenylpropane units cross-linked via a variety of chemically stable bonds, which makes it particularly difficult for microbial degradation. In addition to the recalcitrance, lignin as a biorefinery waste is also highly heterogeneous in terms of chemical linkage abundance, functional group distribution, molecular weight and uniformity. In order to achieve better lignin bioconversion, it is critical to understand how lignin chemistry impacts the conversion performance. Throughout the project, we have clearly found that ‘not all lignin are created equal’, in particular, for the lignin out of different biochemical conversion processes. In addition, we have clearly shown that lignin chemistry defines its functionality and processiblity for bioconversion. In particular, the S/G/H ratio, functional groups, chemical linkages and molecular weight defines the reactivity of lignin during bioconversion. Throughout the project, we have advanced the understanding for the impact of lignin chemistry on reactivity, developed a series of lignin 3 fractionation technologies to improve lignin processibility, and ultimately advanced biorefinery design holistically for simultaneous optimization of carbohydrate release and lignin reactivity. In addition to lignin reactivity, another important consideration for lignin bioconversion is the available microbial strains for bioconversion. Natural biomass utilization systems evolved mechanisms for lignin degradation, often with at least three sequential steps: lignin depolymerization into aromatic compounds, subsequent degradation of aromatic compounds into common intermediates like acetyl-CoA, and eventually the conversion of central metabolites into target bioproducts like lipid and PHA. Extensive work has gone into understanding the mechanisms of lignin depolymerization in model systems such as white rot fungi and termites. Despite interest, no fungus-based commercial process for lignin depolymerization has been developed, due in large part to the practical challenges of fungal genetic manipulations. In addition to white rot fungi, a few members of the Actinobacteria, α-Proteobacteria, and γ-Proteobacteria (including Rhodococcus putida) have lignin degradation activity in vitro. However, most of lignindegrading bacteria has very low lignin conversion efficiency, in particular, for lignin depolymerization. The project focused on addressing the issue with systems biology-guided microbial design to improve lignin bioconversion. The project first established that laccase enzyme can synergize with bacteria to promote lignin depolymerization, and further utilized synthetic biology approach to designed functional modules for lignin depolymerization, aromatic compound utilization, and bioproduct accumulation. The integration of biorefinery design, synthetic biology, and fermentation optimization eventually allowed us to achieve the target yield of 10g/L lipid production from lignin-containing biorefinery waste. In this summary, we hereby briefly review the major discoveries from the project with reference to different sections of this report.

09 BIOMASS FUELS↗
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↗
NASA NTRSDate not supplied

Exploration Exercise System (EES) Physiology

Exploration class missions will be required to have an exercise device that is lightweight, has a small footprint, and is capable of providing enough physical stimulus and exercise variability to be an effective countermeasure against muscle, bone, aerobic fitness, and sensorimotor loss that results from the microgravity environment. Ground studies, inflight studies, and results analyses are required to inform Moon2Mars Design Reference Missions (DRMs) and vehicle designs. This activity encompasses two related projects, an evaluation study called Zero T2 and a requirements development for an exploration treadmill (ET). Exploration class missions’ mass/power/volume restrictions have resulted in the development of exercise devices that are motorized and flywheel-based to provide both aerobic and resistive training on one platform. These devices provide a variety of full body resistance exercise options as well as rowing and cycling for aerobic exercise, but do not provide ambulatory exercise via a treadmill. Because a treadmill has been available for use on the ISS since 2001, we do not understand the efficacy of exploration exercise modalities on muscle performance, aerobic fitness, bone health, or sensorimotor performance. A retrospective study was performed to quantify the association between total exercise (including specific contributions of treadmill exercise) and functional performance upon landing. Due to the relatively small variation in treadmill usage inflight, there is a need to do a controlled study where inflight crewmembers will not use the ISS treadmill (called T2) for the duration of their flights. There will be 3 arms to the Zero T2 study: 1) participants use all currently available exercise platforms aboard ISS (control group), 2) participants use only ARED and CEVIS, and 3) participants use only European Enhanced Exploration Device (E4D). (Please see the TechPort Entry for Exploration Exercise System Development for further detail on the E4D.) This study started in Oct 2020 and will continue through Sep 2027. In a parallel effort, a team will be developing requirements for an exploration treadmill. This project will be broken into two phases. Phase 1 will develop and implement energy prediction models based on available in-flight exercise data to quantify energy expenditure during 0g treadmill running. Phase 2 will then use the available biomechanical and metabolic evidence from Phase 1 to determine the hardware specification requirements essential to providing the desired exercise stimulus. Subject Matter Experts (SMEs) will then leverage the Standardized Process for Evaluating Exercise Devices (SPEED) developed by the team to 1) assess available technologies for candidate treadmills which meet defined requirements, 2) perform hands-on evaluations of top candidates, 3) provide evaluation outcomes and recommendations to developers. This project will begin in 2025 and continue through Sep 2030.

Kent Lawrence Kalogera↗
DOE OSTI2015

Aircraft measurements of BrO, IO, glyoxal, NO 2 , H 2 O, O 2 –O 2 and aerosol extinction profiles in the tropics: comparison with aircraft-/ship-based in situ and lidar measurements

Here, tropospheric chemistry of halogens and organic carbon over tropical oceans modifies ozone and atmospheric aerosols, yet atmospheric models remain largely untested for lack of vertically resolved measurements of bromine monoxide (BrO), iodine monoxide (IO) and small oxygenated hydrocarbons like glyoxal (CHOCHO) in the tropical troposphere. BrO, IO, glyoxal, nitrogen dioxide (NO 2 ), water vapor (H 2 O) and O 2 –O 2 collision complexes (O 4 ) were measured by the University of Colorado Airborne Multi-AXis Differential Optical Absorption Spectroscopy (CU AMAX-DOAS) instrument, aerosol extinction by high spectral resolution lidar (HSRL), in situ aerosol size distributions by an ultra high sensitivity aerosol spectrometer (UHSAS) and in situ H 2 O by vertical-cavity surface-emitting laser (VCSEL) hygrometer. Data are presented from two research flights (RF12, RF17) aboard the National Science Foundation/National Center for Atmospheric Research Gulfstream V aircraft over the tropical Eastern Pacific Ocean (tEPO) as part of the "Tropical Ocean tRoposphere Exchange of Reactive halogens and Oxygenated hydrocarbons" (TORERO) project (January/February 2012). We assess the accuracy of O 4 slant column density (SCD) measurements in the presence and absence of aerosols. Our O 4 -inferred aerosol extinction profiles at 477 nm agree within 6% with HSRL in the boundary layer and closely resemble the renormalized profile shape of Mie calculations constrained by UHSAS at low (sub-Rayleigh) aerosol extinction in the free troposphere. CU AMAX-DOAS provides a flexible choice of geometry, which we exploit to minimize the SCD in the reference spectrum (SCD REF , maximize signal-to-noise ratio) and to test the robustness of BrO, IO and glyoxal differential SCDs. The RF12 case study was conducted in pristine marine and free tropospheric air. The RF17 case study was conducted above the NOAA RV Ka'imimoana (TORERO cruise, KA-12-01) and provides independent validation data from ship-based in situ cavity-enhanced DOAS and MAX-DOAS. Inside the marine boundary layer (MBL) no BrO was detected (smaller than 0.5 pptv), and 0.2–0.55 pptv IO and 32–36 pptv glyoxal were observed. The near-surface concentrations agree within 30% (IO) and 10% (glyoxal) between ship and aircraft. The BrO concentration strongly increased with altitude to 3.0 pptv at 14.5 km (RF12, 9.1 to 8.6° N; 101.2 to 97.4° W). At 14.5 km, 5–10 pptv NO 2 agree with model predictions and demonstrate good control over separating tropospheric from stratospheric absorbers (NO 2 and BrO). Our profile retrievals have 12–20 degrees of freedom (DoF) and up to 500 m vertical resolution. The tropospheric BrO vertical column density (VCD) was 1.5 × 10 13 molec cm –2 (RF12) and at least 0.5 × 10 13 molec cm –2 (RF17, 0–10 km, lower limit). Tropospheric IO VCDs correspond to 2.1 × 10 12 molec cm –2 (RF12) and 2.5 × 10 12 molec cm –2 (RF17) and glyoxal VCDs of 2.6 × 10 14 molec cm –2 (RF12) and 2.7 × 10 14 molec cm –2 (RF17). Surprisingly, essentially all BrO as well as the dominant IO and glyoxal VCD fraction was located above 2 km (IO: 58 ± 5%, 0.1–0.2 pptv; glyoxal: 52 ± 5%, 3–20 pptv). To our knowledge there are no previous vertically resolved measurements of BrO and glyoxal from aircraft in the tropical free troposphere. The atmospheric implications are briefly discussed. Future studies are necessary to better understand the sources and impacts of free tropospheric halogens and oxygenated hydrocarbons on tropospheric ozone, aerosols, mercury oxidation and the oxidation capacity of the atmosphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Developing Methods for Exercise System Kinematic Tracking

BACKGROUND How to quantify the load and forces produced by exercise equipment and their Vibration Isolation and Stabilization (VIS) platforms in-flight is an active area of investigation. Kinematic tracking paired with system modeling can provide insights as well as verification and validation of simulations used for system design and development. Traditional motion capture methods can require significant cost in equipment procurement and crew-time, but newer lessons learned can be leveraged [1]. The VIS systems of current and future exercise hardware on the International Space Station (ISS) such as the Cycle Ergometer with Vibration Isolation System (CEVIS) and the European Enhanced Exploration Exercise Device (E4D) are not currently outfitted with IMUs or similar measurement devices. Video-based methods would enable use of multi-purpose, crew-familiar flight equipment. An initial exploration of video-based solutions was performed utilizing 2-camera video from crew cycling on Teal-CEVIS on the ISS. METHODS AND RESULTS Our group has scoped a variety of video-based object tracking methods. To date, we have primarily investigated computer vision toolkits such as open CV. Techniques explored include key-point detection, background subtraction, region-of interest tracking, color-based tracking, tag masking and tracking, and corner detection. Although object-tracking and 6D pose estimation is a rich field, space applications are a unique problem that are challenging for existing software and toolkits. The majority of the existing object-tracking applications involve vehicles/pedestrians and household objects with simple backgrounds. We have identified the following features which pose particular challenges for on-station exercise equipment tracking: 1. Busy and visually cluttered background 2. Low-textured tracking object with relatively small motions 3. Occlusions and motion by human subject and loose, floating objects 4. Limited number of video cameras with no fixed global references 5. Limited ability to add tags, markers, or visual references to the tracking object 6. Lack of training data for Machine Learning (ML) algorithms CONCLUSION We will summarize the efficacy of techniques tested for a ground mock-trial and the on-station exercise trial. It is likely that human-in-loop feedback or a conglomerate of methods is required. ML-based methods, like those implemented for human body tracking [2], may still be a viable option, but more training data and validation is needed.

L Nilsson↗
DOE OSTI2018

Numerical Modelling of Radionuclide Migration in the Context of a Near-Surface LILW Disposal Facility - 18618

The reference solution for the management of low and intermediate level short-lived waste (LILW) in France consists in their emplacement in concrete vaults built at the near-surface. The corresponding disposal facility under operation by the French waste management organization (Andra), is located in North-eastern France and planned to be closed in 2060. The safety case developed by Andra has to show the evidence that, amongst other aspects, the radionuclide transfer through the confinement barriers will occur in such a way that the human and environmental exposure will be as low as reasonably achievable (ALARA) and in accordance with the dose constraints. For this purpose, modelling is used to integrate uncertainties about the characteristics of concrete components, waste and water table and their evolution in time and space. At last, modelling is finally used to investigate the radiological impact of specific scenarios in which the failure of safety barriers occurs. In the frame of its technical review of the safety assessment, the French Institute for Radiological Protection and Nuclear Safety (IRSN) also performs such kind of numerical modelling of radionuclide fate from the waste to the hydrogeological outlets, in order to check whether the assessment performed by Andra takes into account sufficiently conservative assumptions. The present work aims to evaluate the role of components of the vaults on radionuclide transport during the post-closure period over a few thousand years. For this purpose, IRSN carries out numerical calculations with MELODIE software, which simulates water flow and transport of species, in order to simulate radionuclide migration from waste packages emplaced in a concrete vault to the surrounding groundwater. The developed transport model focuses on 3-H radionuclide and integrates Andra's main assumptions (release model from the waste packages, chemical and hydraulic properties of the vault components, shrinkage and cracks in concrete components). An adaptive mesh method, based on an a posteriori error estimator, was used to model the cracks and the consideration of a strong heterogeneity in the hydraulic and transport parameters values in the different components of the model. This method shows an improvement of results of about 40 times, whereas the mesh size was increased by only 2%. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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↗
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↗
NASA NTRS1999

Three-Dimensional, Transgenic Cell Models to Quantify Space Genotoxic Effects

The space environment contains radiation and chemical agents known to be mutagenic and carcinogenic to humans. Additionally, microgravity is a complicating factor that may modify or synergize induced genotoxic effects. Most in vitro models fail to use human cells (making risk extrapolation to humans more difficult), overlook the dynamic effect of tissue intercellular interactions on genotoxic damage, and lack the sensitivity required to measure low-dose effects. Currently a need exists for a model test system that simulates cellular interactions present in tissue, and can be used to quantify genotoxic damage induced by low levels of radiation and chemicals, and extrapolate assessed risk to humans. A state-of-the-art, three-dimensional, multicellular tissue equivalent cell culture model will be presented. It consists of mammalian cells genetically engineered to contain multiple copies of defined target genes for genotoxic assessment,. NASA-designed bioreactors were used to coculture mammalian cells into spheroids, The cells used were human mammary epithelial cells (H184135) and Stratagene's (Austin, Texas) Big Blue(TM) Rat 2 lambda fibroblasts. The fibroblasts were genetically engineered to contain -a high-density target gene for mutagenesis (60 copies of lacl/LacZ per cell). Tissue equivalent spheroids were routinely produced by inoculation of 2 to 7 X 10(exp 5) fibroblasts with Cytodex 3 beads (150 micrometers in diameter). at a 20:1 cell:bead ratio, into 50-ml HARV bioreactors (Synthecon, Inc.). Fibroblasts were cultured for 5 days, an equivalent number of epithelial cells added, and the fibroblast/epithelial cell coculture continued for 21 days. Three-dimensional spheroids with diameters ranging from 400 to 600 micrometers were obtained. Histological and immunohistochemical Characterization revealed i) both cell types present in the spheroids, with fibroblasts located primarily in the center, surrounded by epithelial cells; ii) synthesis of extracellular matrix; and iii,, mitotic cells located throughout the spheroids. Spheroidal integrity and cell viability were retained for the 30-day test period after removal of spheroids from the bioreactor. Potential utility of this three-dimensional, transgenic model for genotoxicity was initially assessed by exposure of spheroids to 0-2 Gy neon at dose rates of 0.3 to 1.5 Gy/min (National Institute of Radiological Sciences, Chiba, Japan). Quantification of mutation at the lacl gene revealed a linear dose response for mutation induction. Limited sequencing analysis of mutant clones revealed higher frequencies of deletions and multiple base sequence changes with increasing dose. These results suggest that our three-dimensional, transgenic model is applicable to a wide variety of studies involving the quantification, identification, and characterization of genotoxicity incurred in space and on Earth. This model uniquely allows investigation of the interaction of relevant factors, namely cell-to-cell interactions and the mechanistic interaction of microgravity with radiation insults and DNA repair. Using this three-dimensional model will allow us to obtain dual genotoxic information (i.e., mutation rate plus chromosome aberration data) from the same system so that one endpoint can be used to reference the other, thereby increasing the fidelity of the data set. Moreover, the tissue-equivalent nature of the three-dimensional model provides high confidence for relevance of risk assessment, i.e., the establishment of quality factors directly applicable to the microgravity environment.

Gonda, S. R.↗
NASA NTRS2011

Z-Pinch Magneto-Inertial Fusion Propulsion Engine Design Concept

Fusion-based nuclear propulsion has the potential to enable fast interplanetary transportation. Due to the great distances between the planets of our solar system and the harmful radiation environment of interplanetary space, high specific impulse (Isp) propulsion in vehicles with high payload mass fractions must be developed to provide practical and safe vehicles for human spaceflight missions. Magneto-Inertial Fusion (MIF) is an approach which has been shown to potentially lead to a low cost, small fusion reactor/engine assembly (1). The Z-Pinch dense plasma focus method is an MIF concept in which a column of gas is compressed to thermonuclear conditions by an estimated axial current of approximately 100 MA. Recent advancements in experiments and the theoretical understanding of this concept suggest favorable scaling of fusion power output yield as I(sup 4) (2). The magnetic field resulting from the large current compresses the plasma to fusion conditions, and this is repeated over short timescales (10(exp -6) sec). This plasma formation is widely used in the field of Nuclear Weapons Effects (NWE) testing in the defense industry, as well as in fusion energy research. There is a wealth of literature characterizing Z-Pinch physics and existing models (3-5). In order to be useful in engineering analysis, a simplified Z-Pinch fusion thermodynamic model was developed to determine the quantity of plasma, plasma temperature, rate of expansion, energy production, etc. to calculate the parameters that characterize a propulsion system. The amount of nuclear fuel per pulse, mixture ratio of the D-T and nozzle liner propellant, and assumptions about the efficiency of the engine, enabled the sizing of the propulsion system and resulted in an estimate of the thrust and Isp of a Z-Pinch fusion propulsion system for the concept vehicle. MIF requires a magnetic nozzle to contain and direct the nuclear pulses, as well as a robust structure and radiation shielding. The structure, configuration, and materials of the nozzle must meet many severe requirements. The configuration would focus, in a conical manner, the Deuterium-Tritium (D-T) fuel and Lithium-6/7 liner fluid to meet at a specific point that acts as a cathode so the Li-6 can serve as a current return path to complete the circuit. In addition to serving as a current return path, the Li liner also serves as a radiation shield. The advantage to this configuration is the reaction between neutrons and Li-6 results in the production of additional Tritium, thus adding further fuel to the fusion reaction and boosting the energy output. To understand the applicability of Z-Pinch propulsion to interplanetary travel, it is necessary to design a concept vehicle that uses it. The propulsion system significantly impacts the design of the electrical, thermal control, avionics, radiation shielding, and structural subsystems of a vehicle. The design reference mission is the transport of crew and cargo to Mars and back, with the intention that the vehicle be reused for other missions. Several aspects of this vehicle are based on a previous crewed fusion vehicle study called Human Outer Planet Exploration (HOPE), which employed a Magnetized Target Fusion (MTF) propulsion concept. Analysis of this propulsion system concludes that a 40-fold increase of Isp over chemical propulsion is predicted. This along with a greater than 30% predicted payload mass fraction certainly warrants further development of enabling technologies. The vehicle is designed for multiple interplanetary missions and conceivably may be suited for an automated one-way interstellar voyage.

Miernik, Janie H.↗
NASA NTRS2016

Building Climate Resilience at NASA Ames

NASA Ames Research Center, located at the southern end of the San Francisco Bay (SFB) estuary has identified three primary vulnerabilities to changes in climate. The Ames Climate Adaptation Science Investigator (CASI) workgroup has studied each of these challenges to operations and the potential exposure of infrastructure and employees to an increased frequency of hazards. Sea level rise inundation scenarios for the SFB Area generally refer to projected scenarios in mean sea level rather than changes in extreme tides that could occur during future storm conditions. In the Summer of 2014, high resolution 3-D mapping of the low lying portion of Ames was performed. Those data are integrated with improved sea level inundation scenarios to identify the buildings, basements and drainage systems potentially affected. We will also identify the impacts of sea level and storm surge effects on transportation to and from the Center. This information will help Center Management develop future Master Plans. Climate change will also lead to changes in temperature, storm frequency and intensity. These changes have potential impacts on localized floods and ecosystems, as well as on electricity and water availability. Over the coming decades, these changes are going to be imposed on top of ongoing land use and land cover changes, especially those deriving from continued urbanization and increase in impervious surface areas. These coupled changes have the potential to create a series of cascading impacts on ecosystems, including changes in primary productivity and disturbance of hydrological properties and increased flood risk.The majority of the electricity used at Ames is supplied by hydroelectric dams, which will be influenced by reductions in precipitation or changes in the timing or phase of precipitation which reduces snow pack. Coupled with increased demand for summertime air conditioning and other cooling needs, NASA Ames is at risk for electricity shortfalls. To assess the anticipated energy usage as climate changes, the Ames CASI team is collecting historical energy usage data from Ames facilities, historical weather data, and projected future weather parameters from the CASI Climate subgroup. This data will be incorporated into the RETScreen model to predict how energy usage at Ames will change over the coming century.

climate↗
DOE OSTI2016

Cooling thermal parameters and microstructure features of directionally solidified ternary Sn–Bi–(Cu,Ag) solder alloys

Low temperature soldering technology encompasses Sn–Bi based alloys as reference materials for joints since such alloys may be molten at temperatures less than 180 °C. Despite the relatively high strength of these alloys, segregation problems and low ductility are recognized as potential disadvantages. Thus, for low-temperature applications, Bi–Sn eutectic or near-eutectic compositions with or without additions of alloying elements are considered interesting possibilities. In this context, additions of third elements such as Cu and Ag may be an alternative in order to reach sounder solder joints. The length scale of the phases and their proportions are known to be the most important factors affecting the final wear, mechanical and corrosions properties of ternary Sn–Bi–(Cu,Ag) alloys. In spite of this promising outlook, studies emphasizing interrelations of microstructure features and solidification thermal parameters regarding these multicomponent alloys are rare in the literature. In the present investigation Sn–Bi–(Cu,Ag) alloys were directionally solidified (DS) under transient heat flow conditions. A complete characterization is performed including experimental cooling thermal parameters, segregation (XRF), optical and scanning electron microscopies, X-ray diffraction (XRD) and length scale of the microstructural phases. Experimental growth laws relating dendritic spacings to solidification thermal parameters have been proposed with emphasis on the effects of Ag and Cu. The theoretical predictions of the Rappaz-Boettinger model are shown to be slightly above the experimental scatter of secondary dendritic arm spacings for both ternary Sn–Bi–Cu and Sn–Bi–Ag alloys examined. - Highlights: • Dendritic growth prevailed for the ternary Sn–Bi–Cu and Sn–Bi–Ag solder alloys. • Bi precipitates within Sn-rich dendrites were shown to be unevenly distributed. • Morphology and preferential region for the Ag{sub 3}Sn growth depend on Ag content and Ṫ{sub L}. • Rappaz-Boettinger model reasonably estimated the experimental scatter of λ{sub 2}.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigating the Electronic Structure of Prospective Water-Splitting Oxide BaCe 0.25 Mn 0.75 O 3-δ before and after Thermal Reduction

BaCe 0.25 Mn 0.75 O 3-δ (BCM), a non-stoichiometric oxide with a layered perovskite-like crystal structure, has recently emerged as a prospective contender for application in renewable energy harvesting by solar thermochemical hydrogen generation. Using solar-thermal energy and a reducing environment, oxygen vacancies can be created in high-temperature BCM, and the reduced crystal so obtained can, in turn, produce H 2 by stripping oxygen from H 2 O. Therefore, a first step toward understanding the working mechanism and optimizing the performance of BCM is a thorough and comparative analysis of the electronic structure of the pristine and the reduced material. In this paper, we probe the electronic structure of BCM using the combined effort of first-principles calculations and experimental O K-edge X-ray absorption spectroscopy (XAS). The computed projected density of states (PDOS) and orbital plots are used to propose a simplified model for orbital mixing between the oxygen and metal atoms. With the help of state-of-the-art simulations, we are able to find the origins of the XAS peaks and categorize them on the basis of contribution from Ce and Mn. For the reduced crystal, the calculations show that the change in electron density resulting from the reduction is strongly localized around the oxygen vacancy. Experimental measurements reveal a marked lowering of the first O K-edge peak in the reduced crystal. Using theoretical analysis, this is shown to result from lifting of spin degeneracy in the absorption peaks as well as from a diminished O 2p contribution to the frontier unoccupied orbitals, in accordance with the tight binding scheme. The simulated results serve as a reference for the extent of spectral change as a function of the percentage of oxygen vacancies in the reduced crystal. Here, our study paves the way for the investigation of the working mechanism of BCM and for computational and experimental efforts aimed at design and discovery of efficient water-splitting oxides.

08 HYDROGEN↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m-1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave's traveltime, but also the path of a wave, which can be observed in deviations of the wave's incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope's regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗