Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiphysics modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

MOOSE-Workbench integration and MOOSE meshing capability enhancements to facilitate inputs and outputs for multiphysics modeling

The Multiphysics Object-Oriented Simulation Environment (MOOSE) is an open-source framework that supports many of the US Department of Energy’s (DOE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) technical areas (TA). These TAs develop and use NEAMS physics and coupling modules in multiple ways to enable the research and development of complex physics models. In addition to the MOOSE framework, the NEAMS Workbench user interface provides a common analysis environment with user-interaction accelerators that streamline the tasks of model creation, review, execution, and out put inspection. In FY 2024, objectives were realized in the MOOSE framework application development support and user-oriented improvements. Application development improvements support both developers and users with an expanded Reactor Module and Mesh System, stateful material property support for mortar contact, and customizable convergence criteria. Additionally, new user-oriented features were implemented in the MOOSE framework language server, including autocompletion snippets, definition from source and find reference navigations, and syntax overrides. Lastly, improvements were made to the input interpreter necessary to support the MOOSE language server and the NEAMS Workbench so that they can interact with syntactically incomplete user inputs. These improvements and more were intended to address stakeholder feedback and improve developer and user ability to conduct advanced nuclear energy modeling and simulation in support of DOE and industry needs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

MOOSE–Workbench Integration and MOOSE Meshing Capability Enhancements to Facilitate Inputs and Outputs for Multiphysics Modeling

The Multiphysics Object-Oriented Simulation Environment (MOOSE) is an open-source framework that supports many of the US Department of Energy’s (DOE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) technical areas (TA). These TAs develop and use NEAMS physics and coupling modules in multiple ways to enable the research and development of complex physics models. In addition to the MOOSE framework, the NEAMS Workbench user interface provides a common analysis environment with user-interaction accelerators that streamline the tasks of model creation, review, execution, and output inspection. In FY 2024, objectives were realized in the MOOSE framework application development support and user-oriented improvements. Application development improvements support both developers and users with an expanded Reactor Module and Mesh System, stateful material property support for mortar contact, and customizable convergence criteria. Additionally, new user-oriented features were implemented in the MOOSE framework language server, including autocompletion snippets, definition from source and find reference navigations, and syntax overrides. Lastly, improvements were made to the input interpreter necessary to support the MOOSE language server and the NEAMS Workbench so that they can interact with syntactically incomplete user inputs. These improvements and more were intended to address stakeholder feedback and improve developer and user ability to conduct advanced nuclear energy modeling and simulation in support of DOE and industry needs.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Block-Structured Operator Inference for Coupled Multiphysics Model Reduction

This work presents a block-structured formulation of Operator Inference as a way to learn structured reduced-order models for multiphysics systems. The approach specifies the governing equation structure for each physics component and the structure of the coupling terms. Once the multiphysics structure is specified, the reduced-order model is learned from snapshot data following the nonintrusive Operator Inference methodology. In addition to preserving physical system structure, which in turn permits preservation of system properties such as stability and second-order structure, the block-structured approach has the advantages of reducing the overall dimensionality of the learning problem and admitting tailored regularization for each physics component. The numerical advantages of the block-structured formulation over a monolithic Operator Inference formulation are demonstrated for aeroelastic analysis, which couples aerodynamic and structural models. For the benchmark test case of the AGARD 445.6 wing, block-structured Operator Inference provides an average 20% online prediction speedup over monolithic Operator Inference across subsonic and supersonic flow conditions in both the stable and fluttering parameter regimes while preserving the accuracy achieved with monolithic Operator Inference.

42 ENGINEERING↗
DOE OSTI2024

High-Fidelity Multiphysics Modeling of a Heat Pipe Microreactor Using BlueCrab

