Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compiler optimization”

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

Toward Performance Portable Programming for Heterogeneous System-on-Chips: Case Study with Qualcomm Snapdragon SoC

Future heterogeneous Domain-Specific System-on-Chips (DSSoC) will be extraordinarily complex in terms of processors, memory hierarchies, and interconnection networks.To manage this complexity, architects, system software designers, and application developers need programming technologies that are flexible, accurate, efficient, and productive. These technologies will need to be as independent of any one specific architecture as is practical, because the sheer dimensionality and scale of the complexity will not allow porting and optimizing applications foreach given DSSoC. To address these issues, we are developing Cosmic Castle, a performance portable programming toolchain for streaming applications on heterogeneous architectures. The primary focus of Cosmic Castle is on enabling efficient and performant code generation through the smart compiler and intelligent runtime system. This paper presents the preliminary evaluation of our ongoing work toward Cosmic Castle. Specifically, we detail our code porting efforts and evaluate various benchmarks on the Qualcomm Snapdragon SoC using tools developed through Cosmic Castle.

Cabrera, Anthony↗
DOE OSTI2022

Scientific Computational Imaging Code (SCICO)

Scientific Computational Imaging Code (SCICO) is a Python package for solving the inverse problems that arise in scientific imaging applications. Its primary focus is providing methods for solving ill-posed inverse problems by using an appropriate prior model of the reconstruction space. SCICO includes a growing suite of operators, cost functionals, regularizers, and optimization routines that may be combined to solve a wide range of problems, and is designed so that it is easy to add new building blocks. SCICO is built on top of JAX rather than NumPy, enabling GPU/TPU acceleration, just-in-time compilation, and automatic gradient functionality, which is used to automatically compute the adjoints of linear operators. An example of how to solve a multi-channel tomography problem with SCICO is shown in Figure 1. The SCICO source code is available from GitHub, and pre-built packages are available from PyPI. It has extensive online documentation, including API documentation and usage examples, which can be run online at Google Colab and binder.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Evaluating integration and performance of containerized climate applications on a Hewlett Packard Enterprise Cray system

Containers have taken over large swaths of cloud computing as the most convenient way of packaging and deploying applications. The features that containers offer for packaging and deploying applications translate to high performance computing (HPC) as well. At The National Oceanic and Atmospheric Administration, containers provide an easy way to build and distribute complex HPC applications, allowing faster collaboration, portability, and experiment computer environment reproducibility amongst the scientific community. The challenge arises when applications rely on message passing interface (MPI). This necessitates investigation into how to properly run these applications with their own unique requirements and produce performance on par with native runs. We investigate the MPI performance for benchmarks and containerized climate models for various containers covering selection of compiler and MPI library combinations from the Cray provided programming environments on the Cray XC supercomputer GAEA. Performance from the benchmarks and the climate models shows that for the most part containerized applications perform on par with the natively built applications when the system optimized Cray MPICH libraries are bound into the container, and the hybrid model containers have poor performance in comparison. We also describe several challenges and our solutions in running these containers, particularly challenges with heterogeneous jobs for the containerized model runs.

Abraham, Subil↗
DOE OSTI2022

Developing an ELM Ecosystem Dynamics Model on GPU with OpenACC

Porting a complex scientific code, such as the E3SM land model (ELM), onto a new computing architecture is challenging. The paper presents design strategies and technical approaches to develop an ELM ecosystem dynamics model with compiler directives (OpenACC) on NVIDIA GPUs. The code has been refactored with advanced OpenACC features (such as deepcopy and routine directives) to reduce memory consumption and to increase the levels of parallelism through parallel loop reconstruction and new data structures. As a result, the optimized parallel implementation achieved more than a 140-time speedup (50 ms vs 7600 ms), compared to a naive implementation that uses OpenACC routine directive and parallelizes the code across existing loops on a single NVIDIA V100. On a fully loaded computing node with 44 CPUs and 6 GPUs, the code achieved over a 3.0-times speedup, compared to the original code on the CPU. Furthermore, the memory footprint of the optimized parallel implementation is 300 MB, which is around 15% of the 2.15 GB of memory consumed by a naive implementation. This study is the first effort to develop the ELM component on GPUs efficiently to support ultra-high-resolution land simulations at continental scales.

