Search NASA⌕ Search

SEARCH · Search NASA

Results for “reference cell model”

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.

405 records · Page 23
All sourcesNASA NTRSDOE OSTI
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 OSTI2012

A Preliminary Analysis of Dose Rates Associated with ITER CVCS Equipment/Area Location

A preliminary analysis of the ITER Chemical and Volume Control System (CVCS) Area was performed to assess dose rates outside the walls and ceiling of the facility after 1.5 years of operation at shutdown, 2 days, and 10 days after shutdown. For this purpose a simplified Monte Carlo computer model was developed using the MCNP (MCNP5 Ver. 1.51) code. Two components are included: the smaller filter tank and the larger ion exchanger. These pieces of equipment are associated with the Integrated Blanket ELM Divertor Primary Heat Transfer System, which will have the largest dose rates associated with activated corrosion products during operation in comparison with other systems. The ion exchanger contained two source regions, a 1.2-m-thick resin bed above a 0.55 m-thick skirt, and a 0.8-m-thick water region. The filter constituted an additional source. Thus the model consisted of three sources (filter, resin, water), homogeneously distributed within the appropriate source regions. However, much of the results (that address individual isotopes) are presented with the two sources in the ion exchanger combined. In these cases the sources are referred to as the 'ion exchanger source' and the 'filter source.' Dimensions for the facility and components, as well as source isotopes and strengths, and material densities, were supplied by US ITER. Because of its simplification, the model does not contain pipes. Consequently, radiation streaming through pipe penetrations, radiation emanating from the pipes, and shielding from the pipes were not considered in this analysis. Dose rates on the outside of two walls and the ceiling were calculated. The two walls are labeled as the 'long' wall (aligned with the X-axis) and the 'short' wall (aligned with the Y-axis). These walls and ceiling were nominally set to 30-cm-thick concrete. In the original analysis, standard concrete (2.3 g/cc density) was used. In addition to the shielding walls/ceiling, a floor and an additional wall opposite the long wall were added for photon scattering contributions. These were both 10-cm-thick, standard concrete structures. Other components (tanks, pipes, etc.), that were not included in the model, would potentially add additional scattering and shielding. Possibly these additional effects will be addressed in a later, more detailed analysis. The room was 29.6 m in length (X-axis), but was limited to 15 m in the model. The inside width (Y axis) and height (Z axis) were 4 m and 3.4 m, respectively. The origin for the model was located inside the room at the corner opposite the long wall and adjacent to the short wall at the floor level. The room was filled with air at standard temperature and pressure. The stainless steel (SS304) wall thicknesses for the ion exchanger and filter were 2.2 cm and 0.8 cm, respectively. The axial center of the filter was located 140 cm from the short wall and 100 cm from the long wall (outer surface). The axial center of the ion exchanger was located 440 cm from the short wall and 250 cm from the long wall (inner surface). The resin was assumed to be a homogeneous mixture of equal atom density fractions of hydrogen and carbon* at a specified density of 1.136 g/cc. The filter material was assumed to be homogeneous carbon at a specified density of 1.8 g/cc. If the filter media were stainless steel and the accumulated activity were the same, the dose rate outside the filter would be lower, provided the density of the stainless steel sintered filter material is significantly higher than 1.8g/cc. The densities of the water and air were assumed to be 1.0 g/cc and 1.096E-3 g/cc, respectively. The model included 10 small volumes placed outside the model adjacent to the wall at locations where the dose rates were expected to be highest and which showed to a degree the fall-off of the dose rate with distance along a particular wall or ceiling. These 'tally cells' are shown in Fig. 2 and are described in Table 1. Each cell had dimensions of 50 cm x 50 cm x 1 cm and was oriented so that the 1-cm-thick dimension was perpendicular to the wall against which the cell was located. Also, smaller tally volumes were placed at 5 cm (updated model), 30 cm, 50 cm, and 100 cm from each source, both radially from the outside of each component at the height of each source mid-plane, and vertically, along the axial centerline, above the ion-exchanger and below the filter. (For the ion exchanger, the vertical distance was limited to 75 cm because of interference from the ceiling.) These tally regions can be seen in Fig. 3. Each cell had dimensions of 5 cm x 5 cm x 1 cm, oriented so that the 1-cm-thick dimension was perpendicular to the X, Y or Z axis.

08 HYDROGEN↗
DOE OSTI2012

Effect of Curing Temperature and Composition on the Performance Properties of Saltstone - 12473

