Search NASA⌕ Search

SEARCH · Search NASA

Results for “manage simulation input data”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Hydrogeologic modeling and the development of preliminary remediation goals for the Guterl Specialty Steel Site, New York - 15431

The Former Guterl Specialty Steel Corporation Site (Guterl Site) is located 32 kilometers (20 miles) northeast of Buffalo, New York, in Lockport, Niagara County, New York. Between 1948 and 1952, up to 15,875 metric tons (35 million pounds) of natural uranium metal (U) were processed at the Guterl Site. The resulting milling dust, shavings, thermal scale, and associated on-site land disposal contaminated both the facility and soils. These combined manufacturing and disposal processes promoted the development of a large-scale uranium plume in the groundwater underlying the site. Site soils are anthropogenic fill and re-worked, glacially-derived sediments that vary between 0.6 meters (m) and 3 meters (2.0 to 10.0 feet [ft]) in thickness. These soils blanket a highly permeable weathered and fractured zone in the upper Lockport Dolostone bedrock. Aerial recharge through the thin soils flushes soluble constituents to groundwater, as exemplified by variations in specific conductivity and groundwater levels in the bedrock during precipitation events. Seasonally higher groundwater levels also promote contact with uranium contamination in deeper soil, thus complicating the transport characteristics of the site (e.g., oxidizing metallic uranium in the soil that increases leach rates). The Seasonal Soil (SESOIL) compartment model was used to simulate a series of site soil conditions and predict the leaching of uranium to groundwater. The SESOIL input included physical soil characteristics, uranium profiles common in five soil areas, climate data, and an assumed oxidation-state (or chemical valence) for metallic uranium in the soil. The iterative modeling process evaluated site-specific variables to achieve a final deterministic set of input that was applied to the five soil areas. The leaching simulations estimated time-dependent influx rates for uranium, which were input to a numerical groundwater flow (MODFLOW) and transport (MT3DMS) model. The combined modeling effort produced spatially and temporally variable transport conditions that were simulated over a 1,000-year period. The leaching models predicted a preliminary remediation goal (PRG) for uranium of 11 milligrams per kilogram (mg/kg) that is protective of groundwater. However, this low uranium PRG has field implementation concerns, so a modified direct-exposure PRG of 70 mg/kg was developed and evaluated in the flow and transport model. The primary difference in the PRGs is the longevity and location of the uranium plume, which will affect the remedial costs and potential stakeholder acceptance. (authors)

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

Glass Formulation and Testing with RPP-WTP HLW Simulants (Final Report)

This report presents the results from HLW glass formulation work that was conducted at the Vitreous State Laboratory of The Catholic University of America during Part B1 of the RPP-WTP Project (formally known as TWRS-P) in support of the Privatization contractor, BNFL, Inc. Glass formulation development is required to provide data that are essential for design of the vitrification facility in the near term and for safe and reliable operation of that facility in the longer-term. These data also form the bases from which contractually and technically acceptable immobilized high-level waste (IHLW) products will be devised. Glass formulation therefore impacts many aspects of the vitrification facility design, overall process economics, and contractual compliance. As an integral part of the overall system design, development of glass formulations requires input from many sources including tank waste characterization, pretreatment process definition, flow sheet development, and processing schedules. Accordingly, the glass formulation effort must be responsive to new information as the facility design basis evolves.

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

A matrix-inversion method for gamma-source mapping from gamma-count data

In a previous paper it was proposed that a simple matrix inversion method could be used to extract source distributions from gamma-count maps, using simple models to calculate the response matrix. The method was tested using numerically generated count maps. In the present work a 100 kBq Co{sup 60} source has been placed on a gridded surface and the count rate measured using a NaI scintillation detector. The resulting map of gamma counts was used as input to the matrix inversion procedure and the source position recovered. A multi-source array was simulated by superposition of several single-source count maps and the source distribution was again recovered using matrix inversion. The measurements were performed for several detector heights. The effects of uncertainties in source-detector distances on the matrix inversion method are also examined. The results from this work give confidence in the application of the method to practical applications, such as the segregation of highly active objects amongst fuel-element debris. (authors)

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

