Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphysics 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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Improved internal short circuit models for thermal runaway simulations in lithium-ion batteries

Thermal runaway (TR) modeling is one of the primary tools that can be used to overcome challenges associated with lithium-ion battery (LIB) safety. Among all LIB accidents that have occurred over the past decade, Internal Short Circuit (ISC) remains the most common trigger mechanism. Many available models in the literature either use a simplified approach to simulate ISC or completely ignore its contribution. The aim of this study is to understand the nature of the heat released for different types of ISC scenarios, including aluminum-anode, anode–cathode, and copper-cathode ISC. We study ISC behavior using a coupled electrochemical–thermal model with an integrated TR chemical kinetics solver built in the COMSOL Multiphysics framework. The time duration of heat release and the magnitude of the peak ISC current are studied as functions of parameters such as the size of the penetrating filament and the capacity of the cell. The numerical results are used to build an empirical model validated against the published experimental TR propagation data. Our model can be successfully used as a viable low-cost substitute in lower order (lumped) TR simulations to enable TR prevention and mitigation.

Singh, Bakhshish Preet (ORCID:0000000264751992)↗
DOE OSTI2023

The Kimberlina synthetic multiphysics dataset for CO 2 monitoring investigations

Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

58 GEOSCIENCES↗
DOE OSTI2022

Status of the NEAMS and ARC fast reactor tools integration to the NEAMS Workbench

The Workbench initiative was launched in FY-2017 within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to facilitate the transition from conventional tools to high-fidelity tools. The NEAMS Workbench provides a common user interface for model creation, real-time validation, execution, output processing, and visualization for integrated codes. The integration of the Argonne Reactor Computation (ARC) suite of codes into the NEAMS Workbench through the PyARC module was initiated in FY-2017. The ARC codes, which focus on fast reactor multiphysics analysis, contain both legacy codes like DIF3D and REBUS-3 that were developed with over 30 years of experience, and newer NEAMS additions like MC 2 -3, PERSENT, and PROTEUS. Recent work expended this integration to other tools to support the U.S. fast reactor community such as DASSH for sub-channel thermal-hydraulics, Griffin for high-fidelity deterministic transport calculations, and OpenMC for Monte Carlo simulations (with Shift integration planned for FY-2023). The integration of the “extended ARC” suite of codes into the NEAMS Workbench interface relies on the PyARC and PyGriffin modules to handles the pre- and post-processing of these codes input, and the runtime environment. The PyARC module together with the NEAMS Workbench interface are both released under Open Source Software licenses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Berkeley Lab Finite Element Framework (BELFEM) v0.1

The software program, referred to as BELFEM, is a specialized finite element code designed for the magnetodynamic modeling of high-temperature superconducting (HTS) tapes. It incorporates novel mixed finite element formulations, particularly the h-ϕ-formulation with thin-shell simplification, to efficiently simulate larger geometries. This methodology is extremely promising for predicting the electrodynamic performance of HTS tapes used in superconducting cables and magnets, offering the benefit of reduced computational cost. Compared to similar technologies like COMSOL Multiphysics and GetDP, BELFEM's performance benchmarking indicates superior efficiency in its thin-shell implementation. In the future, it will also support features like thermal coupling and inter-tape current sharing, enhancing its utility in research and development, particularly in nuclear fusion applications. The intent is to develop BELFEM as a robust and efficient tool for the HTS community, contributing to the analysis and design of superconducting cables and magnets.

Messe, Christian↗
DOE OSTI2022

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

IMPACT: Design of Integrated Multiphysics Producible Additive Components for Turbomachinery

The overall objective of the IMPACT program was to enable a dramatic reduction in design maturation time for an additive hot-section turbomachinery component through the following: • A fast crack-risk producibility surrogate model generated from machine learning applied to additive process simulation data generated via exascale computing, • Linking this surrogate model to multi-physics topology optimization (TO) to enable the creation of producible, near-optimal structural/thermal designs for additive hot-section components, • Maturing this toolset to reduce hot-section component design-for-manufacturing iterations by a large fraction, and eventually, • Using these tools to develop more efficient gas turbines in much shorter design cycle times.

42 ENGINEERING↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Modeling and experimental validation of a W f /W-fabrication by chemical vapor deposition and infiltration

