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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Wind Farm Simulation and Layout Optimization in Complex Terrain: Preprint

This work reports on incorporating complex terrain into wind farm simulations for the purpose of layout optimization. Adding complex terrain boundary conditions to NREL's medium fidelity computational fluid dynamics model, WindSE, produces significant separation, flow curvature, and speedup effects that would otherwise be difficult to capture with lower-fidelity models or a flat-terrain assumption. These flow features, in turn, can significantly impact the optimal turbine array layout. We demonstrate the impact of complex terrain on flow in both an idealized and real-world setting, and discuss modifications to the code that enable gradient-based optimization using terrain-aware adjoint gradients. Through several optimization case studies, we show that the layout optimization process takes advantage of speedup effects on terrain high points, and leverages flow curvature effects that modify wake trajectories. This yields substantial power improvements over gridded layouts, and hints at future research directions in simulation and optimization for wake trajectories in complex terrain.

17 WIND ENERGY↗
DOE OSTI2024

Transformational Regional-Scale Earthquake Simulations with the DOE EarthQuake SIMulation Exascale Framework

Earthquakes present worldwide risk to economic and human safety. The 2023 earthquakes in Turkiye provided a reminder of the potential for catastrophic consequences with 50,700 deaths and 15.7 million people affected. The ability to predict ground motions and infrastructure damage for earthquakes continues to be a challenging problem for scientists and engineers. Until now, estimates of ground motions have been performed empirically by looking at sparse data from past earthquakes. This approach can provide statistical information on intensity amplitudes but cannot inform site-specific ground motions essential to developing the most effective resilience. Interest has grown in large-scale computational models to simulate earthquakes at regional scale. The U.S. Department of Energy EarthQuake SIMulation (EQSIM) framework was developed for regional-scale earthquake simulations at unprecedented fidelity, taking advantage of emerging GPU-accelerated systems. This article describes the EQSIM workflow and demonstrates regional-scale simulations with the new computational capability available to scientists in their quest to mitigate future disasters.

58 GEOSCIENCES↗
DOE OSTI2021

Predicting the Dissolution Rate of Borosilicate Glasses using QSPR analysis based on molecular dynamics simulations

Quantitative Structure Property Relationship (QSPR) analysis based on molecular dynamics (MD) simulations is a promising approach for establishing the composition-property relationships of glass and other materials with complex structures. A series of 20 borosilicate and boroaluminosilicate glasses have been modeled by using MD simulations with recently developed effective potentials. Short- and medium-range structures of these glasses were analyzed and, based on these structural information, QSPR analysis of the initial dissolution rate (r0) was made and compared with measured r0 at 90°C and pH 9 using various structural descriptors such as percentage of bridging oxygen species, network connectivity and average ring size. The structural descriptors, Fnet, containing energetic information such as single bond strength and other structural information were also used. It was found that overall network connectivity, average ring size and Fnet give reasonable predictions of the r0 of studied glasses, given the conditions that the glasses are homogeneous and dissolve congruently. Modifying glass compositions to account preferential release of modifiers gives a better prediction for incongruently dissolving glasses. The results were compared with our recent work of predicting glass dissolution behavior from compositions using the topological-constraints-based models.

Du, Jincheng↗
DOE OSTI2024

The Pele Simulation Suite for Reacting Flows at Exascale

In this work, we present the Pele suite of software tools for compressible and incompressible reacting flows. The Pele suite leverages several different libraries, notably AMReX and SUNDIALS, to achieve performance portability on heterogeneous computing architectures across the supercomputing landscape. The Pele suite is comprised of PeleC, a compressible reacting flow block-structured adaptive mesh refinement solver, PeleLMeX, a low-Mach number reacting flow block-structured adaptive mesh refinement solver, Pele-Physics, a library for transport, thermodynamics, finite rate chemistry, soot, spray and radiation physics. The objective of this paper is (i) to present the code development efforts necessary to achieve highly effective and scalable applications for exascale machines and (ii) to detail the performance results of the Combustion-Pele project applications on Oak Ridge National Laboratory's Frontier. We show good weak and strong scaling results for both PeleC and PeleLMeX up to more than 50 billion cells on more than 4096 Frontier graphics processing unit nodes. We also present a capability demonstration simulation of a dual-fuel pulse compression ignition engine (six adaptive mesh refinement levels, and 60 billion cells or 2.1 trillion degrees of freedom) on Frontier, to date one of the largest simulations performed on the first exascale-class supercomputer.

