Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation framework”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Status Report on the Component Models Developed in the Modelica Framework: Reverse Osmosis Desalination Plant & Thermal Energy Storage

This report has been prepared as part of an effort to design and build a Modeling and Simulation (M&S) framework to assess the economic viability of a Nuclear-Renewable Hybrid Energy System (N-R HES). In order to facilitate dynamic M&S of such an integrated system, research groups in multiple national laboratories and universities have been developing various subsystems as dynamic physics-based components using the Modelica programming language. In Fiscal Years (FYs) 2015, Idaho National Laboratory (INL) performed a dynamic analysis of two region-specific N-R HES configurations, including the gas-to-liquid (natural gas to Fischer-Tropsch synthetic fuel) and brackish water Reverse Osmosis (RO) desalination plants as industrial processes. In FYs 2016–2017, INL developed two additional subsystems in the Modelica framework: (1) a high-temperature steam electrolysis plant as a high priority industrial plant to be integrated with a light water reactor within an N-R HES and (2) a gas turbine power plant as a secondary energy supply. In FY 2018, the RO desalination system model developed in FY 2015 has been updated such that the model is compatible with the most recent version of the ThermoPower library. Special attention has been given to the controller settings based on process models, aiming to improve process dynamics and controllability. A dynamic performance analysis of the updated RO desalination plant was carried out to evaluate the technical feasibility (load-following capability) of such a system operating under highly variable conditions requiring flexible output. Simulation results involving several case studies show that the suggested control scheme could maintain the controlled variables (including the variable electrical load and RO feed pressure) within desired limits under various plant operating conditions. The results also indicate that the proposed RO plant could provide operational flexibility to participate in energy management at the utility scale by dynamically optimizing the use of excess plant capacity within an N-R HES. For a small-scale energy storage system, a sensible Thermal Energy Storage (TES) model has been developed in the Modelica Framework in FY 2018.

25 ENERGY STORAGE↗
DOE OSTI2025

Learning to Simulate Aerosol Dynamics with Graph Neural Networks

Aerosol effects on climate, weather, and air quality depend on characteristics of individual particles, which are tremendously diverse and change in time. Particle-resolved models are the only models able to capture this diversity in particle physiochemical properties, and these models are computationally expensive. As a strategy for accelerating particle-resolved microphysics models, we introduce Graph-based Learning of Aerosol Dynamics (GLAD) and use this model to train a surrogate of the particle-resolved model PartMC-MOSAIC. GLAD implements a Graph Network-based Simulator (GNS), a machine learning framework that has been used to simulate particle-based fluid dynamics models. In GLAD, each particle is represented as a node in a graph, and the evolution of the particle population over time is simulated through learned message passing. Here, we demonstrate our GNS approach on a simple aerosol system that includes condensation of sulfuric acid onto particles composed of sulfate, black carbon, organic carbon, and water. A graph with particles as nodes is constructed, and a graph neural network (GNN) is then trained using the model output from PartMC-MOSAIC. The trained GNN can then be used for simulating and predicting aerosol dynamics over time. Results demonstrate the framework's ability to accurately learn chemical dynamics and generalize across different scenarios, achieving efficient training and prediction times. We evaluate the performance across four scenarios, highlighting the framework's robustness and adaptability in modeling aerosol microphysics and chemistry.

aerosol chemistry dynamics↗
DOE OSTI2021

Signatures of electronic and nuclear coherences in ultrafast molecular x-ray and electron diffraction

Femtosecond x-ray and electron diffraction hold promise to image the evolving structures of single molecules. We present a unified quantum-electrodynamical formulation of diffraction signals, based on the exact many-body nuclear + electronic wavefunction that can be extracted from quantum chemistry simulations. This gives a framework for analyzing various approximate molecular dynamics simulations. We show that the complete description of ultrafast diffraction signals contains interesting contributions involving mixed elastic and inelastic scattered photons that are usually masked by other larger contributions and are neglected. These terms include overlaps of nuclear wavepackets between different electronic states that provide an electronic decoherence mechanism and are important for the time-resolved imaging of conical intersections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

