SEARCH · Engineering Papers

Results for “mass diffusivity”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗

DOE OSTI2024

Sequential Monte Carlo for Cut-Bayesian Posterior Computation

We propose a sequential Monte Carlo (SMC) method to efficiently and accurately compute cut-Bayesian posterior quantities of interest, variations of standard Bayesian approaches constructed primarily to account for model misspecification. We prove finite sample concentration bounds for estimators derived from the proposed method along with a linear tempering extension and apply these results to a realistic setting where a computer model is misspecified. We then illustrate the SMC method for inference in a modular chemical reactor example that includes submodels for reaction kinetics, turbulence, mass transfer, and diffusion. The samples obtained are commensurate with a direct-sampling approach that consists of running multiple Markov chains, with computational efficiency gains using the SMC method. Overall, the SMC method presented yields a novel, rigorous approach to computing with cut-Bayesian posterior distributions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Efficient Reformulation of Linear and Nonlinear Solid-Phase Diffusion in Lithium-ion Battery Models using Symmetric Polynomials: Mass Conservation and Computational Efficiency

Lithium-ion batteries are typically modeled using porous electrode theory coupled with various transport and reaction mechanisms, along with suitable discretization or approximations for the solid-phase diffusion equation. The solid-phase diffusion equation represents the main computational burden for typical pseudo-2-dimensional (p2D) models since these equations in the pseudo r -dimension must be solved at each point in the computational grid. This substantially increases the complexity of the model as well as the computational time. Traditional approaches towards simplifying solid-phase diffusion possess certain significant limitations, especially in modeling emerging electrode materials which involve phase changes and variable diffusivities. A computationally efficient representation for solid-phase diffusion is discussed in this paper based on symmetric polynomials using Orthogonal Collocation and Galerkin formulation (weak form). A systematic approach is provided to increase the accuracy of the approximation (p form in finite element methods) to enable efficient simulation with a minimal number of semi-discretized equations, ensuring mass conservation even for non-linear diffusion problems involving variable diffusivities. These methods are then demonstrated by incorporation into the full p2D model, illustrating their advantages in simulating high C-rates and short-time dynamic operation of Lithium-ion batteries.

25 ENERGY STORAGE↗

DOE OSTI2024

DRIFT: Diffusivity Regulation of Isotopes by nanoconFinementT

Spontaneous isotope fractionation has been reported under nanoconfinement conditions in naturally occurring systems, but the origin of this phenomena is currently unknown. Two existing hypotheses have been proposed, one based on changes in the solvation environment of the isotopes that reduces the non-mass dependent hydrodynamics contribution to diffusion. The other is that isotopes have mass-dependent surface adsorption, varying their total diffusion through nanoconfined channels. To investigate these hypotheses, benchtop experiments, nuclear magnetic resonance (NMR) spectroscopy, and molecule scale modeling were applied. Classical molecular dynamics simulations identified that the Na + and Cl - hydration shells across the three different salt solutions ( 22 Na 35 Cl, 23 Na 35 Cl, 24 Na 35 Cl) did not vary as a function of the Na + isotope, but that there was a significant pore size effect, with larger hydration shells at larger pore sizes. Additionally, while total adsorption times did not vary as a function of the Na + isotope or pore size, the free ion concentration, or those adsorbed on the surface for <5% of the simulation time did exhibit isotope dependence. Experimentally, challenges occurred developing a repeatable experiment, but NMR characterization of water diffusion rates through ordered alumina membranes was able to identify the existence of two distinct water environments associated with water inside and outside the pore. Further NMR studies could be used to confirm variation in hydration shells and diffusion rates of dissolved ions in water. Ultimately, mass-dependence adsorption is a primary driver of variations in isotope diffusion rates, rather than variation in hydration shells that occur under nanoconfinement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Estimating viscoelastic compliance of desiccating cementitious materials using drying prism tests

Highlights: • Determining relative humidity and degree of reaction of desiccating cementitious materials via an integrated method • Extracting moisture diffusion coefficient of cement mortar using the mass loss data of the prism specimens • Estimating uniaxial viscoelastic compliance of cement mortar using the shrinkage strain data of the prism specimens An integrated modeling approach to determining internal relative humidity and degree of reaction of desiccating cementitious materials has been developed. By simultaneously accounting for both self-desiccation and external drying, the model is capable of predicting the evolution of material properties of desiccating cementitious materials. Using the mass loss data of the prism specimens, the moisture diffusion coefficient of a nonlinear model for cement mortar was determined. Additionally, a previous poroviscoelastic shrinkage model was modified and extended. Through comparison of predicted and measured axial strain of the prism specimens, the viscoelastic compliance of cement mortar was extracted. The present work revealed that it is possible to use a simple drying prism test to fit the moisture diffusion coefficient and creep compliance of cementitious materials, though the comparison of the extracted viscoelastic compliance with measured compliance was only reasonable if the outer layer of the drying mortar was presumed to have microcracking.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Temperature dependence of heavy quark diffusion from (2+1)-flavor lattice QCD

