Search NASA⌕ Search

SEARCH · Search NASA

Results for “physics-based simulations”

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.

188 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Characterizing IHE Response to Multiple Shock Loading

The response of high explosives to shock loading is traditionally measured with a steady loading pressure. In many accident scenarios involving fragment impact, however, a loading duration that is shorter than the build up to detonation may occur. Fragments passing through multiple materials before reaching a high explosive charge may produce loading that is comprised of more than one shock wave. Additionally, the build up to detonation in high explosive corner turning loads the explosive a short duration pressure pulse, since rarefactions can often rapidly overtake the reactive wave. For these reasons, we have studied the response of the insensitive high explosive (IHE) materials PBX 9502 and LX-17 to complex loadings of varied intensity and duration. We refer to a single loading of limited duration as a “thin pulse”, whereas more complex scenarios were studied with an impactor that produces a double shock in the explosive. The following report presents experimental data and analyses of thin pulse shock initiation and double shock experiments designed to guide development of models of Insensitive High Explosives (IHEs) under controlled one-dimensional conditions relevant to accident scenarios and corner turning. Thin pulse shock initiation data on PBX 9502 and LX-17 were obtained under varied pulse duration, pressed density, and temperature conditions in order to probe various parameters essential for the development of a physics-based Cheetah reactive flow hotspot model. In situ pressure gauges provide insight into the degree of reaction in the explosive that are not obtainable with optical PDV measurements or distance measurements such as run to detonation. Double shock data was obtained to inform a Composition Aware Cheetah model which can be applied to any TATB-based IHEs. This model supports efforts to find a new IHE formulation and potentially incorporate new binders into IHE formulations. Simulations of each experiment are included to demonstrate the utility of these focused experiments to developing models of HE behavior. One-dimensional gas gun experiments are essential for characterizing shocked HE behavior and informing HE models.

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

A physics-based ensemble machine-learning approach to identifying a relationship between lightning indices and binary lightning hazard

To convert lightning indices generated by numerical weather prediction experiments into binary lightning hazard, a machine-learning tool was developed. This tool, consisting of parallel multilayer perceptron classifiers, was trained on an ensemble of planetary boundary layer schemes and microphysics parameterizations that generated four different lightning indices over 1 week. In a subsequent week, the multi-physics ensemble was applied and the machine-learning tool was used to evaluate the accuracy. Unintuitively, the machine-learning tool performed better on the testing dataset than the training dataset. Much of the error may be attributed to mischaracterizing the convection. The combination of the machine learning model and simulations could not differentiate between cloud-to-cloud lightning and cloud-to-ground lightning, despite being trained on cloud-to-ground lightning. It was found that the simulation most representative of the local operational model was the most accurate simulation tested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
DOE OSTI2025

Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1

Abstract. Groundwater plays a key role in meeting water demands, supplying over 40 % of irrigation water globally, with this role likely to grow as water demands and surface water variability increase. A better understanding of the future role of groundwater in meeting sectoral demands requires an integrated hydro-economic evaluation of its cost and availability. Yet substantial gaps remain in our knowledge and modeling capabilities related to groundwater availability, recharge, feasible locations for extraction, extractable volumes, and associated extraction costs, which are essential for large-scale analyses of integrated human–water system scenarios, particularly at the global scale. To address these needs, we developed Superwell, a physics-based groundwater extraction and cost accounting model that operates at sub-annual temporal and at the coarsest 0.5° (≈50 km × 50 km) gridded spatial resolution with global coverage. The model produces location-specific groundwater supply–cost curves that provide the levelized cost to access different quantities of available groundwater. The inputs to Superwell include recent high-resolution hydrogeologic datasets of permeability, porosity, aquifer thickness, depth to water table, recharge, and hydrogeological complexity zones. It also accounts for well capital and maintenance costs, as well as the energy costs required to lift water to the surface. The model employs a Theis-based scheme coupled with an amortization-based cost accounting formulation to simulate groundwater extraction and quantify the cost of groundwater pumping. The result is a spatiotemporally flexible, physically realistic, economics-based model that produces groundwater supply–cost curves. We show examples of these supply–cost curves and the insights that can be derived from them across a set of scenarios designed to explore model outcomes. The supply–cost curves produced by the model show that most (90 %) nonrenewable groundwater in storage globally is extractable at costs lower than USD 0.57 m−3, while half of the volume remains extractable at under USD 0.108 m−3. The global unit cost is estimated to range from a minimum of USD 0.004 m−3 to a maximum of USD 3.971 m−3. We also demonstrate and discuss examples of how these cost curves could be used by linking Superwell's outputs with other models to explore coupled human–environmental system challenges, such as water resources planning and management, or broader analyses of multisectoral feedbacks.

Global Change Analysis Model (GCAM)↗
DOE OSTI2024