SU-F-T-507: Modeling Cerenkov Emissions From Medical Linear Accelerators: A Monte Carlo Study

Purpose: Cerenkov emissions are a natural byproduct of MV radiotherapy but are typically ignored as inconsequential. However, Cerenkov photons may be useful for activation of drugs such as psoralen. Here, we investigate Cerenkov radiation from common radiotherapy beams using Monte Carlo simulations. Methods: GAMOS, a GEANT4-based framework for Monte Carlo simulations, was used to model 6 and 18MV photon beams from a Varian medical linac. Simulations were run to track Cerenkov production from these beams when irradiating a 50cm radius sphere of water. Electron contamination was neglected. 2 million primary photon histories were run for each energy, and values scored included integral dose and total track length of Cerenkov photons between 100 and 400 nm wavelength. By lowering process energy thresholds, simulations included low energy Bremsstrahlung photons to ensure comprehensive evaluation of UV production in the medium. Results: For the same number of primary photons, UV Cerenkov production for 18MV was greater than 6MV by a factor of 3.72 as determined by total track length. The total integral dose was a factor of 2.31 greater for the 18MV beam. Bremsstrahlung photons were a negligibly small component of photons in the wavelength range of interest, comprising 0.02% of such photons. Conclusion: Cerenkov emissions in water are 1.6x greater for 18MV than 6MV for the same integral dose. Future work will expand the analysis to include optical properties of tissues, and to investigate strategies to maximize Cerenkov emission per unit dose for MV radiotherapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Numerical simulations of yield-based sooting tendencies of aromatic fuels using ReaxFF molecular dynamics

In this work, we present the first ReaxFF Molecular Dynamics (MD) simulations to quantitatively predict the sooting tendencies of various fuels. The specific sooting tendency metric used in this work is the Yield Sooting Index (YSI), which quantifies the effects of fuel molecular structure on soot yield. YSI has been experimentally measured and numerically simulated using computational fluid dynamics for a large range of fuels, but there is no existing reactive MD framework for YSI simulations. To adopt the experimental YSI concept, a multi-stage simulation procedure is designed using ReaxFF. As a proof-of-concept, toluene and phenol are selected as test fuels, since both have relatively well-understood reaction pathways. The ReaxFF YSI simulations are shown to capture key reaction events for both fuels selected that are consistent with existing chemical kinetic understanding. Toluene is shown to mostly retain its original aromatic ring structure and directly grow to larger aromatic compounds with multiple rings. On the other hand, the aromatic growth process from phenol is accompanied by carbon-loss reactions with CO release. In addition, a quantitative YSI formulation is also derived in ReaxFF and the ReaxFF-predicted YSI values are compared with measurement data and reasonably good agreement is achieved. The results reported in this work demonstrates that the ReaxFF-based framework can potentially be used to quantitatively predict the relative sooting tendencies, especially for fuels with unknown or poorly-known chemistry, to understand their sooting properties and in search for soot relevant reaction pathways from these fuels.

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

A general kinetic framework for dislocation mobility derived from probabilistic cellular automaton simulations of the kink-pair mechanism

Dislocation mobility laws are essential components of dislocation-density-based crystal plasticity models. For dislocations governed by the kink-pair mechanism, however, existing formulations are typically limited to specific regimes due to the com plex interplay between stochastic kink-pair nucleation and lateral kink migration. In this work, we develop a general kinetic framework that expresses the average dislocation velocity as a function of mechanism-level variables: positive/negative kink pair nucleation rates, kink migration velocity, dislocation segment length, critical kink-pair width, and kink height. Probabilis tic cellular automaton simulations are used to capture the behavior of conceptual dislocation segments between the limiting conditions of migration outpacing nucleation on the one end and nucleation outpacing migration on the other. An elemen tary functional form that captures the system dynamics is suggested and fitted against the simulation results. This framework remains valid for arbitrary combinations of the six variables and is, therefore, compatible with any admissible constitutive re lations that describe their stress and temperature dependence. Comparisons with established approaches and experimental results confirm the robustness and physical consistency of the formulation, making it broadly applicable to material systems in which dislocation motion is governed by the kink-pair mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hybrid-RL-MPC4CLR (Hybird-Reinforcement-Learning-Model-Predictive-Control-for-Reserve-Policy-Assisted-Critical-Load-Restoration-in-Distribution-Grids)