Modeling Streamflow and Sediment Transport Under Current and Future Climates at the West Valley Site - 20231

Performance Assessment models often require inputs representing specific hydrologic processes such as streamflow and sediment flux. This paper details the development of a process-level hydrologic model and its use in a Probabilistic Performance Assessment (PPA) system-level model for the West Valley Site in western New York. For this PPA model, a watershed model was developed using SWAT (Soil and Water Assessment Tool) and used to evaluate water and sediment fluxes in nearby streams draining the Site. A particular challenge of probabilistic modeling over long periods of time is scaling parameters appropriately so that outcomes are not grossly over- or underestimated. This challenge is especially important when the spatial and temporal scales of process and system level models differ. In this instance, the SWAT watershed model and GoldSim PPA model have similar spatial scales of stream reaches and sub-catchments within a single watershed (approximately hundreds to tens of thousands of square meters, but the temporal scales differ considerably. The SWAT model calculates streamflow and sediment transport rates at the daily scale, while the PPA model requires long-term average values that are applied over hundreds to thousands of year-long periods. As such, distributions of daily climatic parameters such as average precipitation rate, temperature, wind speed, and relative humidity were developed to force a transient SWAT model that captures the dynamic behavior of the watershed in response to daily changes in the weather. Output from the model, which includes stream flow (m{sup 3}/y) and sediment transport rate (Mg/y) is then averaged over a 100-year period as an approximation for the long-term average values which can be applied in GoldSim. A single SWAT model was run 5000 times with varying climatic inputs and deterministic soil, land use, and elevation data, and the range of outputs were compiled into a distribution which can be implemented stochastically into the PPA Model. Additionally, this paper includes a high-level discussion of how changes in future climate affect weather input distributions, and a summary of how those changes in weather affect the outcomes from the SWAT model and the inputs into the PPA model. In this analysis, future climate distributions were developed based on climate research documenting the expected changes to precipitation and temperature in western New York. The 5,000 SWAT simulations were then repeated with these future climate distributions, and the resulting output was implemented into the PPA Model as hydrologic conditions under a future climate. In the current version of the PPA Model, hydrologic parameters change linearly from the 'current climate' conditions to the 'future climate' conditions over a period of 100 years, after which they remain constant. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2009

Group Capability Model

The Group Capability Model (GCM) is a software tool that allows an organization, from first line management to senior executive, to monitor and track the health (capability) of various groups in performing their contractual obligations. GCM calculates a Group Capability Index (GCI) by comparing actual head counts, certifications, and/or skills within a group. The model can also be used to simulate the effects of employee usage, training, and attrition on the GCI. A universal tool and common method was required due to the high risk of losing skills necessary to complete the Space Shuttle Program and meet the needs of the Constellation Program. During this transition from one space vehicle to another, the uncertainty among the critical skilled workforce is high and attrition has the potential to be unmanageable. GCM allows managers to establish requirements for their group in the form of head counts, certification requirements, or skills requirements. GCM then calculates a Group Capability Index (GCI), where a score of 1 indicates that the group is at the appropriate level; anything less than 1 indicates a potential for improvement. This shows the health of a group, both currently and over time. GCM accepts as input head count, certification needs, critical needs, competency needs, and competency critical needs. In addition, team members are categorized by years of experience, percentage of contribution, ex-members and their skills, availability, function, and in-work requirements. Outputs are several reports, including actual vs. required head count, actual vs. required certificates, CGI change over time (by month), and more. The program stores historical data for summary and historical reporting, which is done via an Excel spreadsheet that is color-coded to show health statistics at a glance. GCM has provided the Shuttle Ground Processing team with a quantifiable, repeatable approach to assessing and managing the skills in their organization. They now have a common frame of reference across NASA/contractor lines to communicate and mitigate any critical skills concerns.

Olejarski, Michael↗
DOE OSTI2019

Solid Secondary Waste Immobilization in Cementitious Waste Forms at the Hanford Site - 19081

During operations of the Hanford Tank Waste Treatment and Immobilization Plant, solid secondary wastes (SSW) will be generated. These materials include HEPA filters, granular activated carbon (GAC), and silver mordenite (AgM). The planned disposal pathway for SSW is through solidification (blending) or encapsulation in a cementitious waste form and then disposal on site at the Hanford Site's Integrated Disposal Facility. Performance assessment evaluations of SSW forms have been performed; however, limited experimental data on simulated SSW is available as input. This paper will cover efforts to fill this gap in experimental data through testing of simulated SSW materials. Two grout formulations were used in the tests: 1) Hanford grout mix 5, a mix of 25 wt% ordinary Portland cement (OPC) and 75 wt% fly ash (FA); and 2) the Cast Stone dry blend of 47 wt% blast furnace slag, 45 wt% FA, and 8 wt% OPC. The sorption (Kd) of potential contaminants (iodide, iodate, technetium) to the grouts and SSW materials was investigated in simulated grout pore water. Waste forms containing GAC and AgM were prepared using the two grout formulations and characterized using scanning electron microscopy energy dispersive X-ray spectroscopy and X-ray diffraction. Acquiring this information will further support the candidacy of cementitious waste forms to be used for disposal of SSW at Hanford. (authors)

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

A deep learning-based workflow for fast prediction of 3D state variables in geological carbon storage: A dimension reduction approach

Deep learning (DL) models are extensively used as surrogate models for high-fidelity simulations of multiphase fluid flow in porous media at large scales, enabling fast forecasts of the spatial–temporal evolution of three-dimensional (3D) state variables in geological carbon storage (GCS). However, training these models in high-dimensional space remains computationally demanding and prone to overfitting because of limited training data. This paper presents a novel workflow to address these challenges by integrating dimension reduction (DR) methods. Here, the proposed workflow employed pre-trained DR models to extract the latent variables of geological models and state variables and utilized the multi-layer perceptron (MLP) for constructing mapping functions between the input and output variables in latent spaces. Subsequently, the pre-trained reconstruction models converted the MLP-predicted latent state variables to their original high-dimensional form. Furthermore, we proposed a novel strategy for the DR and reconstruction of 3D saturation fields to account for the unique data characteristics of sparsity, nonuniformity, and discontinuity. The proposed strategy applied PCA and inverse PCA for 2D average saturation fields and developed a DL-based 3D reconstruction model, leveraging three 2D average saturation fields as input to produce a 3D saturation field as output. The pre-training of DR and reconstruction models and training of MLP models were conducted on 84 Gulf of Mexico (GoM) simulations and evaluated on 12 testing simulations. Each simulation contained 720 monthly time steps, with the first 360 months as the injection period and the rest as the post-injection period. The proposed workflow, incorporating DR and DL models, accurately predicts the normalized 3D pressure fields, achieving mean square error (MSE) of 2.92 × 10 -7 compared to the ground truth obtained from a full-physics simulator. Furthermore, the proposed strategy outperformed PCA and convolutional autoencoder (CAE) models on 3D saturation fields, resulting in minor workflow prediction errors with an MSE of 2.93 × 10 -5 . The results suggest the proposed workflow provides sufficient predictive fidelity across temporal and spatial scales, and enables a speedup of 160 times compared to the full-physics simulator, facilitating improved decision-making and risk assessment for large-scale GCS management in real-time scenarios.

3D reconstruction model↗
DOE OSTI2022

Supplementary Data for "Evaluation of the Economic Implications of Varied Pressure Drawdown Strategies Generated Using a Real-time, Rapid Predictive, Multi-fidelity Model for Unconventional Oil and Gas Wells" by Bello, K., Vikara, D., Sheriff, A., Viswanathan, H., Carr, T., Sweeney, M., O'Malley, D., Marquis, M., Vactor, R.T., and Cunha, L.

The Bello et al. study evaluates the impact of contrasting pressure drawdown on gas productivity and the resulting economics of a well in the Marcellus Shale of the Appalachian Basin. This research applies a techno-economic analysis approach to help identify potential ways pressure management strategies can be used to improve cumulative recovery of hydraulically fractured horizontal wells while maintaining project profitability. Gas production forecast outlook scenarios of the Marcellus Shale Energy and Environment Laboratory Laboratory's MIP-3H well were generated under varying pressure drawdowns using two approaches: 1) a novel physics-informed machine learning (PIML) workflow and 2) via traditional reservoir simulation in Computer Modeling Group’s (CMG) GEM Compositional & Unconventional Simulator. Cash flow and other economic metrics of interest were compiled on the production outlook using the U.S. Department of Energy's (DOE) National Energy Technology Laboratory (NETL) Unconventional Shale Well Economic Model (UShWEM).The sheets within this Microsoft ExcelTM workbook provide the economic metric outputs for the baseline condition and the one-at-a-time (OAT) sensitivity analysis of UShWEM's input parameters for each of the production scenarios evaluated.

Fracture Network Model↗
DOE OSTI2016

Secondary Waste Form Development and Data Package for Hanford's Integrated Disposal Facility Performance Assessment - 16479

The Hanford Low-Activity Waste (LAW) Vitrification facility is planned to immobilize Hanford tank waste for on-site, near-surface disposal at Hanford's Integrated Disposal Facility (IDF). During this process, liquid and solid secondary wastes will be generated which are also slated for disposal at the IDF. In order to receive authorization to dispose of the waste in the IDF, a Performance Assessment (PA) must be completed in accordance with DOE Order 435.1, Radioactive Waste Management, and its accompanying Manual. The PA must include calculations of potential doses to a member of the public for a 1,000-year period to provide a reasonable expectation that performance objectives will not be exceeded and that the waste will be managed in a manner that is protective of worker and public health and safety, and the environment. Liquid and solid secondary wastes are planned to be immobilized within cementitious waste forms similar to those used elsewhere in the United States and around the world for permanent disposal of low-level radioactive waste (LLW). Data packages and supporting documents are being developed to provide the technical underpinning for the analyses of contaminant releases to the environment from these waste forms for the IDF PA. Traditionally, releases from cementitious waste forms have utilized a diffusion-limited release model using contaminant-specific diffusivities measured in standard leach tests as inputs. More mechanistic modeling techniques are under development to increase the rigor of the PA analyses. In addition, recent testing has shown diffusivities of key contaminants can be reduced by several orders-of-magnitude through optimized formulation of the waste form and use of additives to bind specific contaminants. (authors)

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

Thermo-Hydraulic-Mechanical Modelling of the Full-Scale Emplacement (FE) Experiment at the Mont Terri Rock Laboratory - 16490

The Full-Scale Emplacement Experiment at the Mont Terri underground research laboratory is a heater test carried out in a gallery constructed in Opalinus Clay. The aim of the experiment is both to develop emplacement procedures for radioactive waste canisters and their surrounding backfill as well as understanding the thermo-hydro-mechanical effects of disposing of radioactive waste on the engineering materials and the host rock. Coupled thermo-hydro-mechanical models have been used to aid planning of the locations of sensors, to inform decisions on the schedule for turning on the heaters as well as to make predictions of the temperatures and water pressures in the bentonite and the Opalinus Clay. The FE Experiment was preceded by a number of other experiments at the Mont Terri laboratory at smaller scales, but with similar materials, which provided crucial input to the models and gave confidence in the general process understanding. Blind predictions were made prior to turning on the heaters, of the temperature and relative humidity in the bentonite surrounding the heaters. The comparison between the blind predictions and early data from the experiment show that knowledge gained from smaller scale experiments on similar materials has enabled good blind predictions to be made and gives confidence that temperature in the experiment is unlikely to rise to a level that could damage the sensors. It also demonstrates that the thermal properties of the disposal system are understood well enough to produce reliable estimates of temperature in the backfill. (authors)

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

Liquid secondary waste: Waste form formulation and qualification

The Hanford Site Effluent Treatment Facility (ETF) currently treats aqueous waste streams generated during site cleanup activities. When the Hanford Tank Waste Treatment and Immobilization Plant (WTP) begins operations, including Direct Feed Low Activity Waste (DFLAW) vitrification, a liquid secondary waste (LSW) stream from the WTP will need to be treated. The volume of effluent for treatment at the ETF will increase significantly. The powdered salt waste form produced by the ETF will be replaced by a stabilized solidified waste form for disposal in Hanford’s Integrated Disposal Facility (IDF). Washington River Protection Solutions is implementing a Secondary Liquid Waste Immobilization Technology Development Plan to address the technology needs for a waste form and solidification process to treat the increased volume of waste planned for disposal at the IDF. Waste form testing to support this plan is composed of work in the near term to provide data as input to a performance assessment (PA) for Hanford’s IDF. In 2015, three Hanford Liquid Secondary Waste simulants were developed based on existing and projected waste streams. Using these waste simulants, fourteen mixes of Hanford Liquid Secondary Waste were prepared and tested varying the waste simulant, the water-to-dry materials ratio, and the dry materials blend composition.1 In FY16, testing was performed using a simulant of the EMF process condensate blended with the caustic scrubber—from the Low Activity Waste (LAW) melter—, processed through the ETF. The initial EMF-16 simulant will be based on modeling efforts performed to determine the mass balance of the ETF for the DFLAW.2 The compressive strength of all of the mixes exceeded the target of 3.4 MPa (500 psi) to meet the requirements identified as potential IDF Waste Acceptance Criteria in Table 1 of the Secondary Liquid Waste Immobilization Technology Development Plan.3 The hydraulic properties of the waste forms tested (hydraulic conductivity and water characteristic curves) were comparable to the properties measured on the Savannah River Site (SRS) Saltstone waste form. Future testing should include efforts to first; 1) determine the rate and amount of ammonia released during each unit operation of the treatment process to determine if additional ammonia management is required, then; 2) reduce the ammonia content of the ETF concentrated brine prior to solidification, making the waste more amenable to grouting, or 3) manage the release of ammonia during production and ongoing release during storage of the waste form, or 4) develop a lower pH process/waste form thereby precluding ammonia release.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2009

