Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗
DOE OSTI2026

Fuel Performance Modeling Internship Final Presentation

This study examines the performance of U-Zr and U-Pd-Zr annular metallic fuel rodlets and details the current status of modeling efforts regarding U-Pu-Zr solid 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 metallic 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↗
DOE OSTI2020

Enhancement of Photoemission on p-Type GaAs Using Surface Acoustic Waves

We demonstrate that photoemission properties of p-type GaAs can be altered by surface acoustic waves (SAWs) generated on the GaAs surface due to dynamical piezoelectric fields of SAWs. Multiphysics simulations indicate that charge-carrier recombination is greatly reduced, and electron effective lifetime in p-doped GaAs may increase by a factor of 10× to 20×. It implies a significant increase, by a factor of 2× to 3×, of quantum efficiency (QE) for GaAs photoemission applications, like GaAs photocathodes. Conditions of different SAW wavelengths, swept SAW intensities, and varied incident photon energies were investigated. Essential steps in SAW device fabrication on a GaAs substrate are demonstrated, including deposition of an additional layer of ZnO for piezoelectric effect enhancement, measurements of current–voltage (I–V) characteristics of the SAW device, and ability to survive high-temperature annealing. Results obtained and reported in this study provide the potential and basis for future studies on building SAW-enhanced photocathodes, as well as other GaAs photoelectric applications.

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

thornado-hydro: A Discontinuous Galerkin Method for Supernova Hydrodynamics with Nuclear Equations of State

This paper describes algorithms for non-relativistic hydrodynamics in the toolkit for high-order neutrino radiation hydrodynamics (thornado), which is being developed for multiphysics simulations of core-collapse supernovae (CCSNe) and related problems with Runge–Kutta discontinuous Galerkin (RKDG) methods. More specifically, thornado employs a spectral type nodal collocation approximation, and we have extended limiters — a slope limiter to prevent non-physical oscillations and a bound-enforcing limiter to prevent non-physical states — from the standard RKDG framework to be able to accommodate a tabulated nuclear equation of state (EoS). To demonstrate the efficacy of the algorithms with a nuclear EoS, we first present numerical results from basic test problems in idealized settings in one and two spatial dimensions, employing Cartesian, spherical-polar, and cylindrical coordinates. Then, we apply the RKDG method to the problem of adiabatic collapse, shock formation, and shock propagation in spherical symmetry, initiated with a 15 M ⊙ progenitor. Herein, we find that the extended limiters improve the fidelity and robustness of the RKDG method in idealized settings. The bound-enforcing limiter improves robustness of the RKDG method in the adiabatic collapse application, while we find that slope limiting in characteristic fields is vulnerable to structures in the EoS — more specifically, in the phase transition from nuclei and nucleons to bulk nuclear matter. The success of these applications marks an important step toward applying RKDG methods to more realistic CCSN simulations with thornado in the future.

Nuclear astrophysics (1129)↗
DOE OSTI2026

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2026

Dielectric Resonator Design for Low Power and Low Temperature Microwave Plasma

Waveguide-based microwave plasmas generally operate at high temperatures (2000 - 6000K)[1], making it difficult to directly interface solid materials with the plasma without significant thermal damage. Dielectric microwave resonators (DMRs), long studied for wave-based manipulation of electromagnetic radiation for telecom and optics, can focus radiation to extremely small mode volumes, creating intense localized fields with low-power input.[2] This phenomenon can be used for applications ranging from efficient plasma electronics to near-ambient plasma-materials interactions. Such DMR-based plasmas have been demonstrated a handful of times in the literature, but the majority of research towards this utilize the lowest frequency resonance mode.[3], [4], [5] By carefully controlling the geometry of cylindrical resonators, a variety of electromagnetic modes can be excited. In this work, COMSOL Multiphysics simulations are used to study the electric field enhancement and absorption properties of CaTiO3 DMRs as a function of geometry and excitation frequency. Whereas previous studies have utilized the HEM111 resonance frequency to drive low power plasma excitation, we find that higher order resonance frequencies are more effective at field enhancement and result in less power loss within the dielectric material, hence less wasted heating. The effectiveness of these modes is also geometry dependent and can be computationally optimized for plasma generation. Complementing these computational efforts, we demonstrate a new closed-system reactor design built in a WR-650 waveguide and experimentally demonstrate the formation of atmospheric argon microwave plasma using < 30 W input power on DMR dimers. We observe a shifting resonance frequency as the DMRs heat in response to microwave excitation and develop a Python-based lock-in mechanism to effectively track the DMR resonance over time, leading to stable plasma operation. We use infrared thermal imaging to monitor the temperature of the DMR dimers and surrounding quartz chamber, demonstrating thermal temperatures < 60 degreesC. Finally, we utilize optical emission spectroscopy (OES) to probe the plasma properties as a function of the resonance mode.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2020