Hybrid-RL-MPC4CLR was developed as a hybrid controller for active distribution grid critical load restoration, combining deep reinforcement learning (RL) and model predictive control (MPC) aiming at maximizing total restored load following an extreme event. The RL determines a policy for quantifying operating reserve requirements, thereby hedging against uncertainty, while the MPC models grid operations incorporating the RL policy actions (i.e., reserve requirements), renewable (wind and solar) power predictions, and load demand forecasts. The developers formulated the reserve requirement determination problem as a sequential decision-making problem based on the Markov Decision Process (MDP) and design an RL learning environment based on the OpenAI Gym framework and MPC simulation. The RL agent reward and MPC objective function aim to maximize and monotonically increase total restored load and minimize load shedding and renewable power curtailment. The software is developed using various software packages in Python. The MPC's optimal power flow (OPF) model is implemented using the Pyomo package, the RL simulation environment is implemented using the MPC simulation with various scenarios of renewable energy and load demand profiles and power outage beginning times, based on the OpenAI Gym framework. The RL agent training is performed using the RLlib Ray package. The RL algorithm is trained offline using historical forecasts of renewable generation and load demand profiles. Simulation analysis and performance tests are conducted using a modified IEEE 13-bus distribution test feeder containing wind turbine, photovoltaic, microturbine, and battery.

Eseye, Abinet Tesfaye↗
DOE OSTI2024

A parallel and performance portable implementation of a full-field crystal plasticity model

We have developed a parallel implementation of an Elasto-Viscoplastic Fast Fourier Transform-based (EVPFFT) micromechanical solver to enable computationally efficient crystal plasticity modeling for polycrystalline materials. Our primary focus lies in achieving performance portability, allowing a single EVPFFT implementation to run optimally on various homogeneous architectures, including multi-core Central Processing Units (CPUs), as well as on heterogeneous computer architectures comprising multi-core CPUs and Graphics Processing Units (GPUs) from different vendors. To accomplish this goal, we have leveraged MATAR, a C++ software library that simplifies the creation and utilization of multidimensional dense or sparse matrix and array data structures. These data structures are designed to be portable across diverse architectures through the use of Kokkos, a performance-portable library. Additionally, we have employed the Message Passing Interface (MPI) to efficiently distribute the computational workload among processors. The heFFTe (Highly Efficient FFT for Exascale) library is used to facilitate the performance portability of the fast Fourier transforms (FFTs) computation. The computational performance of EVPFFT is evaluated and presented in terms of parallel scalability and simulation runtime on different high-performance computing (HPC) architectures. As a result, the utility of the developed framework to efficiently simulate the micro-mechanical fields in polycrystalline microstructures in engineering applications is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of hydrogen on plasticity of $α$-Fe: A multi-scale assessment

A multi-scale study was carried out to quantify the effect of interstitial hydrogen concentration on plasticity in $α$-Fe. Here, in this work, the influence of hydrogen on the screw dislocation glide behavior was examined across several length-scales. The insights obtained were integrated to provide an accurate continuum description for the effect of hydrogen on the dislocation based plasticity in polycrystalline $α$-Fe. At the outset of this work, a new Fe—H interatomic potential was formulated that enhanced the atomistic estimation of the variation in dislocation glide behavior in presence of hydrogen. Next, the dislocation core reconstruction observed due to the addition of hydrogen using atomistic simulations was validated with the help of large-scale DFT calculations based on the DFT-FE framework. Several atomistic simulations were carried out to comprehensively quantify the effect of hydrogen on the non-Schmid behavior exhibited during the dislocation glide in $α$-Fe. Finally, crystal plasticity simulations were carried out to understand the effect of hydrogen on the meso-scale deformation behavior of polycrystalline $α$-Fe.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Data-centric framework for crystal structure identification in atomistic simulations using machine learning

