Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiphysics modelling”

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 37 records · Page 2

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.↗

Multiphysics Modeling of Hydrokinetic Turbine Energy Conversion System

Hydrokinetic turbines (HKTs) hold great promise as a renewable energy source, but high maintenance costs and limited energy output hinder their widespread adoption. The lack of comprehensive research on HKT drivetrain designs creates a knowledge gap in enhancing generation efficiency and cost reduction. A model of the HKT system with a focus on electric drivetrains and power converters is required to address this knowledge gap. This paper first introduces a MATLAB-averaged model integrating electrical-mechanical-thermal domains and aging behaviors within multi-time frames. Then a PLECS model, which incorporates maximum power point tracking and dq reference framed control for AC-DC-AC power converters, is enriched by a dynamic thermal model to predict fast transients accurately. These models optimize the design at the component level, resulting in improved integrated system performance. Furthermore, the validation of the models is carried out through hardware experiments for the averaged model and hardware-in-the-loop testing for the dynamic model.

ADVANCED PROPULSION SYSTEMS,HYDRO ENERGY↗

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↗

Isopod (inverse Optimization And Design): Moose-based Application For Performing Inverse Optimization On Multiphysics Modeling And Simulation

ISOPOD (InverSe OPtimizatiOn and Design) will be a MOOSE (Multiphysics Object Oriented Simulation Environment)-based framework which will iteratively estimate the parameters by minimizing the difference between simulated and experimental observables. While inverse optimization has long been used for parameter estimation, the unique feature of ISOPOD is rapid solution of new inverse problems, leading to accelerated Research, Development and Demonstration (RD&D) for complex parameter estimation problems. This will be achieved by extending the software structure of MOOSE to inverse optimization, as well as building on the unique automatic differentiation feature of MOOSE.

Munday, LynnB↗

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↗

Multiphysics modeling of printed surface acoustic wave thermometers

Surface acoustic wave (SAW) transducers are a class of sensors and actuators that operate on the fundamental principle of piezoelectricity. Aerosol jet printing and other additive manufacturing techniques have resulted in the low-cost fabrication of low-power and small-footprint SAW devices that are suitable for sensing in high-temperature and radioactive environments. In this work, we developed a series of temperature-dependent finite element models for a SAW transducer consisting of printed silver interdigitated transducers (IDTs) deposited onto piezoelectric lithium niobate. Modeling accuracy was evaluated experimentally from room temperature to 200 °C using an aerosol-jet-printed SAW thermometer. A time-domain study enabled visualization of the wave propagation and successfully guided the denoising of the scattering parameter measurement. Additionally, frequency-domain models using traditional modal analysis or the unique port boundary condition feature in COMSOL Multiphysics accurately predicted the temperature-driven natural frequency drift in the SAW thermometer. The finite element models developed in this study serve to facilitate the computer-aided design of future SAW transducers for applications in harsh environments.

47 OTHER INSTRUMENTATION↗

Multiphysics Modeling of Impurity Transport for FNSF Startup Scenario with ERO2.0

This paper focuses on performing a multiphysics study using the ERO2.0 and UEDGE codes for two standard double null configurations for the Fusion Nuclear Science Facility: (a) 100% recycling and (b) 99% recycling. Results show that the main contributor to tungsten erosion along the divertor plates is impurities from the midplane waveguides. In addition, the standard high-recycling case (100% recycling) shows a significantly higher buildup of impurities along the divertor tiles during the startup phase, which can lead to a higher increase of energy loss in the plasma during steady-state operation. Last, for high recycling, anomalous diffusion can dominate over parallel field diffusion. The work performed in this study can be iteratively applied to a full operation scenario with additional physics such as those from neutrals, wall shaping, and additional external fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphysics modelling in PyLith: poroelasticity

SUMMARY PyLith, a community, open-source code for modelling quasi-static and dynamic crustal deformation with an emphasis on earthquake faulting, has recently been updated with a flexible multiphysics implementation. We demonstrate the versatility of the multiphysics implementation by extending the code to model fully coupled continuum poromechanics. We verify the newly incorporated physics using standard benchmarks for a porous medium saturated with a slightly compressible fluid. The benchmarks include the 1-D consolidation problem as outlined by Terzaghi, Mandel’s problem for the 2-D case, and Cryer’s problem for the 3-D case. All three benchmarks have been added to the PyLith continuous integration test suite. We compare the closed form analytical solution for each benchmark against solutions generated by our updated code, and lastly, demonstrate that the poroelastic material formulation may be used alongside the existing fault implementation in PyLith.

