Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiphysics 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.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Machine learning–aided real-time detection of keyhole pore generation in laser powder bed fusion

Porosity defects are currently a major factor that hinders the widespread adoption of laser-based metal additive manufacturing technologies. One common porosity occurs when an unstable vapor depression zone (keyhole) forms because of excess laser energy input. With simultaneous high-speed synchrotron x-ray imaging and thermal imaging, coupled with multiphysics simulations, we discovered two types of keyhole oscillation in laser powder bed fusion of Ti-6Al-4V. Amplifying this understanding with machine learning, we developed an approach for detecting the stochastic keyhole porosity generation events with submillisecond temporal resolution and near-perfect prediction rate. Finally, the highly accurate data labeling enabled by operando x-ray imaging allowed us to demonstrate a facile and practical way to adopt our approach in commercial systems.

42 ENGINEERING↗
DOE OSTI2026

Two-Level Sketching Alternating Anderson Acceleration for Complex Physics Applications

We present a novel two-level sketching extension of the Alternating Anderson–Picard (AAP) method for accelerating fixed-point iterations in challenging single- and multiphysics simulations governed by discretized PDEs. Our approach combines a static, physics-based projection that reduces the least-squares (LS) problem to the most informative field (e.g., via Schur-complement insight) with a dynamic, algebraic sketching stage driven by a backward stability analysis under Lipschitz continuity. We introduce inexpensive estimators for stability thresholds and cache-aware randomized selection strategies to balance computational cost against memory access overhead. The resulting algorithm solves reduced LS systems in place, minimizes memory footprints, and seamlessly alternates between low-cost Picard updates and Anderson mixing. Implemented in Julia, our two-level sketching AAP achieves up to 50% time-to-solution reductions compared to standard Anderson acceleration—without degrading convergence rates—on benchmark problems including Stokes, 𝑝-Laplacian, bidomain, and Navier–Stokes formulations at varying problem sizes. These results demonstrate the method’s robustness, scalability, and potential for integration into high-performance scientific computing frameworks. Our implementation is available open source in the AAP.jl library.

Barnafi, Nicolas [University of Chile, Santiago]↗
DOE OSTI2026

Discovery of Probabilistic Dirichlet-to-Neumann Maps on Graphs

Dirichlet-to-Neumann maps enable the coupling of multiphysics simulations across computational subdomains by ensuring continuity of state variables and fluxes at artificial interfaces. We present a novel method for learning Dirichlet-to-Neumann maps on graphs using Gaussian processes, specifically for problems where the data obey a conservation law arising from an underlying partial differential equation. Our approach combines discrete exterior calculus and nonlinear optimal recovery to infer relationships between vertex and edge values. This framework yields data-driven predictions with uncertainty quantification across the entire graph, even when observations are limited to a subset of vertices and edges. By minimizing the reproducing kernel Hilbert space norm while penalizing kernel complexity through maximum likelihood estimation, our method ensures that the resulting surrogate strictly enforces conservation laws without overfitting. We demonstrate our method on two representative applications: subsurface flow in fracture networks and arterial blood flow. Finally, the results demonstrate that the method maintains high accuracy and well-calibrated uncertainty estimates even under severe data scarcity, highlighting its potential for scientific applications where limited data and reliable uncertainty quantification are critical.

Dirichlet-to-Neumann map↗
DOE OSTI2026

Trilinos: Enabling Scientific Computing across Diverse Hardware Architectures at Scale

Trilinos is a community-developed, open-source software framework that facilitates building large-scale, complex, multiscale, multiphysics simulation code bases for scientific and engineering problems. Since the Trilinos framework has undergone substantial changes to support new applications and new hardware architectures, this document is an update to “An Overview of the Trilinos project” by Heroux et al. (ACM Transactions on Mathematical Software, 31(3):397–423, 2005). It describes the design of Trilinos, introduces its new organization in product areas, and highlights established and new features available in Trilinos. Particular focus is put on the modernized software stack based on the Kokkos ecosystem to deliver performance portability across heterogeneous hardware architectures. This article also outlines the organization of the Trilinos community and the contribution model to help onboard interested users and contributors.