Watershed and Hydrodynamic Modeling for Evaluating the Impact of Land Use Change on Submerged Aquatic Vegetation and Seagrasses in Mobile Bay

There is a continued need to understand how human activities along the northern Gulf of Mexico coast are impacting the natural ecosystems. The gulf coast is experiencing rapid population growth and associated land cover/land use change. Mobile Bay, AL is a designated pilot region of the Gulf of Mexico Alliance (GOMA) and is the focus area of many current NASA and NOAA studies, for example. This is a critical region, both ecologically and economically to the entire United States because it has the fourth largest freshwater inflow in the continental USA, is a vital nursery habitat for commercially and recreational important fisheries, and houses a working waterfront and port that is expanding. Watershed and hydrodynamic modeling has been performed for Mobile Bay to evaluate the impact of land use change in Mobile and Baldwin counties on the aquatic ecosystem. Watershed modeling using the Loading Simulation Package in C++ (LSPC) was performed for all watersheds contiguous to Mobile Bay for land use Scenarios in 1948, 1992, 2001, and 2030. The Prescott Spatial Growth Model was used to project the 2030 land use scenario based on observed trends. All land use scenarios were developed to a common land classification system developed by merging the 1992 and 2001 National Land Cover Data (NLCD). The LSPC model output provides changes in flow, temperature, sediments and general water quality for 22 discharge points into the Bay. These results were inputted in the Environmental Fluid Dynamics Computer Code (EFDC) hydrodynamic model to generate data on changes in temperature, salinity, and sediment concentrations on a grid with four vertical profiles throughout the Bay s aquatic ecosystems. The models were calibrated using in-situ data collected at sampling stations in and around Mobile bay. This phase of the project has focused on sediment modeling because of its significant influence on light attenuation which is a critical factor in the health of submerged aquatic vegetation. The impact of land use change on sediment concentrations was evaluated by analyzing the LSPC and EFDC sediment simulations for the four land use scenarios. Such analysis was also performed for storm and non-storm periods. In- situ data of total suspended sediments (TSS) and light attenuation were used to develop a regression model to estimate light attenuation from TSS. This regression model was used to derive marine light attenuation estimates throughout Mobile bay using the EFDC TSS outputs. The changes in sediment concentrations and associated impact on light attenuation in the aquatic ecosystem were used to perform an ecological analysis to evaluate the impact on seagreasses and Submerged Aquatic Vegetation (SAV) habitat. This is the key product benefiting the Mobile Bay coastal environmental managers that integrates the influences of sediments due to land use driven flow changes with the restoration potential of SAVs.