Geochemistry & Geophysics↗

Recent Advanced Reactor Multiphysics Model Highlights in the Virtual Test Bed (VTB)

The Virtual Test Bed (VTB) host over 30 distinct simulations that showcase state-of-the art capabilities across the national lab complex. An update on the status of models on the VTB is summarized here, along with a more detailed overview of select recent new capabilities to showcase. All of the major advanced reactor types are represented in the VTB. The first example consists of a multiphysics simulation to track the transport of species in Molten Salt Reactors using depletion, advection, and thermochemical calculations. The second consists of a coupled neutronic and thermal hydraulic simulation to validate a gas cooled reactor. The third consist of pebble-bed equilibrium model for a fluoride high-temperature reactor. The fourth is a high-fidelity neutronic and thermal hydraulic model of a liquid metal reactor assembly. And lastly the fifth consists of transient multiphysics simulations of heat pipe microreactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimentally validated multiphysics modeling of fracture induced by thermal shocks in sintered UO 2 pellets

Uranium Dioxide (UO 2 ) fuel powers almost all commercial Nuclear Power Plants (NPPs) worldwide, generating carbon-free energy and contributing to the fight against climate change. UO 2 fuel incurs damage and fractures due to large thermal gradients that develop across the fuel pellet during normal and transient operating conditions. A comprehensive understanding of the underlying mechanisms by which these processes take place is still lacking. A combined experimental and computational approach is utilized here to quantify the behavior of UO 2 fuel fracture induced by thermal shock. Here, this work introduces both (1) an experimental study to understand the fuel fracturing behavior of sintered UO 2 pellets when exposed to thermal shock, and (2) a Multiphysics phase-field fracture model capable of simulating this process. Parametric studies were conducted to evaluate the effects of uncertainties in fracture properties on the fracture behavior of UO 2 due to thermal shocking. A set of energy release rate (or equivalently fracture toughness) and contract area (the part of the fuel pellet in direct contact with the cold bath) were able to capture the overall fracture trends of the corresponding experimental data. Our combined approach presents a new method for accounting for the effects of microstructure and sample size on the energy release rate/fracture toughness. The experimental data were collected from multiple experiments that exposed UO 2 pellets to high-temperature conditions (589–676 °C) followed by a quench in sub-zero water. This work demonstrates that joint experimental and computational efforts are able to advance the understanding of thermal fracture in the primary fuel source for existing and future NPPs.

36 MATERIALS SCIENCE↗

Release of a High Temperature Engineering Test Reactor (HTTR) Steady-State Multiphysics Model to the Virtual Test Bed

The National Reactor Innovation Center (NRIC) accelerates the demonstration and deployment of advanced nuclear energy through inspiring stakeholders and the public, empowering innovators, and delivering successful outcomes. The Virtual Test Bed (VTB) supports this mission by providing an open-source repository of advanced test reactor models for industry, academia, and the public to utilize. During my time at Idaho National Laboratory (INL), I: • Verified the input files of the High Temperature Engineering Test Reactor • Generated computational results • Created user documentation for the HTTR • Uploaded a steady-state HTTR model to the Virtual Test Bed • Authored a conference publication on the model

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Release of a High Temperature Engineering Test Reactor (HTTR) Steady State Multiphysics Model to the Virtual Test Bed