Heterogeneous Hardware Architectures↗
DOE OSTI2023

Electrochemical Residence Time Distribution as a Diagnostic Tool for Redox Flow Batteries

The fluid dynamic and electrochemical performance of redox flow batteries (RFBs) stems from the relationship between the flow field and the porous electrode, whose interplay determines how active species move and react during device operation. While characterization techniques, such as residence time distribution, offer insights into species mobility within a reactive volume for a traditional chemical reactor, electrochemical reactors also enable simultaneous measurement of the redox reactions, unlocking another dimension of analysis. Herein, we demonstrate how potentiodynamic measurements, using injections of electrolyte examined through moment analysis, can provide electrode-specific performance scaling relationships across a matrix of carbon paper and cloth electrodes with flow through and interdigitated flow fields. We further combine experimental campaigns with multiphysics simulations to demonstrate how electrode surface area can be estimated with this technique, which we then validate with activated and unactivated commercial carbon cloth electrodes. These studies reveal the multiscale observations that potentiodynamic measurements afford, augmenting existing electrochemical techniques for holistic electrochemical reactor diagnostics.

Electrochemistry↗
DOE OSTI2023

NED_Asteroid_Energy_Deposition

In the event of a potentially catastrophic asteroid impact, with sufficient warning time, deploying a nuclear device remains a powerful option for planetary defense if a kinetic impactor proves insufficient. Predicting the effectiveness of a potential nuclear deflection or disruption mission depends on accurate multiphysics simulations of the device's x-ray energy deposition into the asteroid and the resulting material ablation. These simulations span many orders of magnitude, require a variety of different complex physics packages, and are computationally expensive. Having an efficient and accurate way of modeling this system is necessary for exploring a mission's sensitivity to the asteroid's range of physical properties. To expedite future simulations, we present a completed x-ray energy deposition model developed using the radiation-hydrodynamics code Kull which can be used to initiate a nuclear mitigation mission hydrocode. The model spans a wide variety of possible mission initial conditions: four different asteroid-like materials (Silicon Dioxide, Forsterite, Iron, and Ice), two different source spectra (1 and 2 keV blackbodies), and then a broad range of radiation fluences (0.0001 to 1 kt per square meter), source durations (10 to 100 ns), and asteroid porosities (0 to 80 percent). Using blowoff momentum as the primary metric, the modelinitiated simulation results match the full radiation-hydrodynamics results to within 10 percent. Please reference the journal article: Burkey et al., X-Ray Energy Deposition Model for Simulating Asteroid Response to a Nuclear Planetary Defense Mitigation Mission, Planetary Science Journal, (2023) for more information.

Managan, RobertA↗
DOE OSTI2021

Griffin-m

Griffin-m is a MOOSE-based reactor multiphysics application that streamlines the analysis of a variety of nuclear multiphysics simulations, including steady-state and transient radiation transport, core performance, fuel depletion, criticality and decay heat calculations, reprocessing and post-irradiation examination. This streamlining is accomplished via enhanced flexibility of the tools, uniform syntax in the MOOSE framework, dynamic linking of all relevant physics and a single point of execution. The design for flexible multi-physics, multi-radiation, multi-scheme tasks demands and ultimately makes Griffin-m a highly extendable code system. A software quality assurance (SQA) procedure is enforced during Griffin-m development. The Griffin-m version will be the main production version and it includes the MCC3 code, which was originally developed at Argonne National Laboratory to prepare cross sections for fast reactors.

