Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

High-Enthalpy Testing to Validate Simulation of an Aerosol Capture Probe

AERACEPT (Aerosol Rapid Analysis Combined Entry Probe/sonde Technology) is a developing technology enabling in-situ aerosol particle sampling and analysis in a small spacecraft mission envelope. It integrates a passive aerosol sample collection system into a probe’s thermal protection system (TPS) to remove the need for heat shield separation and active descent control (parachutes, gliders, etc.). AERACEPT is being validated against the requirements of the Nephele mission concept, which targets the middle and lower Venus cloud layers. The proposed technology employs 3D Carbon-Carbon (3D-CC) at the probe’s nose to withstand the extreme environments encountered in a Venus entry without producing pyrolysis gases that would contaminate the mission sample collection. Any shape change of this 3D-CC sample inlet will affect airflow through the sampling system, thereby impacting aerosol sample collection bias and efficiency. Accurate prediction of the nose’s material response during entry is there-fore of fundamental importance to the development and practical use of this technology. As a result, the AERACEPT project is planning a test campaign in the newly built PlasmatronX facility with a novel “Open Iso-Q” test article to validate state-of-the-art material response tools for a simulated Venus entry. The PlasmatronX is a 350 kW inductively-coupled plasma facility developed and run by the Center for Hypersonics and Entry Systems Studies at the University of Illinois at Urbana-Champaign. The facility can support ground testing for a variety of planetary destinations by simulating entry conditions in Nitrogen, Air, and Carbon Dioxide. To best simulate the Venus atmosphere and match flight-like recession, the planned AERACEPT test campaign will use Carbon Dioxide as its test gas. AERACEPT has designed two test articles for this campaign: 1) a standard test article whose curved surface approximates a constant applied heat flux (“Closed Iso-Q” model) and 2) a novel inlet test article with a through-hole at the stagnation point (“Open Iso-Q” model). Both article types consist of a 3D-CC sample bonded to a graphite fixture. All articles will be laser scanned before and after testing to estimate shape change and recession, and the “Closed Iso-Q” articles will be instrumented with thermocouples to provide temperature histories. The test article temperatures and inlet shape change will be compared against simulations run with the Porous material Analysis Toolbox based on OpenFOAM (PATO) to validate the AERACEPT 3D-CC material response model. Due to the high thermal conductivity of 3D-CC, multidimensional effects have a large influence on the test article temperature and the standard 1-D material response tools are not sufficient for test planning. The test campaign features (e.g. article design, run duration, heat fluxes) are therefore informed by a series of 3-D simulations using state-of-the-art material response tools such as PATO. The test campaign is scheduled for spring of 2024. Expected results include 1) Key drivers of selected test conditions (e.g. matching mission flight atmosphere & recession), 2) Simulation results that informed test article design, 3) Acquired test data such as recession measurements, temperature histories, and photos, and 4) Comparisons between test results and simulation predictions.

AERACEPT↗
NASA NTRSDate not supplied

Simulated Multipath Using Software Generated GPS Signals

