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

Cybersecurity for Grid Connected eXtreme Fast Charging (XFC) Station (CyberX) (Final Scientific/Technical Report)

This report summarizes the activities conducted under the DOE VTO funded project DE- EE0008451, where ABB Inc. (ABB), in collaboration with Idaho National Laboratory (INL), APS Global (APS), and XOS Trucks (XOS) pursued the development of a cyber-resilient extreme fast charging (XFC) management system. This project entitled Cybersecurity for Grid Connected eXtreme Fast Charging (XFC) Station (CyberX) focuses on a resilient architecture for smart charging EV Supply Equipment (EVSE) device control and Coordinated Anomaly Detection System (CADS) features that can be added at the charging site depot level to increase cybersecurity. The project was split into two budget periods focused first on developing the threat model and resilient control concepts and second on testing, improving, and validating those developed resilient control algorithms and features with a focus on key vulnerabilities identified during the threat assessment portion of the project. During the first budget period of the CyberX project, the ABB led team focused on activities to identify, model, and quantitatively prioritize high-impact attack scenarios with potential cyber-physical effects while also modeling and developing concepts for a resilient control system that could securely address integration of DERs and other resources with EV charging. Development of the security focused XFC management system (XMS) was accomplished first by offline simulation using a developed XFC station or depot with 480V input level and simulating measurement inputs to monitoring and control systems in concept development. A representative distribution grid model was developed supporting an EV charging site model with BESS and 6 general EV charging models. These EV charging models allowed multiple configurations of charging level, multiple connected protection and measurement devices, and simulation function to show general compromise of EV, BESS, and protection features based on parallel threat analysis. During the second budget period, the EV site and supporting systems model was developed in more detail and converted from offline model to real-time to real-time with EV charging hardware in the loop (HIL). The resilient control architecture developed as concept in the first part of the project was further tested and validated for integration of local energy resources and XFC charging station site equipment while maintaining cybersecure operating principles. The proposed resilient architecture for smart charging and cybersecurity features consists of two main concepts developed and tested within the project. The first concept is an XFC management system (XMS) consisting of a hardware gateway, software platform, and Supervisory Control and Data Acquisition (SCADA) or Distribution Management System integration components. The second concept is a Coordinated Anomaly Detection System (CADS) which forms a primarily software-related subsystem of the total CyberX solution focused on monitoring system measurements, estimation of measurement states, and predicting current at the utility point of interaction based on machine learning for anomaly detection.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2012

Fate and transport of phenol in a packed bed reactor containing simulated solid waste

Highlights: Black-Right-Pointing-Pointer Anaerobic column experiments were conducted at 37 Degree-Sign C using a simulated waste mixture. Black-Right-Pointing-Pointer Sorption and biodegradation model parameters were determined from batch tests. Black-Right-Pointing-Pointer HYDRUS simulated well the fate and transport of phenol in a fully saturated waste column. Black-Right-Pointing-Pointer The batch biodegradation rate and the rate obtained by inverse modeling differed by a factor of {approx}2. Black-Right-Pointing-Pointer Tracer tests showed the importance of hydrodynamic parameters to improve model estimates. - Abstract: An assessment of the risk to human health and the environment associated with the presence of organic contaminants (OCs) in landfills necessitates reliable predictive models. The overall objectives of this study were to (1) conduct column experiments to measure the fate and transport of an OC in a simulated solid waste mixture, (2) compare the results of column experiments to model predictions using HYDRUS-1D (version 4.13), a contaminant fate and transport model that can be parameterized to simulate the laboratory experimental system, and (3) determine model input parameters from independently conducted batch experiments. Experiments were conducted in which sorption only and sorption plus biodegradation influenced OC transport. HYDRUS-1D can reasonably simulate the fate and transport of phenol in an anaerobic and fully saturated waste column in which biodegradation and sorption are the prevailing fate processes. The agreement between model predictions and column data was imperfect (i.e., within a factor of two) for the sorption plus biodegradation test and the error almost certainly lies in the difficulty of measuring a biodegradation rate that is applicable to the column conditions. Nevertheless, a biodegradation rate estimate that is within a factor of two or even five may be adequate in the context of a landfill, given the extended retention time and the fact that leachate release will be controlled by the infiltration rate which can be minimized by engineering controls.

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

