Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation and modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Sequential Bayesian Experimental Design for Calibration of Expensive Simulation Models

Simulation models of critical systems often have parameters that need to be calibrated using observed data. For expensive simulation models, calibration is done using an emulator of the simulation model built on simulation output at different parameter settings. Using intelligent and adaptive selection of parameters to build the emulator can drastically improve the efficiency of the calibration process. The article proposes a sequential framework with a novel criterion for parameter selection that targets learning the posterior density of the parameters. The emergent behavior from this criterion is that exploration happens by selecting parameters in uncertain posterior regions while simultaneously exploitation happens by selecting parameters in regions of high posterior density. Furthermore, the advantages of the proposed method are illustrated using several simulation experiments and a nuclear physics reaction model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Biomass burning emission analysis based on MODIS aerosol optical depth and AeroCom multi-model simulations: implications for model constraints and emission inventories

We assessed the biomass burning (BB) smoke aerosol optical depth (AOD) simulations of 11 global models that participated in the AeroCom phase III BB emission experiment. By comparing multi-model simulations and satellite observations in the vicinity of fires over 13 regions globally, we (1) assess model-simulated BB AOD performance as an indication of smoke source–strength, (2) identify regions where the common emission dataset used by the models might underestimate or overestimate smoke sources, and (3) assess model diversity and identify underlying causes as much as possible. Using satellite-derived AOD snapshots to constrain source strength works best where BB smoke from active sources dominates background non-BB aerosol, such as in boreal forest regions and over South America and southern hemispheric Africa. The comparison is inconclusive where the total AOD is low, as in many agricultural burning areas, and where the background is high, such as parts of India and China. Many inter-model BB AOD differences can be traced to differences in values for the mass ratio of organic aerosol to organic carbon, the BB aerosol mass extinction efficiency, and the aerosol loss rate from each model. The results point to a need for increased numbers of available BB cases for study in some regions and especially to a need for more extensive regional-to-global-scale measurements of aerosol loss rates and of detailed particle microphysical and optical properties; this would both better constrain models and help distinguish BB from other aerosol types in satellite retrievals. More generally, there is the need for additional efforts at constraining aerosol source strength and other model attributes with multi-platform observations.

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

A hybrid calorimetry-simulation model of mixing enthalpy for molten salt

Calorimetric determination of enthalpies of mixing (ΔH mix ) in multicomponent molten salts is often interpreted using empirical models that lack physically meaningful parameters. However, for improving pyrochemical separation of spent nuclear fuel, where lanthanides are major fission products and critical elements, a deeper thermodynamic understanding of the link between excess thermodynamic properties and solvation structure is critically needed. In this work, we implement a hybrid and physics-informed framework, MIVM+Calorimetry+AIMD, which integrates experimentally measured ΔH mix (via high temperature drop calorimetry) with solvation structures from ab initio molecular dynamics (AIMD). This approach is demonstrated using LaCl 3 mixed with eutectic LiCl-KCl (58 mol% – 42 mol%) at 873 K and 1133 K. MIVM-derived parameters enable extrapolation of excess Gibbs energy and La 3+ activity across compositions. In contrast, direct ΔH mix predictions from AIMD and polarizable ion model simulations deviate significantly. By incorporating experimentally benchmarked solvation structures into an interpretable thermodynamic model, the MIVM+Calorimetry+AIMD formalism achieves higher accuracy and generalizable method for studying molten salts, offering a robust path for understanding and optimizing molten salt chemistry relevant to nuclear fuel cycles and separation science.