Tungsten (W) has a unique combination of excellent thermal properties, low sputter yield, low hydrogen retention, and acceptable activation. Therefore, W is presently the main candidate for the first wall material in future fusion devices. However, its intrinsic brittleness and its further embrittlement during operation bears the risk of a sudden and catastrophic component failure. As a countermeasure, tungsten fiber-reinforced tungsten (W f /W) with extrinsic toughening is being developed. A possible synthesis route is chemical vapor deposition (CVD) using heated W fabrics as substrate. The challenge is that the growing CVD-W can isolate domains from precursor access leading to strength-reducing pores. To deepen the process understanding and to optimize the CVD parameters, models were developed with COMSOL Multiphysics and validated experimentally. W deposition rate equations as function of the temperature and the partial pressures of the precursors H 2 and WF 6 were experimentally validated in previous work. In the present article, these equations are applied to obtain partial pressures within the CVD reactor. The results are taken as input for transient simulations in the microscale, in which W coatings, growing onto multiple adjacent W fibers, were simulated via mesh deformation and remeshing. The surface-to-surface contact of the W coatings and the corresponding potential pore formation were simulated by implementing sophisticated deposition rate stop conditions. Within the measuring uncertainties of ≃ ±1%, the models are validated successfully by experimental comparison regarding the deposition rate, pore structure, and relative densities ranging from 0.6 to 0.9.

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

Modeling of flash sintering of ionic ceramics

A fundamental understanding of the influence of defects in ionic ceramics at the atomic, microstructural, and macroscopic levels, before, during, and after the flash sintering event is key to the development of ceramic processing operations that lead to fast, low cost, and environmentally safe fabrication of materials. The observed phenomenology of the flash process encompasses multiple time and length scales and have resulted on a wide variety of what sometimes appears to be contradictory explanations. Here, this article summarizes the latest developments on the modeling and simulation of flash sintering, specifically those related to the understanding of the equilibrium and kinetic properties and the corresponding microstructural evolution of ionic ceramics. Challenges and opportunities in the development of theoretical analyses that include unidentified multiphysical effects are discussed, as they pertain to the processing of technologically relevant ceramic materials for advanced structures and devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Smoothed particle hydrodynamics modeling and analysis of oxide reduction process for uranium oxides

A common kinetic feature for oxide reduction chemical/electrochemical processes is oxygen transport via a porous metallic layer, which has been considered as a rate-determining step for reducing uranium oxides to metallic uranium. Accounting this kinetic behavior must involve the resolution of the moving reactive interface between shrinking oxide and expanding metal phases. This study presents a numerical model using smoothed particles hydrodynamics (SPH) to effectively deal with the evolution of the shrinking core reaction interface and oxygen transport via mass transfer of lithium oxide (Li 2 O) species in multiple mass transfer domains. We successfully validated the proposed model against a theoretical derivation for the oxide reduction process handling a shrinking oxide core with molten salt and metal ash medium on a simple planar geometry. Armed with successful validation results, the model examined a realistic reactant geometry to extend the arguments beyond the one-dimensional analyses, allowing the proposed model to apply to general application scenarios with multi-dimensional geometry. Here, this study demonstrated that the proposed model could simulate and evaluate an arbitrary reaction basket design without iterative experimental trials, which is prohibitive for a scaled high-temperature molten salt study in an inert environment. The construct potentially provides not only deep insights on multiphysics behaviors governing the process dynamics but also a robust framework for evaluating and screening candidate basket designs in the most cost-effective manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment FY2020: Part 3

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: Task2c: Explicit modeling of pebble transient temperature response. In this simulation, the 400 MWth Pebble-Bed Modular Reactor (PBMR) design, PBMR- 400, experiences a 20-second power ramp from 100% to 150% power. This is followed by a similar reduction in the power back to 100%. Several multiscale pebble coupling approaches are tested with one pebble per mesh element in the active core region. The results show good conservation behavior and the stability of the coupling.; Extended scope part 1: An assessment of the computational efficiency of the Discontinuous Finite Element Method (DFEM) heat transfer solver shows good scalability. The DFEM solver is a factor of 4 more expensive in solution time than the Finite Element Method (FEM) solver for heat transfer problems due to the increased number of degrees of freedom. Nonetheless, the DFEM approach provides the user with the flexibility to model gap heat transfer problems.; Extended scope part 2: The GapHeatTransferInterfaceMaterial was improved to give the user increased flexibility with the modeling of heat transfer through gaps with the DFEM solver. A number of gap parameters can now be coupled both through functions and variables.; Extended scope part 3: Demonstration of how the gap width between hexagonal fuel cells can be calculated during a heat-up transient and used in the GapHeatTransferInterface model. A full-domain DFEM model with gap expansion is coupled to a SubApp that models the thermal expansion of the base plate. The results show the expected physical behavior, although have not been fully bench-marked at this point in time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment (FY2020--Part 3)