adaptive mesh refinement↗
DOE OSTI2025

Hydrologic Model Data for the East Fork Poplar Creek Watershed Simulated with the Advanced Terrestrial Simulator (ATS): Streamflow and Network Expansion–Contraction Dynamics

This dataset supports hydrologic modeling and stream network expansion–contraction analysis for the East Fork Poplar Creek (EFPC) Watershed in Tennessee. It includes a Jupyter notebook for model setup, model configuration files, simulation outputs, and derived products used to evaluate model performance and investigate stream dynamics under varying hydrologic conditions. The dataset was generated using the Watershed Workflow Python package and the Advanced Terrestrial Simulator (ATS), enabling integrated surface–subsurface hydrologic simulations using a stream-aligned mesh. Outputs include high-resolution time series of streamflow, active network length, water table depth, and related hydrologic variables. Also included are spatially explicit stream persistency indices and classifications of reaches as perennial or non-perennial. These data facilitate reproducibility and support further research on stream intermittency and variability in network extent.The model data archive is organized in following directories:1) model_setup_inputsContains the Watershed Workflow Jupyter notebooks (accessed through any open source code editor), selected input datasets, and resulting ATS input files, including XML files (access through any open source code editor), computational mesh (.exo files can be viewed using Paraview), and meteorological forcing files (.h5 files can be accessed through h5py python package and HDFView open source software). 2) model_outputsIncludes ATS simulation outputs relevant to this study. Time series of spatially integrated or averaged variables (e.g., streamflow, water table depth) are provided as CSV files. Select spatial fields (e.g., ponded depth and water table depth) are saved as pickled Python objects to reduce file size, and can be accessed through pickle package in Python. Key geometry objects from Watershed Workflow—such as the surface mesh and river tree—are also included to support analysis of streamflow persistency and expansion–contraction dynamics. These files can also be accessed through Watershed Workflow Python package.3) model_evaluationProvides observed streamflow time series and field survey-based flow regime classifications used to evaluate model performance. Jupyter notebooks for processing ATS outputs and comparing model predictions with observations to build confidence in the model prior to scientific analysis are also included.4) Q_L_relationshipsContains workflows for generating time series of discharge, active network length, and related hydrologic variables used in the stream network expansion–contraction analysis. Includes routines for delineating baseflow-dominated periods. For each catchment, notebooks and processed data (as pickled DataFrames accessed through Pandas Python package) are provided. 5) figure_scriptsProvides the Jupyter notebooks used to generate the figures presented in the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Monte Carlo simulation analysis for radiation damage on Glidcop Al-15 caused by 17-20 MeV/u heavy ions

The Facility for Rare Isotope Beams (FRIB) SRF linear heavy-ion accelerator is designed to accelerate all ions up to uranium to a maximum beam power of 400 kW. Several beam intercepting devices (BIDs) are essential to the successful operation of the accelerator, including a low power charge selector (LPCS) made of copper containing 0.15% precipitated aluminum oxide by weight called Glidcop Al-15. As FRIB ramps up the primary beam power beyond the current 20 kW level, the charge selector must withstand higher radiation damage rates, typically measured in displacements per atom (dpa). Significant beam induced radiation damage including significant deformation, like swelling, blistering, and cracking has been observed on a recently removed LPCS. We also observed physical features up to about 5.5 mm wide and appear to be deeper than the projected range of any ions. This paper presents the results of Monte Carlo simulations carried out using the Particle and Heavy Ion Transport code System (PHITS), quantifying the total damage dose and ion concentration, etc. Simulating an accurate irradiation history is essential to determining the scope of post irradiation examination (PIE) work.

Fletcher, Ethan [Michigan State U.]↗
DOE OSTI2021

Impact of Rheology Models on Horizontal Well Polymer Flooding in a Heavy Oil Reservoir on Alaska North Slope: A Simulation Study