When freshly mixed Saltstone is placed into a vault at the Savannah River Site, the heat of hydration produced during curing causes the temperature of the grout to increase. It is not uncommon for the temperature of the grout to increase by 30 to 50 deg. C above the starting temperature depending on the vault size and location, environmental conditions and pour schedule. For example, temperatures measured by thermocouples in one cell have exceeded 70 deg. C. The elevated temperatures, once reached, can extend for weeks or even months. Recent studies at Savannah River National Laboratory (SRNL) have shown that curing simulated Saltstone mixes at elevated temperatures can lead to reduced compressive strength and dynamic Young's modulus as well as increased total porosity when compared to the same grout mixes cured at 20 deg. C. Samples cured at 60 deg. C can have hydraulic conductivities approximately 600 times greater than samples cured at 20 deg. C. Previous studies have shown that the 60 deg. C cured mix had a dynamic Young's modulus that was half that of the mix cured at 20 deg. C. However, high temperature curing may not always have negative effects on the performance properties of the grout. These results indicate that higher curing temperatures can significantly alter the performance properties that are inputs for the Saltstone Performance Assessment. The Saltstone formulation (mix design) must produce a grout waste form that meets both placement and performance properties. A statistically designed set of mixes was developed to determine key process and compositional factors that affect the performance properties of Saltstone grout. A total of 27 mixes were batched and tested. There are eight baseline mixes that contain high and low values of aluminate concentration, varying w/p ratio, and varying fly ash content in the premix. Each of the eight mixes was cured at 22, 40 and 60 deg. C. The three additional reference mixes batched at the beginning, middle and end of testing showed good reproducibility throughout the study. The results of varying the operational and compositional factors on the performance properties of Saltstone grout are: - The Young's modulus is inversely affected by the curing temperature of the grout. - The aluminate content of the salt solution increases the Young's modulus of the Saltstone for samples cured at 20 and 40 deg. C. - All grouts cured at 60 deg. C have low Young's moduli and the composition of the grout mix did not have any positive effects on the performance properties after curing at this temperature. - The fly ash content in the premix had not discernible effect on the Young's modulus. - The water to premix ratio does not affect the Young's modulus. - Grouts made with the high aluminate salt solution have lower hydraulic conductivities than those made with the low aluminate simulant. - Temperature has an inverse effect on the hydraulic conductivity. - In general, increasing the water to premix ratio and the fly ash content will increase the hydraulic conductivity. - The porosity and cured density of the Saltstone grout are primarily a factor of the water to premix ratio of the grout formulation. The correlation between hydraulic conductivity and Young's modulus investigated in previous studies was evaluated further in this report. The current data as well as data from previous studies was fitted by a second degree polynomial function with a 77.1 % R-squared value. This model can be used to estimate the hydraulic conductivity of a sample based its Young's modulus for a limited range (E values up to 7 GPa). More data points should be included to further define the relationship between these performance properties are high Young's moduli and low hydraulic conductivities. The Young's modulus should still be measured on grouts formulated in future studies to be used as screening method for samples to be sent for hydraulic conductivity measurements. Saltstone grout has a complicated microstructure that is affected by the composition of the salt solution, the dry feed mixture composition, and the conditions under which it cures. The performance properties investigated are interrelated, which can make it difficult to determine individual relationships between various factors. Porosity is a large factor in determining both the Young's modulus and the hydraulic conductivity of the samples. Although the cured density is a straightforward measurement, it is not a significant performance property that needs further investigation. The porosity measurement provides similar and more meaningful data in relation to other performance properties. Based on the results of this study, the aluminate content of the salt solution does have an effect on the Young's modulus and the hydraulic conductivity and its effects should continue to be investigated. However, a new projection of aluminate contents should be used since projections change over time. The variability of fly ash in the premix did not prove to be a contributing factor that effects the performance properties of Saltstone and does not need to be included in future studies unless significant deviations from the baseline mix are proposed. The water to premix ratio has a significant effect on the performance properties of interest and the effects of varying this factor need further investigation. It is recommended that smaller ranges of w/p are investigated in order to provide data that is more relevant to the Saltstone being placed in the vaults. (authors)

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

High Temperature Steam Electrolysis Process Performance and Cost Estimates