Depending on the environment, multipath can be one of the largest error sources contributing to degradation in Global Navigation Satellite System (GNSS) (e.g., GPS) performance. Multipath is a phenomenon that occurs as radio signals reflect off of surfaces, such as buildings, producing multiple copies of the original signal. When this occurs with GPS signals, it results in one or more delayed signals arriving at the receiver with or without the on-time/direct GPS signal. The receiver measures the composite of these signals which, depending on the severity of the multipath, can substantially degrade the accuracy of the receiver's calculated position. Multipath is commonly experienced in cities due to tall buildings and its mitigation is an ongoing area of study. This research demonstrates a novel approach for simulating GPS multipath through the modification of an open-source tool, GPS-SDR-SIM. The resulting additional testing capability could allow for improved development of multipath mitigating technologies. Currently, open-source tools for simulating GPS signals are available and can be used in the testing and evaluation of GPS receiver equipment. These tools can generate GPS signals that, when used by a GPS receiver, result in computation of a position solution that was pre-determined at the time of signal generation. That is, the signals produced are properly formed for the pre-determined location and result in the receiver reporting that position. This allows for a GPS receiver under test to be exposed to various simulated locations and conditions without having to be physically subjected to them. Additionally, while these signals are generated by a software simulation, they can be processed by real or software defined GPS receivers. This work utilizes the GPS-SDR-SIM software tool for GPS signal generation and while this tool does implement some sources of error that are inherent to GPS, it cannot inject multipath. GPS-SDR-SIM was modified in this effort to produce additional copies of signals with pre-determined delays. These additional delayed signals mimic multipath and represent what happens to GPS signals in the real world as they reflect off of surfaces and arrive at a receiver in place of or alongside the direct GPS signal. A successful proof of concept was prototyped and demonstrated using this modified version of GPS-SDR-SIM to produce simulated GPS signals as well as additional simulated multipath signals. The generated data was processed using a software defined GPS receiver and it was found that the introduction of simulated multipath signals successfully produced the expected characteristics of a composite multipath signal. Further maturation of this work could allow for the development of a GPS receiver testing and evaluation framework and aid in the development of multipath mitigating technologies.

GPS↗
NASA NTRSDate not supplied

High-Enthalpy Testing to Validate Simulation of an Aerosol Capture Probe

AERACEPT (Aerosol Rapid Analysis Combined Entry Probe/sonde Technology) is a developing technology enabling in-situ aerosol particle sampling and analysis in a small spacecraft mission envelope. It integrates a passive aerosol sample collection system into a probe’s thermal protection system (TPS) to remove the need for heat shield separation and active descent control (parachutes, gliders, etc.). AERACEPT is being validated against the requirements of the Nephele mission concept, which targets the middle and lower Venus cloud layers. The proposed technology employs 3D Carbon-Carbon (3D-CC) at the probe’s nose to withstand the extreme environments encountered in a Venus entry without producing pyrolysis gases that would contaminate the mission sample collection. Any shape change of this 3D-CC sample inlet will affect airflow through the sampling system, thereby impacting aerosol sample collection bias and efficiency. Accurate prediction of the nose’s material response during entry is there-fore of fundamental importance to the development and practical use of this technology. As a result, the AERACEPT project is planning a test campaign in the newly built PlasmatronX facility with a novel “Open Iso-Q” test article to validate state-of-the-art material response tools for a simulated Venus entry. The PlasmatronX is a 350 kW inductively-coupled plasma facility developed and run by the Center for Hypersonics and Entry Systems Studies at the University of Illinois at Urbana-Champaign. The facility can support ground testing for a variety of planetary destinations by simulating entry conditions in Nitrogen, Air, and Carbon Dioxide. To best simulate the Venus atmosphere and match flight-like recession, the planned AERACEPT test campaign will use Carbon Dioxide as its test gas. AERACEPT has designed two test articles for this campaign: 1) a standard test article whose curved surface approximates a constant applied heat flux (“Closed Iso-Q” model) and 2) a novel inlet test article with a through-hole at the stagnation point (“Open Iso-Q” model). Both article types consist of a 3D-CC sample bonded to a graphite fixture. All articles will be laser scanned before and after testing to estimate shape change and recession, and the “Closed Iso-Q” articles will be instrumented with thermocouples to provide temperature histories. The test article temperatures and inlet shape change will be compared against simulations run with the Porous material Analysis Toolbox based on OpenFOAM (PATO) to validate the AERACEPT 3D-CC material response model. Due to the high thermal conductivity of 3D-CC, multidimensional effects have a large influence on the test article temperature and the standard 1-D material response tools are not sufficient for test planning. The test campaign features (e.g. article design, run duration, heat fluxes) are therefore informed by a series of 3-D simulations using state-of-the-art material response tools such as PATO. The test campaign is scheduled for spring of 2024. Expected results include 1) Key drivers of selected test conditions (e.g. matching mission flight atmosphere & recession), 2) Simulation results that informed test article design, 3) Acquired test data such as recession measurements, temperature histories, and photos, and 4) Comparisons between test results and simulation predictions.

