Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Modeling & 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 55 records · Page 3

Plateau to River Model Predictive Simulations for All Ensemble Realizations to Support Modeling Work in Fiscal Year 2025

The purpose of this environmental calculation file (ECF) is to document predictions of flow and hydraulic head on the Central Plateau of the Hanford Site using the Plateau-to-River (P2R) Model (CP-57037, Model Package Report for the Plateau-to-River Model: Version 9.1). This calculation documents the simulation of the groundwater for the parent model domain of the P2R Model as a basis for use in other applications of the P2R Model. This application is unique from the standpoint that it will simulate all ensemble member models of the P2R Model whereas other applications may only utilize specific ensemble members. These simulations provide results that can be used in the process of selecting an appropriate subset of ensemble members for other applications.

54 ENVIRONMENTAL SCIENCES

Opportunities and challenges in process modeling and simulation of electrochemical systems

Electrochemical technologies have garnered intense interest in both academic research and industrial applications due to their potential to increase energy efficiency and reduce carbon footprint. However, electrochemical process fundamentals have been absent from the chemical process simulators available to millions of chemical engineers worldwide. To expedite process research and development of electrochemical technologies in the chemical industry, it is imperative that essential electrochemical process fundamentals be incorporated into process simulators to support modeling and simulation of electrochemical processes. Here, this study examines three process fundamentals key to the research and development of electrochemical processes: electrochemical reaction kinetics, electrolyte thermodynamics, and heat and mass transfer. It further illustrates application of these process fundamentals with a case study modeling an electrochemical process for the conversion of acrylonitrile to adiponitrile. The modeling example highlights the roles of applied voltage on reaction rates, and electron flow rates on the performance of electrochemical processes. It suggests that the applied voltage and electron flow rates are two unique concepts that should be included in simulators for electrochemical processes.

09 BIOMASS FUELS

Comparing computational times for simulations when using PBPK model template and stand-alone implementations of PBPK models

Introduction We previously developed a PBPK model template that consists of a single model “superstructure” with equations and logic found in many physiologically based pharmacokinetic (PBPK) models. Using the template, one can implement PBPK models with different combinations of structures and features. Methods To identify factors that influence computational time required for PBPK model simulations, we conducted timing experiments using various implementations of PBPK models for dichloromethane and chloroform, including template and stand-alone implementations, and simulating four different exposure scenarios. For each experiment, we measured the required computational time and evaluated the impacts of including various model features (e.g., number of output variables calculated) and incorporating various design choices (e.g., different methods for estimating blood concentrations). Results We observed that model implementations that treat body weight and dependent quantities as constant (fixed) parameters can result in a 30% time savings compared with options that treat body weight and dependent quantities as time-varying. We also observed that decreasing the number of state variables by 36% in our PBPK model template led to a decrease of 20–35% in computational time. Other factors, such as the number of output variables, the method for implementing conditional statements, and the method for estimating blood concentrations, did not have large impacts on simulation time. In general, simulations with PBPK model template implementations of models required more time than simulations with stand-alone implementations, but the flexibility and (human) time savings in preparing and reviewing a model implemented using the PBPK model template may justify the increases in computational time requirements. Conclusion Our findings concerning how PBPK model design and implementation decisions impact computational speed can benefit anyone seeking to develop, improve, or apply a PBPK model, with or without the PBPK model template.

Bernstein, Amanda S.

Trusted Simulation: Considering Model Quality in the Context of User Trust

A high‐quality simulation model should help its users to easily and appropriately calibrate their trust in the model. Traditional evaluation metrics such as validation and robustness are necessary but insufficient for this task. Trust calibration depends on factors like the model's transparency, applicability to intended use, usability, reputation, and consideration of potential bias. This article proposes a framework for designing and evaluating system dynamics models by considering factors that contribute to the proper calibration of user trust. This framework takes inspiration from trusted artificial intelligence, broadening our traditional concept of model quality and explicitly focusing on what users need to consider a model trustworthy and to understand the model's relevance to its intended purpose. The trusted simulation framework can improve our integration of model quality activities throughout the modeling process, leading to more impactful and better‐targeted model design, development, and evaluation.