On the influence of the power plant operational history on the inventory and the uncertainties of radionuclides relevant for the final disposal of PWR spent fuel - 15149

A mandatory requisite for official granting safe final disposal approval of SNF is the comprehensive specification and declaration of the nuclear inventory in SNF by the NPP-operators. The whole radiochemical characterization of each spent fuel assembly hides many difficulties such as very high radiation burden, destruction method, the inhomogeneous burn-up of fuel assembly, very high cost. Thus, non-destructive characterization is usually applied on industrial scale. In this case only several nuclides can be measured, and the rest is calculated using simulations and typical irradiation history. Like any other German rad-waste acceptance requirements all these declarations will have to be checked and verified by an independent inspector. In this verification process not only values of radionuclide (RN) activities are important, but also their uncertainties. The simple simulation using burn-up calculations based on typical and generic reactor operational parameters do not answer the question about possible uncertainties observed in real reactor operations. At the same time, the details of irradiation history are often unknown, which complicates the assessment of declared RN inventories. Therefore, the development of dedicated generic proof tools, which can be applied to different types of waste, often of unknown origin and history is worth a serious consideration. Here, we have compiled a set of burn-up calculations considering the operational history of 339 published or anonymized real PWR fuel assemblies. This gives us the possibility to calculate the realistic variation or spectrum of RN inventories. These histories were used as a basis for 'SRP analysis', to provide information about the range of the values of the associated secondary reactor parameters (SRP's), i.e. initial enrichment (IE), specific power (SP), irradiation time (IT), downtime (DT), boric acid concentration (BA), moderator and fuel density (MD and FD) and temperatures (MT and MT) used as input parameters for the simulation models. Then, burn-up calculation codes provide the anticipated results for the entire nuclear inventory such as characteristic activity values and the corresponding minimum and maximum limits as a function of the fuel burn-up. These results are being validated with the experimental data from various publications and the online database - SFCOMPO and subsequently transferred and summarized into our own database. Ultimately, this database will provide the basis of proof tools for the operator independent verification of the declared nuclear inventories for the quality and safety assessment of a radioactive waste compound and bounding values for certain RN as a function of burn-up. The RN inventory can hence be verified without taking credit of the whole of the operational history of the SNF. (authors)

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

Case Study: Nuclear-Renewable-Water Integration in Arizona

This document reports the application of the N-R HES software framework to a case study for Arizona Public Service (APS); the manager and part owner of the Palo Verde (PV) nuclear power plant. The case study is a work in progress of which this report presents a detailed description of the current model input data, assumptions and the corresponding results produced by the developed software framework. The goal of the report is, that with this information APS (together with INL) should be able to identify where more detailed data and more accurate model assumptions are necessary to reproduce all driving physical and economic phenomena in order to capture enough complexity of the reality to allow APS to use the results for their strategic decisions. The N-R HES software framework started being developed at INL two years ago [ref]. The framework has reached some level of maturity so that it can be applied to more than simple demonstration cases, e.g. real industry problems. Nevertheless, more capabilities are constantly added to accommodate the special needs of these challenging real life problems. The N-R HES framework is built on top of the RAVEN code [ref] and uses it as driver and workflow manager for all calculations. The framework has specifically been developed for the economic assessment of N-R HES systems. There are four main cornerstones of the N-R HES simulation framework: 1) generation of stochastic time series, 2) a probabilistic analysis and optimization set of algorithms available in RAVEN, 3) a set of models for representation of the physical behavior of N-R HES, and 4) a RAVEN plug-in that maps physical performance into economic performance called CashFlow. Within this framework, a broad spectrum of questions related to N-R HES can be addressed. One of the challenges currently of high interest is that the increasing penetration of variable renewables is altering the profile of the net demand (demand after removing all non curtailable renewable energy sources), with which the other generators on the grid have to cope. The N-R HES software framework is capable to analyze the potential feasibility of mitigating the resultant volatility in the net electricity demand. One possible solution to that problem currently intensively studied by the energy industry is adding stabilizing loads to the grid. These loads can be external industrial processes that are able to ramp production up and down quickly or adding variable internal loads that will lift the base load of the power plant. The latter is a possible solution currently studied by APS. APS if currently anticipating several challenges that will likely arise in the near future. One of them is, as mentioned, coping with the rapid growth of variable renewable energy (VRE) sources put on the grid. Although APS can sell its electricity for a fixed retail price to cover the local demand first, it also sells electricity in excess of the local demand at the Palo Verde energy hub. Due to the growing VRE penetration, the Palo Verde hub electricity spot price is negative more and more frequently. To mitigate this price volatility, APS is seeking to add more base load, so that the quantity of excess energy to be sold for prices potentially less than the internal retail price or even negative prices at the hub is reduced. A second challenge APS is facing is that their cooling water acquisition contract with the Sub Regional Operating Group (SROG) will expire soon and a renewal can only be done for a significantly higher price of the water. Therefore, APS is also seeking for alternative sources for their cooling water. An opportunity currently under investigation is to pump a limited quantity of brackish water from the regional ground water. Although much cheaper than the water from the new SROG contract, the salinity of the brackish water is so high that a blend of brackish and SROG water would need additional treatment to improve its quality for use in the cooling towers. An on-site reverse osmosis (RO) desalination plant is envisaged to reduce the salinity to an acceptable level. This report investigates the economic impact of such an RO plant. As one can see, the RO plant would help with both problems APS is facing: First, it would rise the APS base load helping to mitigate VRE induced hub price volatility. From this perspective, the RO plant can be seen as a stabilizing load on the grid as would be any other industrial process. Second, the RO’s goods, e.g. the clean water produced can be used by APS itself for a (hopefully) lower cost than 100% SROG water would be. The report will first describe the APS water procurement strategy in detail and lay out the possible alternative scenarios studied in this document. This is covered in Chapter 2. Next, in Chapter 3, the APS case model developed within the N-R HES software framework is described. This chapter includes also a description of the needed framework developments to accommodate the APS study.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Guidelines for Publicly Archiving Terrestrial Model Data to Enhance Usability, Intercomparison, and Synthesis

