Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Use Case”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Numerical Modeling Framework for Flocculation and Cohesive Sediment Transport in the Wave Bottom Boundary Layer

Flocculation, a critical process in coastal and estuarine systems, plays a significant role in sediment transport, nutrient cycling, and ecological health. This study develops a cohesive sediment transport modeling framework tailored to the wave bottom boundary layer under dilute and equilibrium conditions, explicitly incorporating flocculation effects via a Population Balance Equation (PBE). Using Direct Numerical Simulation, six baseline cases, each with a distinct sediment concentration profile resulting from a constant settling velocity and critical erosion shear stress, are generated to drive the PBE flocculation model for given floc yield strength and stickiness. Results reveal that flocculation significantly influences sediment concentration profiles promoting three distinct stages, well‐mixed, transition to lutocline, and well‐developed lutocline. At low concentrations with well‐mixed profiles, cohesive floc properties are less significant, and turbulence is a main flocculation driver. In contrast, as concentration increases, cohesive floc properties become crucial, facilitating lutocline formation. Here, the analysis also highlights limitations of depth‐averaged settling velocity as a parameterization. It is suitable for well‐mixed and transitional profiles but fails in well‐developed lutoclines, where empirical formulations that explicitly incorporate turbulent shear rate and sediment concentration better capture variability. This study underscores the necessity of incorporating flocculation effects into sediment transport models to enhance predictions of sediment dynamics in wave bottom boundary layers.

Penaloza‐Giraldo, Jorge A. [Oak Ridge National Lab↗
DOE OSTI2025

2D reactive transport model of shale chemical weathering and biogeochemical fluxes along a mountainous hillslope, East River Watershed, Colorado: Input files and simulation results

This data package contains input files and simulation results for a two-dimensional (2D) reactive transport model used to quantitatively analyze the coupled hydrological and biogeochemical processes governing shale weathering and associated biogeochemical fluxes under realistic environmental conditions in the high-elevation East River Watershed. These data support the conclusions presented in Stolze et al. (Water Resources Research, under review), "Model-based interpretation of solute exports and carbon partitioning during shale weathering in a mountainous hillslope". The model simulates atmospheric-subsurface gas exchange, subsurface water flow, and shale weathering processes under dynamic, year-scale conditions along a shale-underlain hillslope located in the East River watershed. The simulations were performed using the PFLOTRAN flow and reactive transport code and executed on the Perlmutter supercomputer to leverage its large-scale parallel computing capabilities. The data package contains two zipped folders, "model_input_files" and "simulation_results", and one readme.txt file. "model_input_files" contains the necessary input files to run the calibrated base-base model presented in Stolze et al. (Water Resources Research, under review). "simulation_results" contains a single hdf5 file ("Output_2D_hillslope_model.h5") which includes the results of simulation performed using the base-case model. This file can be opened with HDFView 3.1.4, Python, or MATLAB. "readme.txt" contains relevant information about the base-case model and provides guidelines on how to run the associated input files provided in the folder "model_input_files". Furthermore, readme.txt provides information regarding the model results provided in "Output_2D_hillslope_model.h5" such as matrix dimensionality and output units. Field datasets used to evaluate model performance were collected at three monitoring wells located along a hillslope transect (PLM1, PLM2, and PLM3). Dissolved ion concentration data were collected from November 2016 to October 2021 for Ca, Mg, DIC, Na, K, SO4 (Dong et al., 2025 - dic_npoc_data_2014_2024.zip - DOI:10.15485/1660459; Williams et al., 2025 - anion_data_2014_2024.zip - DOI:10.15485/1668054; Dong et al., 2025 - cation_data_2014_2024.zip - DOI:10.15485/1668055). Note that we used the files named er_PLM1_xx_yy, er_PLM2_xx_yy, and er_PLM3_xx_yy where xx stands for the name of the aqueous species and yy stands for the depth where the measurements were performed. Soil water content ([0 - 1] m) and water table depth were collected from November 2016 to October 2021 (Wan et al., 2024 - Dynamic_water_table__depthsFig2b.csv and Soil_water_content_Fig4e.csv - DOI:10.15485/2322567). Gaseous CO2 concentration were collected from October 2020 to December 2021(Wan et al., 2024 - Soil_CO2_concentrations_Fig4h.csv - DOI:10.15485/2322567) Gaseous CO2 flux from the subsurface to the atmosphere were collected in the vicinity of PLM2 from October 2019 to May 2022 (Wu et al., 2025). Soil microbial biomass concentration was measured from August 2016 to June 2017 (Sorensen et al., 2019 - 2017_East_River_Pumphouse_Microbial_Biomass__1_.csv - DOI:10.15485/1577267) All field data are published as CSV files compatible with Microsoft Excel, MATLAB, and Python, or as text files. The coordinates of the monitoring wells and the CO2(g) flux sensor in the coordinate system WGS84 are: -PLM1: [38.9197710 ; -106.9492750] -PLM2: [38.9201580 ; -106.9487170] -PLM3: [38.9207843 ; -106.9483668] -PLM4: 38.9210060 ; -106.9479528] -CO2(g) flux sensor: [38.9199180 ; -106.9489906] ------------------------------------------------------------------------------------------- This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a Department of Energy User Facility using NERSC award BER-ERCAP 23980, BER-ERCAP 28550, and BER-ERCAP 33789.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Designing SOC electrodes for better lifetime performance