Schwartz, Peter↗
DOE OSTI2023

A semi-agnostic ansatz with variable structure for variational quantum algorithms

Quantum machine learning—and specifically Variational Quantum Algorithms (VQAs)—offers a powerful, flexible paradigm for programming near-term quantum computers, with applications in chemistry, metrology, materials science, data science, and mathematics. Here, one trains an ansatz, in the form of a parameterized quantum circuit, to accomplish a task of interest. However, challenges have recently emerged suggesting that deep ansatzes are difficult to train, due to flat training landscapes caused by randomness or by hardware noise. This motivates our work, where we present a variable structure approach to build ansatzes for VQAs. Our approach, called VAns (Variable Ansatz), applies a set of rules to both grow and (crucially) remove quantum gates in an informed manner during the optimization. Consequently, VAns is ideally suited to mitigate trainability and noise-related issues by keeping the ansatz shallow. We employ VAns in the variational quantum eigensolver for condensed matter and quantum chemistry applications, in the quantum autoencoder for data compression and in unitary compilation problems showing successful results in all cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Data driven investigation to understand the influence of total solids on biological biogas upgrading

In situ biogas upgrading achieves CO 2 conversion to CH 4 via hydrogenotrophic methanogenesis; however, gas-liquid mass transfer constraints limit the upgrading performance. Recognizing that optimization studies often underrepresent the effects of total solids (TS) and organic loading rate (OLR), this study undertook a holistic, statistics driven assessment of operating conditions for in situ H 2 assisted biogas upgrading, centering the analysis on TS and OLR. A dataset of 31 studies was compiled and comprised 99 observations. A rigorous analytical framework was employed, combining data standardization, fixed- and random-effects (REML) weighted regressions with cluster-robust errors, stratified analyses, and machine learning. Mixed-effects meta regression indicated that TS was the main factor explaining differences of methane fraction (CH 4 %) when considering the between studies heterogeneity. Focusing on a near-stoichiometric subset (H 2 /CO 2 ≈ 4:1), TS remained significant. Stratified results showed a stronger negative relationship between TS and CH 4 % in UASB reactors than in CSTRs, with a negative effect under mesophilic conditions and no significant effect under thermophilic conditions. A Random Forest model corroborated the statistical findings, consistently ranking H 2 /CO 2 ratio, OLR, TS, and hydrogen injection rate (HIR) as the most influential predictors. These findings delineate trends across increasing TS levels, particularly between 1% and 10%, and provide preliminary insights for TS above 15% in in situ biogas upgrading. They further provide insights for the influence of TS by reactor type and temperature, thereby advancing the evidence base for implementing biological CO 2 conversion to CH 4 in practice.

In situ biogas upgrading↗
DOE OSTI2023

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). 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 8. 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 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 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 4. 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 5. 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 6. 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 7. 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 8. OUTCARL: 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 9. 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 10. 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 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 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.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tailoring MoS 2 for Small-Molecule Electroreduction: The Role of Metal Doping and Heterostructures

The electrification of chemical transformations central to sustainable fuel production and waste valorization, such as overall water splitting (OWS), hydrogen evolution reaction (HER), and electrochemical reduction of CO 2 (CO 2 R), presents a powerful opportunity to advance carbon-neutral energy technologies. Transition metal dichalcogenides (TMDs), particularly MoS 2 , have emerged as promising electrocatalyst candidates, owing to their abundance, tunable active sites, and defect-rich structures. This review highlights recent progress in leveraging metal doping and heterostructure engineering of MoS 2 to enhance the electrocatalytic activity and selectivity. By compiling insights from experimental studies and density functional theory (DFT) predictions, we examine how defect creation, electronic structure modification, and interface design contribute to improved charge transport and catalytic efficiency. Particular emphasis is placed on rational design principles, synthetic strategies, and operando characterization methods that provide a pathway to understanding and optimizing MoS 2 -based materials. We also discuss the challenges of stability, mechanistic ambiguity, and scaling while outlining opportunities to bridge theory and experiment. Collectively, this review underscores how defect and heterostructure engineering of MoS 2 can accelerate the development of efficient, sustainable electrocatalysts for both fuel generation and waste-to-value generation.