We present a lattice determination of the heavy-quark diffusion coefficient in (2+1)-flavor QCD with almost physical quark masses. The momentum and spatial diffusion coefficients are extracted for a wide temperature range, from T = 163 MeV to 10 GeV. The results are in agreement with previous works from the HotQCD collaboration, and show fast thermalization of the heavy quark inside the QGP. Near the chiral crossover temperature T c ≃ 150 MeV, our results are close to the AdS/CFT estimation computed at strong coupling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Conversion Kinetics and Ionic Conductivity in Na-β”-Alumina + YSZ (Naβ”AY) Sodium Solid Electrolyte via Vapor Phase Conversion Process

Sodium ion batteries have been receiving increasing attention and may see potential revival in the near future, particularly in large-scale grid energy storage coupling with wind and solar power generation, due to the abundant sodium resources, low cost, and sufficiently high energy density. Among the known sodium ion conductors, the Na-β”-alumina electrolyte remains highly attractive because of its high ionic conductivity. This study focuses on the vapor phase synthesis of a Na-β”-Alumina + YSZ (Naβ”AY) composite sodium electrolyte, which has higher mechanical strength and stability than conventional single phase β”-Alumina. The objectives are the measurement of conversion kinetics through a newly developed weight-gain based model and the determination of sodium ionic conductivity in the composite electrolyte. Starting samples contained ~70 vol% α-Alumina and ~30 vol% YSZ (3 mol% Y 2 O 3 stabilized Zirconia) with and without a thin alumina surface layer made by sintering in air at 1600 °C. The sintered samples were placed in a powder of Na-β”-alumina and heat-treated at 1250 °C for various periods. Sample dimensions and weight were measured as a function of heat treatment time. The conversion of α-Alumina in the α-Alumina + YSZ composite into Naβ”AY occurred by coupled diffusion of sodium ions through Na-β”-alumina and of oxygen ions through YSZ, effectively diffusing Na2O. From the analysis of the time dependence of sample mass and dimensions, the effective diffusion coefficient of Na 2 O through the sample, D eff , was estimated to be 1.74 x 10 -7 cm2 s -1 , and the effective interface transfer parameter, k eff , was estimated as 2.33 x 10 -6 cm s -1 . By depositing a thin alumina coating layer on top of the bulk composite, the chemical diffusion coefficient of oxygen through single phase Na-β”-alumina was estimated as 4.35 x 10 -10 cm 2 s -1 . An AC impedance measurement was performed on a fully converted Naβ”AY composite, and the conductivity of the composite electrolyte was 1.3 x 10 -1 S cm -1 at 300 °C and 1.6 x 10 -3 S cm -1 at 25 °C, indicating promising applications in solid state or molten salt batteries at low to intermediate temperatures.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Graphene oxide precursor effects on 3D-printed carbon scaffolds

Manganese oxide (MnO 2 ), an earth-abundant material, is a promising component for energy storage devices, with uses in both pseudocapacitors and batteries. However, high MnO 2 loading often leads to reduced performance due to poor ion diffusion. 3D printing, particularly using the direct ink writing (DIW) technique, offers a solution by enabling the fabrication of electrodes with hierarchical porous structures and open channels that enhance mass transport and ion diffusion. Previous work demonstrated that 3D-printed graphene aerogels with MnO 2 coatings exhibited excellent electrochemical performance, even with thick electrodes, due to their optimized structure. Building on this work, the current study investigates the performance differences between aerogels developed using graphene oxide (GO) and reduced graphene oxide (rGO) as carbon precursors. Both materials were incorporated into thixotropic inks, 3D-printed into lattice structures, and carbonized. Despite expected similarities between the final graphene aerogel, rGO-based aerogels exhibited superior areal capacitance, compared to GO-based aerogels. These differences are attributed to the lower oxygen content and defect density of rGO, which influence its interaction with cellulose viscosifiers in the ink formulation. Brunauer–Emmett–Teller (BET) surface area analysis revealed that rGO aerogels exhibit a larger surface area and mesoporous structure, further enhancing their performance. When coated with MnO 2 , rGO-based aerogels maintained their superior capacitive behavior over GO-based aerogels. This study highlights the effect of carbon precursor on the end performance of graphene aerogels.