Atomic-level modeling performed at large scales enables the investigation of mesoscale materials properties with atom-by-atom resolution. The spatial complexity of such cross-scale simulations renders them unsuitable for simple human visual inspection. Instead, specialized structure characterization techniques are required to aid interpretation. These have historically been challenging to construct, requiring significant intuition and effort. Here we propose an alternative framework for a fundamental structural characterization task: classifying atoms according to the crystal structure to which they belong. Our approach is data-centric and favors the employment of Machine Learning over heuristic rules of classification. A group of data-science tools and simple local descriptors of atomic structure are employed together with an efficient synthetic training set. We also introduce the first standard and publicly available benchmark data set for evaluation of algorithms for crystal-structure classification. Further, it is demonstrated that our data-centric framework outperforms all of the most popular heuristic methods—especially at high temperatures when lattices are the most distorted—while introducing a systematic route for generalization to new crystal structures. Moreover, through the use of outlier detection algorithms our approach is capable of discerning between amorphous atomic motifs (i.e., noncrystalline phases) and unknown crystal structures, making it uniquely suited for exploratory materials synthesis simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modelling Beam Loss Within the LANSCE Proton Storage Ring

Several upgrades are being considered for the proton storage ring (PSR) at the Los Alamos Neutron Science Center (LANSCE) to reduce beam loss and thereby reduce the cooldown period of the PSR following a beam run. First, we have considered an increased beam pipe diameter would reduce beam loss due to beam scraping caused by misalignments and tuning errors. However, this would require increased pole-to-pole gap height within the dipoles and quadrupoles, which would change their effective length and alter their fringe fields. The effect of different magnet gaps on the beam optical parameters and on beam loss was studied using the simulation codes MAD-X and PyORBIT. Second, we are developing a detailed particle tracking model within the framework of the simulation package General Particle Tracer (GPT) for the PSR H - stripping system. This model will be used to study the effect of the stripper foil parameters (position, composition, areal density, depth, etc.) on first-turn losses, where most of the observed beam loss and emittance growth occurs due to foil scattering, foil stripping, and Lorentz stripping. The model implements C++ custom elements created to model foil scattering, foil stripping, and Lorentz stripping, as these are not built-in features of GPT. The preliminary results of the MAD-X and PyORBIT simulations, as well as the GPT simulation model, is described in this tech note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Using geometric simulation software ‘GASP’ to model conformational flexibility in a family of zinc metal–organic frameworks

Here, a new tripodal tricarboxylic acid ligand, 4,4'-(4'-(4'-carboxy-[1,1'-biphenyl]-4-yl)-[2,2':6',2''-terpyridine]-5,5''-diyl)-dibenzoic acid (H 3 cbt ), was synthesised using a three-step convergent strategy. Subsequent reactions with zinc( II ) nitrate hexahydrate yielded three metal–organic frameworks (MOFs). The three MOFs, [Zn(H cbt )]·4DMF ( 1 ), [Zn(H cbt )]·4DMSO·1.5H 2 O·DMF ( 2 ), and [Zn(H cbt )]·2DMF·3H 2 O ( 3 ), each adopt flexible interdigitated 2D net topologies. Framework 1 has DMF-filled channels that retain porosity upon desolvation, with a measured BET surface area of 248 m 2 g -1 . Framework 2 possesses larger DMSO-containing channels that collapse upon desolvation, resulting in near-equivalent porosity values to framework 1 . In silico calculations and topological considerations determined using the geometric simulation software GASP dictate that framework 2 can feasibly alter conformation to approximate 1 , but cannot perfectly replicate the interdigitated motif. Framework 3 formed when wet solvents were used to synthesise 1 . Interestingly, the interdigitated structure of 3 contains a unique carboxylate binding mode that precludes its subsequent adoption by either 1 or 2 upon their exposure to water. This diverse array of structural considerations recommends this MOF family for modelling using GASP. Interrogating frameworks 1–3 using this software provided insights that justified experimentally observed conformational trends, as well as barriers to interconversion between members of this MOF family. In a broader sense, this work demonstrates the wider applicability of GASP software to modelling structural changes within flexible MOF materials.

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