This report details progress and activities of Idaho National Laboratory (INL) on the NRC project ”Development and Modeling Support for Advanced Non-Light Water Reactors.” The tasks completed for this report are as follows: First, Task 2c: Explicit modeling of pebble transient temperature response. In this simulation, the PBMR-400 reactor experiences a 20 second power ramp from 100% to 150% power. This is followed by a similar reduction in the power back to 100%. Several multiscale pebble coupling approaches are tested with one pebble per mesh element in the active core region. The results show good conservation behavior and the stability of the coupling. Next, Extended scope part 1: An assessment of the computational efficiency of the DFEM heat transfer solver shows good scalability. The DFEM solver is a factor of 4 more expensive in solution time than the FEM solver for heat transfer problems due to the increased number of degrees of freedom. Nonetheless, the DFEM approach provides the user with the flexibility to model gap heat transfer problems. Then Extended scope part 2: the GapHeatTransferInterfaceMaterial was improved to give the user increased flexibility with the modeling of heat transfer through gaps with the DFEM solver. A number of gap parameters can now be coupled both through functions and variables. Finally, Extended scope part 3: demonstration of how the gap width between hexagonal fuel cells can be calculated during a heat up transient and used in the GapHeat- TransferInterface model. A full domain DFEM model with gap expansion is coupled to a SubApp that models the thermal expansion of the base plate. The results show the expected physical behavior, although have not been fully bench- marked at this point in time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

LDRD-DR: Accelerating Nuclear Fuel Qualification through Integrated Multiscale and Multiphysics Models [Slides]

Nuclear fuel qualification implies ensuring that a new fuel can safely operate in a reactor under both normal and accident conditions. Historically, extensive irradiation campaigns probe fuel performance relationships and interconnections in an integral sense, which is very time-consuming and costly. The presentation focus is on introducing advanced modeling and simulation into the qualification process in order to reduce the qualification timeline by targeting fewer integral tests

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

PyJMAK: An Open-Source Python Toolkit for Modeling Solid-State Metallurgical Phase Transformations

Accurate prediction of metallurgical phase transformations is an essential basis for autonomous optimization and rapid part qualification. Several methods can be used to estimate the evolution of phase fractions such as JMAK kinetics-based models, phase-field models, thermodynamic models, and data-driven machine learning models. Thermodynamic and phase-field-based methodologies solve multiphysics equations requiring numerous calibration parameters and significant computational resources. As a result, the computation domain is limited to a point or on order of micron-meters. The data-driven models rely on large datasets from experiments and simulations. While the JMAK model only provides information about phase fraction evolution, it can predict this evolution in near real-time using thermal history and thermodynamic data without restriction on the domain. JMAK models have been popularly used by researchers to model phase transformations occuring during additive manufacturing or over arbitrary temperature profiles. Commercial proprietary software such as Abaqus and Ansys or closed-source in-house implementations offer the ability to model JMAK based kinetics to predict phase transformation. However, these software packages are not open-source or freely available for use and development in conjunction with manufacturing machines, sensors, and machine learning algorithms. In addition, the use of the model is restricted by a license token. In contrast, given temperature profiles at multiple points in the domain, this Python-based PyJMAK model can compute phase evolution in parallel due to its stand-alone modular, voxel-based structure, and it can be executed on high-performance computing resources without any license restrictions.

Prabhune, Bhagya [Oak Ridge National Laboratory (O↗
DOE OSTI2021

Tissue scale agent-based simulation of premalignant progressions in Barrett’s esophagus

Barrett’s esophagus (BE) is a benign condition of the distal esophagus that initiates a multistage pathway to esophageal adenocarcinoma (EAC). Short of frequent intrusive (and costly) surveillance, effective screening for neoplasia in BE populations is yet to be established since progressors are rare and virtually undetectable without routine biopsies, which often sample only a small portion of the BE tissue. As a result, reliable estimation of the true prevalence of dysplasia in a BE population and evidence-based optimization of screening for at-risk individuals is challenging. Data-driven microsimulations, i.e., model-generated instances of disease history in a predefined virtual population, have found utility in the EAC screening literature as low-overhead alternatives to real-world hypothesis testing of optimal interventions for dysplasia. Despite the successes, computational limitations, paucity of knowledge and data on Barrett’s dysplasia, and the complexities of disease progression as a multiscale multiphysics process have hindered the treatment of disease progression in BE as a spatial process. Agent-based modeling of nucleation and proliferation processes in dysplasia warrants exploration in this context as an approximation that operates at a trade-off between computational tractability and precise representation of the composition and physics of the substrate (tissue). In this study, we describe spatially resolved simulations of premalignant progression toward EAC in a coarse-grained model of Barrett’s tissue that resolves the metaplastic tissue at a length scale of 0.42 mm (~3300 crypts/mm 2 ). Finally, the model is calibrated to reproduce historical high-grade dysplasia prevalence when model-generated patients are screened using the Seattle protocol.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Safety Applications Project FY20-21 Highlights, FY2022 Goals [Slides]

Collaborations leverage the investment and expertise of others. Code capabilities determine the questions that can be addressed using SAP tools. Infrastructure acts as a force multiplier for software developer skills. Improved code performance allows simulations to finish faster or with greater fidelity. Code performance has dramatically improved over the past ~4 years.

42 ENGINEERING↗