Scientific communities are increasingly publishing data to evaluate, accredit, and build on published research. However, guidelines for curating data for publication are sparse for model-related research, limiting the usability of archived simulation data. In particular, there are no established guidelines for archiving data related to terrestrial models that simulate land processes and their coupled interactions with climate. Terrestrial modelers have a unique set of challenges when publishing data due to the diversity of scientific domains, research questions, and the types and scales of simulations. Researchers in the U.S. Department of Energy’s (DOE) projects use a variety of multiscale models to advance robust predictions of terrestrial and subsurface ecosystem processes. Here, we synthesize archiving needs for data associated with different DOE models, and provide guidelines for publishing terrestrial model data components following FAIR (Findable, Accessible, Interoperable, Reusable) principles. The guidelines recommend archiving model inputs and testing data used in final simulation runs along with associated codes, workflow scripts, and metadata in public repositories. Researchers should consider archiving model outputs if they are within the storage limits of the repository. We also provide considerations for how to bundle files into different data publications with citable digital object identifiers. Finally, we identify repository features and tools that would enable storage and reuse of model data. Given the diversity of DOE terrestrial models, these guidelines are transferable to other model types and will enable efficient reuse of simulation data for purposes such as model intercomparisons, initialization, benchmarking, synthesis, and comparisons with field observations.

58 GEOSCIENCES↗
DOE OSTI2018

Spectroscopic and Computational Analysis of Chemical Pollutants Interactions on Steel - 18409