Polymer rheology can have either a positive or a negative effect on polymer flooding performance under varied circumstances. Many researchers have studied the effect of polymer rheology in a vertical well, but no field scale studies have been conducted to investigate whether polymer rheology is beneficial to polymer flooding in heavy oil reservoirs developed by horizontal wells. In this paper, we conducted a numerical simulation study to examine the effect of HPAM polymer rheology on a polymer flooding pilot, which is the first-ever project conducted on a heavy oil reservoir from Alaska North Slope (ANS) developed by horizontal wells. Three rheology types were considered in the study including the apparent viscosity measured during coreflooding of using a HPAM polymer, the bulk viscosity measured with a viscometer, and a Newtonian flow model. The results suggest that using the bulk viscosity in simulation underestimates the conformance control and the water-oil-ratio reduction capability of the HPAM polymer solution. When the apparent viscosity is used, the incremental oil and sweep were largely increased, and the optimal recovery period of polymer flooding was extended greatly, especially for the heterogeneous formations. Therefore, the rheology type of polymer plays a significant role in the incremental oil recovery and injection profile of the horizontal well system given the pilot test conditions. This study has provided practical guidance to field operators for the ongoing polymer flooding pilot on ANS and will also provide valuable information for other polymer projects conducted in similar conditions.

Leng, Jianqiao↗
DOE OSTI2021

Identifying the sources of uncertainty in climate model simulations of solar radiation modification with the G6sulfur and G6solar Geoengineering Model Intercomparison Project (GeoMIP) simulations

We present here results from the Geoengineering Model Intercomparison Project (GeoMIP) simulations for the experiments G6sulfur and G6solar for six Earth system models participating in the Climate Model Intercomparison Project (CMIP) Phase 6. The aim of the experiments is to reduce the warming that results from a high-tier emission scenario (Shared Socioeconomic Pathways SSP5-8.5) to that resulting from a medium-tier emission scenario (SSP2-4.5). These simulations aim to analyze the response of climate models to a reduction in incoming surface radiation as a means to reduce global surface temperatures, and they do so either by simulating a stratospheric sulfate aerosol layer or, in a more idealized way, through a uniform reduction in the solar constant in the model. We find that over the final two decades of this century there are considerable inter-model spreads in the needed injection amounts of sulfate (29 ± 9 Tg-SO 2 /yr between 2081 and 2100), in the latitudinal distribution of the aerosol cloud and in the stratospheric temperature changes resulting from the added aerosol layer. Even in the simpler G6solar experiment, there is a spread in the needed solar dimming to achieve the same global temperature target (1.91 ± 0.44 %). The analyzed models already show significant differences in the response to the increasing CO 2 concentrations for global mean temperatures and global mean precipitation (2.05 K ± 0.42 K and 2.28 ± 0.80 %, respectively, for SSP5-8.5 minus SSP2-4.5 averaged over 2081–2100). With aerosol injection, the differences in how the aerosols spread further change some of the underlying uncertainties, such as the global mean precipitation response (–3.79 ± 0.76 % for G6sulfur compared to –2.07 ± 0.40 % for G6solar against SSP2-4.5 between 2081 and 2100). These differences in the behavior of the aerosols also result in a larger uncertainty in the regional surface temperature response among models in the case of the G6sulfur simulations, suggesting the need to devise various, more specific experiments to single out and resolve particular sources of uncertainty. The spread in the modeled response suggests that a degree of caution is necessary when using these results for assessing specific impacts of geoengineering in various aspects of the Earth system. However, all models agree that compared to a scenario with unmitigated warming, stratospheric aerosol geoengineering has the potential to both globally and locally reduce the increase in surface temperatures.

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

Kinetic Plasma Simulation in the MOOSE Framework: Verification of Electrostatic Particle In Cell Capabilities

In magnetic confinement nuclear fusion reactors, the interaction between the plasma edge and plasma facing components is extremely important. At the plasma edge, a kinetic representation such as particle-in-cell (rather than a fluid representation) is required to accurately capture the plasma behavior. General purpose particle-in-cell plasma simulation capabilities have been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. This new capability is a part of the development of a new MOOSE-based framework for modeling plasma facing components, the Fusion ENergy Integrated multiphys-X (FENIX) framework. In this work, the verification of foundational particle-in-cell capabilities in FENIX is presented. This new plasma simulation capability has three main components: moving particles in discrete steps on the finite element mesh, mapping charge density from the particle's location to the finite element mesh, and solving for the electrostatic potential based on the charge density mapped from particles to the mesh. In this paper, simple verification problems demonstrating each of these new capabilities are presented, and future work includes electromagnetic capabilities and Monte Carlo collisions with neutral gas particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Bridging particle and reactor scales in the simulation of biomass fast pyrolysis by coupling particle resolved simulation and coarse grained CFD-DEM