Naugle, Asmeret Bier [Sandia National Laboratories

A hybrid Penman-Monteith and machine learning model for simulating evapotranspiration and its components

Integrating physical processes with machine learning has advanced evapotranspiration (ET) simulation, yet most hybrid models fail to partition total ET into its components: soil evaporation (E) and vegetation transpiration (T). This study introduces Residual Neural Network–Penman–Monteith (RNN-PM), a novel hybrid dual-source ET model designed to overcome this limitation. The model synergizes the physically-based Penman–Monteith framework with three specialized residual neural networks trained to estimate key conductance parameters (canopy conductance, soil surface conductance, and aerodynamic conductance). Furthermore this explicit parameterization allows for the direct partitioning of total ET. Validation at National Ecological Observatory Network (NEON) flux sites using high-frequency partitioned E and T shows that RNN-PM reliably reproduces ET and the transpiration fraction (T/ET). For ET, the model achieves an average Kling–Gupta efficiency (KGE) of 0.89 and a root-mean-square error (RMSE) of 0.55 mm/day; for T/ET, the KGE is 0.87 with an RMSE of 0.06. Furthermore, RNN-PM demonstrates robust generalization, accurately simulating ET and its components well beyond the initial training dataset, even under extreme climatic conditions. This study extended the analysis by comparing the RNN-PM model with seven established dual-source ET models. The results indicate that RNN-PM outperforms both conventional machine learning models and purely physical process-based models in simulating ET components in most cases. Among the purely physical process-based dual-source models, those based on surface temperature decomposition showed improved performance as the leaf area index (LAI) decreased when evaluated against high-frequency ET component datasets. In contrast, the performance of conductance-based dual-source models declined with decreasing LAI. Although purely machine learning-based models can produce relatively accurate simulations of ET components, they often exhibit limited generalization capability, an issue that the RNN-PM model effectively overcomes. Ultimately, the RNN-PM model represents a significant advance in simulating ET components, offering a novel and scalable approach for improving the representation of land–atmosphere interactions in Earth system models.

54 ENVIRONMENTAL SCIENCES

Comparing multi-model ensemble simulations with observations and decadal projections of upper atmospheric variations following the Hunga eruption

The Hunga Tonga–Hunga Ha'apai Model–Observation Comparison (HTHH–MOC) project aims to comprehensively investigate the evolution of volcanic water vapor and sulfur emissions and their subsequent atmospheric impacts and underlying response mechanisms using state-of-the-art global climate models. This study evaluates multi-model ensemble simulations participating in the HTHH–MOC free-run experiment with climate projections for 10 years (2022–2032). Model results are evaluated against satellite observations to assess their ability to reproduce the observed evolution of stratospheric water vapor, aerosols, temperature, and ozone from 2022 to 2024. The participating models accurately capture the observed distribution patterns and associated upper atmospheric responses, providing confidence for their future projections. Model simulations suggest that the Hunga eruption-induced stratospheric water vapor anomaly lasts 4–7 years, with a water vapor e-folding time of 31–43 months. This prolonged water vapor perturbation leads to significant stratospheric and mesospheric cooling, resulting in significant ozone loss in the upper stratosphere and lower mesosphere for 7–10 years. Comparisons between simulations with both SO 2 and H 2 O emissions and those with H 2 O-only emissions indicate that the pronounced dipole response with upper-stratospheric cooling and lower-stratospheric warming is driven by the combined effects of SO 2 and H 2 O injections. These results highlight the prolonged atmospheric impacts of the Hunga eruption and the potential critical role of stratospheric water vapor in modulating long-term atmospheric chemistry and dynamics.

Zhuo, Zhihong [Univ. of Quebec, Montreal, QC (Cana

Upscaling Methods Applied to a Fine-Scale Reservoir Model

This study was conducted as part of the Southwest Regional Partnership on Carbon Sequestration (SWP) project to evaluate how upscaling fine-scale simulation models to coarse-scale simulation models impacted the results. The focus was on the Farnsworth Unit (FWU) and its Morrow' B' Sandstone reservoir, specifically the west half of the field. Due to data limitations and the geologic characteristics of the surrounding area, the upscaling was limited to the west half of the FWU rather than a broader basinscale model. The primary aim was to explore how upscaling impacts numerical simulation models, particularly regarding CO 2 -enhanced oil recovery (EOR) and storage capacity predictions. Upscaling was necessary to reduce computational demands when transitioning from high-resolution geological models to coarser grids, as large-scale simulations with finer grids can be computationally prohibitive. This study expands on previous work by the SWP to understand how additional upscaling, applied to already fine-scale numerical simulation models, affects reservoir performance simulations (Ampomah, Balch, & Grigg, 2015). This is key to understanding how loss of resolution can affect coarsescale model results that may be used for large sensitivity analyses, uncertainty quantifications, and training data for machine learning applications.

02 PETROLEUM

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING

Assessment of Physical Security Modeling and Simulation in the Vulnerability Assessment Process

This report provides a comprehensive assessment of physical security modeling and simulation tools available for use in the vulnerability assessment (VA) process for nuclear facilities. It outlines the historical evolution of VA methodologies, emphasizing the transition from traditional layer-based approaches to a more holistic framework that integrates detection probabilities directly into combat simulations. The document details the critical components of the VA process, including the characterization of targets, threats, and protective measures, as well as the development of adversary scenarios that reflect both insider and outsider threats. It highlights the importance of performance assurance programs, emphasizing the need for continuous evaluation and testing of security systems to ensure their effectiveness against evolving threats. Additionally, the report discusses the significance of utilizing accredited modeling and simulation tools in accredited areas to accurately represent adversary actions and the corresponding responses of protective forces. By establishing a systematic approach to VA, this document aims to enhance the overall security posture of nuclear facilities, ensuring compliance with regulatory standards while effectively mitigating risks associated with potential adversarial actions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF

Nuclear Modeling and Simulation Using MOOSE-based Tools

Overview presentation of the variety of modeling and simulation activities for nuclear applications to be presented to visitors (professors and their students) from the United Kingdom. Presentation compiled by listed author with contributions from others at the laboratory. Others are listed within the acknowledgements of the presentation rather than a long list of authors on the title slide.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS

A protocol for model intercomparison of impacts of marine cloud brightening climate intervention

A modeling protocol (defined by a series of climate model simulations with specified model output) is introduced. Studies using these simulations are designed to improve the understanding of climate impacts using a strategy for climate intervention (CI) known as marine cloud brightening (MCB) in specific regions; therefore, the protocol is called MCB-REG (where REG stands for region). The model simulations are not intended to assess consequences of a realistic MCB deployment intended to achieve specific climate targets but instead to expose responses to interventions in six regions with pervasive cloud systems that are often considered candidates for such a deployment. A calibration step involving simulations with fixed sea surface temperatures (SSTs) is first used to identify a common forcing, and then coupled simulations with forcing in individual regions and combinations of regions are used to examine climate impacts. Synthetic estimates constructed by superposing responses from simulations with forcing in individual regions are considered a means of approximating the climate impacts produced when MCB interventions are introduced in multiple regions. A few results comparing simulations from three modern climate models (CESM2, E3SMv2, and UKESM1) are used to illustrate the similarities and differences between model behavior and the utility of estimates of MCB climate responses that were synthesized by summing responses introduced in individual regions. Cloud responses to aerosol injections differ substantially between models (CESM2 clouds appear much more susceptible to aerosol emissions than the other models), but patterns in precipitation and surface temperature responses were similar when forcing is imposed with similar amplitudes in the same regions. A previously identified La Niña-like response to forcing introduced in the Southeast Pacific is evident in this study, but the amplitude of the response was shown to markedly differ across the three models. Other common response patterns were also found and are discussed. Forcing in the Southeast Atlantic consistently (across all three models) produces weaker global cooling than that in other regions, and the Southeast Pacific and South Pacific show the strongest cooling. This indicates that the efficiency of a given intervention depends on not only the susceptibility of the clouds to aerosol perturbations, but also the strength of the underlying radiative feedbacks and ocean responses operating within each region. These responses were generally robust across models, but more studies and an examination of responses with ensembles would be beneficial.

54 ENVIRONMENTAL SCIENCES

Data for Zheng et al. (2025), "AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry"

Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystems and soil function. Salinity affects soil carbon cycling by directly impacting microbial activity and indirectly altering soil physicochemical properties, but current models inadequately represent these complexities. This dataset contains the observational and modeling data from Zheng et al. (2025), which described a process-based modeling framework that couples soil solution chemistry with microbial carbon cycling reactions to study the impacts of soil salinization. This conceptual model is implemented numerically into the open-source geochemical program PHREEQC 3.0 (Parkhurst and Appelo, 2013). This dataset consists of: - Figure2_AquaMEND_salinity_buffer: Contains model simulation outputs to assess the impact of three different cation exchange and surface complexation processes on salinity buffering (Fig. 2 from Zheng et al. 2025). - Figure3_Salinity_function: Contains salinity function fitting for literature data (Fig. 3 from Zheng et al. 2025). - Figure4_AquaMEND_microbial_mechanisms: Contains model simulation outputs for testing various microbial process-based hypotheses related to soil salinization, including microbial mortality, carbon use efficiency (CUE), extracellular enzyme activity, and other microbial mechanisms (Fig. 4 from Zheng et al. 2025). - Figure5_AquaMEND_Redox: Contains on model simulation outputs to evaluate shifts among key redox processes, such as aerobic respiration, sulfate reduction, and methanogenesis (Fig.5 from Zheng et al. 2025). - Figure6_AquaMEND_sorption: Contains on model simulation outputs for investigating the effects of salinity on dissolved organic matter (DOM) sorption and desorption processes (Fig. 6 from Zheng et al. 2025). - Figure7_AquaMEND_process_couple: Contains on model simulation outputs for exploring coupled biotic-abiotic processes and their interactions (Fig. 7 from Zheng et al. 2025). - data: Includes datasets used to develop salinity response functions and evaluate salinity buffering capacity. Datasets for MEND model calibration. - database: Contains the `.dat` file required by PHREEQC for model execution. - README.md: A Markdown plain text file describing the computational tools and directories. Files are a mixture of plain text CSV (comma-separated value) and plain text *.dat files written by the model; no special software is required to read them.

EARTH SCIENCE > AGRICULTURE > SOILS > SOIL SALINIT

Enhancing Modeling and Simulation for Effective Protection Strategies

This report was created by Sandia National Laboratories (SNL) to document the principles and methodology of performance data collection and integration with modeling and simulation tools to better facilitate the performance evaluation of physical protection systems (PPS). Results and conclusions from the use of modeling and simulation tools are only as good as the data employed by the tools when conducting analysis. Acquiring the performance testing data necessary to ensure effective evaluation can be a complex and sometimes daunting process. It is the desire of the organization to provide guidance that eases the burdens associated with pursuit of these objectives. This document draws heavily upon longstanding principles of systems engineering that have been developed and employed by SNL in the discipline of security since the 1970s. The scope of this document is constrained to the testing of system components and integration of data that is applicable within the context of PPS performance analysis using two tools that have been developed by SNL, PathTrace© and Scribe3D©. For guidance related to testing and evaluation more broadly, the manuals and reports referenced by this document can be consulted.

42 ENGINEERING

CHARMM-GUI Bicelle Builder : An Extension of Membrane Builder for Modeling and Simulation of Bicelle Systems

Membrane mimetics, such as detergent micelles, nanodiscs, and amphipol complexes, which can provide membrane-like environments while retaining small and soluble features, have been utilized to study membrane proteins. A bicelle, composed of varying lipids and detergents, is a useful membrane mimetic because the lipid-to-detergent ratio, the q-value, can be adjusted to alter the properties of the aggregate, including the thickness and size of the bicelle. However, building a bicelle model for modeling and simulation studies requires nontrivial efforts, even for experts. We introduce CHARMM-GUI Bicelle Builder, a web-based platform that can generate various all-atom bicelle systems via a graphical user interface with all available lipids and detergents in Membrane Builder. To illustrate and validate Bicelle Builder with practical systems, we have modeled and simulated pure bicelles consisting of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipids with 1,2-dihexanoyl-sn-glycero-3-phosphocholine (C6DHPC) detergents and protein–bicelle complexes, composed of DMPC with C6DHPC, foscholine-10 (FOS10), and lysophosphatidylcholine-12 (LPC12) detergents. Our simulation results indicate that Bicelle Builder can generate reliable and robust bicelle models with and without proteins that retain DMPC bilayer characteristics. Bicelle Builder is expected to help researchers better understand not only bicelles themselves but also atomistic-level structures of protein–bicelle complexes that are often difficult to access through experimental approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Finding Real Uncertainties From Physical Simulations

Modeling strong gravitational lenses is computationally expensive for the complex data from modern and next-generation cosmic surveys. Deep learning has emerged as a promising approach for finding lenses and predicting lensing parameters, such as the Einstein radius. Mean-variance Estimators (MVEs) are a common approach for obtaining aleatoric (data) uncertainties from a neural network prediction. However, neural networks have not been demonstrated to perform well on out-of-domain target data successfully - e.g., when trained on simulated data and applied to real, observational data. In this work, we perform the first study of the efficacy of MVEs in combination with unsupervised domain adaptation (UDA) on strong lensing data. The source domain data is noiseless, and the target domain data has noise mimicking modern cosmology surveys. We find that adding UDA to MVE increases the accuracy on the target data by a factor of about two over an MVE model without UDA. Including UDA also permits much more well-calibrated aleatoric uncertainty predictions. Advancements in this approach may enable future applications of MVE models to real observational data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system