CO2 reduction↗
DOE OSTI2022

6.26 Low Cost Basalt Fiber for Automotive Applications

Vehicle lightweighting is an essential component to the automotive industry to improve fuel economy of internal combustion engine (ICE) vehicles to meet ever improving emission standards and to improve the range of electric vehicles (EV). Composite materials offer high specific modulus and specific strength, which makes them appealing for these light weighting efforts. Sheet molding compounds (SMC) are particularly interesting from an automotive perspective because of the relatively low cost and high volume of producing SMC parts. Traditionally, composite materials for automotive application are glass fiber reinforced because of the attractive price - performance ratio, but basalt fibers are a cost and recycling competitive reinforcement alternative in this market. The aim of this project was to examine the feasibility of utilizing basalt fiber for automotive applications. More specifically, an effort was made to examine different fiber sizings on basalt fiber combined with vinyl-ester (VE) resin, and their performance as part of an SMC process. In addition to offering vehicle lightweighting with fiber reinforced polymer composites, basalt fiber is a fully recyclable material and thus supports the IACMI technical goal of: Demonstrate that the technology is capable at a sufficient scale for >80% recyclability or reuse of fiber reinforced polymer composites in five years into useful components with projected cost and quality at commercial scale competitive with virgin materials on a pathway to 95% recyclability or reuse starting in ten years. Three different fiber sizings were applied to a continuous roving of basalt fiber and compared to a traditional Electrical/Chemical Resistance (E-CR) glass fiber that is typically used in these types of applications. Fiber tows were examined for Loss on Ignition percentage (LOI%), Tex, and tow strength. Some sizings clearly outperformed others, and the ability to process these fibers on a pilot scale SMC line was demonstrated. A test plan for the manufacturing and mechanical testing of SMC panels was developed. This work continues outside the time frame allocated for this project. When this work is completed, it will be added to this report and posted as Appendix C. Glass fiber reinforced SMC materials have already proven feasible as a light weighting method for traditionally steel parts like the Volkswagen (VW) Atlas lift-gate (Figure 1); this project team is seeking the feasibility of basalt fiber as a drop-in replacement for glass fiber reinforced SMC. Sizing development for basalt fibers has proven that the mechanical properties are better than E-glass and closer to S-glass, which makes it an interesting material for SMC applications. Better mechanical properties translate to less material needed to achieve load case requirement for target applications. The business case has already been demonstrated for 100,000 parts per year of glass fiber reinforced SMC Atlas lift-gates compared to traditional steel manufacturing processes. Reduced overhead and assembly costs are offset by glass fiber SMC higher cost per kg beyond 100,000 parts per year, which is still a relatively low volume for the automotive industry. For basalt fiber reinforced SMC to become feasible for automotive applications, the price-performance ratio has to be precisely determined. Based on the mechanical performance it is possible to establish a range of applications and technical solutions in which the potential of basalt SMC can be utilized, while the price of the material can be used to compile the business case for such applications. Based on these business cases and the sustainability indicators, glass fiber reinforcement (or other) materials can be directly substituted. Volkswagen’s commitment to reducing carbon emissions cannot be understated. Basalt fiber shows promise of reducing the carbon footprint in SMC materials, especially if sizing optimizations can be made with thermoplastic based SMC. To fully realize the value of basalt fiber reinforced materials, a lifecycle cost analysis should be performed on basalt’s production and recycling, and then compared against E-glass. From this assessment, a true judgement can be made on the commercialization potential of this material. Figure 1. Example of Fiber Reinforced Polymer Composite Liftgate As a conclusion, we can state that Mafic basalt fiber is not a direct replacement for E-glass or E-CR glass based on price, but should be considered a technical solution when E-glass does not provide adequate performance in a composite design and S-glass, aramid and carbon fibers are too costly. Mafic basalt fiber can be placed on the high-performance fiber spectrum next to S-glass for performance but at one third the price. It should be considered for more technically challenging structural designs wherein the performance can demonstrate 20-25% performance enhancement over E-glass to elicit more strength or a weight reduction. Both Michelman and Mafic produce thermoplastic sizings which, in combination with Nylon and polypropylene resin and fibers, can further advance high speed composite implementations while maintaining an eco-friendly manufacturing process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

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