Technology readiness levels (TRLs) of electrolysis systems have dramatically increased in recent years as the interest in clean hydrogen production and decarbonization of transportation, industrial and other sectors increases across the globe. This is especially true of high temperature steam electrolysis (HTSE) / solid oxide electrolysis cell (SOEC) systems which show promise of much higher system efficiencies than other more developed electrolysis technologies. This possibility of higher efficiencies of HTSE / SOEC systems has been previously assumed to be theoretically possible but in recent years it has become less theoretical and more realistic as an increasing amount of suppliers complete lab and pilot tests showing very promising results. Research in the areas of manufacturing techniques, material selection, electrode and electrolyte compositions, and balance of plant size and integration continues at a fast pace as an increasing number of suppliers both internationally and domestically become involved. The advantages of HTSE become more pronounced when HTSE is coupled with nuclear power plants (NPPs). This is because thermal energy produced by the nuclear reactor can be used in a series of heat transfer loops and heat exchangers to vaporize HTSE feedwater, which drastically improves the economics of the process. Idaho National Laboratory (INL) has been very involved in the research and modeling of HTSE systems for a number of years, in collaboration with other national laboratories, academia, and industry stakeholders both on the hydrogen production as well as the hydrogen demand side. The modeling completed over the years on a large variety of projects has led to a wealth of knowledge at INL including in the area of the technoeconomic assessment (TEA) of HTSE systems. TEAs include process modeling of the HTSE systems to calculate system energy requirements and equipment sizing, followed by estimation of capital and operating costs to enable calculation of the levelized cost of hydrogen (LCOH). The TEA work performed has produced incremental improvements and tuning of the methods, assumptions, models, and results of the analyses as well as providing some opportunities for validating these results. The purpose of this document is to record the current baseline HTSE analyses led by INL to show the current status of assumptions and costs of these systems. Given the rapid development of this technology, the variety of suppliers entering the space, and the increasing attention government and industry are giving to such systems, this document may be updated on a periodic basis with updated analysis and assumptions. This document compiles various analyses results and approaches completed over a period of years into a single document to be used as a baseline going forward. It represents what the INL HTSE analysis group assumes to be the internal best estimate of the current operation, costs, and landscape of the HTSE industry state of the art capability for current SOEC technology in an Nth-of-a-Kind (NOAK) plant, which in this study is defined as existence of the manufacturing capacity to support previous deployment of N = 100 count of 25 MWe modular HTSE blocks (with modular equipment component cost reductions specified as following a 95% learning curve). That said, This is a public document and as such so no proprietary data was used or included in this report. There may be HTSE suppliers that have performance specifications, and cost estimates, and test data that differ from the analysis presented in this document. This document is meant to be a best conservative estimate of the technology and not an absolute reference.

08 HYDROGEN↗
DOE OSTI2022

High Temperature Steam Electrolysis Process Performance and Cost Estimates - DOE Hydrogen Program AMR Presentation

Technology readiness levels (TRLs) of electrolysis systems have dramatically increased in recent years as the interest in clean hydrogen production and decarbonization of transportation, industrial and other sectors increases across the globe. This is especially true of high temperature steam electrolysis (HTSE) / solid oxide electrolysis cell (SOEC) systems which show promise of much higher system efficiencies than other more developed electrolysis technologies. This possibility of higher efficiencies of HTSE / SOEC systems has been previously assumed to be theoretically possible but in recent years it has become less theoretical and more realistic as an increasing amount of suppliers complete lab and pilot tests showing very promising results. Research in the areas of manufacturing techniques, material selection, electrode and electrolyte compositions, and balance of plant size and integration continues at a fast pace as an increasing number of suppliers both internationally and domestically become involved. The advantages of HTSE become more pronounced when HTSE is coupled with nuclear power plants (NPPs). This is because thermal energy produced by the nuclear reactor can be used in a series of heat transfer loops and heat exchangers to vaporize HTSE feedwater, which drastically improves the economics of the process. Idaho National Laboratory (INL) has been very involved in the research and modeling of HTSE systems for a number of years, in collaboration with other national laboratories, academia, and industry stakeholders both on the hydrogen production as well as the hydrogen demand side. The modeling completed over the years on a large variety of projects has led to a wealth of knowledge at INL including in the area of the technoeconomic assessment (TEA) of HTSE systems. TEAs include process modeling of the HTSE systems to calculate system energy requirements and equipment sizing, followed by estimation of capital and operating costs to enable calculation of the levelized cost of hydrogen (LCOH). The TEA work performed has produced incremental improvements and tuning of the methods, assumptions, models, and results of the analyses as well as providing some opportunities for validating these results. The purpose of this document is to record the current baseline HTSE analyses led by INL to show the current status of assumptions and costs of these systems. Given the rapid development of this technology, the variety of suppliers entering the space, and the increasing attention government and industry are giving to such systems, this document may be updated on a periodic basis with updated analysis and assumptions. This document compiles various analyses results and approaches completed over a period of years into a single document to be used as a baseline going forward. It represents what the INL HTSE analysis group assumes to be the internal best estimate of the current operation, costs, and landscape of the HTSE industry state of the art capability for current SOEC technology in an Nth-of-a-Kind (NOAK) plant, which in this study is defined as existence of the manufacturing capacity to support previous deployment of N = 100 count of 25 MWe modular HTSE blocks (with modular equipment component cost reductions specified as following a 95% learning curve). That said, this is a public document and as such so no proprietary data was used or included in this report. There may be HTSE suppliers that have performance specifications, and cost estimates, and test data that differ from the analysis presented in this document. This document is meant to be a best conservative estimate of the technology and not an absolute reference.