Presentation for talk at the 2024 Pacific Rim Meeting of the Electrochemical Society October 6-11, 2024. The talk covers a review of NETL's solid oxide cell (SOC) performance degradation modeling framework and present case studies on how simulations are used to design better SOC electrodes.

Abernathy, Harry↗
DOE OSTI2025

Impact of Urbanization on Convection, Lightning, and Precipitation Over the Houston Metropolitan Area: Case Study Simulation From the TRACER Campaign

This study investigates the effects of urbanization, specifically land use change and anthropogenic emissions (AE), on convection, lightning, and surface precipitation for a case of summertime sea‐breeze convection observed over the Houston metropolitan area. The unique capabilities of the NASA‐Unified Weather Research and Forecasting model allows us to conduct a series of sensitivity experiments with complex configurations, in particular including multi‐year land model spin‐up simulations, treatment of aerosols and their precursors, and explicit cloud charging and lightning. The simulation results show that urban land use primarily alters the temporal evolution of convection, lightning, and surface precipitation, leading to late afternoon thunderstorm development. The decrease in latent heat flux from the land surface caused by urbanization weakens convection in the early afternoon, while a condition suitable for convection development is maintained in the late afternoon due to less stabilization of the lower troposphere by the weaker convection development and high sensible heat flux from the surface. On the other hand, anthropogenic aerosols directly enhance convection, lightning, and surface precipitation by increasing convective updrafts due to the aerosol‐induced convective invigoration. The combined effects of urban land use and AE lead to even stronger thunderstorms in the late afternoon, mostly consistent with observations. These results indicate that urbanization increases the probability of late afternoon thunderstorms over the Houston area during the summer season. Advanced weather forecasting models that incorporate these urbanization effects might support sustainable urban planning to better mitigate the impacts of urbanization on local weather and public safety.