Goncharov, Vitaliy G. [Washington State Univ., Pul↗
DOE OSTI2024

Observationally constrained analysis of sulfur cycle in the marine atmosphere with NASA ATom measurements and AeroCom model simulations

The atmospheric sulfur cycle plays a key role in air quality, climate, and ecosystems, such as pollution, radiative forcing, new particle formation, and acid rain. In this study, we compare the spatially and temporally resolved measurements from the NASA Atmospheric Tomography (ATom) mission with simulations from five AeroCom III models for four sulfur species (dimethyl sulfide (DMS), sulfur dioxide (SO 2 ), particulate methanesulfonate (MSA), and particulate sulfate (SO 4 )). We focus on remote regions over the Pacific, Atlantic, and Southern oceans from near the surface to ~12 km altitude range covering all four seasons. In general, the differences among model results can be greater than 1 order of magnitude. Comparing with observations, model-simulated SO 2 is generally low, whereas SO 4 is generally high. Simulated DMS concentrations near the sea surface exceed observed levels by a factor of 5 in most cases, suggesting potential overestimation of DMS emissions in all models. With GEOS model simulations of tagging emission from anthropogenic, biomass burning, volcanic, and oceanic sources, we find that anthropogenic emissions are the dominant source of sulfate aerosol (40 %–60 % of the total amount) in the ATom measurements at almost all altitudes, followed by volcanic emissions (18 %–32 %) and oceanic sources (16 %–32 %). Similar source contributions can also be derived at broad ocean basins and on monthly scales, indicating the representativeness of ATom measurements for global ocean. Our work presents the first assessment of AeroCom sulfur study using ATom measurements, providing directions for improving sulfate simulations, which remain the largest uncertainty in radiative forcing estimates in aerosol climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Model simulations of Plum Island Ecosystems LTER low marsh site using ELM-PFLOTRAN

Model simulations using the E3SM Land Model (ELM) coupled to the PFLOTRAN reactive transport model via the Alquimia interface. The simulations were conducted for a tidal salt marsh at the Plum Island Ecosystems LTER near Rowley, Massachusetts, USA. Model simulations were forced using site-specific tidal cycles and salinity, and the simulations used a biogeochemical reaction network including aerobic decomposition, sulfate reduction, iron reduction, and methanogenesis. Model outputs include simulated carbon stocks, carbon dioxide and methane fluxes, and porewater concentrations of key solutes related to sulfur, iron, and carbon cycling. The model simulations included a saline simulation (with tidal sulfate inputs), a fresh simulation (with low salinity and low sulfate inputs), and a saline simulation with lower vegetation productivity to represent the effect of salinity on vegetation. These simulations were conducted to demonstrate that a new model framework incorporating subsurface redox and biogeochemical interactions into a land surface model could reproduce measured surface greenhouse gas fluxes and biogeochemical dynamics in tidal marsh ecosystems, and to test whether including redox interactions in a land surface model would allow the model to resolve contrasts in biogeochemical cycling and greenhouse gas production between saline and freshwater wetlands.The data package includes gzipped tar archives (which can be expanded using standard tar and gzip utilities) of model outputs from three model configurations: saline subsurface and reduced vegetation productivity related to salinity; saline subsurface with vegetation productivity not reduced; and freshwater. Also included are code for the modified E3SM model, Alquimia interface, and PFLOTRAN reactive transport simulator in gzipped tar format; plain text parameter and configuration files; python code files for visualizing model output and defining model configurations; and model output, tide and salinity forcing, and configuration files in netCDF format. See the README.md file in the data package for a detailed description of all files contained in the package. All files are in netCDF (.nc), gzipped tar archive (.tar.gz or .tgz), or text (all other files).Updated: May 13, 2024. Model output, E3SM code, PFLOTRAN input files, and python codes for visualizing results were updated to reflect changes made for the manuscript revision. The updated archive reflects the code and model output from the final accepted manuscript. Changes included updated reaction parameters reflecting improved parameterization and additional comparisons with field measurements. E3SM code changes included better support for multiple grid cells and improved flow and transport parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

ELM model simulations of Plum Island Ecosystems LTER low marsh site 2018-2020

Model simulations using the Department of Energy's Energy Exascale Earth System Model (E3SM) land model (ELM) with improved capabilities to represent vegetation response to salinity and inundation. The simulations were conducted for a tidal salt marsh at Plum Island Ecosystems Long Term Ecological Research (LTER) site near Rowley, Massachusetts, USA; the site is a low marsh dominated by Spartina alterniflora. The model was forced with site-specific meteorology, salinity and tidal cycles from 2018-2020. Four sets of model simulations are included and described below:1. Parameterization of the salinity response function. These simulations tested different combinations of values for optimal salinity and salinity tolerance.2. Model evaluation. This comparison conducted simulations using the default model, the salinity function only, the submergence function only, and both the salinity and submergence functions. 3. Salinity scenarios. These simulations used the 2018 salinity input data varied by -5 to +10 ppt salinity.4. Water level scenarios. These simulations used the tide height varied by -10 to +50 cm. These simulations were used to demonstrate how the salinity and submergence functions better represent carbon uptake by tidal salt marshes.The data package includes netCDF files used as forcing files for tide height and salinity, one for each year 2018-2020 at observed salinity concentrations, and an additional three forcing files in which salinity concentrations were varied 5 ppt lower, 5 ppt higher, and 10 ppt higher than the measured 2018 time series. Also included are python scripts for creating forcing files, plain text parameter and command files for running simulations, model outputs in netCDF format, and python scripts for visualizing outputs. Code for the modified E3SM model is archived in Sulman et al 2023 at doi:10.15485/1991625. More detail about files is provided in the README.md file.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Co-Simulation Model for Optimal Wind-Hydro Coordination Using Wind Farm Control Dynamics

The growing share of Variable Renewable Energy sources (VRES) in power systems presents challenges for regula- tors, grid operators and energy producers. The VRES’ operation has limited flexibility in their operations, as they are highly dependent on ambient environments. To address these challenges, decision-makers must consider multiple objectives, among these are revenue, power system services and mechanical load on wind turbines. Coordinated operation of power plants and different wind farm control strategies are examples of measures that can benefit these objectives. This study proposes a Multi-Objective Linear Programming (MOLP) model to simulate the optimum operation of wind and hydropower plants that share limited transmission capacity. Further, wind farm control dynamics are included to obtain realistic output power and accumulated damage. From this, a case study based on a relevant location in Norway is presented to analyze the improved effect of wind- hydro coordination and wind farm control in achieving the objectives of accumulated wind turbine damage and total revenue of the hybrid power system. In addition, the study considers the potential advantages of adding a variable-speed pump to the hydropower plant. The results demonstrate that by considering multiple objectives in the optimization, one may achieve better overall performance of the objectives. By utilizing the flexibility of hydro storage, the decision maker may adjust to obtain the most desired outcome. Moreover, the added flexibility of utilizing a pump for hydro storage shows great improvements for the combined revenue of the power plants and reduced curtailment of the wind farm. However, less impact is observed from using a variable speed pump compared to a fixed speed pump.

13 HYDRO ENERGY↗
DOE OSTI2024

Simulations, Modeling and Data Analysis of Parity Violating Electron Scattering Experiments

In the Standard Model (SM) of nuclear and particle physics, parity violation is incorporated through the representation of the weak interaction as a chiral gauge interaction. Only the left-handed components of particles and right-handed components of antiparticles participate in weak interactions in the Standard Model. This implies that parity is asymmetric for the weak interaction. Parity violating electron scattering (PVES) experiments are designed to probe the physics parameters related to the SM, with the possibility to discover physics beyond the SM (BSM) by measuring the parity violating asymmetry ¿¿¿ of longitudinally polarized electrons scattered off unpolarized targets with high precision. This dissertation will be focused on two PVES experiments, the next 208Pb Lead Radius Experiment (PREX-II), and the Measurement of a Lepton-Lepton Electroweak Reaction (MOLLER) experiment, as well as in some small sections, the Calcium Radius Experiment (CREX) and P2 experiment which are also PVES experiments). PREX-II and CREX experiments, performed in Hall A at the Thomas Jefferson National Accelerator Facility (Jefferson Lab), measured ¿¿¿ in the elastic scattering of longitudinally polarized electrons from 208Pb and 48Ca targets to provide a precise model independent determination of the neutron skin thickness of 208Pb and 48Ca nuclei, respectively. The MOLLER experiment, proposed to start in 2027 and also to be performed in Hall A at Jefferson Lab, is to measure ¿¿¿ of longitudinally polarized electrons scattered off unpolarized electrons (Møller scattering) to determine the weak charge of electrons ¿¿¿ and the weak mixing angle ¿¿ with high precision. As for the P2 experiment, which will be performed at the upcoming MESA accelerator in Mainz Germany, it is to measure the weak charge of proton ¿¿¿ using ¿¿¿ in the elastic electron-proton scattering of polarized electrons off unpolarized protons. The final results from the PREX-II experiment are presented as ¿¿¿=550±16 (¿¿¿¿)±8 (¿¿¿¿) parts-per-billion (ppb). Combining the PREX-I and PREX-II results, the neutron skin thickness from PREX experiments is determined as ¿¿-¿¿=0.283±0.071 ¿¿ in 208Pb. This thesis lists the software and computational contribution of the author to these PVES experiments, including writing scripts and software to help with the PREX-II/CREX experiments, analyzing data to provide useful information and systematic uncertainty for the PREX-II experiment, modeling and simulations for the MOLLER, and providing an alternative design of an electronic equipment for the P2 experiment.