Shape Optimization of Vortex Generators to Control Mack Mode Amplification

This paper demonstrates the potential to use shape optimization for the design of vortex generators in an axisymmetric boundary layer. This shape optimization increases the amplitude of stationary streaks created by streamwise vorticity, with the goal of reducing the amplification of Mack mode instabilities. The test case under consideration matches a trajectory point during the ascent phase of the HIFiRE-1 flight experiment. Wall-mounted vortex generators are added to the cone and their shape is optimized to control the amplification of Mack mode instabilities that are known to initiate laminar-turbulent transition in this flow. An empirical objective function is developed based on the previous studies of optimal streaks. A constraint on the maximum streak amplitude is added to avoid the initiation of bypass transition. The shape optimization is conducted using SU2, an open-source suite for multiphysics simulation and design. A significant improvement is observed in an integral metric of the streak amplitude, while the maximum streak amplitude is maintained close to the baseline level. The qualitative features of the optimized geometry are discussed.

Vortex generators↗
NASA NTRSDate not supplied

A Multiphysics Study to Improve Specific Energy of Primary Batteries for Low Temperature Operation for Deep Space Missions

Several lander missions on the outer planets such as Europa, Enceladus, and Titan require electrical power to operate scientific and communication equipment. The traditional power generation methods, such as a photovoltaic array, are not feasible as their efficiency drops significantly at these vast distances. The novel radioisotope power systems are not practical today based on current lander designs and the effectiveness of these systems. To perform in situ science on distant planets, a high specific energy battery (>700 Wh/kg) needs to operate for about 480 hours under cold temperatures (-40C or 0C) [1]. While a primary battery such as Li-CFx can provide high specific energy at room temperature, its specific capacity decreases significantly at low temperatures. One of the causes for this drop is low ion and electrical conductivity, and slower reaction kinetics. Slower transport and facile kinetics lead to an increase in the battery’s resistance and higher voltage drops during the cell operation, thus reducing specific capacity. Both the transport and kinetics show a strong dependence on temperature. Thus, a small temperature rise can lead to an increase in the reaction rate and ion conductivity; since the temperature, cell resistance, and specific capacity are interdependent. A conventional battery model accounts for ohmic, thermodynamic, and, electrochemical, and chemical decomposition heating. The ohmic heating can be controlled by designing a resistive microstructure and varying the ratios of the active materials [2]. The kinetics can be improved by increasing the surface area, reducing the particle size, or adding a catalyst. These parameters are often optimized to achieve high specific energy at room temperatures. A similar optimization study is not available at low temperatures and for a primary (high specific energy) battery. For this presentation, we will explore the effect of geometrical, microstructural, and material properties on optimal specific capacity at low temperatures through multiphysics simulations. The ion transport resistance depends on the porosity and the tortuosity of an electrode and the separator.

M. Mehta↗
NASA NTRSDate not supplied

Optimizing the Electroforming Process to Enhance the Thickness Uniformity of Full Shell X-Ray Optics