AERACEPT↗
NASA NTRSDate not supplied

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗
NASA NTRS2024

Expansion of Check-Cases for 6DOF Simulation

This effort expands upon a previous NASA activity that developed flight simulation benchmark check-cases to include new check-cases for the Cislunar domain, comparing multiple NASA simulation tools. The results of this effort describe the benefits of standardizing inputs, simulation comparisons and describe an interactive website that enables comparison of externally provided simulation data. Participating simulations improved their software and identified implementation errors. This activity elevated simulation credibility and provided a measure of validation for the simulations actively in use for NASA’s Human Landing Systems (HLS).

Modeling↗
DOE OSTI2017

Optimizing for Large Planar Fractures in Multistage Horizontal Wells in Enhanced Geothermal Systems Using a Coupled Fluid and Geomechanics Simulator

Enhanced Geothermal Systems (EGS) could potentially use technological advancements in coupled implementation of horizontal drilling and multistage hydraulic fracturing techniques in tight oil and shale gas reservoirs along with improvements in reservoir simulation techniques to design and create EGS reservoirs. In this study, a commercial hydraulic fracture simulation package, Mangrove by Schlumberger, was used in an EGS model with largely distributed pre-existing natural fractures to model fracture propagation during the creation of a complex fracture network. The main goal of this study is to investigate optimum treatment parameters in creating multiple large, planar fractures to hydraulically connect a horizontal injection well and a horizontal production well that are 10,000 ft. deep and spaced 500 ft. apart from each other. A matrix of simulations for this study was carried out to determine the influence of reservoir and treatment parameters on preventing (or aiding) the creation of large planar fractures. The reservoir parameters investigated during the matrix simulations include the in-situ stress state and properties of the natural fracture set such as the primary and secondary fracture orientation, average fracture length, and average fracture spacing. The treatment parameters investigated during the simulations were fluid viscosity, proppant concentration, pump rate, and pump volume. A final simulation with optimized design parameters was performed. The optimized design simulation indicated that high fluid viscosity, high proppant concentration, large pump volume and pump rate tend to minimize the complexity of the created fracture network. Additionally, a reservoir with 'friendly' formation characteristics such as large stress anisotropy, natural fractures set parallel to the maximum horizontal principal stress (SHmax), and large natural fracture spacing also promote the creation of large planar fractures while minimizing fracture complexity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Fast Learning for Immersive Engagement in Energy Simulations

The fast computation which is critical for immersive engagement with and learning from energy simulations would be furthered by developing a general method for creating rapidly computed simplified versions of NREL's computation-intensive energy simulations. Created using machine learning techniques, these 'reduced form' simulations can provide statistically sound estimates of the results of the full simulations at a fraction of the computational cost with response times - typically less than one minute of wall-clock time - suitable for real-time human-in-the-loop design and analysis. Additionally, uncertainty quantification techniques can document the accuracy of the approximate models and their domain of validity. Approximation methods are applicable to a wide range of computational models, including supply-chain models, electric power grid simulations, and building models. These reduced-form representations cannot replace or re-implement existing simulations, but instead supplement them by enabling rapid scenario design and quality assurance for large sets of simulations. We present an overview of the framework and methods we have implemented for developing these reduced-form representations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Multi-physics Object Oriented Simulation Environment (MOOSE)

This two-day workshop is a hands-on introduction to the Multi-Physics Object Oriented Simulation Environment (MOOSE). MOOSE is a simulation framework built from PETSc that enables efficient simulation of complex physical systems with minimal programming background. Attendees are encouraged to visit the MOOSE website at www.mooseframework.org and walk through the step by step process of installing the framework on their computer. Support through this process is provided through the MOOSE Google Group message board. The workshop itself will consist of interactive modules where instructors from the MOOSE development team will work with attendees to build a simple plasma model. Each module will build on previous modules to introduce fundamental MOOSE features that will combine to demonstrate simulation capabilities for low temperature plasma modeling. The intention is not to provide an exhaustive overview of simulation capabilities but to instead provide attendees with the necessary tools to build unique simulation capabilities for their specific research challenges. Users should attend the workshop with laptop computers with MOOSE installed to allow for interactive instruction. At the completion of the workshop attendees will be able to install and run MOOSE and associated applications as well as carry out plasma simulations studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