Researchers who are actively developing nuclear microreactors are planning to employ innovative designs and features using traditional commercial modeling tools that may be inadequate for their design and licensing activities. The codes developed under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program provide flexibility in terms of geometry modeling and multiphysics coupling and are particularly well suited for modeling novel microreactor concepts. To test the maturity of these codes, this paper introduces a conceptual heat pipe microreactor (HP-MR) designed to gather various technologies of interest to microreactor developers such as control drums, heat pipes, and hydride moderators. Here, the objective of this effort is to demonstrate NEAMS tools capability to perform high-fidelity multiphysics simulations, using coupled neutronics (via the Griffin code), heat conduction (via the BISON code), heat pipe modeling (via the Sockeye code), and hydrogen redistribution in hydride metal moderator (via the SWIFT code). Codes are coupled in-memory through the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, which permits flexible multiphysics data transfer schemes. The analysis confirmed two key aspects of the HP-MR concept: (1) its ability to follow the power load requested from the heat pipe and (2) its ability to avoid heat pipe cascading failure unless designed with high power close to operating failure limits of its heat pipes. The developed computational model was distributed publicly on the Virtual Test Bed for training purposes to accelerate adoption by industry and to provide a high-fidelity multiphysics solution for benchmarking against other tools. Additional multiphysics analyses including other transients and coupled physics were identified as necessary future work, together with a focus on validating multiphysics behavior against experiments.

Microreactor↗
NASA NTRSDate not supplied

Comprehensive 3D Multiphysics Model on Electrochemical Recovery of O 2 from Metabolic CO 2 at the International Space Station (ISS)

The International Space Station (ISS) is presently equipped with an elaborate, heavy, and high-power consuming system that recovers approximately 50% of O 2 from metabolic CO 2 as part of the atmospheric revitalization (AR) at the ISS habitat. Future long-duration missions will require a sustainable and efficient system capable of yielding a minimum of 75% O 2 recovery to reach the self-sufficiency required for long space missions beyond earth’s low orbit. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to not only increase significantly current O 2 recovery efficiency, improving self-sufficiency on AR at the ISS habitat and future long missions, but also reduce the complexity of the system. The authors have developed and deployed a comprehensive 3D multiphysics model that thoroughly replicates the actual configuration and fluid/material domains of the MFECR. The coupled physics in this multiphysics model include multicomponent-multiphase electrochemical-driven reactions, non-ideal mass transport mechanism, free and porous flow, heat transfer, CO 2 solubility on alkaline electrolyte, water condensation on porous medium, and DC electrical current generation along with Joule heating effect. This model is aimed to conduct quantitive benchmark on three different MFECR’s layouts, one without serpentine paths (plain) and two with serpentines leading to four and twelve paths respectively. Once experimental data is generated via a test matrix of 200 tests, the model will be validated to conduct MFECR’s process optimization and revalidate the quantitive benchmark on three different MFECR’s layouts.

Jesus A. Dominguez↗
DOE OSTI2025

Summary of Multiphysics Modeling for Sublimation Mass Transfer

This report summarizes the effort to develop a multiphysics modeling framework for sublimation mass transfer processes. The aim of the project was to develop a multiphysics code capable of predicting phase change between solid and gas in a closed container, and the movement of material within such a container when exposed to various exterior environmental conditions. Two modeling frameworks were developed towards this aim: one with a high-fidelity computational fluid dynamics (CFD) structure, and the other as a fast reduced-order model. Ultimately, both attempts were unsuccessful due to instabilities in the code and physical processes for which there is no adequate numerical representation. Both modeling attempts are briefly detailed before providing a brief survey of recent updates in the literature, and a recommendation for future work on this subject.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

3D Multiphysics Model for Large Scale Planar Cell: Inductance Investigation in Impedance Analysis

Developed a comprehensive 3D multiphysics model to analyze impedance behavior in large-scale planar Solid Oxide Fuel Cells (SOFCs). The model investigates how impedance varies across operating conditions and cell components, with a particular focus on inductive loops or negative capacitance at the air electrode. It was found that the Inductive loop disappears in the low-frequency zone when the activation overpotential or faradic current is not temperature-dependent. These findings offer new insights into SOFC impedance behavior.

impedence analysis↗
DOE OSTI2024

Multiphysics modeling of accelerators through code integration

This work aims to improve the ability of particle accelerator researchers to develop high-performance accelerator cavity designs by creating an overall multiphysics framework that integrates and couples existing application codes. This framework will allow accelerator researchers to build multiphysics models that will optimize cavity design, improve understanding of whole-device performance, and reduce the development and fabrication costs of accelerator research. We utilize the open-source VizSchema data standard as an intermediate data structure interface layer to standardize interfaces between individual application codes. VizScema is extensively documented online, and plugins for VizSchema are available for popular visualization packages, including VisIt and ParaView. Currently, the work focuses on coupling the EM field solver COMSOL and the electron gun code MICHELLE to allow COMSOL field-solve results to be seamlessly used by MICHELLE for particle-solve. Later work will extend this integration to include other fields, particles, and thermodynamics simulation codes.