There is limited understanding of the scientific interaction of the chemical contaminants to DOE facility materials which can alter the stability of the materials, affecting the general safety. The resulting studies will have a direct impact towards developing an enhanced cost-effective decontamination method and understanding the negative effects caused by the chemical waste on the structural integrity of the solid surface in contact. First order kinematic models provide a better fit to the experimental data, implying that one layer of Hg adsorbs to the steel surface. Energy dispersive X-ray analysis confirmed that stainless steel samples are mainly composed of Fe, Ni, and Cr, and have a low carbon content. Objective: To determine the physiochemical interactions of Mercury (Hg) and Cesium (Cs) to steel and concrete using spectroscopic, adsorption and computational studies. To develop a multiscale simulation that shows the effects of adsorption and diffusion of alkali metals on the mechanical performance of steel alloys used in nuclear power plants. Approach: Molecular Dynamics (MD) simulations were used to compute the stress gradient generated in steel due to Hg adsorption. Finite element models were developed for simulating the effects of temperature and pressure in a diffusion analysis. Based on experimental data, one layer of Hg was used for the MD simulations. Material properties at the atomic level were determined from the MD simulations The stress gradient from the MD simulations is used as an input parameter in the FE simulation. Conclusions and Future Work: Computational methods such as MD and FE simulations were used to determine the effects of mercury adsorption in stainless steel samples. The adsorption process induces a stress gradient in the facility materials, which can affect the mechanical performance of facility materials under service. Tensile stresses help the diffusion process while compressive have the opposite effect. Temperature has a positive effect in the mass diffusion process. Future work will evaluate the effects of temperature on the material properties at the molecular level using both experimental and MD methods. Future work will evaluate the mechanical performance of facility materials at a wide range of environmental variables.

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

Fast and self-learning indoor airflow simulation based on in situ adaptive tabulation

Fast simulation for stratified indoor airflow distributions is desired for various applications, such as design of advanced indoor environments, emergency management, and coupled annual energy simulation for buildings with stratified air distributions. Reduced order models trained by pre-computed computational fluid dynamics results are fast, but their prediction may be inaccurate when applied for conditions outside the training domain. To overcome this limitation, we propose in this paper a fast and self-learning model based on an in situ adaptive tabulation (ISAT) algorithm, which is trained by a fast fluid dynamics (FFD) model as an example. The idea is that the ISAT will retrieve the solutions from an existing data set if the estimated prediction error is within a pre-defined tolerance. Otherwise, the ISAT will execute the FFD simulation, which is accelerated by running in parallel on a graphics processing unit, for a full-scale simulation. This paper systematically investigates the feasibility of the ISAT for indoor airflow simulations by presenting the ISAT-FFD implementation alongside results related to its overall performance. Using a stratified indoor airflow as an example, we evaluated how the training time of ISAT was impacted by four factors (training methods, error tolerances, number of inputs, and number of outputs). Then we demonstrated that a trained ISAT model can predict the key information for inputs both inside and outside the training domain. The ISAT was able to answer query points both inside and close to training domain using retrieve actions within a time less than 0.001 s for each query. In conclusion, we provided suggestions for using the ISAT for building applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Position Papers for the ASCR Workshop on Cybersecurity and Privacy for Scientific Computing Ecosystems