Status on the Component Models Developed in the Modelica Framework: High-Temperature Steam Electrolysis Plant & Gas Turbine Power Plant

This report has been prepared as part of an effort to design and build a modeling and simulation (M&S) framework to assess the economic viability of a nuclear-renewable hybrid energy system (N-R HES). In order to facilitate dynamic M&S of such an integrated system, research groups in multiple national laboratories have been developing various subsystems as dynamic physics-based components using the Modelica programming language. In fiscal year (FY) 2015, Idaho National Laboratory (INL) performed a dynamic analysis of two region-specific N-R HES configurations, including the gas-to-liquid (natural gas to Fischer-Tropsch synthetic fuel) and brackish water reverse osmosis desalination plants as industrial processes. In FY 2016, INL has developed two additional subsystems in the Modelica framework: a high-temperature steam electrolysis (HTSE) plant and a gas turbine power plant (GTPP). HTSE has been proposed as a high priority industrial process to be integrated with a light water reactor (LWR) in an N-R HES. This integrated energy system would be capable of dynamically apportioning thermal and electrical energy (1) to provide responsive generation to the power grid and (2) to produce alternative industrial products (i.e., hydrogen and oxygen) without generating any greenhouse gases. A dynamic performance analysis of the LWR/HTSE integration case was carried out to evaluate the technical feasibility (load-following capability) and safety of such a system operating under highly variable conditions requiring flexible output. To support the dynamic analysis, the detailed dynamic model and control design of the HTSE process, which employs solid oxide electrolysis cells, have been developed to predict the process behavior over a large range of operating conditions. As first-generation N-R HES technology will be based on LWRs, which provide thermal energy at a relatively low temperature, complementary temperature-boosting technology was suggested for integration with the HTSE process that requires higher temperature input. Simulation results involving several case studies show that the suggested control scheme could maintain the controlled variables (including the steam utilization factor, cathode stream inlet composition, and temperatures of the process streams at various locations) within desired limits under various plant operating conditions. The results also indicate that the proposed HTSE plant could provide operational flexibility to participate in energy management at the utility scale by dynamically optimizing the use of excess plant capacity within an N-R HES. A natural-gas fired GTPP has been proposed as a secondary energy supply to be included in an N-R HES. This auxiliary generator could be used to cover rapid dynamics in grid demand that cannot be met by the remainder of the N-R HES. To evaluate the operability and controllability of the proposed process during transients between load (demand) levels, the dynamic model and control design were developed. Special attention was given to the design of feedback controllers to regulate the power frequency, and exhaust gas and turbine firing temperatures. Several case studies were performed to investigate the system responses to the major disturbance (power load demand) in such a control system. The simulation results show that the performance of the proposed control strategies was satisfactory under each test when the GTPP experienced high rapid variations in the load.

08 HYDROGEN↗
DOE OSTI2016

SU-F-J-14: Kilovoltage Cone-Beam CT Dose Estimation of Varian On-Board Imager Using GMctdospp Monte Carlo Framework