Materials science↗

DOE OSTI2024

Ionic Diffusion in Slurry Electrolytes for Redox Flow Batteries

Slurry electrodes have been proposed as a method to decouple the storage and power capacities of hybrid redox flow batteries by allowing the reduced metal to adhere to a flowing dispersion of electrically conductive particles rather than accumulate in the flow cell. In this work, the ionic mass transport through these slurry electrolytes is investigated via voltammetry at a rotating disc electrode and in a parallel channel flow cell. A modified Levich equation, i l = 0.6 2 nFC ω 1 / 2 ν − 1 / 6 ( D + r 2 Φ ω m ) 2 / 3 , is developed to describe limiting currents to a rotating disc electrode when using a particle filled electrolyte. Mass transfer is demonstrated to be enhanced by greater particle loading so long as there are sufficient particles such that their respective flow velocity boundary layers interact. In the flow cell, faradaic current in a carbon black slurry electrolyte was observed to decrease over time. This suggests that the carbon particles adhere and accumulate on the stationary electrode, impeding mass transport by imposing diffusion path tortuosity. This accumulation is undesirable in hybrid redox flow batteries as it encourages metal deposition in the cell, limiting the scalability of the technology.

Electrochemistry↗

DOE OSTI2026

Diffusion and Solvation Dynamics of Ions in Water: Beyond the Brownian Approximation

The coupled dynamics of ions and water molecules in their first hydration shell impact a variety of processes including ion diffusion, selective ion transport in water-filled nanopores, and the kinetics of ion-pairing, ion adsorption, and metal-ligand binding reactions. In this work, we study these coupled dynamics for alkali metals (Li, Na, K, Rb, Cs), alkaline Earth metals (Mg, Ca, Sr, Ba), and chloride through the lens of their dependence on ion isotopic mass. Results are validated against previous measurements of the isotopic mass-dependence of ion diffusion coefficients in water and previous ab initio calculations of ion high-frequency dynamics in water. We find that the vibrational power spectra of ions in water consistently exhibit either two or three peaks, i.e., ions have several rattling frequencies within their solvations shells as previously reported for a subset of the species examined here. These frequencies have different sensititivies to isotopic mass that may serve as signatures of ion solvation processes (such as the tendency of ions to orient their first-shell water molecules) and that also may relate to Hofmeister-like effects including the relative affinity of different metals for ribonucleic acid (RNA).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Sulfur-enhanced dynamics of coinage metal(111) surfaces: Step edges versus terraces as locations for metal-sulfur complex formation

The propensity of trace amounts of sulfur adsorbed on coinage metal(111) surfaces to dramatically enhance surface dynamics has been demonstrated by STM observations of accelerated 2D island decay for Cu and Ag. It is generally accepted that this enhancement is due to the formation of adsorbed metal-sulfur complexes, which facilitate surface mass transport of the metal. These complexes were originally proposed to form on terraces following the extraction of metal atoms from step edges and subsequent combination with sulfur on the terraces. However, even when thermodynamically feasible, this mechanism may not be kinetically viable for some complexes due to limited coupling of the complex concentration to the surface diffusion flux of metal atoms. Focusing on the case of Cu, we assess various scenarios where complexes are formed either on terraces or instead directly at step edges, the latter being a new paradigm. In this work, a new pathway is proposed for the formation on terraces. A rich variety of structures incorporating S at step edges exist, which could provide a viable source for complexes, at least from a thermodynamic perspective. However, it is necessary to also assess the activation barrier for complex formation and detachment from step edges. This is facilitated by the nudged-elastic-band analysis of the minimum energy path for this process utilizing machine-learning derived potentials based on density functional theory energetics for the metal-sulfur system.

36 MATERIALS SCIENCE↗

DOE OSTI2022

New Particle Formation and Growth Dynamics for α-Pinene Ozonolysis in a Smog Chamber and Implications for Ambient Environments

Secondary organic aerosol (SOA) plays an important role in the Earth’s energy balance and air quality. Understanding the formation of SOA accompanying the oxidation of volatile organic compounds (VOCs) remains incomplete because of the complex reactions and partitioning mechanisms. In this study, we investigated new particle formation (NPF) of SOA derived from α-pinene ozonolysis in the Harvard Environmental Chamber using the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the condensing species as two products, one of lower volatility (LVOC, 0.01 μg m –3 ) and one of higher volatility (SVOC, 1.0 μg m –3 ). A module for the classical nucleation theory was included. The surface tension (σ), bulk diffusivity (D b ), and mass accommodation coefficient (α) were varied until convergence between the simulated and observed particle number size distributions. α of LVOC was constrained as 0.03–0.1 with the reported σ of α-pinene-ozonolysis-derived SOA and Tolman surface correction (~23.0–27.5 mN m –1 ) and Db from 10 –15.5 to >10 –11 cm 2 s –1 . The observed particle number size distribution of α-pinene ozonolysis was well-simulated to derive the physical parameters, including σ, Db, and α, suggesting the possible application to group the reported complex oxidation products and further application for other atmospheric-related oxidations for quantifying NPF.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Effects of diffusion barriers on reaction wave stability in Co/Al reactive multilayers

