SEARCH · Search NASA

Results for “Fluid System”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Symmetry Determining Equations of the Euler Equations

This purpose of this document is to provide a clear derivation of the symmetry determining equations (SDEs) of the Euler equations to supplement the lecture document. The Euler equations are a coupled system of three first-order quasilinear hyperbolic partial differential equations (PDEs). The Euler equations are a system of partial differential equations (PDEs) for fluid velocity v(t, r), fluid density ρ(t, r), and fluid pressure P(t, r), where t is time and r is the spatial coordinate.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

CFD modeling and simulation for corrosive wear of refractory in molten slag

The development of high-wear resistant refractories having minimal production costs is facilitated by characterizing the wear mechanisms associated with their corrosive wear. Static cup testing is a commonly used method for comparing the corrosion resistance performance of two or more refractory materials. Although the static cup test conditions are not as severe as dynamic tests, this study shows that the thermal gradient present within the system during heating and cooling stages serves to generate movement of the slag leading to mechanical wear. The thermal gradient within the refractory, and between the slag and the refractory, occurs during the ramp stage of the test and lasts until the soaking stage is reached bringing the system to a thermal equilibrium. Using computational fluid dynamics (CFD) capabilities embedded within ANSYS software, this study modelled and quantified the convection currents within the slag and associated shear stresses generated on the refractory walls due to the thermal gradient. A traditional ladle furnace was employed as a case study to verify the results of the studied CFD model. The corrosion rate of the refractory lining was found to depend on the mass transfer coefficient of the refractory dissolution into the slag, and a velocity term which governs the extent of corrosion at any given location. This velocity term is a function of slag viscosity, as well as the concentration gradient and/or temperature gradient at the triple points. In this study, wall shear stress was used as a reliable proxy for identifying high-velocity regions prone to excessive corrosive wear. Elevated wall shear stress near the slag/air and slag/molten steel interfaces align with observed corrosion grooves, which reflects the intensified corrosive wear at these locations.