Purpose: Although cone-beam CT (CBCT) imaging became popular in radiation oncology, its imaging dose estimation is still challenging. The goal of this study is to assess the kilovoltage CBCT doses using GMctdospp - an EGSnrc based Monte Carlo (MC) framework. Methods: Two Varian OBI x-ray tube models were implemented in the GMctpdospp framework of EGSnrc MC System. The x-ray spectrum of 125 kVp CBCT beam was acquired from an EGSnrc/BEAMnrc simulation and validated with IPEM report 78. Then, the spectrum was utilized as an input spectrum in GMctdospp dose calculations. Both full and half bowtie pre-filters of the OBI system were created by using egs-prism module. The x-ray tube MC models were verified by comparing calculated dosimetric profiles (lateral and depth) to ion chamber measurements for a static x-ray beam irradiation to a cuboid water phantom. An abdominal CBCT imaging doses was simulated in GMctdospp framework using a 5-year-old anthropomorphic phantom. The organ doses and effective dose (ED) from the framework were assessed and compared to the MOSFET measurements and convolution/superposition dose calculations. Results: The lateral and depth dose profiles in the water cuboid phantom were well matched within 6% except a few areas - left shoulder of the half bowtie lateral profile and surface of water phantom. The organ doses and ED from the MC framework were found to be closer to MOSFET measurements and CS calculations within 2 cGy and 5 mSv respectively. Conclusion: This study implemented and validated the Varian OBI x-ray tube models in the GMctdospp MC framework using a cuboid water phantom and CBCT imaging doses were also evaluated in a 5-year-old anthropomorphic phantom. In future study, various CBCT imaging protocols will be implemented and validated and consequently patient CT images will be used to estimate the CBCT imaging doses in patients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Gene expression of functionally-related genes coevolves across fungal species: detecting coevolution of gene expression using phylogenetic comparative methods

Researchers often measure changes in gene expression across conditions to better understand the shared functional roles and regulatory mechanisms of different genes. Analogous to this is comparing gene expression across species, which can improve our understanding of the evolutionary processes shaping the evolution of both individual genes and functional pathways. One area of interest is determining genes showing signals of coevolution, which can also indicate potential functional similarity, analogous to co-expression analysis often performed across conditions for a single species. However, as with any trait, comparing gene expression across species can be confounded by the non-independence of species due to shared ancestry, making standard hypothesis testing inappropriate. We compared RNA-Seq data across 18 fungal species using a multivariate Brownian Motion phylogenetic comparative method (PCM), which allowed us to quantify coevolution between protein pairs while directly accounting for the shared ancestry of the species. Our work indicates proteins which physically-interact show stronger signals of coevolution than randomly-generated pairs. Interactions with stronger empirical and computational evidence also showing stronger signals of coevolution. We examined the effects of number of protein interactions and gene expression levels on coevolution, finding both factors are overall poor predictors of the strength of coevolution between a protein pair. Simulations further demonstrate the potential issues of analyzing gene expression coevolution without accounting for shared ancestry in a standard hypothesis testing framework. Furthermore, our simulations indicate the use of a randomly-generated null distribution as a means of determining statistical significance for detecting coevolving genes with phylogenetically-uncorrected correlations, as has previously been done, is less accurate than PCMs, although is a significant improvement over standard hypothesis testing. These methods are further improved by using a phylogenetically-corrected correlation metric. Our work highlights potential benefits of using PCMs to detect gene expression coevolution from high-throughput omics scale data. This framework can be built upon to investigate other evolutionary hypotheses, such as changes in transcription regulatory mechanisms across species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Estimation of Fission Product Transport Parameters for Cesium in the AGR-3/4 TRISO Fuel Experiment