In response to climate change, global governments and private industry have established a common goal of achieving net-zero emissions by 2050 \cite{osti_1865910}. This goal requires a reassessment of current energy demands and production methods. Reducing emissions at an affordable cost while maintaining grid reliability requires a nationwide collaborative effort among government and industry in the United States. Nuclear power is the leading low-carbon electricity generation method. In the past 50 years, the use of nuclear power has reduced carbon dioxide emissions by over 60 gigatons and has played a crucial role in the security of energy supply~\cite{IEA}. In the U.S., nuclear power accounts for 20\% of the electrical supply and provides energy reliably. Advanced reactors will operate at higher temperatures, operate more efficiently, utilize more energy stored within fuel, and reduce the amount of waste produced \cite{osti_1616270}. To face these challenges and goals, the U.S. Department of Energy has created an initiative to focus on the modeling and simulation tools to support future nuclear power plant design, licensing, and operations. The Virtual Test Bed (VTB)~\cite{vtb2023} was launched by the National Reactor Innovation Center (NRIC) in collaboration with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to support the advanced nuclear reactor community. The VTB involves teams from both Idaho National Laboratory and Argonne National Laboratory and aims to provide example models for a broad range of both current and future advanced reactor designs. A feature of the VTB is the automatic testing of these models to ensure continued functionality as simulation tools are further developed. The VTB and the advanced reactor models documented there are important resources for this initiative. This work describes the inclusion of a new model on the VTB---a High Temperature Engineering Test Reactor (HTTR) steady-state model \cite{LABOURE2023109838}.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multiphysics modeling of SOFC performance degradation caused by interface delamination and active layer cracking

In solid oxide device operation, delamination and cracking cause more catastrophic damage to cell life than other degradation mechanisms, leading to more rapid cell failure. The electrochemical performance degradation caused by electrode/electrolyte interface delamination and active layer cracking were investigated through in-house developed three-dimensional multiphysics simulations. Here the multiphysics simulations, fully calibrated against experimental datasets for button cells in a previous study, were extended to quantitatively investigate the effects of electrode/electrolyte interface delamination and active layer cracking. The simulated polarization curves and impedance behavior with various levels of delamination and cracking provided insights of the performance degradation. Furthermore, the local current distributions within the cell were analyzed for three different operating modes: constant voltage, constant current, and constant power. Finally, a practical strategy was proposed to mitigate the performance degradation in the realistic cell operations, extending overall cell lifetime.

08 HYDROGEN↗

Multiphysics Modeling Framework for Soft PVC Gel Sensors with Experimental Comparisons

Polyvinyl chloride (PVC) gels have recently been found to exhibit mechanoelectrical transduction or sensing capabilities under compressive loading applications. This phenomenon is not wholly understood but has been characterized as an adsorption-like phenomena under varying amounts and types of plasticizers. A different polymer lattice structure has also been tested, thermoplastic polyurethane, which showed similar sensing characteristics. This study examines mechanical and electrical properties of these gel sensors and proposes a mathematical framework of the underlying mechanisms of mechanoelectrical transduction. COMSOL Multiphysics is used to show solid mechanics characteristics, electrostatic properties, and transport of interstitial plasticizer under compressive loading applications. The solid mechanics takes a continuum mechanics approach and includes a highly compressive Storakers material model for compressive loading applications. The electrostatics and transport properties include charge conservation and a Langmuir adsorption migration model with variable diffusion properties based on plasticizer properties. Results show both plasticizer concentration gradient as well as expected voltage response under varying amounts and types of plasticizers. Experimental work is also completed to show agreeance with the modeling results.

42 ENGINEERING↗

Multiphysics modeling of a critical dual-purpose canister in a saturated geological repository

Dual purpose dry cask storage canisters for spent nuclear fuel are designed for storage and transportation, but are not licensed for permanent disposal in a geological repository. If dual purpose canisters were to be used to dispose of spent nuclear fuel in a geological repository, they would be expected to eventually breach and be flooded with groundwater, and it is shown that some fraction of these canisters will achieve criticality. In order to evaluate the consequences of canisters going critical in a repository, an initial capability has been developed for estimating the quasi-static power level of a critical canister using loosely coupled multiphysics simulations. The low power level in a critical canister enables coupling through precomputed physics proxies. This calculated power level is then used to compute the change in the critical canister’s isotopic inventory as a function of time. Three as-loaded canisters are evaluated and two were found to have power levels below 4 kW, with a modest effect on the radiological inventory over time. This effort also shows that although some DPCs will have extremely peaked power shapes, the relatively low power and long-time scales result in relatively homogenized thermohydraulic properties in the water within the DPC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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↗