Bimetallic, reactive multilayers are uniformly structured materials composed of alternating sputter-deposited layers that may be ignited to produce self-propagating mixing and formation reactions. These nanolaminates are most commonly used as rapid-release heat sources. The specific chemical composition at each metal/metal interface determines the rate of mass transport in a mixing and formation reaction. The inclusion of engineered diffusion barriers at each interface will not only inhibit solid-state mixing but also may impede the self-propagating reactions by introducing instabilities to wavefront morphology. Here, this work examines the effect of adding diffusion barriers on the propagation of reaction waves in Co/Al multilayers. The Co/Al system has been shown to exhibit a reaction propagation instability that is dependent on the bilayer thickness, which allows for the occurrence of unstable modes in otherwise stable designs from the inclusion of diffusion barriers. Based on the known stability criteria in the Co/Al multilayer system, the way in which the inclusion of diffusion barriers changes a multilayer's heat of reaction, thermal conductivity, and material mixing mechanisms can be determined. These factors, in aggregate, lead to changes in the wavefront velocity and stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Experimental and computational design tools for industrial drying processes: Challenges in process‐limit prediction

Abstract The coating and drying of inks and slurries are important steps to manufacture a plethora of products. Drying processes, particularly, comprise energy‐intensive steps that affect product cost and quality. Prior work has highlighted failures of various multicomponent diffusivity models to conserve mass in dryer modeling and challenges in predicting process limits given variability in published values of key thermodynamic parameters. Herein, we develop a computational model and benchtop drying experiments to investigate these concerns for drying polymer‐laden coatings. Model predictions of process limits in a single‐zone drying oven demonstrate that published variability in Flory–Huggins parameter yields large variations in predicted operating temperatures above which blistering occurs. This indicates that caution should be exercised when choosing approaches to obtain or predict the Flory–Huggins parameter, and that both benchtop drying experiments and a set of additional experiments, such as sorption experiments, are needed to fully characterize and optimize a given drying process.

Parrish, Chance↗

DOE OSTI2026

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗

DOE OSTI2023

Improved melt model for power flow

Accelerators that drive z-pinch experiments transport current densities in excess of 1 MA/cm 2 in order to melt or ionize the target and implode it on axis. These high current densities stress the transmission lines upstream from the target, where rapid electrode heating causes plasma formation, melt, and possibly vaporization. These plasmas negatively impact accelerator efficiency by diverting some portion of the current away from the target, referred to as “current loss”. Simulations that are able to reproduce this behavior may be applied to improving the efficiency of existing accelerators and to designing systems operating at ever higher current densities. The relativistic particle-in-cell code CHICAGO ® is the primary code for modeling power flow on Sandia National Laboratories’ Z accelerator. We report here on new algorithms that incorporate vaporization and melt into the standard power-flow simulation framework. Taking a hybrid approach, the CHICAGO® kinetic/multi-fluid treatment has been expanded to include vaporization while the quasi-neutral equation-of-motion has been updated for melt at high current-densities. For vaporization, a new one-dimensional substrate model provides a more accurate calculation of electrode thermal, mass, and magnetic field diffusion as well as a means of emitting absorbed contaminants and vaporized metal ions. A quasi-fluid model has been implemented expressly to mimic the motion of imploding liners for accurate inductance histories. For melt, a multi-ion Hall-MHD option has been implemented and benchmarked against Alegra MHD. This new model is described with sufficient detail to reproduce these algorithms in any hybrid kinetic code. Physics results from the new code are also presented. A CHICAGO ® Hall-MHD simulation of a radial transmission line demonstrates that Hall physics, not included in Alegra, has no significant impact on the diffusion of electrode material. When surface contaminant desorption is mocked in as a hydrogen surface plasma, both the surface and bulk-material plasmas largely compress under the influence of the j × B force. Similar results are seen in Alegra, which also shows magnetic and material diffusion scaling with peak current. Test vaporization simulations using MagLIF and a power-flow experimental geometry show Fe + ions diffuse only a few hundred µm from the electrodes, so present models of Z power flow remain valid.

43 PARTICLE ACCELERATORS↗

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↗