Ortensi, Javier [Idaho National Laboratory (INL), ↗
DOE OSTI2021

Godiva Experiments for the Nuclear Criticality Safety Program (NCSP)

Godiva IV is a fast burst critical assembly constructed of approximately 65 kg of highly enriched uranium (HEU) fuel alloyed with 1.5 percent molybdenum for strength. Godiva is one of the last such critical assemblies in the United States, and can be used for studies of super-prompt critical behavior as well as irradiations and demonstrations. A demonstration of a Godiva burst is usually performed as a highlight of the hands-on portion of the Criticality Safety Training Classes taught at the National Criticality Experiment Research Center (NCERC). The Godiva burst is used to demonstrate the concept of super-prompt critical and the time-scale of a criticality accident. In addition, several NCSP projects have been conducted on Godiva IV over the past two years. One experiment focused on collecting data to support multiphysics simulations using Photo-Doppler Velocimetry (PDV) to measure surface movement and gamma detectors to measure the burst output as the burst develops from background to peak over ten orders of magnitude. Another experiment was performed to demonstrate the functionality of the Criticality Accident Alarm System (CAAS) system developed for installation in the Y-12 Uranium Processing Facility (UPF). The system must not only respond to a criticality event and alarm, but must also be shown to operate in a high dose environment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Modeling and Validation of Irradiation Damage in Ni-based Alloys for Long-Term LWR Applications

Grizzly is a multiphysics simulation code under development at Idaho National Laboratory to simulate aging mechanisms and their effects on the integrity of critical LWR components, including Ni-Cr alloy components. Predictions of the LWR component integrity rely on the ductile-to-brittle transition temperature (DBTT) curve to determine the toughness of the alloy after a period of exposure to in-service reactor conditions. The macroscopic DBTT curve is governed by the microstructure evolution of the alloy, thus, a capability to predict the microstructure evolution in response to thermal ageing is being implemented into Grizzly. Within this report section, the implementation of the Avrami equation to estimate the phase fraction of the aged long range ordered precipitates is discussed. The initial evaluation of current crystal plasticity models within Grizzly to model these Ni-Cr alloys is also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Report on Aging and Irradiation Theoretical Model Development for Alloy 690 in Grizzly

This report describes the current efforts to implement a model to predict the engineering scale behavior of irradiated and aged Ni 2 Cr model alloy and Alloy 690 into the Grizzly software. Grizzly is a multiphysics simulation code under development at Idaho National Laboratory to simulate aging mechanisms and their effects on the integrity of critical LWR components. Engineering scale behavior of materials are strongly impacted by the underlying microstructure: the distribution and interaction of dislocations and defects within the individual grains dictates how the material will respond to loading. The evolution of these dislocations and defects determines how the engineering scale material will perform under different loading and environmental conditions. Mesoscale models can function as a bridge to connect microstructure evolution models to engineering models of the ductile-brittle transition temperature (DBTT) curve and fracture. The toughness of a material, which governs the engineering scale fracture behavior, is a function of both elastic deformation and plastic deformation. While the yield stress measurements capture the elastic contribution to a material’s toughness, information about the dislocation movement is required to understand the plastic component of toughness. Crystal plasticity models predict dislocation density evolution and thereby connect microstructure evolution to engineering scale behavior predictions. These models calculate the plastic strain within a metal as a function of dislocations and other crystal defects, including precipitates and interstitial loops, among others, by tracking dislocation movement and the interaction of these dislocations with the crystal defects. Within NiCr alloys, long-range ordered (LRO) precipitates form under aging conditions. These precipitates act as barriers to dislocation motion and are either sheared by dislocations or pin dislocations. Additional defects, such as stacking fault tetrahedra (SFTs) from irradiation damage, also impede dislocation motion in these alloys. A crystal plasticity mode for Ni 2 Cr model alloy and Alloy 690 is being currently implemented in the Grizzly code at INL. This report includes a description of the theoretical models incorporated into the Grizzly crystal plasticity code. These theoretical models include mechanisms for dislocation glade and twinning as well as dislocation-defect interaction relationship for both LRO precipitates and SFTs. We have selected these dislocation models based on the deformation mechanisms observed in experimental micropillar compression tests conducted on the Ni 2 Cr model alloy. The models are then extended to Alloy 690 with ongoing parameter fitting.

36 MATERIALS SCIENCE↗
DOE OSTI2019

DRACO and Nonequilibrium Statistical Mechanics of Aging [Slides]

Topics covered include: Types of problems which DRACO (Multiphysics simulation) was designed to address; Tools from nonequilibrium statistical mechanics; Rare events (Extreme value statistics); Transients (Diverging time scales); and, Thresholds (Absorbing phase transitions).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Charged Wellbore Casing Controlled Source Electromagnetics (CWC-CSEM) for Reservoir Imaging and Monitoring (Final report)

This project addresses the needs of the U.S. Department of Energy (DOE) to develop advanced monitoring technologies and protocols to track the fate of subsurface carbon dioxide (CO2) plumes for carbon storage. Specifically, the project seeks to develop and test a unique and novel system of technologies consisting of electromagnetic data acquisition, coupled multiphysics imaging, and reservoir model enhancement to understand the migration and long-term distribution of CO2 in the subsurface. The overarching objective is to develop an integrated approach for long term monitoring of carbon storage. The two main components of the project include the methodology development and the test of the method at a field site. The methodology component consists of 1) developing the field procedure and protocol for collecting time-lapse controlled-source electromagnetic (CSEM) data with source electric current injected into the subsurface through wellbore casings; 2) building of background 3D electrical conductivity utilizing multiple sources of data such as supplemental surface transient EM (TEM) surveys, well-logs, and seismic structural information, for enhancing CSEM signal from reservoir depths; and 3) coupled multiphysics simulations and inversion of CSEM data constrained by production data and by structural information from seismic imaging of the reservoir and overlying formations. The testing component used the field site of Bell Creek Oil Field, which served both as a field laboratory for the method development as well as a test site to evaluate the CSEM signal strengths and the methodology developed in this research project. We have accomplished all the proposed tasks and developed the methodology as planned. These include the procedure for time-lapse CSEM data acquisition, data processing techniques, integration with 3D conductivity model building, fast reservoir simulation for history matching using machine learning, and interpreting CSEM data with coupling to the reservoir modeling. Collectively, the outcome of these tasks form a coherent workflow that can be applied to monitor dedicated carbon storage in saline reservoirs. The testing component evaluated the applicability and limitations of the method, and concluded that the method would be ideal for monitoring dedicated carbon storage sites utilizing saline reservoirs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Expanded verification and validation studies of hypersonic aerodynamics with multiple physics-fidelity models