VALIDATION AND DETERMINATION OF SIGNIFICANT SIMULATION PARAMETERS USING THE SMOOTHED PARTICLE HYDRODYNAMIC CODE NEUTRINO

In recent years there has been an increase in flooding concern to nuclear power plants; thus, tools to reduce the cost of mitigation efforts are valuable to industry. Simulating flooding events is achieved through different methods and at varied levels of detail. However, these simulations must be able to match application conditions, which requires code validation as well as knowledge of which parameters affect the accuracy of the simulation. To validate and determine the significant parameters for these tools, a sloshing experiment was conducted at the George Washington University and then simulated with the smoothed particle hydrodynamic code Neutrino. The simulation setup was constructed to match the experiment as closely as possible, and the end-wall pressure results were analysed. The 90% and 50% end-wall pressure bounds were compared along with the pressure impulse. The simulation bounds fell mostly within the experimental bounds, with exceptions in insignificant low-pressure areas, such as at the beginning (~40% difference) and the end of the cycle (~40% difference). Another exception was seen at a point in the peak pressure (~75% difference). The simulated pressure impulse matched the experiments well and was within 10%. Additionally, seven parameters were qualitatively assessed and then investigated to determine their significance. The parameters’ significance was determined by sampling five values across the parameters’ range using the Risk Analysis Virtual Environment, then comparing the average pressure for each run. All seven parameters were identified as having some significance. However, the particle size, interaction-radius to particle-size ratio, and fluid settling seem to have a greater effect due to the larger fluctuation of the results. Now that the significant parameters are identified, further research is being done to quantify the significance of each parameter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

SimNet: Accurate and High-Performance Computer Architecture Simulation using Deep Learning

While cycle-accurate simulators are essential tools for architecture research, design, and development, their practicality is limited by an extremely long time-to-solution for realistic applications under investigation. This work describes a concerted effort, where machine learning (ML) is used to accelerate microarchitecture simulation. First, an ML-based instruction latency prediction framework that accounts for both static instruction properties and dynamic processor states is constructed. Then, a GPU-accelerated parallel simulator is implemented based on the proposed instruction latency predictor, and its simulation accuracy and throughput are validated and evaluated against a state-of-the-art simulator. Leveraging modern GPUs, the ML-based simulator outperforms traditional CPU-based simulators significantly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Roadmap for IES Modeling and Simulation Activities

Since inception of the DOE-NE Integrated Energy Systems (IES) program, several program and modeling visions and roadmaps have been published. In 2017, the last modeling and simulation capability development plan was published. In 2020, a comprehensive roadmap for the DOE-NE Integrated Energy Systems program has been published. The current report is meant to update the 2017 modeling and simulation capability roadmap and complement the IES overall program roadmap published in 2020. The role of modeling and simulation within IES is to support the demonstration of new coupled integrated technologies along every step of the technology maturation from the strategic analysis of preferred system architecture to preliminary design, to laboratory testing, up to full commercial testing and integration. To achieve this role, modeling and simulation must be able to assess the technical performance and the economic viability of potential IES. Also, modeling and simulation must provide support for experimental evaluations, i.e., support component design and real time operations of experimental demonstration systems. Finally, modeling and simulation can help scale-up experimental systems to the final commercial systems. This report details the current state of the modeling and simulation efforts within IES for the above-mentioned areas, provides a gap analysis and proposes next steps for a time horizon over the next 4 to 5 years.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Simulating the Effect of Fracture Connectivity on Repository Performance with GDSA Framework - 18589