43 PARTICLE ACCELERATORS↗
NASA NTRS2010

Multiphysics Modeling for Dimensional Analysis of a Self-Heated Molten Regolith Electrolysis Reactor for Oxygen and Metals Production on the Moon and Mars

The technology of direct electrolysis of molten lunar regolith to produce oxygen and molten metal alloys has progressed greatly in the last few years. The development of long-lasting inert anodes and cathode designs as well as techniques for the removal of molten products from the reactor has been demonstrated. The containment of chemically aggressive oxide and metal melts is very difficult at the operating temperatures ca 1600 C. Containing the molten oxides in a regolith shell can solve this technical issue and can be achieved by designing a self-heating reactor in which the electrolytic currents generate enough Joule heat to create a molten bath. In a first phase, a thermal analysis model was built to study the formation of a melt of lunar basaltic regolith irradiated by a focused solar beam This mode of heating was selected because it relies on radiative heat transfer, which is the dominant mode of transfer of energy in melts at 1600 C. Knowing and setting the Gaussian-type heat flux from the concentrated solar beam and the phase and temperature dependent thermal properties, the model predicts the dimensions and temperature profile of the melt. A validation of the model is presented in this paper through the experimental formation of a spherical cap melt realized by others. The Orbitec/PSI experimental setup uses an 3.6-cm diameter concentrated solar beam to create a hemispheric melt in a bed of lunar regolith simulant contained in a large pot. Upon cooling, the dimensions of the vitrified melt are measured to validate the thermal model. In a second phase, the model is augmented by multiphysics components to compute the passage of electrical currents between electrodes inserted in the molten regolith. The current through the melt generates Joule heating due to the high resistivity of the medium and this energy is transferred into the melt by conduction, convection and primarily by radiation. The model faces challenges in two major areas, the change of phase as temperature increases, and the dominance of radiative heat flux as heat transfer mechanism within the melt the change of phase concerns the regolith itself which is present in states ranging from a fine grain regolith with low thermal conductivity and low density to a vitrified melt with much higher thermal conductivity, and higher density. As the regolith is heated, it starts to soften around 1300 C the melt iS very viscous and evolving gas bubbles out in thick, lava-like fashion. By 1600 C the regolith is completely melted and the viscosity is low The second challenge resides in the proper modeling of the radiative heat flux requiring the addition of the computing-demanding radiative-heat-transfer function to the general heat transfer equation. The model Includes temperature-dependent properties (density, thermal conductivity, heat capacity, and viscosity, and absorption coefficients) and solves the radiative heat flux equation assuming gray (fine grains) and semi-transparent (melt) media and using an absorption coefficient spectral found in the literature for terrestrial minerals similar in composition to those of lunar regolith simulant

Dominguez, Jesus A.↗
NASA NTRS2023

Multiphysics Modeling of Electrode Dissolution under Hydrodynamic Conditions in Low Conductivity Water

Efficient plug-flow electrochemical reactors are characterized by their ability to effectively convert species with significant ionic strength and/or with the assistance of a well-behaved electrolyte. Because of the classic electroconductivity requirement, electrochemical reactors are not particularly attractive options for applications with low-conductivity conditions, especially for low-pressure feeds. However, these reactors can still find application in dedicated processes such as ultra-clean water conditioning and biological-coupled separation systems. Consequently, a multi-physics model was developed using COMSOL Multiphysics to characterize different reactor-design conditions based on the Tertiary Current Distribution which accounts for the effects of expected variations in electrolyte composition and ionic strength on the electrochemical process, as well as solution resistance and electrode kinetics. The modelling framework presented here employs single-phase laminar fluid flow, the Nernst-Planck equation, the water-based electroneutrality condition, and concentration-dependent overpotentials. This investigation particularly explores the electrolysis process of “anodic dissolution,” in which elemental species dissociate from the solid matrix of the electrode and enter the liquid phase of the electrolyte as soluble ions. The reactor configuration of interest has a rectangular parallel-electrode design and is operated galvanostatically. A series of parametric studies are carried out to inspect the response of the system when hydrodynamic, kinetic, and geometric conditions are changed. The parametric sweep resulted in variable effluent concentrations and system voltage, which reveal underlaying optimization patterns for electrode dissolution in low conductivity water.