At the request of the Department of Energy's (DOE) Office of Advanced Scientific Computing Research (ASCR), this program committee has been tasked with organizing a workshop to identify basic research needs in cybersecurity and privacy to better support DOE's science and energy mission. As part of the process, the program committee is soliciting community input in the form of position papers to help identify significant use cases, facility issues, and other barriers to enabling verifiably trustworthy computational science while preserving data confidentiality as appropriate for scientific workflows of interest to DOE. The program committee will review these position papers and based on the fit of their area of expertise and interest, selected contributors will have the opportunity to participate in the workshop currently planned as a virtual event November 3-5th, 2021. The thrust areas that will be explored by this workshop are the following: (1) Algorithms for secure, scalable, privacy-enhancing technologies and frameworks, including: Federated AI/ML, Differential privacy, Randomized algorithms, Adversarial modeling & simulation, Graph algorithms, and Formal methods; (2) Platforms to support the entire scientific-computing ecosystem, including edge computing for large-scale experiments, focusing on heterogeneous systems and distributed systems, including: Heterogeneous computing systems, Distributed computing systems, and Secure data architectures; and (3) Data workflows to allow agile use of data while preserving integrity and privacy, making the important properties verifiable either at runtime or post-computation, including: Integrity and provenance and Data management infrastructure. Topics that are out-of-scope for the workshop include discussing specific proposed solutions or areas that are clearly out of DOE's fundamental and applied-sciences mission scope, e.g., cryptography, enterprise security, and general-operations technology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Integration and Demonstration of Monitoring, Modeling, and Prediction of DV-1 Amendment Performance at the Bench Scale: DV-1 Amendment Demonstration

During fiscal years 2024 and 2025, the U.S. Department of Energy’s Hanford Field Office commissioned Pacific Northwest National Laboratory to conduct applied research aimed at reducing the cost, time, and uncertainty associated with in situ treatment of vadose zone contaminants at the Hanford Site. This report outlines the integration of three key research efforts into a meso-scale demonstration designed to advance field-scale solutions that aim to (1) optimize the delivery of chemical amendments to contaminated soils, (2) reduce uncertainty in amendment delivery performance assessment using advanced monitoring techniques, and (3) provide real-time insights into when and where amendment-induced precipitation reactions occur in the subsurface. To achieve these objectives, the tank-scale (~ 1 cubic meter) Geophysical Imaging of Flow and Transport (GIFT) system was developed. GIFT enables experimental testing of amendment delivery while incorporating automated multi-modal monitoring approaches, including pressure measurements, direct fluid sampling, and remote time-lapse geophysical imaging. The data generated from these monitoring techniques will serve as inputs for a generative artificial-intelligence-driven digital twin – a numerical simulation model designed to honor observed data while quantifying uncertainty in simulation accuracy. Using this simulator, researchers will refine an amendment injection strategy to maximize delivery efficiency within a low-permeability soil zone. Monitoring data will be interpreted through simulated outputs to enhance understanding of the injection process. The efficacy of this integrated approach will be evaluated through direct sampling at the conclusion of the experiment.

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

User Guidelines and Best Practices for CASL VUQ Analysis Using Dakota