Iguchi, T. [Univ. of Maryland, College Park, MD (U↗
DOE OSTI2026

Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran↗
DOE OSTI2024

Simulation-based inference for parameter estimation of complex watershed simulators

High-resolution, spatially distributed process-based (PB) simulators are widely employed in the study of complex catchment processes and their responses to a changing climate. However, calibrating these PB simulators using observed data remains a significant challenge due to several persistent issues, including the following: (1) intractability stemming from the computational demands and complex responses of simulators, which renders infeasible calculation of the conditional probability of parameters and data, and (2) uncertainty stemming from the choice of simplified representations of complex natural hydrologic processes. Here, we demonstrate how simulation-based inference (SBI) can help address both of these challenges with respect to parameter estimation. SBI uses a learned mapping between the parameter space and observed data to estimate parameters for the generation of calibrated simulations. To demonstrate the potential of SBI in hydrologic modeling, we conduct a set of synthetic experiments to infer two common physical parameters – Manning's coefficient and hydraulic conductivity – using a representation of a snowmelt-dominated catchment in Colorado, USA. We introduce novel deep-learning (DL) components to the SBI approach, including an “emulator” as a surrogate for the PB simulator to rapidly explore parameter responses. We also employ a density-based neural network to represent the joint probability of parameters and data without strong assumptions about its functional form. While addressing intractability, we also show that, if the simulator does not represent the system under study well enough, SBI can yield unreliable parameter estimates. Approaches to adopting the SBI framework for cases in which multiple simulator(s) may be adequate are introduced using a performance-weighting approach. The synthetic experiments presented here test the performance of SBI, using the relationship between the surrogate and PB simulators as a proxy for the real case.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

MPACT Safeguards Modeling: FY25 Update

Sandia National Laboratories develops and maintains several open-source software packages to support material accountancy analyses. This includes the Material Accountancy Performance Indicator Toolkit (MAPIT), the Fissile Facility Flow Modeler (F3M) and the Separation and Safeguards Performance Model Library (SSPM-L). MAPIT is responsible for performing statistical safeguards analyses on bulk and itemized data from nuclear fuel cycle facilities and can operate on real or synthetic data. MAPIT is the only open-source software for such analyses. F3M is a library of modules, built in MATLAB Simulink, that contain pre made blocks to represent different generic fuel cycle processes. These blocks can be used together in a modular fashion to represent and simulate nuclear fuel cycle processes with the goal of improving facility-level accountancy during the design phase. F3M is also an open-source library. Finally, the SSPM-L library is a series of completed models built from F3M. The library includes facility models such as a generic PUREX facility and a fuel fabrication facility. The SSPM-L library is not open source, but is available to collaborators with a relevant use case. These tools include modeling and simulation pipelines to simulate nuclear fuel cycle facilities and the underlying software needed to simulate measurement uncertainty and perform statistical analyses. Together, these tools can perform end-to-end nuclear material accountancy analyses. This report documents the various improvements made to these tools in FY25. Specifically, we added new statistical test, new statistical modeling capabilities, new fuel cycle facility models, and launched a new open-source model component library.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Factors Controlling Precipitation Onset and Maintenance: Inferences From Large Eddy Simulations of Two Afternoon Deep Convective Regimes Over Amazon

To investigate the environmental factors controlling the onset and maintenance of afternoon precipitation over tropical rainforests, two contrasting cases are created for large eddy simulations (LES) using GoAmazon observations: one features a single-pulse rain dissipating quickly in earlier afternoon, while the other shows double pulses lasting until evening. For these two specific cases under consideration, LES confirms that early-morning relative humidity dominates afternoon rain patterns. The attributed impacts are distinct: moisture determines the rain onset timing while temperature affects the peak rain intensity. Single-pulse day observes one round of strong precipitation and cold pools, which further suppresses convection. On double-pulse day, the first precipitation peak results from an intermediate development of congestus clouds, whose detrainment leads to a gradual moistening of the lower-to-middle troposphere. This behavior of convection favors a second pulse of stronger precipitation with more convective organization, whose development decouples from surface fluxes and sustains until evening.

Amazon↗
DOE OSTI2026

T&D LIBRA: Simulation Framework for Transmission and Distribution Load and IBR Assessment

This paper introduces T&D LIBRA, a HELICS-based dynamic co-simulation framework designed for transmission and distribution (T&D) load and inverter-based resource assessment. This framework is designed to model interactions between single-phase induction motors (i.e., Motor-D), distributed energy resources (DERs) and the transmission network, emphasizing system stability during fault conditions. The primary objective of this research is to study the fault-induced delayed voltage recovery (FIDVR) behavior of the T&D networks using the co-simulation framework and to compare the results with conventional PSS/E only simulations. Simulations are validated using T&D co-simulations with a single distribution node and the IEEE 123-bus distribution test case. The co-simulation results demonstrate the value of using T&D co-simulations to capture motor stalling, DER tripping and other sources of voltage diversity on distribution networks.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Influence of Surfactants on Noble Metal Liquid-Gas Interface Mass Transfer in Molten Salt Reactors

Noble metals, a group of insoluble fission products (FPs), can accumulate in molten salt reactors (MSRs) and influence system performance and safety through deposition. Upon generation in fission, noble metals are transported by circulating salt and may either deposit on structural surfaces or enter the circulating gas bubble interfaces. While sparging was primarily designed for xenon and krypton removal, the mass transfer of noble metals into bubble interfaces highlights a potential pathway for their transport and partial extraction. The properties of the liquid-gas interface therefore play an important role in noble metal deposition and removal behaviors. Observations from earlier studies suggest that noble metals accumulated at the bubble interface can act as surfactants, making the bubble interface partially immobile. This behavior reduces the mass transfer capacity of circulating bubbles and affects both the overall distribution of noble metals in the loop and their removal efficiency. To investigate the influence of surfactants on noble metal transport, deposition, and extraction behaviors in MSRs, a previously developed and validated species transport model incorporating noble metal interphase mass transfer was used to simulate these behaviors. Two designed cases were studied to predict the distribution of 132Te along the Molten Salt Reactor Experiment (MSRE) loop under different bubble interface conditions. The results show that the mass transfer coefficients at the liquid-gas interface, determined by bubble interface characteristics, significantly influence noble metal distribution within the reactor loop. These findings emphasize the importance of continuously introducing fresh circulating bubbles into the molten salt to improve the removal efficiency of insoluble FPs from the MSR primary loop.

MSR↗
DOE OSTI2024

Development of an Optimal Variable-Pitch Controller for Floating Axial-Flow Marine Hydrokinetic Turbines

This article discusses the development of an optimal variable-pitch controller for floating, axial-flow marine turbines. Recently, OpenFAST, an open-source wind turbine modeling tool, has been extended to model marine turbines. A controller is necessary to simulate marine turbines for different load cases using OpenFAST, which greatly impacts the performance of the energy system. Previous studies have designed controllers using a linearized model of the marine turbine, which can be time-consuming and require the expertise of a control engineer. In this study, we use an automated approach that uses generic models of the marine turbine to identify the controller gains, which can expedite the process of designing a controller. Using an optimizer to identify the control system parameters can additionally improve the controller's performance. The optimal controller tuned using such an approach results in a 20% reduction in the tower-base damage equivalent loading and better tracking of the rated generator speed and power.

axial flow↗
DOE OSTI2024

Development of an Optimal Variable-Pitch Controller for Floating Axial-Flow Marine Hydrokinetic Turbines: Preprint

This article discusses the development of an optimal variable-pitch controller for floating, axial-flow marine turbines. Recently, OpenFAST, an open-source wind turbine modeling tool, has been extended to model marine turbines. A controller is necessary to simulate marine turbines for different load cases using OpenFAST, which greatly impacts the performance of the energy system. Previous studies have designed controllers using a linearized model of the marine turbine, which can be timeconsuming and require the expertise of a control engineer. In this study, we use an automated approach that uses generic models of the marine turbine to identify the controller gains, which can expedite the process of designing a controller. Using an optimizer to identify the control system parameters can additionally improve the controller's performance. The optimal controller tuned using such an approach results in a 20% reduction in the towerbase damage equivalent loading and better tracking of the rated generator speed and power.

closed-loop control↗
DOE OSTI2025

Data-model files associated with the manuscript "Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California"

This package contains the data, simulation setups, notebooks and figures used in “Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California” (Xu et al., 2025). In this study, we selected Elkhorn Slough, a tidal estuary, in California, to investigate the impact of wetland restoration and sea level rise on coastal hydrology using the process-based coastal hydrologic model, Advanced Terrestrial Simulator (ATS), informed by site-specific data. We designed a novel modeling workflow for incorporating wetland restoration features into land cover and soil properties for the model parameterization. The validation results demonstrate a strong agreement between modeled and observed data. We studied the characteristics of coastal watershed hydrology, then focused on the surface water dynamics at two wetland sites within Elkhorn Slough, a reference site and a restored site. Our simulation results indicate that the restored site successfully maintains surface elevation, resulting in reduced surface inundation. We also examined the impact of wetland restoration under expected sea level rise over the next few decades. The low-lying Yampah Marsh, the reference site, is likely to be inundated due to future sea level rise when highest tides arrive; while a higher percentage of Hester Marsh, the restored site, would retain marsh vegetation in coming decades, regardless of tidal conditions. Our study provides important information for examining the outcome of restoration practices that include surface elevation in tidal wetlands under climate changes.Several files can be found from this data package.1. README.md: This file describes the title, journal, co-authors, abstract, repository structure and model version.2. Simulation_Setups.zip: The file contains the model configuration files (XML format) for ATS. 3. Notebooks.zip: The file contains the Jupyter notebooks for generating the pre- and post-restoration meshes and the meshes of future scenarios. 4. Figures.zip: The file contains the figures used in the manuscript.5. Data.zip: The file contains the data used to drive the model simulations, including watershed and wetlands boundaries, mesh files and references to additional datasets (e.g., meteorological forcing, tidal dataset, DEMs, land cover, soil properties). Also, it contains water level observations at the restored wetland.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Enhancement of mid-/high-Z impurity transport by continuous Li-granule dropping in a stellarator plasma

An enhancement of core impurity transport is observed in high-density plasmas of the stellarator large helical device, heated by neutral beam injection, when continuous lithium (Li) granule dropping is performed. In the reported experiments, in which the tracer-encapsulated solid pellet is employed to inject trace amounts of titanium (Ti) and molybdenum (Mo) into the plasma core, confinement times for these impurities are seen to reduce significantly when Li dropping is applied, this reduction being more notable for Mo. To gain some initial insight into these observations, simulations are performed using the drift-kinetic transport code SFINCS for the Mo case. These simulations indicate that, while neoclassical transport prevails for the main plasma components (electrons, majority ions, and low-Z impurities), the classical contribution appears dominant for transporting Mo impurities. In summary, this work reports the first experimental observation of the degradation of mid-Z and high-Z impurity confinement induced by the continuous dropping of Li granules into a high-density stellarator plasma. In the case of the Mo impurity, simulations suggest that classical transport is the key mechanism underlying the enhanced impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Simulation budgeting for hybrid effective field theories

In this work, we forecast the number of, and requirements on, N-body simulations needed to train hybrid effective field theory (HEFT) emulators for a range of use cases, using a hybrid of HMcode and perturbation theory as a surrogate model. Our accuracy goals, determined with careful consideration of statistical and systematic uncertainties, are 1% accurate in the high-likelihood range of cosmological parameters, and 2% accurate over a broader parameter space volume for k < 1 h Mpc -1 and z < 3. Focusing in part on the 8-parameter w 0 w a CDM+ m ν cosmological model, we find that < 225 simulations are required to meet our error goals over our wide parameter space, including models with rapidly evolving dark energy, given our simulation and emulator recommendations. For a more restricted parameter space volume, as few as 80 simulations are sufficient. We additionally present simulation forecasts for example use cases, and make the code used in our analyses publicly available. These results offer practical guidance for efficient emulator design and simulation budgeting in future cosmological analyses.

cosmological parameters from LSS↗
DOE OSTI2025

A Digital Twin of Scalable Quantum Clouds

Quantum computing has emerged as a transformative technology capable of solving complex problems beyond the limit of classical systems. The rapid development of quantum processors has led to the proliferation of cloud-based quantum computing services offered by platforms such as IBM, Google, and Amazon. These platforms introduce unique challenges in resource allocation, job scheduling, and multi-device orchestration as quantum workloads become increasingly complex. In this work, we present a digital twin of quantum cloud infrastructures: a framework designed to model and simulate the behavior of real quantum cloud systems. Developed in Python using the SimPy discrete-event simulation library, the framework replicates key aspects of quantum cloud environments, including detailed quantum device modeling, job lifecycle management, and job fidelity. It incorporates noise-aware fidelity estimation, making it the first of its kind to simulate superconducting gate-based quantum cloud systems at an administrative level with job fidelity. We present use cases as proof of concept, demonstrating that our quantum cloud simulation framework can act as a digital twin of a quantum cloud and support the modeling and implementation of practical systems.

Luo, Waylon [Kent State University]↗
DOE OSTI2025

Hybrid Simulations of FRC Merging and Compression

An improved understanding of field-reversed configuration (FRC) merging and stability in high acceleration and compression magnetic fields is needed to speed up the development of the pulsed fusion concept developed at Helion Energy. All previous theoretical and simulation work on FRC merging and compression was performed using two-dimensional (2D) magnetohydrodynamic (MHD) models. The results of novel 2D hybrid simulations (fluid electrons and full-orbit kinetic ions) of FRC merging and compression are presented. Results of kinetic and MHD simulations, computed using the HYM code, are compared and analyzed. In cases without axial magnetic compression, both the MHD and hybrid simulations show a high sensitivity to the initial parameters (i.e. FRC separation, velocity, normalized separatrix radius, and plasma viscosity), showing that FRCs with large elongation and separatrix radius either do not merge or merge partially, forming a doublet FRC. In conclusion, application of a mirror coil field at the FRC ends with increasing strength is shown to lead to fast and complete merging of the FRCs in MHD and kinetic simulations.

FRC↗