Chen, Yufan↗
DOE OSTI2024

Simulations, Modeling and Data Analysis of Parity Violating Electron Scattering Experiments

In the Standard Model (SM) of nuclear and particle physics, parity violation is incorporated through the representation of the weak interaction as a chiral gauge interaction. Only the left-handed components of particles and right-handed components of antiparticles participate in weak interactions in the Standard Model. This implies that parity is asymmetric for the weak interaction. Parity violating electron scattering (PVES) experiments are designed to probe the physics parameters related to the SM, with the possibility to discover physics beyond the SM (BSM) by measuring the parity violating asymmetry ¿¿¿ of longitudinally polarized electrons scattered off unpolarized targets with high precision. This dissertation will be focused on two PVES experiments, the next 208Pb Lead Radius Experiment (PREX-II), and the Measurement of a Lepton-Lepton Electroweak Reaction (MOLLER) experiment, as well as in some small sections, the Calcium Radius Experiment (CREX) and P2 experiment which are also PVES experiments). PREX-II and CREX experiments, performed in Hall A at the Thomas Jefferson National Accelerator Facility (Jefferson Lab), measured ¿¿¿ in the elastic scattering of longitudinally polarized electrons from 208Pb and 48Ca targets to provide a precise model independent determination of the neutron skin thickness of 208Pb and 48Ca nuclei, respectively. The MOLLER experiment, proposed to start in 2027 and also to be performed in Hall A at Jefferson Lab, is to measure ¿¿¿ of longitudinally polarized electrons scattered off unpolarized electrons (Møller scattering) to determine the weak charge of electrons ¿¿¿ and the weak mixing angle ¿¿ with high precision. As for the P2 experiment, which will be performed at the upcoming MESA accelerator in Mainz Germany, it is to measure the weak charge of proton ¿¿¿ using ¿¿¿ in the elastic electron-proton scattering of polarized electrons off unpolarized protons. The final results from the PREX-II experiment are presented as ¿¿¿=550±16 (¿¿¿¿)±8 (¿¿¿¿) parts-per-billion (ppb). Combining the PREX-I and PREX-II results, the neutron skin thickness from PREX experiments is determined as ¿¿-¿¿=0.283±0.071 ¿¿ in 208Pb. This thesis lists the software and computational contribution of the author to these PVES experiments, including writing scripts and software to help with the PREX-II/CREX experiments, analyzing data to provide useful information and systematic uncertainty for the PREX-II experiment, modeling and simulations for the MOLLER, and providing an alternative design of an electronic equipment for the P2 experiment.