08 HYDROGEN↗
DOE OSTI2015

A Dual Stage Sodium Thermo-Electro-Chemical Converter

The need for reliable and distributable energy sources is evident given the current forecasting of the world’s energy demand [1]. One prominent obstacle is the lack of an efficient direct heat-to-electricity conversion technology. The sodium thermo-electro-chemical converter (Na-TECC) is a highly efficient thermally regenerative electrochemical system (TRES) that can fill this void. The advantage of the Na-TECC over other direct conversion technologies is that it is a closed device (unlike fuel cells) requiring nothing more than a heat input (similar to thermoelectrics). It also has a high power density (> 1 W/cm2) and specific power up to 0.5 kW/kg , with the potential to achieve > 45% heat to electricity conversion efficiencies [2]. Furthermore, the converter rejects heat at a temperature hot enough to make it amenable for cogeneration, potentially under air-cooled operation. This electrochemical technology uses a special ion-selective solid-electrolyte called beta-alumina (Na2O•11Al2O3), that is highly conductive to sodium cations. A more conducting variation of this compound is beta”-alumina (Na2O•MgO•5Al2O3), where the magnesium acts as a stabilizer [3]. In the device, liquid sodium enters an evaporator where it is heated isobarically and vaporizes. The vapor is then oxidized at the interface between the anode and the solid-electrolyte. The Na-TECC generates electric power by allowing the high pressure oxidized sodium vapor ions to expand isothermally through the solid-electrolyte from the anode to the cathode. The liberated electrons are simultaneously conducted through an external load and recombine with the cations at the cathode. The low pressure sodium vapor is subsequently condensed and a wick pumps the liquid into the high pressure plenum to complete the cycle. Previous Na-TECCs have achieved efficiencies of 18% due to practical losses, whereas the potential efficiencies are above 45% [4]. Thermodynamic limitations and design constraints have limited previous devices from attaining these record efficiencies. Specifically, the maximum potential voltage in these previous devices was substantially reduced by the kinetic limitations of the charge transfer in the cathode (referred to as the overpotential) and by Joule heating losses due to the solid-electrolyte resistance, electrode resistance, and current collector resistance. Moreover, previous devices experienced a steady loss in performance throughout their lifetime. The drop in performance has been attributed to the electrolyte undergoing thermal and chemical breakdown due to crack propagation, loss of sodium oxide, and the formation of molten dendrites that create electrical shorts [5]. Finally, previous designs have also had a substantial problem with proper thermal management, with some showing thermal losses as great as 50% of the heat input. To improve the efficiency, a dual stage conversion process is being explored. Rather than using one isothermal expansion step, the dual cycle will employ two separate electrolytes to allow two separate expansion steps: one at the high temperature, and another at an intermediate temperature. This new cycle is thermodynamically favorable as it allows for regeneration from the isobaric cooling of the first expansion step, and reheat to complete the expansion in the second step. This allows the engine to reach a lower condenser temperature than its single stage counterpart, thus increasing the overall efficiency. The temperature drop across each electrolyte in the dual stage device (e.g., ~ 400 K) is smaller than the temperature drop across a comparable single stage device (e.g., ~ 800 K). This is amenable for better thermal management and allows for less heat loss. This will also mitigate the thermal shock across the very thin (< 1 mm) solid electrolytes, reduce crack propagation and limit the percolation of the molten dendrites, and ultimately increase the lifetime of the device. As an additional benefit the dual stage can be designed using planer geometries, rather than tubular geometries used in previous single stage models, which can be manufactured less expensively. This modified geometry can be tailored to improve the thermal management by reducing the radiation losses, which has been the biggest source of thermal losses [2]. A dual stage concept for the Na-TECC has never been demonstrated in literature. In this presentation, the operation of a single stage device will be revisited and, by adopting concepts used for similar converters, the thermodynamic operation of the dual stage will be developed [6]. Using the dual stage approach, the efficiency and power equations have been redefined to include regeneration and reheat. Questions such as what current densities optimize the power output or which intermediate temperature maximizes the efficiency will also be answered. Initial calculations of the dual stage show that the thermal management is the most sensitive parameter. Although the ideal voltage produced by the dual stage can be higher, there are potential drawbacks. For example, by using a lower condenser temperature the dual stage device will have a higher open circuit voltage, but it will also have a larger overpotential due to the limited cathodic kinetics that occur at a lower temperature. In short, there are several trade-offs in this device that need to be balanced. It will be demonstrated that the dual stage device must have lower thermal losses than a single stage device if it is to have a higher efficiency.