Promising approaches for the assembly of the catalytically active, recombinant Desulfomicrobium baculatum hydrogenase with substitutions at the active site

Background: Hydrogenases (H2ases) are metalloenzymes capable of the reversible conversion of protons and electrons to molecular hydrogen. Exploiting the unique enzymatic activity of H2ases can lead to advancements in the process of biohydrogen evolution and green energy production. Results: Here we created of a functional, optimized operon for rapid and robust production of recombinant [NiFe] Desulfomicrobium baculatum hydrogenase (Dmb H2ase). The conversion of the [NiFeSe] Dmb H2ase to [NiFe] type was performed on genetic level by site-directed mutagenesis. The native dmb operon includes two structural H2ase genes, coding for large and small subunits, and an additional gene, encoding a specific maturase (protease) that is essential for the proper maturation of the enzyme. Dmb, like all H2ases, needs intricate bio-production machinery to incorporate its crucial inorganic ligands and cofactors. Strictly anaerobic, sulfate reducer D. baculatum bacteria are distinct, in terms of their biology, from E. coli. Thus, we introduced a series of alterations within the native dmb genes. As a result, more than 100 elements, further compiled into 32 operon variants, were constructed. The initial requirement for a specific maturase was omitted by the artificial truncation of the large Dmb subunit. The assembly of the produced H2ase subunit variants was investigated both, in vitro and in vivo. This approach resulted in 4 recombinant [NiFe] Dmb enzyme variants, capable of H 2 evolution. The aim of this study was to overcome the gene expression, protein biosynthesis, maturation and ligand loading bottlenecks for the easy, fast, and cost-effective delivery of recombinant [NiFe] H2ase, using a commonly available E. coli strains. Conclusion: The optimized genetic constructs together with the developed growth and purification procedures appear to be a promising platform for further studies toward fully-active and O 2 tolerant, recombinant [NiFeSe] Dmb H2ase, resembling the native Dmb enzyme. It could likely be achieved by selective cysteine to selenocysteine substitution within the active site of the [NiFe] Dmb variant.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Efficient Trainability of Linear Optical Modules in Quantum Optical Neural Networks

The existence of “barren plateau landscapes” for generic discrete-variable quantum neural networks, which obstructs efficient gradient-based optimization of cost functions defined by global measurements, would be surprising in the case of generic linear optical modules in quantum optical neural networks due to the tunability of the intensity of continuous variable states and the relevant unitary group having exponentially smaller dimension. We demonstrate that coherent light in m modes can be generically compiled efficiently if the total intensity scales sublinearly with m, and extend this result to cost functions based on homodyne, heterodyne, or photon detection measurement statistics, and to noisy cost functions in the presence of attenuation. Here, we further demonstrate efficient trainability of m mode linear optical quantum circuits for variational mean field energy estimation of positive quadratic Hamiltonians for input states that do not have energy exponentially vanishing with m.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Adaptive job and resource management for the growing quantum cloud

As the popularity of quantum computing continues to grow, efficient quantum machine access over the cloud is critical to both academic and industry researchers across the globe. And as cloud quantum computing demands increase exponentially, the analysis of resource consumption and execution characteristics are key to efficient management of jobs and resources at both the vendor-end as well as the client-end. While the analysis and optimization of job / resource consumption and management are popular in the classical HPC domain, it is severely lacking for more nascent technology like quantum computing.This paper proposes optimized adaptive job scheduling to the quantum cloud taking note of primary characteristics such as queuing times and fidelity trends across machines, as well as other characteristics such as quality of service guarantees and machine calibration constraints. Key components of the proposal include a) a prediction model which predicts fidelity trends across machine based on compiled circuit features such as circuit depth and different forms of errors, as well as b) queuing time prediction for each machine based on execution time estimations. Altogether, this proposal is evaluated on simulated IBM machines across a diverse set of quantum applications and system loading scenarios, and is able to reduce wait times by over 3x and improve fidelity by over 40% on specific usecases, when compared to traditional job schedulers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