Chen, Yufan↗
DOE OSTI2024

How Well Does the DOE Global Storm Resolving Model Simulate Clouds and Precipitation Over the Amazon?

This study assesses a 40-day 3.25-km global simulation of the Simple Cloud-Resolving E3SM Model (SCREAMv0) using high-resolution ground-based observations from the Atmospheric Radiation Measurement (ARM) Green Ocean Amazon (GoAmazon) field campaign. SCREAMv0 reasonably captures the diurnal timing of boundary layer clouds yet underestimates the boundary layer cloud fraction and mid-level congestus. SCREAMv0 well replicates the precipitation diurnal cycle, however it exhibits biases in the precipitation cluster size distribution compared to scanning radar observations. Specifically, SCREAMv0 overproduces clusters smaller than 128 km, and does not form enough large clusters. Such biases suggest an inhibition of convective upscale growth, preventing isolated deep convective clusters from evolving into larger mesoscale systems. This model bias is partially attributed to the misrepresentation of land-atmosphere coupling. This study highlights the potential use of high-resolution ground-based observations to diagnose convective processes in global storm resolving model simulations, identify key model deficiencies, and guide future process-oriented model sensitivity tests and detailed analyses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Scientist-in-the-Loop Data Analytics Framework for Intelligent Simulation Model Tuning and Validation

This project developed a scientist-in-the-loop data analytics framework for intelligent simulation model tuning and validation, targeting the Weather Research and Forecasting (WRF) model and its solar energy variant, WRF-Solar-BNL. Domain experts, such as climate scientists, depend on large-scale numerical simulations for knowledge discovery and decision-making, yet the complexity of parameter tuning and the disconnect between automated optimization and domain expertise pose significant challenges. We extended an interactive visual analytics framework that enables domain experts to observe and intervene in the computational steering process by identifying disagreements between the simulation model, surrogate model, and the expert’s domain knowledge. Using Bayesian Optimization with Gaussian Process Regression as the surrogate model, our system allows users to probe parameter relationships, analyze correlation patterns, and adjust tuning parameters in real time. We developed use cases for solar irradiance forecasting through sustained collaboration with Brookhaven National Laboratory, resolving critical model configuration challenges and achieving meaningful reductions in prediction error. The project supported one PhD student, one MS student, and eight undergraduate students across three Data Science Capstone projects, resulting in one master’s thesis.