Sandia’s Dakota software (available at http://dakota.sandia.gov) supports science and engineering transformation through advanced exploration of simulations. Specifically it manages and analyzes ensembles of simulations to provide broader and deeper perspective for analysts and decision makers. This enables them to enhance understanding of risk, improve products, and assess simulation credibility. In its simplest mode, Dakota can automate typical parameter variation studies through a generic interface to a physics-based computational model. This can lend efficiency and rigor to manual parameter perturbation studies already being conducted by analysts. However, Dakota also delivers advanced parametric analysis techniques enabling design exploration, optimization, model calibration, risk analysis, and quantification of margins and uncertainty with such models. It directly supports verification and validation activities. Dakota algorithms enrich complex science and engineering models, enabling an analyst to answer crucial questions of - Sensitivity: Which are the most important input factors or parameters entering the simulation, and how do they influence key outputs?; Uncertainty: What is the uncertainty or variability in simulation output, given uncertainties in input parameters? How safe, reliable, robust, or variable is my system? (Quantification of margins and uncertainty, QMU); Optimization: What parameter values yield the best performing design or operating condition, given constraints? Calibration: What models and/or parameters best match experimental data? In general, Dakota is the Consortium for Advanced Simulation of Light Water Reactors (CASL) delivery vehicle for verification, validation, and uncertainty quantification (VUQ) algorithms. It permits ready application of the VUQ methods described above to simulation codes by CASL researchers, code developers, and application engineers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

E-Area Vault Concrete Material Property And Vault Durability/Degradation Projection Recommendations

Subsequent to the 2008 E-Area Low-Level Waste Facility (ELLWF) Performance Assessment (PA) (WSRC 2008), two additional E-Area vault concrete property testing programs have been conducted (Dixon and Phifer 2010 and SIMCO 2011a) and two additional E-Area vault concrete durability modeling projections have been made (Langton 2009 and SIMCO 2012). All the information/data from these reports has been evaluated and consolidated herein by the Savannah River National Laboratory (SRNL) at the request of Solid Waste Management (SWM) to produce E-Area vault concrete hydraulic and physical property data and vault durability/degradation projection recommendations that are adequately justified for use within associated Special Analyses (SAs) and future PA updates. The Low Activity Waste (LAW) and Intermediate Level (IL) Vaults structural degradation predictions produced by Carey 2006 and Peregoy 2006, respectively, which were used as the basis for the 2008 ELLWF PA, remain valid based upon the results of the E-Area vault concrete durability simulations reported by Langton 2009 and those reported by SIMCO 2012. Therefore revised structural degradation predictions are not required so long as the mean thickness of the closure cap overlying the vaults is no greater than that assumed within Carey 2006 and Peregoy 2006. For the LAW Vault structural degradation prediction (Carey 2006), the mean thickness of the overlying closure cap was taken as nine feet. For the IL Vault structural degradation prediction (Peregoy 2006), the mean thickness of the overlying closure cap was taken as eight feet. The mean closure cap thicknesses as described here for both E-Area Vaults will be included as a key input and assumption (I&A) in the next revision to the closure plan for the ELLWF (Phifer et al. 2009). In addition, it has been identified as new input to the PA model to be assessed in the ongoing update to the new PA Information UDQE (Flach 2013). Once the UDQE is approved, the SWM Key I&A database will be updated with this new information.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Biospheric feedback effects in a synchronously coupled model of human and Earth systems

Fossil fuel combustion and land-use change are the first and second largest contributors to industrial-era increases in atmospheric carbon dioxide concentration, which is itself the largest driver of present-day climate change1. Projections of fossil fuel consumption and land-use change are thus fundamental inputs for coupled Earth system models (ESM) used to estimate the physical and biological consequences of future climate system forcing2,3. While empirical datasets are available to inform historical analyses4,5, assessments of future climate change have relied on projections of energy and land use based on energy economic models, constrained using historical and present-day data and forced with assumptions about future policy, land-use patterns, and socio-economic development trajectories6. Here we show that the influence of biospheric change – the integrated effect of climatic, ecological, and geochemical processes – on land ecosystems has a significant impact on energy, agriculture, and land-use projections for the 21st century. Such feedbacks have been ignored in previous ESM studies of future climate. We find that synchronous exposure of land ecosystem productivity in the economic system to biospheric change as it develops in an ESM results in a 10% reduction of land area used for crop cultivation; increased managed forest area and land carbon; a 15-20% decrease in global crop price; and a 17% reduction in fossil fuel emissions for a low-mid range forcing scenario7. These simulation results demonstrate that biospheric change can significantly alter primary human system forcings to the climate system. This synchronous two-way coupling approach removes inconsistencies in description of climate change between human and biosphere components of the coupled model, mitigating a major source of uncertainty identified in assessments of future climate projections8-10.

climate change↗
DOE OSTI2024

Improved Fuel Cycle Capability of Griffin for Fast Reactor Applications

Griffin is a MOOSE based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, the fuel cycle capability has been significantly extended by improving the assembly shuffling option to allow flexible fuel reloading in the multi-cycle depletion calculation and incorporating decay between cycles. An equilibrium core calculation capability was also implemented to find an equilibrium core. Additionally, an enrichment search capability was added to determine the enrichment condition that allows a core to reach an equilibrium cycle with the end-of-cycle k-effective meeting a user-specified target value. The updated fuel management capability has been extensively tested using the three-dimensional ABTR problem with different batch schemes, exhibiting reasonable solutions in terms of manual shuffling, equilibrium cycle, and enrichment search calculations. The cross section generation workflow capability for fast reactors was further verified to produce microscopic cross sections as well as Griffin core inputs for an ring-heterogeneous configuration. Rigorous verification tests using the ABTR problem demonstrated that the RH core calculations, with cross sections and Griffin inputs generated from the cross section workflow of Griffin, produced accurate solution for fast reactor problems. In addition, an option to convert delay neutron parameter data generated from MC 2 -3 in the DLAYXS format into XML format was added to support transient calculations using MC 2 -3-generated data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

HDF5 in the exascale era: Delivering efficient and scalable parallel I/O for exascale applications

Accurately modeling real-world systems requires scientific applications at exascale to generate massive amounts of data and manage data storage efficiently. However, parallel input and output (I/O) faces challenges due to new application workflows and the state-of-the-art memory, interconnect, and storage architectures considered in exascale designs. The storage hierarchy has expanded with node-local persistent memory, solid-state storage, and traditional disk and tape-based storage, thus requiring efficiency at each layer and much more efficient data movement among these layers. This paper discusses how the ExaHDF5 project improved the I/O performance and data management for exascale architectures by enhancing HDF5, a widely used parallel I/O library. The team developed an Asynchronous I/O Virtual Object Layer (VOL) connector that allowed overlapping I/O with computation. They also created a Cache VOL to complement asynchronous I/O by incorporating fast storage layers, such as burst buffer and node-local storage, into the parallel I/O workflow through caching and staging data. Additionally, the team enabled data aggregation and I/O at the node level by using a Subfiling Virtual File Driver (VFD). To demonstrate superior I/O performance with HDF5 at exascale, the ExaHDF5 team collaborated with several exascale applications. In this paper, we show I/O performance improvements for three applications: Cabana (a particle-based simulation library), EQSIM (a regional earthquake simulation software), and E3SM (a climate system modeling library).

Asynchronous I/Ol↗
DOE OSTI2017

EpiK: A Knowledge Base for Epidemiological Modeling and Analytics of Infectious Diseases

Computational epidemiology seeks to develop computational methods to study the distribution and determinants of health-related states or events (including disease), and the application of this study to the control of diseases and other health problems. Recent advances in computing and data sciences have led to the development of innovative modeling environments to support this important goal. The datasets used to drive the dynamic models as well as the data produced by these models presents unique challenges owing to their size, heterogeneity and diversity. These datasets form the basis of effective and easy to use decision support and analytical environments. As a result, it is important to develop scalable data management systems to store, manage and integrate these datasets. In this paper, we develop EpiK —a knowledge base that facilitates the development of decision support and analytical environments to support epidemic science. An important goal is to develop a framework that links the input as well as output datasets to facilitate effective spatio-temporal and social reasoning that is critical in planning and intervention analysis before and during an epidemic. The data management framework links modeling workflow data and its metadata using a controlled vocabulary. The metadata captures information about storage, the mapping between the linked model and the physical layout, and relationships to support services. EpiK is designed to support agent-based modeling and analytics frameworks—aggregate models can be seen as special cases and are thus supported. We use semantic web technologies to create a representation of the datasets that encapsulates both the location and the schema heterogeneity. The choice of RDF as a representation language is motivated by the diversity and growth of the datasets that need to be integrated. A query bank is developed—the queries capture a broad range of questions that can be posed and answered during a typical case study pertaining to disease outbreaks. The queries are constructed using SPARQL Protocol and RDF Query Language (SPARQL) over the EpiK. EpiK can hide schema and location heterogeneity while efficiently supporting queries that span the computational epidemiology modeling pipeline: from model construction to simulation output. As a result, we show that the performance of benchmark queries varies significantly with respect to the choice of hardware underlying the database and resource description framework (RDF) engine.

60 APPLIED LIFE SCIENCES↗
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↗