High Performance Approximate Computing

This code repository contains the implementation of the "High-Performance Approximate Computing" (HPAC) toolkit. The toolkit allows you to approximate your own C/C++. The developer uses "pragma's" to annotate code regions as approximate. The compiler extensions lower these pragmas to either compiletime approximate techniques or runtime approximation techniques. At execution time, the implemented runtime system decides which annotated regions it should approximate. HPAC also provides a set of script utilities. The utilities perform a grid search within approximation parameters and performance. The user can analyze the raw data to identify optimal approximation techniques for the application.

Parasyris, Konstantinos↗
DOE OSTI2020

Conceptual Model Update of Mercury Sources and Flux at Y-12 and Upper East Fork Poplar Creek, Oak Ridge,Tennessee

The interactions between the key processes that control the fate and transport of mercury (Hg) in the sediments, groundwater, and stream water within the Y-12 National Security Complex (Y-12) in Oak Ridge, Tennessee are complex, and many aspects are not well understood. In contaminated environments such as Y-12, conceptual models are frequently developed to aid in visualizing and understanding the dynamic nature of the hydrologic, geochemical, and physical environment. These conceptual models integrate data in an internally consistent manner to understand processes that control the fate and transport of contaminants. Over the past few decades of environmental investigation at Y-12, a number of conceptual models have been developed to identify and define various technical processes at various scales. The main purpose of the activities described in this report is to update the 2011 conceptual model for Upper East Fork Poplar Creek using the most recent Hg concentration and flux data and informed scientific interpretation. Emphasis was placed on improving the previous model by implementing a more detailed spatial and temporal approach to visualize transport pathways in the watershed and trends in flux and concentration over time, and to compare baseflow and stormflow system dynamics. Detailed descriptions of historical and current Hg sources and transport pathways are also provided for the east and west ends of the facility. This conceptual model will allow the US Department of Energy to evaluate past and present remedial activities and provide a strong technical basis for prioritizing and optimizing remedial responses in a cost-effective and efficient manner. Twenty years ago, Hg flux from Outfall (OF) 200 at the headwaters of East Fork Poplar Creek represented approximately 20% of the overall flux leaving Y-12. By the time of the 2011 conceptual model report, that percentage had increased to 70%–80%. Flux estimates compiled for this report using data from 2009 to 2018 show that OF200 accounted for roughly 60% of the flux leaving Y-12. The relative role of the four storm drain conduits to OF200 appeared similar to 2011, with OF163 being the greatest contributor to downstream flux. An increase in Hg concentrations and flux occurred throughout and downstream of west end storm drains in 2011 in response to a storm drain cleanout. Leading up to the 2011 storm drain cleanout, annual baseflow flux was higher at OF200A6 than at Station 17. From 2011 to 2018, the trend reversed with concentrations and flux being higher at Station 17 than at OF200A6. A spike in concentrations occurred in July 2018, potentially due to a one-time influx of Hg that occurred during COLEX decontamination and decommissioning activities at Y-12. This report includes recommendations for improving our conceptual understanding of Hg sources, transport pathways, and flux at Y-12. Important recommendations include further evaluation of transport pathways, concurrent measurement of Hg concentration and flux under baseflow and stormflow conditions at multiple sites, and research to understand the connections among periphyton and methylmercury concentration, bioavailability, and bioaccumulation. Ongoing facility decontamination and decommissioning efforts, site characterization, remedial actions, and research are facilitating the collection of new Hg data at Y-12. The conceptual model update provided herein takes advantage of the extensive recent Hg sampling to provide an updated assessment of Hg mass balance and support a refined understanding of Hg behavior at or near Y-12. This assessment should assist in future environmental management decisions and in mitigating the impacts of Hg on the surrounding environment.

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

Numerical Simulations of Noisy Variational Quantum Eigensolver Ansatz Circuits