Hypersonic aerothermodynamics is an important domain of modern multiphysics simulation. The Multi-Fidelity Toolkit is a simulation tool being developed at Sandia National Laboratories to predict aerodynamic properties for compressible flows from a range of physics fidelities and computational speeds. These models include the Reynolds-averaged Navier–Stokes (RANS) equations, the Euler equations with momentum-energy integral technique (MEIT), and modified Newtonian aerodynamics with flat-plate boundary layer (MNA+FPBL) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities, as well as model validation with uncertainty considerations. Code verification activities on the MNA+FPBL model build on previous work by focusing on the viscous portion of the model. Viscous quantities of interest are compared against those from an analytical solution for flat-plate, inclined-plate, and cone geometries. The code verification methodology for the MEIT model is also presented. Test setup and results of code verification tests on the laminar and turbulent models within MEIT are shown. Solution-verification activities include grid-refinement studies on simulations that model the HIFiRE-1 wind tunnel experiments. These experiments are used for validation of all model fidelities. A thorough validation comparison with prediction error and uncertainty is also presented. Three additional HIFiRE-1 experimental runs are simulated in this study, and the solution verification and validation work examines the effects of the associated parameter changes on model performance. Finally, a study is presented that compares the computational costs and fidelities from each of the different models.

42 ENGINEERING↗
DOE OSTI2025

Increasing Radon Daughter Plate-out onto Superconducting Qubits