Dasgupta, Aritra [New Jersey Institute of Technolo↗
DOE OSTI2026

A comprehensive framework for validating simulation models of power system equipment using terminal measurements

Accurate simulation of power-plants is essential to the planning and operation of modern power grids. The current methods used to periodically check power-plant simulation models have many open questions about their limitations and accuracy. The research in this project explored using Monte-Carlo Experimentation (MCE) as a means for answering these important questions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Dynamics of streamflow permanence in a headwater network: Insights from catchment-scale model simulations

The hillslope and channel dynamics that govern streamflow permanence in headwater systems have important implications for ecosystem functioning and downstream water quality. Recent advancements in process-based, semi-distributed hydrologic models that build upon empirical studies of streamflow permanence in well-monitored headwater catchments show promise for characterizing the dynamics of streamflow permanence in headwater systems. However, few process-based models consider the continuum of hillslope-stream network connectivity as a control on streamflow permanence in headwater systems. The objective of this study was to expand a process-based, catchment-scale hydrologic model to better understand the spatiotemporal dynamics of headwater streamflow permanence and to identify controls of streamflow expansion and contraction in a headwater network. Further, we aimed to develop an approach that enhanced the fidelity of model simulations, yet required little additional data, with the intent that the model might be later transferred to catchments with limited long-term and spatially explicit measurements. This approach facilitated network-scale estimates of the controls of streamflow expansion and contraction, albeit with higher degrees of uncertainty in individual reaches due to data constraints. Our model simulated that streamflow permanence was highly dynamic in first-order reaches with steep slopes and variable contributing areas. The simulated stream network length ranged from nearly 98±2% of the geomorphic channel extent during wet periods to nearly 50±10% during dry periods. The model identified a discharge threshold of approximately 1 mm d –1 , above which the rate of streamflow expansion decreases by nearly an order of magnitude, indicating a lack of sensitivity of streamflow expansion to hydrologic forcing during high-flow periods. Overall, we demonstrate that process-based, catchment-scale models offer important insights on the controls of streamflow permanence, despite uncertainties and limitations of the model. We encourage researchers to increase data collection efforts and develop benchmarks to better evaluate such models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Fidelity Analysis of Complex Large-Scale Simulation Models of PEs (MARS) in Future Power Grids

With the proliferation of power electronic (PE)-based resources in the electric grid, there have been numerous different kinds of simulation models developed and studied to understand their dynamic behavior. Study of PE-based resources using each simulation model has its own distinct advantages and disadvantages. Therefore it is very important to carefully select the appropriate model for the evaluation of the PE-based resources. In this paper, one such effort is made to analyze the different types of models in the case of a multi-port autonomous reconfigurable solar power plant (MARS). The MARS integrates a hybrid photovoltaic (PV)-energy storage system (ESS) plant to high-voltage direct current (HVdc) system and transmission alternating current (ac) grid. The dynamic models are evaluated under different test cases varying the short-circuit ratio (SCR) and control algorithms. Recommendations are provided on the equipment study needs based on the analysis completed.

Marthi, Phani Ratna Vanamali↗
DOE OSTI2023

Tropical Tree Crop Simulation with a Process-Based, Daily Timestep Simulation Model (ALMANAC): Description of Model Adaptation and Examples with Coffee and Cocoa Simulations

Coffee (Coffea species) and Cocoa (Theobroma cacao) are important cash crops grown in the tropics but traded globally. This study was conducted to apply the ALMANAC model to these crops for the first time, and to test its ability to simulate them under agroforestry management schemes and varying precipitation amounts. To create this simulation, coffee was grown on a site in Kaua’i, Hawai’i, USA, and cocoa was grown on a site in Sefwi Bekwai, Ghana. A stand-in for a tropical overstory tree was created for agroforestry simulations using altered parameters for carob, a common taller tropical tree for these regions. For both crops, ALMANAC was able to realistically simulate yields when compared to the collected total yield data. On Kaua’i, the mean simulated yield was 2% different from the mean measured yield, and in all three years, the simulated values were within 10% of the measured values. For cocoa, the mean simulated yield was 3% different from the mean measured yield and the simulated yield was within 10% of measured yields for all four available years. When precipitation patterns were altered, in Ghana, the wetter site showed lower percent changes in yield than the drier site in Hawai’i. When agroforestry-style management was simulated, a low Leaf Area Index (LAI) of the overstory showed positive or no effect on yields, but when LAI climbed too high, the simulation was able to show the detrimental effect this competition had on crop yields. These simulation results are supported by other literature documenting the effects of agroforestry on tropical crops. This research has applied ALMANAC to new crops and demonstrated its simulation of different management and environmental conditions. The results show promise for ALMANAC’s applicability to these scenarios as well as its potential to be further tested and utilized in new circumstances.

60 APPLIED LIFE SCIENCES↗