This is a case study of the variational quantum eigensolver (VQE) method using numerical simulations to test the influence of noise on the accuracy of the underlying circuit ansatz. We investigate a computational chemistry application of VQE to calculate the electronic ground state and its energy for Sodium Hydride (NaH), a prototypical two-electron problem. Using a one-parameter ansatz derived from unitary coupled cluster (UCC) theory, we simulate the effects of noise on the energy expectation value and variance with respect to the ansatz parameter. These numerical simulations provide insights into the accuracy of the prepared quantum state and the efficiency of the classical optimizer that iteratively refines the ansatz. We conduct a comparative study between analytical results derived for the UCC ansatz in the absence of noise and the noisy numerical simulation results obtained using an isotropic depolarizing noise model for each gate. We also compare the relative increase in noise on logically equivalent UCC ansatz circuits generated by randomized compiling. Notably, we observe that the intrinsic variance in the energy due to the simplicity of the ansatz itself compares with the noise induced by the bare circuit.

Gowrishankar, Meenambika↗
DOE OSTI2021

Development of Energy Storage: Cost models

Energy storage technologies offer a promising solution to electric grid stability issues associated with the integration of variable renewable generators. The capability to match the electrical power output to instantaneous fluctuations in grid demand is crucial to ensure continuity of service. Including energy storage capability in an integrated energy system (IES) can provide the flexibility needed to meet variable electric demand and reduce the load following demands place on the reactor. In this report, economic data collected for energy storage (ES) technologies are described to support the objective of assessing the profitability of ES integration within the IES framework. In particular, there is a growing interest in thermal energy storage (TES) given its unique capability for long-duration storage for improving electricity reliability at a low levelized cost. It is common practice to evaluate the total lifetime cost and profitability before commercializing new technologies. In this report we identify and examine the models needed to better understand the economics of thermal energy storage. We extend the TES cost model in RAVEN in the context of a balance of plant (BOP) that incorporates thermal storage. To focus the discussion, following a general overview of the most promising TES technologies, we consider a use case that involves a sensible heat, two-tank, molten-salt system. Structural and operational details are reported to identify the source of construction capital expenditure and operation and maintenance cost. A detailed description of the different cost items is provided as well as the cost scaling with different storage capacity and power ratings for capacity optimization purpose. In addition, the capital expenditure and the recurring cost of representative two-tank, molten salt coupled with concentrated solar plants are provided for readers’ reference. The ES use case is noteworthy as it is in the pilot stage of commercialization. We identify those areas that would benefit from an increased economic focus to obtain a more complete compilation of cost data. We also describe the thermal coupling issues that arise from integration of the two-tank molten salt thermal energy storage system with a BOP, which is the subject of our current research. Some components of costs will need to be evaluated through dedicated technoeconomic analysis in future modeling activities using modeling procedures proposed in this report. With the data presented and the procedures described in this report, sufficiently accurate models can be implemented for the solution of both the power dispatch and the capacity expansion problems within the RAVEN-based HYBRID framework.

25 ENERGY STORAGE↗
DOE OSTI2023

Preliminary Enhanced LAW Glass Formulation Algorithm

This report summarizes the Preliminary Enhanced Low-Activity Waste Glass Formulation Algorithm (GFA), its background information, and the calculations it performs. The Preliminary Enhanced Low-Activity Waste GFA is a tool developed in MATLAB to formulate glass for a given waste composition while attempting to maximize waste loading. It is intended for use at the Waste Treatment and Immobilization Plant, where nuclear waste will be vitrified into glass. The formulated glass is required to satisfy several processing and product quality constraints. In addition, calculations must account for associated uncertainties in constraint prediction and measurement. The GFA adheres to Pacific Northwest National Laboratory nuclear quality assurance procedures and has been validated and verified. In this second revision to the report, two constraints: the lithium and zirconia target constraints were amended to Table 3.1 and programmed into the algorithm. Additionally, three sections were added. Section 3.2.1 includes the Enhanced LAW Correlation Rules, an addition to the glass optimality criteria. Section 4.0 discusses a stepped process approach taken to ensure the algorithm can adequately react to differences between its recommended output and the actual measurements and transfers made during plant operation, and calculation of specific radionuclide concentrations. Section 5.2 describes changes to the user interaction with the algorithm program including the addition of a graphical user interface and compilation of the software.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