Radon daughter plate-out may pose a fundamental limitation to realizing large-scale superconducting quantum computers. Radon daughters plated-out onto a superconducting qubit chip can decay, releasing energy into the substrate, producing phonons. Phonon induced quasiparticles generate bit flip errors and reduce coherence times. We seek to increase plate-out on chips to enable future study of the impact of Radon on coherence times. In this project, we utilize COMSOL Multiphysics simulations to finalize design for a device to increase radon plate-out.

Kay, Dax [Northwestern U.]↗
DOE OSTI2025

RAG for FLAG: AI Assistance for a Physics Code

Artificial intelligence (AI) has quickly become an important tool in scientific research, where significant efforts are underway to develop tools that will expedite the research process. One area of particular impact is scientific software, which can be particularly complex, and therefore time consuming to learn and use effectively. AI assistants are increasingly helping to streamline the process by performing tasks such as interactively answering user questions or suggesting solutions. Los Alamos National Laboratory (LANL) develops several advanced scientific codes, such as FLAG, which can be used to run multiphysics simulations. With this study, our goal was to develop an AI assistant for FLAG that could help make the process of understanding the software and running physics simulations more efficient. To develop an AI assistant for FLAG, we used a method called retrieval-augmented generation (RAG), which is a technique that uses information from relevant data sources to enhance the accuracy of large language models (LLMs). We used the FLAG user manual and other FLAG documentation as the knowledge base for the RAG system. When a user provides a query, RAG retrieves relevant sections from the knowledge base in response, then uses those excerpts to generate grounded and contextually rich answers. We found that our AI assistant was able to provide context aware answers and source references to user queries. To evaluate performance, we developed a set of 40 benchmark questions and compared the accuracy of the responses to those of two standard LLMs without retrieval. Our AI assistant significantly outperformed the standard LLMs at answering FLAG-related questions, with an 82.5% accuracy rate, compared to 47.5% for both of the standard LLMs. This has the potential to make the process of learning and using FLAG much easier, especially for new users. Ultimately, it supports LANL’s broader mission by empowering scientists and engineers to focus more on discovery and analysis rather than on navigating complex software systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Exact signed distance fields using parallel Fast Sweeping Method

Signed distance fields are often used in multiphysics simulations to track material interfaces. We present a simple methodology based on the fast sweeping method to generate the exact signed distance from triangular meshes and linear paths on Cartesian grids. The methodology propagates the closest primitive to the boundary to the rest of the domain following the characteristics. A local upwind criterion is used to decide between the new and existing closest primitive at each grid point while capturing the correct sign of the global function. The methodology has optimal computational complexity and runs efficiently in distributed-memory architectures. We include 2D and 3D test cases along with a resolution study up to 0.512 trillion zones and 1,000 computer cores. The solution strategy can also be applied to other types of meshes or collections of primitives.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Ab initio equation of state for vanadium

Accurate materials’ equations of state (EOS) are essential for understanding materials properties as well as for use in multiphysics simulations. In particular, hydrodynamics simulations are based on three fundamental conservation laws (mass, momentum and energy) that form an under-determined system of equations. The equation of state serves as an additional closure relation between thermodynamic variables for a given material that enables numerical hydrodynamics simulation. In this report, we focus on the development of an ab initio EOS for the body centered cubic (BCC) phase of Vanadium (V) for eventual integration into a multiphase EOS in the OpenSesame EOS database.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Multiphysics Evaluation of Annular Uranium-Zirconium Metallic Fuels [Poster]

This study examines the performance of U-10Zr annular metallic fuel rodlets which were experimentally evaluated as part of the Advanced Fuels Campaign (AFC). The AFC mission is to develop novel fuel technologies and facilitate the implementation of those technologies by industry partners. A key objective is to improve steady-state and transient performance over current fuel types. The experiments of interest in this study included annular metallic U-Zr fuel rodlets within HT-9 cladding which were placed in SS-316 capsules and inserted in the Advanced Test Reactor (ATR). Certain mechanical and thermal conditions cannot be directly evaluated through experiments and fuel performance modeling is used to shed light on this evolution over time. In this study, BISON Multiphysics simulations are leveraged to investigate the state of the fuel system throughout and after the experimental conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