A one-dimensional (1D) finite-element model of fission product transport in the AGR-3/4 experiment has been developed using the Multiphysics Object Oriented Simulation Environment (MOOSE) framework and implemented in the fuel performance code, BISON. The model resolves capsule-specific geometries, materials, and temperature histories and simulates radial migration of fission products from the fuel compact through the inner ring, outer ring, and into the sink ring. Model parameters governing diffusion and sorption were estimated for key fission products – cesium (Cs), and europium (Eu) – by simultaneously fitting modeled isotopic concentration profiles and total ring inventories to a post-irradiation experimental measurement. These data include gamma scanning, liquid scintillation for Sr-90, radial deconsolidation leach-burn-leach analysis, tomographic reconstructions, and destructive physical sampling. A mortar-based interfacial sorption framework was implemented to enforce physically consistent mass transfer and flux conservation across gas gaps. Two classes of parameter sets were derived: a least-squares best-fit, and a safety-oriented conservative-fit, what applies strong penalties for underprediction of sink inventories. Across all twelve capsules, the model successfully reproduces the dominant radial transport trends for Cs, Sr, with decreasing concentrations from the compact outward through successive rings. Cs behavior is captured most consistently, while strontium predictions reveal systematic trade-offs between compact accuracy and conservative sink-ring bounding. The results demonstrate that sink ring weighted calibration provides conservative, safety relevant bounds on low temperature fission product transport, but at the cost of underpredicting compact inventories for Sr isotopes. These discrepancies highlight the need for additional physics, including fast-slow diffusion model, incorporating trapping mechanism in the transport behavior. Overall, this work establishes a robust, capsule-specific modeling framework for AGR-3/4 fission product transport and provides a defensible basis for parameter selection in source-term and fuel performance analyses for high temperature gas-cooled reactors.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Underground sources of drinking water chemistry changes in response to potential CO 2 leakage

The purpose of this study was to quantify changes to underground sources of drinking water (USDW) quality in response to potential CO 2 leakage from geologic CO 2 sequestration (GCS) reservoirs. Here, the developed a framework of combined laboratory experiments and reactive transport simulations and used this framework to evaluate the Ogallala aquifer overlying the Farnsworth Unit (FWU), an active GCS site, as a case study. Using chemical reaction parameters obtained from laboratory experiments and numerical simulations, site-specific mechanisms of CO 2 -water-sediment interactions at the USDW aquifer were interpreted. Long-term risks of potential CO 2 leakage were then evaluated with field-scale numerical models using the regional hydrogeological characteristics and reaction parameters obtained from our experiments and simulations. Results suggest that carbonate mineral impurity and cation exchange are key mechanisms for interactions between CO 2 and the aquifer sediment. Additionally, for a large leakage rate of 0.1 % injection from one leaky well, the leakage plume might impact an area of 300 m in diameter and significantly affect the local water quality by changing pH and cation concentrations (e.g., Zn, Ba and Sr). After leakage ceases, the zone of impacted fluids would not migrate significantly in subsequent decades due to a low regional groundwater flowrate (for this case study). The relatively small area of impact might not be detected in a monitoring well given the broader spacing in a typical field scenario. Effective early leakage detection may require additional tools, e.g., borehole CO 2 movement, four-dimensional seismicity, CO 2 soil flux, samples from deeper aquifers, etc., to ensure effective leakage detection and long-term safety of GCS projects.

58 GEOSCIENCES↗
DOE OSTI2020

EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers, part II: Regional simulations of building response

The existing observational database of the regional-scale distribution of strong ground motions and measured building response for major earthquakes continues to be quite sparse. As a result, details of the regional variability and spatial distribution of ground motions, and the corresponding distribution of risk to buildings and other infrastructure, are not comprehensively understood. Utilizing high-performance computing platforms, emerging high-resolution, physics-based ground motion simulations can now resolve frequencies of engineering interest and provide detailed synthetic ground motions at high spatial density. This provides an opportunity for new insight into the distribution of infrastructure seismic demands and risk. In the work presented herein, the EQSIM fault-to-structure computational framework described in a companion paper, McCallen et al., is employed to investigate the regional-scale response of buildings to large earthquakes. A representative M = 7.0 strike-slip event is used to explore the distribution and amplitude of building demand, and comparisons are made between building response computed with fault-to-structure simulations and building response computed with existing measured near-fault earthquake records. New information on the distribution and variability of building response from high-performance parallel simulations is described and analyzed, and favorable first comparisons between building response predicted with both fault-to-structure simulations and real ground motions records are presented.

58 GEOSCIENCES↗