Efficient utilization of biomass is challenging due to the complexity of the feedstock at molecular, particle, and reactor scales. Direct coupling of particle scale model and reactor scale simulation is computationally infeasible. In this work, the results from a particle resolved simulation were used to calibrate the heat transfer and reaction kinetics of particles. A hybrid drag model was used to consider the different fluidization properties of sands and biomass particles. The computation cost was reduced using a coarse-grained Discrete Element Method (DEM). The simulated conversion rate, char yield, tar yield, and pyrolysis gas yield compare well with experiment. The residence time and axial distributions of biomass were also analyzed. This novel multi-scale method provides an efficient and accurate tool for the modeling of biomass pyrolysis reactors.

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

SARS-CoV2 Protein-Ligand Simulation Dataset: Layer 1 (Simulation Initial Conditions and Parameters)

A set of 24 protein structures/complexes from the SARS CoV-2 proteome and inputs prepared for simulation using the CHARMM36m forcefield in PDB and gromacs formats. Each system contains a pdb and gromacs top and related input files necessary for running a temperature replica-exchange simulation. In addition, we also include: charmm PSF files (generated from the gromacs topology), a list of temperatures at which replica-exchange simulations were done (tempRamp), example gromacs run input mdp files, and initial minimized structures where available (minimized.pdb).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Multiscale thermal hydraulic coupling methods for boiling water reactor simulation

During the last two years, the Virtual Environment for Reactor Applications (VERA) has been extended to simulate boiling water reactors (BWRs). The thermal hydraulic effects present in BWRs are far more complex than those in pressurized water reactors (PWRs). Therefore, the runtime is significantly increase and convergence behavior worse compared to PWRs. Most of the additional runtime is spent during the first few coupled iterations, when the power shape is still rapidly evolving, dramatically affecting the thermal hydraulics (TH). To alleviate the increased computational expense, a multiscale TH coupling approach was developed in VERA in which a highly efficient, simplified TH model is solved for several coupled iterations until the power shape is partially converged. The simplified solution is used to inform the assembly-wise flow distribution in the high-fidelity TH solver CTF, reducing the amount of work required to properly balance the pressure drop in each channel and reducing runtime for the couple calculation. Because the final TH calculations are performed with CTF, there is ultimately no impact on the accuracy of the converged solution. This paper presents the details of this multiscale TH coupling approach, along with results for a 4x4 array of GE-14 fuel bundles and whole-core coupled simulations of the Hatch core. These two cases show that the multiscale approach dramatically reduces the runtime of coupled BWR simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

An HPC-Based Hydrothermal Finite Element Simulator for Modeling Underground Geothermal Behavior with Example Simulations on The Treasure Island and UC Berkeley Campus

This submission contains the source code of the Hydrothermal Finite Element Simulator used for the Treasure Island and UC Berkeley campus geothermal simulation. It contains a report that summarizes the development and validation of this Hydrothermal Finite Element Simulator, with a case study on Treasure Island site. It also contains a report that investigates the feasibility of upgrading the existing campus energy delivery system at UC Berkeley to a fifth-generation district heating and cooling system that includes geothermal heat/cold storage.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Production Through Simulation: Using Simulation Technologies to Create and Evaluate Nuclear Fuel Fabrication Facility Designs

This project explores how physics-based and discrete-event simulation technologies can be jointly employed to model a suite of nuclear fuel fabrication facility designs with variable radiological environments, optimize their design to lower operational costs subject to a set of design constraints, and evaluate their performance. Initially, three variants of U-20Pu-10Zr metallic nuclear fuel in three specific geometric configurations are modeled and studied across more than 50 physics simulations in nuclear physics software packages SCALE and Monte Carlo N-Particle Transport (MCNP6.2). From these tests, values for effective dose rate and criticality are recorded for each of the nine alloy-geometry configurations. These values are then incorporated into discrete-event models of fuel fabrication in the simulation program ExtendSim Pro 10 as process attributes. This allows one to track the dose to facility personnel as fuel is fabricated, dually enabling us to evaluate the safety of a given facility design in terms of absorbed dose per worker per year and to craft operational guidelines that adhere to federal and local safety regulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