Ramteke, Rajat Rajat Durgesh Ramteke [University o↗

DOE OSTI2025

Numerical assessment of triply periodic minimal surfaces for direct air capture of carbon dioxide

Direct air capture (DAC) systems often consist of packing material wetted by a capture fluid that reacts with CO 2 in the airstream. The efficiency of the contactor is determined by a complex relationship of fluid dynamics, heat and mass transfer, contactor geometry, and chemical properties. The efficiency of the contactor must be balanced with other factors, primarily pressure drop through the system. Triply periodic minimal surfaces (TPMS) are a class of differential surfaces that have been explored in multiple engineering applications and have been shown to exhibit excellent performance when used in heat exchangers. Their tortuous path provides a high surface-to-volume ratio and favorable trade-off between contact area and pressure drop. In this work, a gyroid-type TPMS contactor was evaluated using computational fluid dynamics for a variety of geometric parameters to explore the potential benefit of TPMS shapes for DAC applications. A thin-film model was employed to model the flow and distribution of the capture solvent, allowing efficient simulations of TPMS structures at scale by eliminating the need for a computationally intensive interface capturing method. A liquid-gas mass transfer model was implemented in the commercial software STAR-CCM+ and used to predict the CO 2 capture efficiency and study the trade-off between capture performance and pressure drop through analysis of capture rates, mass transfer coefficients, and other relevant variables. TPMS contactors with a variety of geometric parameters and two capture solvent options were investigated to determine the effect of design choices on the operational performance of DAC systems. In conclusion, results showed that while contactor geometry is the dominant factor in efficiency and pressure drop, the physiochemical properties of the solvent are an important secondary influence on the contactor performance.

CFD↗

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

Thermal Analysis of a Solid Particle Light-Trapping Planar Cavity Receiver Using Computational Fluid Dynamics

Concentrated solar power (CSP) is one of the most effective ways of harnessing solar power to create efficient, durable, and resilient energy systems. This study entails thermal modeling and analysis of a novel central tower receiver configuration. This receiver uses solid particles as the heat transfer fluid (HTF), a promising option for third-generation CSP systems. The configuration considered here is the light-trapping planar cavity receiver (LTPCR) introduced by the National Renewable Energy Laboratory. While heat transfer studies of various LTPCR subsystems have been done, system-level thermal analysis of the LTPCR receiver has not been attempted. This study also presents important sensitivity analyses of the operating parameters of the CSP system, which can help guide the design of future central tower receivers. This study employs Ansys Fluent as a computational fluid dynamics (CFD) tool to model fluid dynamics and heat transfer in the receiver, intending to quantify its thermal performance. The model seamlessly integrates Monte Carlo ray tracing data, which generates absorbed solar flux profiles from the heliostat field design, with the heat transfer characteristics of the fluidized particle bed. This unified model is designed to accurately predict the thermal behavior of the LTPCR. Analysis of preliminary results reveals that the primary loss mechanisms are radiative and natural convective losses, in that order. Based on observations from a baseline case, several strategies are suggested and numerically tested. These solutions include selective cooling of high-temperature regions and manipulation of particle bed parameters. Selective cooling of high-temperature regions reduced the peak temperature by 151 degrees C and decreased thermal losses by 0.9%. Improving the particle-wall heat transfer coefficient (P-W HTC) of the particle bed decreased the thermal losses by 1.7% and decreased the peak temperatures by 57 degrees C. Decreasing the particle inlet temperature (PIT) also reduced thermal losses by 3.5% and decreased peak temperatures by 29 degrees C. Compounding these strategies improved the thermal losses of the receiver from 13.5% in the baseline case to 7.5%. Additionally, the study explores the variation in thermal performance across different locations of the receiver, where a variation of thermal losses from 12.9% to 17.3% is found. This allows a comprehensive evaluation of potential improvements in efficiency and temperature management.

computational fluid dynamics↗

DOE OSTI2026

Pilot Development Progress to Demonstrate Electric Thermal Energy Storage (ETES) Using Low-Cost Particles

The rapid growth of future energy demand increases the need to economically store electrical energy over durations up to several days in long-duration energy storage (LDES) applications. LDES can reduce the grid stress for improving the resilience of the grid. The integration of renewable power and storage has several significant and positive impacts including expanding the renewable energy portion of electricity generation, meeting peak-load demand, and coordinating supply and demand. Several energy storage approaches including mechanical, chemical and electrochemical methods are currently deployed or under development. However, LDES requirements pose unique challenges for scalability, energy capacity, and cost. Thermal energy storage (TES) has the ability to store a large capacity of energy with improved site flexibility compared to incumbent LDES technologies (i.e., pumped hydro) and has attracted significant interest. Our development in TES technology using solid particles as the storage medium has shown feasibility and potential in serving LDES purposes. Particle TES has various advantages over commercialized molten nitrate salt TES and has emerged as a promising storage technology originating from Generation 3 concentrating solar power (CSP) development. Stable, inexpensive silica sand produced in the Midwestern U.S. was identified as a suitable storage medium and has been validated through lab heating and cycling tests. We have successfully developed particle TES for bidirectional electric energy storage to broaden the applications for TES beyond CSP technologies and into standalone electric-thermal energy storage (ETES) systems. Key component designs and performance were verified with prototype testing and model simulations. The ETES system uses high-temperature, low-cost particle thermal energy storage coupled with a unique pressurized fluidized bed heat exchanger that supports a high-efficiency, air-Brayton combined power cycle. At times of low energy demand or high renewable production, an electric particle heater heats large quantities of solid particles to high temperatures in excess of 1100 degrees C. The particles are stored in internally insulated silos. At times of high electricity demand, the hot particles are gravity-fed through the heat exchanger where stored thermal energy is transferred to the power cycle working fluid to drive the gas turbine power generation system, thereby converting the thermal energy back into electricity for the grid. The cold particles leaving the heat exchanger are returned to the storage silos by a well-insulated particle conveyor. This ETES technology using particle TES can be sited anywhere and is predicted to be economic and efficient. Moreover, this ETES technology could be co-located with a retired coal or gas plant to leverage existing power generation and grid infrastructure. Our work has focused on the development of each of the key components (e.g., the electric particle heater, particle storage silos, and the particle-to-gas pressurized fluidized bed heat exchanger) through conceptual designs, prototype testing, and computational modeling. These developments demonstrate technological feasibility and are an important step towards the next phase of development of this promising LDES technology. This presentation will show the progress in developing a pilot demonstration of the ETES system on a commercially relevant scale.

25 ENERGY STORAGE↗

DOE OSTI2026

Filament extension atomizers: A study of integration potential into spray dryers for high-viscosity and non-Newtonian fluids

Filament extension atomizers (FEAs) are an emerging class of spray nozzles designed to atomize high-viscosity and non-Newtonian fluids that are challenging for conventional pressure or two-fluid nozzles. In this study, we investigate the influence of roller geometry, surface velocity, roller material, and fluid rheology on the atomization performance of FEA systems using concentrated whey protein suspensions (50–70 wt. %). Extensional and shear rheology experiments, along with high-speed imaging and particle image velocimetry, reveal that filament breakup dynamics are governed by competition between inertial, capillary, and viscoelastic stresses. High roller rotational velocity leads to narrower spray cones, contradicting rheology experiments and suggests a significant inertial contribution to filament breakup. Smaller rollers operating at the same rotational velocity led to broader spray cones consistent with expectations. An FEA nozzle was integrated into a conventional dryer producing particles in the 100 μm range. Results suggest that FEA technology enables atomization of highly viscous fluids at industrially relevant spray cone angles similar to those generated by pressure nozzles, offering a pathway to improve energy efficiency in spray drying by enabling higher solids loading feedstocks. Furthermore, these insights provide critical guidance for optimizing FEA nozzle designs and process parameters across a range of applications.

Energy efficiency↗

DOE OSTI2024

Exploring the Feasibility of INCONEL® ALLOY 740H® for Power Plant Headers: Integrating Machine Learning with Computational Fluid Dynamics (CFD)

This keynote presentation explores the behavior of headers—essential components of pipeline systems—using ANSYS simulation software and machine learning techniques. The study aims to predict the thermal and mechanical performance of headers under diverse conditions through both steady-state and transient simulations. We investigate critical parameters such as heat transfer coefficient, fluid velocity, and temperature to optimize header design. Conducted as part of a DOE project led by NCAT in collaboration with UNC Charlotte, this research encompasses multiple key topics. The initial section focuses on the behavior of header systems under steady-state conditions using ANSYS simulation. It underscores the importance of headers in industrial infrastructure, especially in the energy sector, and examines the implications of material selection and flow direction on heat transfer dynamics. Methodologically, we employ Computational Fluid Dynamics (CFD) analysis through ANSYS, detailing the development of models, material properties, geometry specifications, boundary conditions, and meshing strategies. Our simulations explore various operational parameters, including temperature and mass flow rates, crucial for predicting heat transfer coefficients and enhancing header design. Results from the study include parametric investigations into mesh sensitivity, viscosity model evaluations, and the effects of heat transfer locations, all validated against theoretical calculations. We conclude with insights on mesh optimization, the suitability of viscosity models, and recommendations for future research aimed at improving header system efficiency and sustainability in industrial applications.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2026

Pumped Thermal Energy Storage Using Low-Cost Particles and a Fluid Bed Heat Exchanger for Maximum Power Efficiency

Robust, efficient, cost-effective long-duration electricity storage (LDES) solutions can enhance grid resiliency, support existing transmission and distribution infrastructure, and enable a greater share of low-cost, variable alternative energy sources to penetrate the market. To meet this need, the project team at the National Renewable Energy Laboratory developed a transformative LDES system based on pumped thermal energy storage (TES) using low-cost particles and a fluid bed heat exchanger for maximum power efficiency (PUMP). The PUMP system is composed of high-temperature, low-cost particle TES coupled with an advanced pressurized fluid bed heat exchanger (PFB HX) that supports a high-efficiency pumped thermal energy storage (PTES) system integrated with concentrating solar thermal power (CSP). The PUMP project developed and de-risked a PFB HX and particle CSP system intended to be integrated with reversible turbomachinery and a modeling tool to assess PTES cost and performance.

14 SOLAR ENERGY↗

DOE OSTI2024

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area: Preprint

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs drilled in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated including a single, double, and multi-lateral design. The subsurface within area is characterized by deep, hot (> 200degreesC) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lies within a target depth of 6 km where temperatures are estimated to approach 300degreesC. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down technoeconomic analysis model. The results show that the performance and cost optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60degreesC/km. By injecting 20 degreesC-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth can be achieved, resulting in an LCOE and LCOH of $136/MWh and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗

DOE OSTI2024

Two-Fluid and Discrete Element Modeling of a Parallel Plate Fluidized Bed Heat Exchanger for Concentrating Solar Power

A novel high-temperature particle solar receiver is developed using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a fluidized bed at the lower section, leaving the system from the bottom at a constant flowrate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built by collaborators at the Colorado School of Mines and their data are used for model validation. In this experimental setup, near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using the discrete element method and a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A new continuum conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2 K.

concentrating solar power↗

DOE OSTI2025

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated, including a single-, double-, and multi-lateral design. The subsurface within the area of interest is characterized by deep, hot (> 200 degrees C) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lie(s) within a target depth of 6 km, where temperatures are estimated to approach 300 degrees C. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot, dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down techno-economic analysis model. The results show that the performance- and cost-optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60 degrees C/km. By injecting 20 degrees C-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth (i.e., 2.2 MWe net plant output) can be achieved, resulting in an LCOE and LCOH of $136/MWhe and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗

DOE OSTI2025

Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems

The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.

Advection-dominated dynamical systems↗

DOE OSTI2025

Bringing solar to agriculture: An interdisciplinary design and analysis of a Concord grape agrivoltaic system

Agrivoltaics presents an opportunity to integrate solar photovoltaics (PV) with agricultural production, but crop-specific challenges and operational constraints remain underexplored. This study develops and evaluates a Concord grape agrivoltaic system in the Lake Erie American Viticulture Area, where vineyards face economic pressures and land use conflicts with solar development. Starting with vertical, tracking, and overhead PV systems, we model photosynthetic photon flux density (PPFD) reduction in grapevines and power generation losses from shading using the Agrivoltaic Radiation Tool (ART). Based on these results, which indicate 0.47 % annual grapevine PPFD loss for vertical designs, 1.6 % for tracking, and up to 25 % for the overhead systems, the vertical and tracking designs are selected for further computational fluid dynamics (CFD) analysis to evaluate airflow interactions. CFD results show that vertical panels do not significantly impact airflow through the grapevine canopy, and that tracking systems in horizontal position may enhance airflow compared to a vine-only scenario. Considering operational constraints for tracking systems, the vertical design is selected for an economic evaluation to reveal key financial outcomes for solar developers (14-year payback period) and growers ($408 reduction in financial losses per acre annually). A sensitivity analysis quantifies uncertainty in power generation (±8.8 %) and PPFD (±5.0 %), ensuring model robustness across different vineyard conditions. Furthermore, these findings provide quantitative evidence for the feasibility of Concord grape agrivoltaics, demonstrating a synergistic opportunity for dual-use solar while preserving cultural heritage in grape-growing regions.

14 SOLAR ENERGY↗

DOE OSTI2025

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗

DOE OSTI2025

Building a new multiphysics workflow in MOOSE: application to tritium migration, trapping and advection in TMAP8

Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Measurement of partial vapor pressures of salt mixtures via combined horizontal transpiration and thermogravimetric analysis

A method combining thermogravimetric analysis (TGA) and horizontal transpiration with elemental analysis via inductively coupled plasma mass spectrometry or ion chromatography enabled calculation of partial pressures of individual salts in molten mixtures. TGA quantified total mass loss, while transpiration identified vapor-phase composition. Furthermore, two chloride (NaCl-MgCl 2 , NaCl-MgCl 2 + UCl 3 ) salts and one mixed halide (LiCl-LiF + Li 2 O) salt were analyzed at 750 °C and 550 °C, respectively. NaCl and MgCl 2 vapor pressures were 2.19–2.61 × 10 -4 atm and 2.47–2.48 × 10 -5 atm (dependent upon the identity of the invesitgated mixture); UCl 3 was 1.42 × 10 -7 atm. LiCl and LiF vapor pressures at 550 °C were 1.53 × 10 -6 and 6.32 × 10 -6 atm, respectively. Additionally, the TGA method was validated against values from the literature for unary LiCl and LiF.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Fast permeability measurement for tight reservoir cores using only initial data of the one chamber pressure pulse decay test

Here, in this study, a mathematical model for fast determination of the permeabilities of tight rocks using measurements taken from the initial period of the One Chamber Pressure Pulse Decay (OC-PPD) test is presented. The model applies to measurements taken both before and after the pressure pulse front has reached the downstream end of the specimen. The analytical solutions for the pressure decay in the upstream chamber are derived based on a parabolic arc approximation of pore pressure distribution along the test specimen. This approximation allows converting the initial–boundary value problem of fluid diffusion in the specimen, governed by partial differential equations, to a system of ordinary differential equations that can be easily solved by explicit formulae. Thus, an explicit formula for the pressure decay rate is obtained, which enables inverse analysis of the initial experimental data to estimate the rock permeability. The proposed method expedites the pulse decay test as it does not require the system to reach equilibrium. The method is validated with three sets of experimental data of the OC-PPD test using helium as the diffusing fluid, for which the relative error of the permeability is found to be less than 6%. This method is particularly useful if the equilibrium time of the pulse decay test for rock specimens with permeabilities in the range of nano-Darcy takes hours or days.

early-time solution↗