Fuel Performance Evaluation of THOR-C Experiments

The Temperature Heatsink Overpower Response Commissioning (THOR-C) and THOR-Metal (THOR-M) experiments will be performed as part of an ongoing project for testing sodium fast reactor fuels with the Japan Atomic Energy Agency (JAEA). The THOR-C experiments consist of fresh metallic fuel pins and have been analyzed using the ABAQUS, Ansys codes and the BISON fuel performance code. THOR-M-Loss of Flow-1 (THOR-M-LOF-1) is designed to test an EBR-II irradiated fuel pin under LOF conditions. Simulation of the THOR-MLOF-1 experiment required first simulating the base irradiation of the fuel pin in EBR-II. MFUEL module of SAS4A/SASSYS-1 [1] is a physics-based metallic fuel performance model applicable to the normal operation, transient scenarios and fuel failure modeling including scenarios with bulk fuel melting. The model has been validated using EBR-II normal operation, separate effect transient tests as well as TREAT M-Series transient tests [2]. In this study, MFUEL models has been utilized together with a new capsule heat transfer model developed in this project. The new heat transfer model was necessary due to (1) significant amount of heat losses that required 2D heat transfer, (2) the presence of a titanium heat sink, rejecting a significant amount of heat, and (3) stagnant coolant conditions, which are inconsistent with SAS4A/SASSYS-1 (SAS) heat transfer model. Updates to SAS4A/SASSYS-1 and MFUEL has been described below, followed by a preliminary validation effort using the results from THOR-C-2 fresh fuel capsule experiment. A previous study for THOR-C-2 analysis using BISON code is also utilized in this study to model this test [3]. [1] D. O’Grady, A. J. Brunett, L. Ibarra, A. Karahan, T. Kim, T. S. Sumner, R. Thomas, T. H. Fanning, “The SAS4A/SASSYS-2 Version 5.7 Safety Analysis Code System,” Argonne National Laboratory,ANL/NSE-SAS/5.7, (2023). [2] A. Karahan, T. Kim, T. Fanning, D. O’Grady, “Validation of MFUEL Metal Fuel Performance Models of SAS4A/SASSYS-1,” Argonne National Laboratory, ANL/NSE-23/11, (2023). [3] M. Mihelish, A. Zabriskie, K. Paaren, P. Medvedev, C. Jensen, “Fuel Performance Predictions for the TREAT THOR-C Experiments,” Idaho National Laboratory, INL/RPT-23-73397, Revision 0, (2023)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Thermo-Fluid Modeling Framework for Supercomputer Digital Twins: Part 1, Demonstration at Exascale

A thermo-fluid modeling framework is being developed for ExaDigiT---an open-source framework for developing comprehensive digital twins of liquid-cooled supercomputers. The work is being conducted in two parts, and discussion is divided into two companion papers. The work documented in this paper focuses on the development of a cooling system library in Dymola for the Frontier supercomputer at Oak Ridge National Laboratory. The second part, outlined in a companion paper, focuses on a templating structure called Auto-CSM for easily creating model-agnostic, physics-based thermo-fluid cooling system models for liquid-cooled supercomputers using a text-based schema. The cooling model is being developed using primarily the open-source Transient Simulation Framework of Reconfigurable Models (TRANSFORM) library. The library follows the templating architecture developed within the TRANSFORM library for modeling subsystems. A full-system validation was performed to validate a very simple model that is integrated with the system controls, and the results are presented herein.

Kumar, Vineet↗
DOE OSTI2025

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗
DOE OSTI2026

A Physics-Informed Reinforcement Learning Framework for Economic-Thermal Co-Optimization of Crypto Mining Data Centers: Preprint

The rapid expansion of cryptocurrency mining has created a new class of high-density data centers characterized by extreme thermal flux and high sensitivity to volatile economic markets. Traditional thermal management strategies, typically reliant on rule-based control, maintain static setpoints that fail to account for fluctuating electricity prices and cryptocurrency values - factors critical to mining profitability. To address this, we present a physics-informed reinforcement learning (PIRL) framework for economic-thermal co-optimization in crypto mining data centers. This framework consists of a proximal policy optimization (PPO) agent, a virtual testbed powered by high-fidelity physics-based models, and an interactive frontend dashboard. The PPO agent is trained using the virtual testbed and strict hardware safety limits. This physics-informed approach allows the agent to learn a stochastic policy that dynamically balances mining revenue against operational costs by co-optimizing HVAC cooling setpoints and IT computational hashrate. The simulation results demonstrate that the integrated framework achieved an 8.62% increase in net operational profit compared to traditional baseline strategies while strictly adhering to safety-critical temperature constraints (coolant supply temperature < 32 degrees C). This work provides a scalable template for the deployment of reinforcement learning in mission critical facilities where economic volatility and physical safety must be managed simultaneously.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