Estes, Maurice G.↗
DOE OSTI2025

Evaluation of Hypotheses for the Fate of Historical Discharges from Hanford Site BC Cribs and Trenches

A preliminary review of declassified reports and historical documents from the 1950s and 1960s indicated they may contain valuable information regarding potential impacts to groundwater from liquid waste discharged to the 216-BC Cribs and Trenches during the late 1950s. Transport of contaminants in the vadose zone beneath the BC Cribs area is simulated by a vadose zone model developed for the Composite Analysis (DOE/RL-2019-52, Composite Analysis for Low-Level Waste Disposal in the Hanford Site Central Plateau ). Results from this model indicate there has been no substantial breakthrough of contaminants (particularly technetium-99) to groundwater from the vadose zone up to the present day. However, there is uncertainty in the anisotropy and other model input parameters as well as acknowledged sparseness of sampling data for this area. In this study, historical documents and other information were reviewed to evaluate the possibility that contaminant breakthrough to groundwater occurred during or shortly after operation of the BC Cribs sites, and whether refinements of the vadose zone model parametrization should be considered.

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

DOE Hybrid and Electric Vehicle Test Platform

Based on the contract NT-42790 to the Department of Energy, “Plug-in Hybrid Ethanol Research Platform”, Advanced Vehicle Research Center (AVRC) Virginia has successfully developed the phase I electric drive train research platform which has been named as Laboratory Rapid Application Testbed (LabRAT). In phase II, LabRAT is to be upgraded into plug-in hybrid research platform, which will be capable of testing power systems for electric vehicles, and plug-in hybrid electric vehicles running on conventional as well as alternative fuels. LabRAT is configured as a rolling testbed with plentiful space for installing various component configurations. Component connections are modularized for flexibility and are easily replaced for testing various mechanisms. LabRAT is designed and built as a full functional vehicle chassis with a steering system, brake system and four wheel suspension. The rear drive axle offers maximum flexibility with a quickly changeable gear ratio final drive to accommodate different motor speed requirements. The electric drive system includes an electric motor which is mechanically connected to the rear axle through an integrated speed/torque sensor. Initially, a 100 kW UQM motor and corresponding UQM motor controller is used which can be easily replaced with another motor/controller combination. A lithium iron phosphate (LiFePO4) battery pack is installed, which consists of 108 cells of 100 AH capacity, giving the total energy capacity of 32.5 kWh. Correspondingly, a fully functional battery management system (BMS) is installed to perform battery cell operation monitoring, cell voltage balancing, and reporting battery real time operating parameters to vehicle controller. An advanced vehicle controller ECU is installed for controlling the drive train. The vehicle controller ECU receives traction or braking torque command from driver through accelerator and brake pedal position sensors and battery operating signals from the BMS through CAN BUS, and then generates motor torque command (traction or braking) to the motor controller based on the control algorithm software embedded in the vehicle controller ECU. The vehicle controller ECU is a re-programmable electronic control unit. Any control algorithm software developed can be easily downloaded to vehicle controller ECU to test any newly developed control strategy. The flexibility of the control system significantly enhances the practical applicability of the LabRAT. A new test methodology has been developed for the LabRAT simulating any vehicles running on road with different weights from compact passenger car to light duty truck on an AC or eddy current dynamometers without much effort for modification of the system. LabRAT is equipped with a fully functional data acquisition system supplied by CyberMetrix. The measurement points along the drive train are DC electric power between battery pack and motor controller input, AC electric power between motor controller and electric motor, mechanical power between motor and rear axle. The data acquisition system is designed with more capability than current requirements in order to meet the requirements for phase II.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS2012