Multiphysics↗
NASA NTRS2014

Multiphysics Modeling of an Annular Linear Induction Pump With Applications to Space Nuclear Power Systems

An annular linear induction pump (ALIP) that could be used for circulating liquid-metal coolant in a fission surface power reactor system is modeled in the present work using the computational COMSOL Multiphysics package. The pump is modeled using a two-dimensional, axisymmetric geometry and solved under conditions similar to those used during experimental pump testing. Real, nonlinear, temperature-dependent material properties can be incorporated into the model for both the electrically-conducting working fluid in the pump (NaK-78) and structural components of the pump. The intricate three-phase coil configuration of the pump is implemented in the model to produce an axially-traveling magnetic wave that is qualitatively similar to the measured magnetic wave. The model qualitatively captures the expected feature of a peak in efficiency as a function of flow rate.

Kilbane, J.↗
NASA NTRS2015

Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials

The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition, after investigating various methods, a Smoothed Particle Hydrodynamics Model (SPH Model) was developed to model wire feeding process. Its computational efficiency and simple architecture makes it more robust and flexible than other models. More research on material properties may be needed to realistically model the AAM processes. A microscale model was developed to investigate heterogeneous nucleation, dendritic grain growth, epitaxial growth of columnar grains, columnar-to-equiaxed transition, grain transport in melt, and other properties. The orientations of the columnar grains were almost perpendicular to the laser motion's direction. Compared to the similar studies in the literature, the multiple grain morphology modeling result is in the same order of magnitude as optical morphologies in the experiment. Experimental work was conducted to validate different models. An infrared camera was incorporated as a process monitoring and validating tool to identify the solidus and mushy zones during deposition. The images were successfully processed to identify these regions. This research project has investigated multiscale and multiphysics of the complex AAM processes thus leading to advanced understanding of these processes. The project has also developed several modeling tools and experimental validation tools that will be very critical in the future of AAM process qualification and certification.

Liou, Frank↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Serpent - Bison - THM Preliminary Multiphysics Modeling of a Nuclear Thermal Propulsion System Fuel Assembly

This work demonstrates the Monte Carlo neutronic and Thermo-Hydraulic coupling scheme using the Serpent code and MOOSE application Bison and Thermo Hydraulic Module. The coupling scheme is then applied to the reference BWX Technologies Nuclear Ther- mal Propulsion system at he fuel assembly level where it’s used to perform an analysis of the isothermal material coefficients and potential material reactivity worth. A method is developed to isolate which feedback effects should be considered for proceeding with reduced order deterministic neutronic modeling where branch off analysis must be con- ducted. The convergence behavior of the coupling scheme is demonstrated where it fol- lows the standard Picard iteration approach. Verification studies for the method of deduc- ing relevant feedback effects is also demonstrated.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computational multiphysics modeling of radioactive aerosol deposition in diverse human respiratory tract geometries

The evaluation of aerosol exposure relies on generic mathematical models that assume uniform particle deposition profiles over the human respiratory tract and do not account for subject-specific characteristics. Here we introduce a hybrid-automated computational workflow that generates personalized particle deposition profiles in 3D reconstructed human airways from computed tomography scans using Computational Fluid and Particle Dynamics simulations. This is the first large-scale study to consider realistic airways variability, where 380 lower and 40 upper human respiratory tract 3D geometries are reconstructed and parameterized. The data is clustered into nine groups using random forest regression. Computational fluid and particle dynamics simulations are conducted on these representative geometries using a realistic heavy-breathing respiratory cycle and radioactive iodine-131 as a source term. Monte Carlo radiation transport simulations are performed to obtain detailed energy deposition maps. Our findings emphasize the importance of personalized studies, as minor respiratory tract variations notably influence deposition patterns rather than global parameters of the lower airways, observing more than 30% variance in the mass deposition fraction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