Electroforming replication technology at Marshall Space Flight Center has a long heritage of producing highquality full-shell X-ray mirrors for various applications. Nickel alloys are electroformed onto a super-polished mandrel in the electroforming process, then separated to form the replicated full-shell optic. Various parameters in the electroplating configuration could result in the nonuniformity of the shell’s thickness. Thickness nonuniformities primarily occur due to non-uniform electric field distributions in the electroforming tank during the deposition. Using COMSOL Multiphysics simulations, we have studied the electric field distributions during the deposition process. Using these studies, we have optimized the electric fields inside the tank using customized shields and insulating gaskets on the mandrel. These efforts reduced thickness non-uniformity from over 20% to under 5% percent. Improving the thickness uniformity of the shell aids in better mounting and alignment of shells in the optics module. Optimization of the electroforming process, in some cases, improved the optical performance of the shells. COMSOL optimizing of the electroforming process and the experimental results validating these simulations are presented in this article.

X-ray optics↗
DOE OSTI2023

Release on the Virtual Test Bed of an MSRE thermal hydraulics model

The nuclear industry is taking leaps in innovations with companies seeking a sustainable energy future through advanced nuclear reactors. The \gls{DOE}’s \gls{neams} program seeks to substantiate and bolster the deployment of advanced reactors through flexible multifidelity, multiphysics simulations of advanced nuclear reactors. Applications like SAM for one-dimensional systems thermal-hydraulics, and Pronghorn for multidimensional coarse mesh thermal-hydraulics, are geared to support innovations in industry by facilitating design, optimization, and licensing of advanced nuclear reactors. Coupling systems thermal-hydraulics and computational fluid dynamics codes can be difficult as the pressure coupling converges slowly; however, it is important to obtain the desired accuracy in each part of the primary loop. The authors of this model created an \gls{odc}~\cite{Mau23} approach to coupling SAM and Pronghorn. Leveraging this coupling technique, a \gls{msre} model was developed and released to the \gls{neams}/\gls{nric} \gls{vtb}. The \gls{msre} was chosen to be modeled because of the wealth of experimental data available and because of the strong physics coupling between the core and primary circuit \cite{doi:10.13182/NT8-2-118}. This document contextualizes the history of the \gls{msre}, describes the thermal hydraulics models used, and detail the implementation of multidimensional thermal-hydraulics and system codes based on the \gls{odc} method~\cite{Penn} for the \gls{msre} model. Finally, this document presents how other modelers could apply the SAM-Pronghorn \gls{odc} for other advanced reactor models. Current item is the set of slides for ANS Winter 23. The release of the model and the ANS summary have already been approved

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Electrochemical Characterization of Molten Salt Chemistry During Atmospheric Ingressions

This report describes the salt chemistry and sedimentation studies that were completed in FY23 as part of the Pyrochemical Fuel Cycles – ANL project. The primary objective of these activities was to employ electrochemical methods in order to quantify and gain insights into the reaction mechanisms associated with the interaction of O2 and moisture impurities with the molten chloride salts used for the pyrochemical processing of nuclear materials. If not monitored and controlled, the ingression of these atmospheric impurities can lead to changes in the salt redox conditions and initiate the formation of oxide particles, potentially resulting in operational and safeguards challenges. By elucidating the mechanisms underlying the formation of oxide particles, we can effectively monitor the conditions of the salt and design systems to mitigate the influence of atmospheric O2 and H2O. To achieve these goals, we designed a test system that allowed for precisely controlled gas ingressions into the molten salt vessel. This system included on-line monitoring provided by electrochemical probes that enable near real-time measurements of the salt conditions throughout the course of the experiment. Using this system, we systematically varied experimental conditions, including O2/Ar flow rates and O2 concentrations, to comprehensively understand their impact on solid particle formation. We also used a particle size and shape analyzer to characterize the particles that were generated. Tests in FY23 concentrated on the formation of CeO2 particles, but we also prepared a separate apparatus targeting UO2 particles for use in FY24. We additionally conducted extensive modeling activities in support of this work. This included multiphysics simulations of the gas ingressions along with the development of a machine learning workflow to enable accelerated molecular dynamics modeling of the molten salt chemistry of LiCl-KCl-UCl3. The combination of experimental and modeling tools has allowed us to begin to get a more complete understanding of the reactions that occur when atmospheric ingressions occur in pyroprocessing systems.

Guo, Jicheng↗