Availability Estimation for Facilities in Extreme Geographical Locations

A value added analysis for the Reliability. Availability and Maintainability of McMurdo Ground Station was developed, which will be a useful tool for system managers in sparing, maintenance planning and determining vital performance metrics needed for readiness assessment of the upgrades to the McMurdo System. Output of this study can also be used as inputs and recommendations for the application of Reliability Centered Maintenance (RCM) for the system. ReliaSoft's BlockSim. a commercial Reliability Analysis software package, has been used to model the availability of the system upgrade to the National Aeronautics and Space Administration (NASA) Near Earth Network (NEN) Ground Station at McMurdo Station in the Antarctica. The logistics challenges due to the closure of access to McMurdo Station during the Antarctic winter was modeled using a weighted composite of four Weibull distributions. one of the possible choices for statistical distributions throughout the software program and usually used to account for failure rates of components supplied by different manufacturers. The inaccessibility of the antenna site on a hill outside McMurdo Station throughout one year due to severe weather was modeled with a Weibull distribution for the repair crew availability. The Weibull distribution is based on an analysis of the available weather data for the antenna site for 2007 in combination with the rules for travel restrictions due to severe weather imposed by the administrating agency, the National Science Foundation (NSF). The simulations resulted in an upper bound for the system availability and allowed for identification of components that would improve availability based on a higher on-site spare count than initially planned.