This work examines the possible migration of radionuclides from a deep geologic nuclear waste repository sited in fractured crystalline host rock. The key safety concern is the potential effect on waste isolation of the inter-connectivity of the fracture network, which is primarily established by the temporal evolution of the temperature and stress fields at the time of original rock intrusion and cooling. Two end members are considered in these simulations, one with a high degree of connectivity to the biosphere, such that advective transport through the fracture and fault network controls radionuclide migration, and the other with a low degree of connectivity, such that slow diffusive transport through the crystalline rock matrix is the controlling process for radionuclide migration to the biosphere. Both end members have several areally extensive, high-transmissivity deformation zones, whose distance to the repository is an important factor for waste isolation capability. Uncertainties in fracture properties (transmissivity, orientation, radius) and fracture distribution also give rise to uncertainty in waste isolation capability (i.e., the overall connectivity of the repository with the surface biosphere). These types of natural system uncertainties are generally present in any geologic site-characterization program and thus represent an important factor in assessing repository performance in hard rock environments. The waste considered in these simulations are spent fuel rod assemblies from the U.S. commercial nuclear reactor fleet. Thermal output of the waste must be considered in the simulations and can influence the rate and timing of waste package failure, waste form degradation, and fluid flux due to thermal expansion around the repository horizon. To achieve a high-fidelity representation of radionuclide transport in fractures and rock matrix, combined with thermal energy transport and fluid flow in fractures and matrix, the mathematical model is solved numerically in a parallel high-performance computing (HPC) environment on a finite volume unstructured grid consisting of approximately 5 million cells, using Geologic Disposal Safety Assessment (GDSA) Framework (https://pa.sandia.gov), an open-source performance assessment tool for deep underground disposal of nuclear waste. GDSA Framework uses PFLOTRAN to solve the balance equations on a three-dimensional grid with heterogeneous properties, using multiple processors in a parallel configuration based on domain decomposition. The fracture networks in these simulations are originally generated as discrete fracture networks (DFNs), which are sets of two-dimensional planes distributed in a three-dimensional domain. The method used in GDSA Framework maps the stochastically generated DFN to an equivalent continuous porous medium (ECPM) domain that allows for the simulation of coupled heat flow, fluid flow, and radionuclide transport, including heat conduction through the matrix of the fractured rock, which is a process not easily modeled in a DFN representation. Computational efficiency is also greatly enhanced using the ECPM method, allowing for a realistic representation and analysis of uncertainties in a multi-realization performance assessment of a deep geologic repository. The effect of fracture connectivity on the waste isolation safety function, as brought to light by these GDSA Framework simulations, points to the importance of including a realistic representation of uncertainties in fracture properties and distribution (effectively, uncertainty in spatial heterogeneity) in repository safety assessment simulations. (authors)

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

The Next Generation System Analysis Model: Capabilities for Simulating a Waste Management System - 19131

The U.S. Department of Energy (DOE) Office of Nuclear Energy Office of Spent Fuel and Waste Disposition's Integrated Waste Management System (IWMS) program is applying integrated waste management system architecture analysis, system engineering, and decision analysis principles to inform potential future decisions regarding potential nuclear waste management system architectures. Architecture analyses of the IWMS are being conducted to support the future deployment of a comprehensive system for managing nuclear waste that considers all major aspects of the back end of the nuclear fuel cycle (i.e., transportation, storage, and disposal). The system architecture analyses consider alternative disposition pathways for commercial spent nuclear fuel (SNF) being stored dry in at-reactor Independent Spent Fuel Storage Installations (ISFSIs) or wet in spent-fuel pools. One method for evaluating these alternatives is simulating the movement of SNF throughout the system components. Successfully simulating a system on the scale of a national waste management system requires - accurate baseline information (e.g., historic fuel discharges from reactors and dry storage loadings at reactor sites), - modeling all logistic interactions between components of the proposed system, and - tracking and reporting all relevant details regarding the simulation's execution. The Next Generation System Analysis Model (NGSAM) is an agent-based simulation software tool designed for the express purpose of modeling the IWMS. It is based on the Performance Analysis Tool, a simulation platform that has been used to model logistics for the Department of Defense and the Federal Emergency Management Agency. Using existing, government-owned software allows NGSAM developers to modify the software at the level of source code and to develop flexible usability and modeling solutions for developing and running simulations of the IWMS. NGSAM imports data from the Oak Ridge National Laboratory (ORNL) Unified Database (e.g., historic assembly information, thermal profiles for assembly heat, at-reactor dry storage loadings) to ensure that the simulation initializes with a realistic representation of the state of commercial SNF in the United States. NGSAM is designed with the flexibility to model and track the behavior of virtually any IWMS architecture assumptions or design considerations. For example, NGSAM can model multiple interim storage facilities or repositories or none at all. Similarly, all assumptions regarding shipment schedules, operating schedules, and reactor operation assumptions can be modeled and adjusted in NGSAM with relative ease. The development of NGSAM was targeted to waste management system analysts as the primary users. With that in mind, analysts from ORNL, Sandia National Laboratories (SNL), and Argonne National Laboratory (ANL) were intimately involved in the development process. That involvement extended beyond merely defining requirements. System analysts have been involved with design decisions and software benchmarking and tests, and they have provided subject-matter expertise where necessary. The result is a high degree of interaction between professional software developers and system analysts, with a multitude of benefits to both groups. This paper explores some of the underlying methodology of the agile software development process used in the creation of NGSAM. It illustrates the capabilities of NGSAM. Some example use cases are discussed as well, with relevant results providing context for how NGSAM is used by system analysts to gain a greater understanding of alternative waste management system architectures. This is a technical paper that does not take into account the contractual limitations under the Standard Contract for Disposal of Spent Nuclear Fuel and/or High-Level Radioactive Waste (Standard Contract) (10 CFR Part 961). Under the provisions of the Standard Contract, DOE does not consider spent nuclear fuel in multi-assembly canisters to be an acceptable waste form, absent a mutually agreed-to contract amendment. To the extent discussions or recommendations in this paper conflict with the provisions of the Standard Contract, the Standard Contract provisions prevail. (authors)

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

Kinetic Monte Carlo Simulations of Chemical Reactions on Surfaces

Lattice based kinetic Monte Carlo (kMC) simulations are a powerful tool for micro-kinetic modeling. The kinetic Monte Carlo technique stochastically simulates chemical reaction processes (e.g. adsorption, desorption, diffusion), where the probability of a given reaction process being selected is weighted by its rate constant. The open source software package, KMOS1, allows fast construction of kMC simulations by automatically implementing a computationally efficient kMC framework, thus allowing the researcher to focus on accurate parameterization and formulation of the kinetic model. From this, the spatiotemporal evolution (e.g. concentration, surface configuration, turnover frequencies) of crystal surfaces or bulk crystalline structures is obtained. KMOS uses the variable step-size method implementation of kMC. At each step, the VSSM updates available events, calculates the corresponding transition rates and simulation time, and selects one of the available events. The lattice based kMC technique has many potential applications within SRNL where chemical kinetics are of interest. By transitioning through time along a trajectory of configuration space, fundamental insights to a system's dominant reactions and parameters can be gained, leading to improved materials and process design. In this work, we have established capabilities for fast implementation of kMC simulations through the installation and testing of the KMOS software on SRNL high performance computing resources. The parameterization of kMC simulations is often performed from first-principles energetic calculations. For the purposes of this work, the test cases use arbitrary species and thus, have arbitrary parameterization. In order to demonstrate the use of KMOS, two primary test cases are shown. Case 1 involves simple adsorption and desorption. Case 2 involves a slightly more complex mechanism. The evolution of the surface concentration and the final surface configuration for simulation Case 2 is shown. Rate-limiting processes and rate-limiting configurations are not always known or understood, but kMC provides information pertaining to both. The capability to implement efficient micro-kinetic models using the kinetic Monte Carlo technique has been established through the testing and installation of KMOS on SRNL HPC resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Regional-scale fault-to-structure earthquake simulations with the EQSIM framework: Workflow maturation and computational performance on GPU-accelerated exascale platforms

Continuous advancements in scientific and engineering understanding of earthquake phenomena, combined with the associated development of representative physics-based models, is providing a foundation for high-performance, fault-to-structure earthquake simulations. However, regional-scale applications of high-performance models have been challenged by the computational requirements at the resolutions required for engineering risk assessments. The EarthQuake SIMulation (EQSIM) framework, a software application development under the US Department of Energy (DOE) Exascale Computing Project, is focused on overcoming the existing computational barriers and enabling routine regional-scale simulations at resolutions relevant to a breadth of engineered systems. This multidisciplinary software development—drawing upon expertise in geophysics, engineering, applied math and computer science—is preparing the advanced computational workflow necessary to fully exploit the DOE’s exaflop computer platforms coming online in the 2023 to 2024 timeframe. Achievement of the computational performance required for high-resolution regional models containing upward of hundreds of billions to trillions of model grid points requires numerical efficiency in every phase of a regional simulation. This includes run time start-up and regional model generation, effective distribution of the computational workload across thousands of computer nodes, efficient coupling of regional geophysics and local engineering models, and application-tailored highly efficient transfer, storage, and interrogation of very large volumes of simulation data. This article summarizes the most recent advancements and refinements incorporated in the workflow design for the EQSIM integrated fault-to-structure framework, which are based on extensive numerical testing across multiple graphics processing unit (GPU)-accelerated platforms, and demonstrates the computational performance achieved on the world’s first exaflop computer platform through representative regional-scale earthquake simulations for the San Francisco Bay Area in California, USA.

58 GEOSCIENCES↗
DOE OSTI2016

Dependence of radiation belt simulations to assumed radial diffusion rates tested for two empirical models of radial transport

Radial diffusion is one of the dominant physical mechanisms that drives acceleration and loss of the radiation belt electrons, which makes it very important for nowcasting and forecasting space weather models. In this paper, we investigate the sensitivity of the two parameterizations of the radial diffusion of Brautigam and Albert (2000) and Ozeke et al. (2014) on long-term radiation belt modeling using the Versatile Electron Radiation Belt (VERB). Following Brautigam and Albert (2000) and Ozeke et al. (2014), we first perform 1-D radial diffusion simulations. Comparison of the simulation results with observations shows that the difference between simulations with either radial diffusion parameterization is small. To take into account effects of local acceleration and loss, we perform 3-D simulations, including pitch angle, energy, and mixed diffusion. We found that the results of 3-D simulations are even less sensitive to the choice of parameterization of radial diffusion rates than the results of 1-D simulations at various energies (from 0.59 to 1.80 MeV). This result demonstrates that the inclusion of local acceleration and pitch angle diffusion can provide a negative feedback effect, such that the result is largely indistinguishable simulations conducted with different radial diffusion parameterizations. Finally, we also perform a number of sensitivity tests by multiplying radial diffusion rates by constant factors and show that such an approach leads to unrealistic predictions of radiation belt dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

The Effect of a Guide Field on Local Energy Conversion During Asymmetric Magnetic Reconnection: Particle-in-Cell Simulations

We use theory and simulations to study how the out-of-plane (guide) magnetic field strength modifies the location where the energy conversion rate between the electric field and the plasma is appreciable during asymmetric magnetic reconnection, motivated by observations (Genestreti et al., 2017). For weak guide fields, energy conversion is maximum on the magnetospheric side of the X line, midway between the X line and electron stagnation point. As the guide field increases, the electron stagnation point gets closer to the X line, and energy conversion occurs closer to the electron stagnation point. We motivate one possible nonrigorous approach to extend the theory of the stagnation point location to include a guidefield. The predictions are compared to two-dimensional particle-in-cell (PIC) simulations with vastly different guide fields. The simulations have upstream parameters corresponding to three events observed with Magnetospheric Multiscale (MMS). The predictions agree reasonably well with the simulation results,capturing trends with the guide field. The theory correctly predicts that the X line and stagnation points approach each other as the guide field increases. The results are compared to MMS observations, Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE) observations of each event,and a global resistive-magnetohydrodynamics simulation of the 16 October 2015 event. The PIC simulation results agree well with the global observations and simulation but differ in the strong electric fields and energy conversion rates found in MMS observations. The observational, theoretical, and numerical resultssuggest that the strong electric fields observed by MMS do not represent a steady global reconnection rate.

79 ASTRONOMY AND ASTROPHYSICS↗