14 SOLAR ENERGY↗
DOE OSTI2019

Size-resolved ice nucleating particle (INP) concentrations from the MOSAiC campaign

Sample collection: Ground-based immersion mode ice nucleation measurements were conducted on samples collected at the U.S. DOE ARM AMF2 onboard the R/V Polarstern P-deck during MOSAiC from 23 Oct 2019 to 01 Oct 2020. Samples were collected using a Davis Rotating-drum Unit for Monitoring (DRUM; DRUMAir, LLC) cascading impactor in 4 stages (i.e., size bins). Instrument information on the model used, the DA400, can be found at the DRUMAir website (http://www.drumair.com/). Stages A, B, C, and D had particle size cuts at 2.96, 1.21, 0.34, and 0.15 &micro;m in particle diameter, respectively. Strips of sample substrate (PFA, perfluoroalkoxy; McMaster) were adhered to each disc in each stage and coated with petrolatum to enable particles to stick on impact. The discs move slowly over time, such that aerosol loading is streaked onto the PFA. Every 24 hours, a blank is created on the PFA to separate daily samples. During MOSAiC, the pump typically pulled 20 &ndash; 35 lpm (average of 32 lpm) at the inlet with all the discs in place in each stage. Daily samples were collected for each of the four stages, equaling 28,800 &ndash; 50,400 total L of air (average of 46,080 L total L of air) per sample. Discs rotated for approximately 24 days before the sampling per strip was complete and changed to a new disc. Discs were stored frozen after collection. Sample preparation for analysis: Following collection and transport from MOSAiC to the Kreidenweis/DeMott laboratory, care of PI Creamean, in Fort Collins, CO, strips were removed from the discs and cut up to separate daily samples. Each daily portion of the strip was stored in separate sterile 15-mL centrifuge tubes and stored in a standard commercial freezer until analysis (approximately 4 &ndash; 12 months after collection). Sample analysis: Drop freezing assays (DFAs) were conducted on one daily sample set every 3 &ndash; 4 days of the study for all four stages using the Colorado State University (CSU) cold plate (CP). Sample set days were chosen to cover the entire year of MOSAiC and avoiding days where the AOS purge system was activated. Immediately after removal from the freezer, 2 mL of 0.1-&micro;m filtered deionized water was added with a single-use sterile syringe directly to the tubes, then mixed using a vortex mixer for at least 30 minutes at 500 rpm to resuspend particles from the PFA into the deionized water. Following sample preparation, another single-use syringe was used to create 100 aliquots on a 3-inch diameter copper plate, then covered with a transparent plastic dome. The drops were approximately the same volume (2.5 &micro;L each), though careful inspection. The plate was cooled at approximately 3 &ndash; 10 K min-1 from room temperature until all drops froze on the plate. Temperature was measured by an Omega thermocouple probes and meter. Drop freezing was visually detected but recorded through monitoring software to provide the time frozen, channel of the meter used, temperature, and cooling rate for each drop. Not all 100 drops were always detected; the total number of rows in each data file equals the number of drops recorded, which typically was > 80 %. Each sample was tested three times (i.e., a new set of 100 drops was created for each test). Data provided are files containing these parameters. From the freezing temperatures recorded, the fraction frozen was and INP concentrations were calculated per L of air using the equation from Vali (1971): [INPs] = &minus;(ln(f)/Valiquot)*(Vsuspension/Vair) where f is the proportion of droplets not frozen, Valiquot is the volume of each aliquot (i.e., 2.5e-6 L), Vsuspension is the volume of the suspension (i.e., 0.002 L), and Vair is the volume of air per sample (i.e., 36000 L). Data usage guidance: There are three files included with this dataset: (1) dumcp_inp_concentrations.csv which includes cumulative INP concentrations binned every 0.5 &deg;C from 0 to &ndash;30 &deg;C for all samples and size ranges (i.e., DRUM stages) processed via DFA using the CSU CP; (2) dumcp_inp_standard_deviation.csv which includes the standard deviations for all samples and temperatures in dumcp_inp_concentrations.csv; and (3) dumcp_blank_concentrations_standard_deviation.csv which includes the INP concentrations and standard deviations for the deionized water and PFA blanks binned every 0.5 &deg;C from 0 to &ndash;30 &deg;C. The first column in each file represents the freezing temperature in &deg;C binned every 0.5 &deg;C from 0 to &ndash;30 &deg;C (first cell in this column = header = &ldquo;temp&rdquo;). The second column in each file represents the standard deviation of each temperature bin in &deg;C (first cell in this column = header = &ldquo;temperr&rdquo;). For file (1) dumcp_inp_concentrations.csv, the remaining columns indicate the INP concentrations per L of air for each temperature bin. For each sample column, the header in the first row indicates the stage (i.e., size bin as &ldquo;A_&rdquo;, &ldquo;B_&rdquo;, &ldquo;C_&rdquo;, or &ldquo;D_&rdquo;) and date (&ldquo;_yyyymmdd) of collection. For file (2) dumcp_inp_standard_deviation.csv, the remaining columns indicate the INP concentration standard deviation per L of air for each temperature bin. For each sample column, the header in the first row indicates the stage (i.e., size bin as &ldquo;A_&rdquo;, &ldquo;B_&rdquo;, &ldquo;C_&rdquo;, or &ldquo;D_&rdquo;) and date (&ldquo;_yyyymmdd) of collection. For file (3) dumcp_blank_concentrations_standard_deviation.csv, the remaining columns indicate the INP concentrations and standard deviations for the deionized water blank (&ldquo;di_0.1&rdquo; and &ldquo;di_0.1sd&rdquo; for the concentrations and standard deviations, respectively) and PFA blank (&ldquo;pfa&rdquo; and &ldquo;pfasd&rdquo; for the concentrations and standard deviations, respectively). Please see Creamean et al. (2019a; 2018a; 2018b; 2019b) for additional details on DRUM sample collection and offline processing. For questions or help with data interpretation, contact Dr. Jessie Creamean (jessie.creamean@colostate.edu). References: Creamean, J. M., Mignani, C., Bukowiecki, N., and Conen, F.: Using spectra characteristics to identify ice nucleating particle populations during winter storms in the Alps, Atmos. Chem. Phys., 19, 8123&ndash;8140, 2019a. Creamean, J. M., Kirpes, R. M., Pratt, K. A., Spada, N. S., Maahn, M., de Boer, G., Schnell, R. C., and China, S.: Marine and terrestrial influences on ice nucleating particles during continuous springtime measurements in an Arctic oilfield location, Atmospheric Chemistry and Physics, 18, 18023-18042, 10.5194/acp-18-18023-2018, 2018a. Creamean, J. M., Primm, K. M., Tolbert, M. A., Hall, E. G., Wendell, J., Jordan, A., Sheridan, P. J., Smith, J., and Schnell, R. C.: HOVERCAT: A novel aerial system for evaluation of aerosol-cloud interactions, Atmospheric Measurement Techniques, 11, 3969-3985, 10.5194/amt-11-3969-2018, 2018b. Creamean, J. M., Cross, J. N., Pickart, R., McRaven, L., Lin, P., Pacini, A., Hanlon, R., Schmale, D. G., Ceniceros, J., Aydell, T., Colombi, N., Bolger, E., and DeMott, P. J.: Ice nucleating particles carried from below a phytoplankton bloom to the Arctic atmosphere, Geophys. Res. Lett., 46, 8572-8581, https://doi.org/10.1029/2019GL083039, 2019b. Vali, G.: Quantitative Evaluation of Experimental Results an the Heterogeneous Freezing Nucleation of Supercooled Liquids, Journal of the Atmospheric Sciences, 28, 402-409, 10.1175/1520-0469(1971)028<0402:qeoera>2.0.co;2, 1971.

4-stage DRUM aerosol collector (DA400),Ice nuclei,↗