Fischer, Gerd M.↗
DOE OSTI2015

T2-weighted four dimensional magnetic resonance imaging with result-driven phase sorting

Purpose: T2-weighted MRI provides excellent tumor-to-tissue contrast for target volume delineation in radiation therapy treatment planning. This study aims at developing a novel T2-weighted retrospective four dimensional magnetic resonance imaging (4D-MRI) phase sorting technique for imaging organ/tumor respiratory motion. Methods: A 2D fast T2-weighted half-Fourier acquisition single-shot turbo spin-echo MR sequence was used for image acquisition of 4D-MRI, with a frame rate of 2–3 frames/s. Respiratory motion was measured using an external breathing monitoring device. A phase sorting method was developed to sort the images by their corresponding respiratory phases. Besides, a result-driven strategy was applied to effectively utilize redundant images in the case when multiple images were allocated to a bin. This strategy, selecting the image with minimal amplitude error, will generate the most representative 4D-MRI. Since we are using a different image acquisition mode for 4D imaging (the sequential image acquisition scheme) with the conventionally used cine or helical image acquisition scheme, the 4D dataset sufficient condition was not obviously and directly predictable. An important challenge of the proposed technique was to determine the number of repeated scans (N{sub R}) required to obtain sufficient phase information at each slice position. To tackle this challenge, the authors first conducted computer simulations using real-time position management respiratory signals of the 29 cancer patients under an IRB-approved retrospective study to derive the relationships between N{sub R} and the following factors: number of slices (N{sub S}), number of 4D-MRI respiratory bins (N{sub B}), and starting phase at image acquisition (P{sub 0}). To validate the authors’ technique, 4D-MRI acquisition and reconstruction were simulated on a 4D digital extended cardiac-torso (XCAT) human phantom using simulation derived parameters. Twelve healthy volunteers were involved in an IRB-approved study to investigate the feasibility of this technique. Results: 4D data acquisition completeness (C{sub p}) increases as NR increases in an inverse-exponential fashion (C{sub p} = 100 − 99 × exp(−0.18 × N{sub R}), when N{sub B} = 6, fitted using 29 patients’ data). The N{sub R} required for 4D-MRI reconstruction (defined as achieving 95% completeness, C{sub p} = 95%, N{sub R} = N{sub R,95}) is proportional to N{sub B} (N{sub R,95} ∼ 2.86 × N{sub B}, r = 1.0), but independent of N{sub S} and P{sub 0}. Simulated XCAT 4D-MRI showed a clear pattern of respiratory motion. Tumor motion trajectories measured on 4D-MRI were comparable to the average input signal, with a mean relative amplitude error of 2.7% ± 2.9%. Reconstructed 4D-MRI for healthy volunteers illustrated clear respiratory motion on three orthogonal planes, with minimal image artifacts. The artifacts were presumably caused by breathing irregularity and incompleteness of data acquisition (95% acquired only). The mean relative amplitude error between critical structure trajectory and average breathing curve for 12 healthy volunteers is 2.5 ± 0.3 mm in superior–inferior direction. Conclusions: A novel T2-weighted retrospective phase sorting 4D-MRI technique has been developed and successfully applied on digital phantom and healthy volunteers.

60 APPLIED LIFE SCIENCES↗