Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse radiation”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Volcano‐like Activity Trends in Au@Pd Catalysts: The Role of Pd Loading and Nanoparticle Size

The addition of palladium (Pd) to preformed gold nanoparticles (Au NPs) enables the formation of core‐shell structures with enhanced catalytic performance in oxidation reactions. However, predicting the precise palladium content required to achieve maximum catalytic activity remains difficult based on current understanding. Herein, Pd was systematically introduced onto titania‐supported Au NPs (2, 6, and 10 nm) to evaluate their performance in benzyl alcohol oxidation. A volcano‐like trend in catalytic activity was observed, where activity increased with Pd addition, peaked, and then declined. The Pd loading required for maximum activity depended on Au NP size: ≈40 at% Pd/Au for 2.6 nm, ≈20 at% Pd/Au for 6.4 nm, and ≈12.5 at% Pd/Au for 10.6 nm. For Au NPs > 6 nm, peak activity aligned with monolayer Pd coverage, while for smaller NPs (2–3 nm), optimal Pd content was below monolayer predictions. X‐ray absorption spectroscopy revealed a core‐shell structure at low Pd content, but higher Pd loadings led to Pd diffusion into the Au core. This structural transformation likely caused activity decline, indicating that AuPd alloying negatively impacts catalysis. These results highlight that core‐shell Au@Pd catalysts outperform AuPd alloys and provide crucial insights for designing highly active bimetallic catalysts.

X-ray absorption spectroscopy↗
DOE OSTI2021

Predicting CaO-(MgO)-Al2O3-SiO2 glass reactivity in alkaline environments from force field molecular dynamics simulations

In this investigation, force field-based molecular dynamics (MD) simulations have been employed to generate detailed structural representations for a range of amorphous quaternary CaO-MgO-Al{sub 2}O{sub 3}-SiO{sub 2} (CMAS) and ternary CaO-Al{sub 2}O{sub 3}-SiO{sub 2} (CAS) glasses. Comparison of the simulation results with select experimental X-ray and neutron total scattering and literature data reveals that the MD-generated structures have captured the key structural features of these CMAS and CAS glasses. Based on the MD-generated structural representations, we have developed two structural descriptors, specifically (i) average metal oxide dissociation energy (AMODE) and (ii) average self-diffusion coefficient (ASDC) of all the atoms at melting. Both structural descriptors are seen to more accurately predict the relative glass reactivity than the commonly used degree of depolymerization parameter, especially for the eight synthetic CAS glasses that span a wide compositional range. Hence these descriptors hold great promise for predicting CMAS and CAS glass reactivity in alkaline environments from compositional information.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Secondary-Source Core Reload Modeling with VERA

The CASL reactor simulation package VERA has been adapted to provide high-fidelity simulation capabilities for modeling source range detector response during subcritical reactor configurations. New features include the activation and shuffling of secondary-source assemblies, use of burned fuel neutron emission data from the ORIGEN depletion solver to the MPACT deterministic neutron transport solver, allowance of user-defined sources in MPACT based on material composition, ability to solve the subcritical source-driven system with neutron multiplication using the MPACT diffusion solver, and transfer of the calculated fission source from MPACT to the continuous-energy Monte Carlo solver Shift for final detector response evaluation using the CADIS methodology for variance reduction. These new capabilities were benchmarked against Watts Bar Unit 1 plant operating data for the first few fuel loading steps and were found to demonstrate excellent agreement with the measured data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Hand Calculations for Nuclear Criticality Safety – Primer Revision

The document “ Hand Calculation Methods for Criticality Safety – A Primer, ” was published by the authors (Los Alamos National Laboratory report LA-14244-M) in December 2006 to provide an overview of the most common hand calculation methods used for criticality safety analysis. Hand calculation methods can be used as a starting point for more advanced calculations, and in many circumstances, they can provide sensitivity and perturbation information quicker than using a criticality code, such as KENO or MCNP. This report has been revised to augment existing basis information and example problems for each of the methods discussed: 1-group and modified 1- group diffusion theory, core-density conversions, buckling conversions, limiting surface density method, density analog method, surface density method, and solid angle method. This primer discusses the applicability of the various methods, illustrates how they are used, and provides an interpretation of the example problem results. The example problems provided are simple, practical problems, and the reference data to solve each example problem is provided. This paper will provide an overview of the Primer revision and how it can be useful for new Nuclear Criticality Safety practitioners.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Atomistic Modeling of Irradiation-Induced Defects and Clusters in Additively-Manufactured Austenitic Stainless Steel

Understanding the progression of irradiation-induced defects in real materials and alloys requires insight into the interaction between point defects, defect clusters, and the existing microstructure. This is particularly true when it comes to additively manufactured (AM) materials with complicated microstructures and residual stresses. In this work, we employ molecular dynamics (MD) to simulate the creation of point defects via primary radiation damage modeling and point defect interactions with extended defects such as voids and dislocations using simplified but representative alloy compositions. We analyze the effect of alloying elements on the formation energies of the defect clusters and measure interactions such as sink strengths and binding energies. This study provides insight into the unique behaviors of irradiation damage in AM materials and information necessary to bridge the gap in length and time scales to connect the atomistic behavior of defects to longer length scale simulations such as cluster dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Frustrated charge density wave and quasi-long-range bond-orientational order in the magnetic kagome FeGe

The intrinsic frustrated nature of a kagome lattice is amenable to the realization of exotic phases of matter, such as quantum spin liquids or spin ices, and the multiple-q charge density waves (CDW) in the kagome metals. Despite intense efforts to understand the mechanism driving the electronic modulations, its origin is still unknown and obscured by competing interactions and intertwined orders. Here, we identify a dimerization-driven 2D hexagonal charge-diffuse precursor in the antiferromagnetic kagome metal FeGe and demonstrate that the fraction of dimerized/undimerized states is the relevant order parameter of the multiple-q CDW of a continuous phase transition. The pretransitional charge fluctuations with propagation vector q = q M at T CDW < T < T*(125 K) are anisotropic, hence holding a quasi-long-range bond-orientational order. The broken translational symmetry emerges from the anisotropic diffuse precursor, akin to the Ising scenario of antiferromagnetic triangular lattices. The temperature and momentum dependence of the critical scattering show parallels to the stacked hexatic B-phases reported in liquid crystals and transient states of CDWs and highlight the key role of the topological defect-mediated melting of the CDW in FeGe.

Subires, D. [Donostia International Physics Center↗
DOE OSTI2021

Electrochromic Properties of Perovskite NdNiO 3 Thin Films for Smart Windows

Semiconductors with electrically tunable band gaps are of great interest in controlling transparency to electromagnetic radiation. Thin films of perovskite nickelate NdNiO 3 (NNO), a class of correlated oxides, were deposited on single-crystal (LaAlO 3 (LAO)) and polycrystalline (fluorine-doped tin oxide-coated glass (FTO)) substrates by magnetron sputtering, chemical solution deposition (CSD), and atomic layer deposition (ALD). Their electrochromic behaviors were investigated using a three-electrode setup in basic (KOH solution, pH = 12) electrolyte. During bleaching/coloration process, the proton intercalation/deintercalation and simultaneous electron compensation in the NNO lattice under electrical bias led to crossover of the material between the pristine-conducting phase (Ni 3+ ) and the strongly correlated insulating phase (Ni 2+ ), which serves as the working principle for electrochromic (tunable opacity in the visible range) behavior. Cyclic voltammetry (CV) scans demonstrate that NNO films are electrochemically stable in basic solutions for all three film deposition methods explored here. CV scans at varying rates enabled the extraction of diffusion coefficient of protons in thin film NNO, which is similar to 10 –7 cm 2 s –1 among all films studied. Large light transmittance modulation by bleaching and coloration was observed on films grown on both LAO and FTO substrates, suggesting its potential as an electrochromic material candidate for smart windows and optical shutter applications. Here, porous NNO films obtained by chemical solution deposition tend to demonstrate stronger electrochromic activity than dense films grown by sputtering or ALD.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Study of interface reaction in a B 4 C/Cr mirror at elevated temperature using soft X-ray reflectivity

Boron carbide is a prominent material for high-brilliance synchrotron optics as it remains stable up to very high temperatures. The present study shows a significant change taking place at 550°C in the buried interface region formed between the Cr adhesive layer and the native oxide layer present on the silicon substrate. An in situ annealing study is carried out at the Indus-1 Reflectivity beamline from room temperature to 550°C (100°C steps). The studied sample is a mirror-like boron carbide thin film of 400 Å thickness deposited with an adhesive layer of 20 Å Cr on a silicon substrate. The corresponding changes in the film structure are recorded using angle-dependent soft X-ray reflectivity measurements carried out in the region of the boron K-edge after each annealing temperature. Analyses performed using the Parratt recursive formalism reveal that the top boron carbide layer remains intact but interface reactions take place in the buried Cr–SiO2 region. After 300°C the Cr layer diffuses towards the substrate. At higher temperatures of 500°C and 550°C the Cr reacts with the native oxide layer and tends to form a low-density compound of chromium oxysilicide (CrSiOx). Depth profiling of Si and Cr distributions obtained from secondary ion mass spectroscopy measurements corroborate the layer model obtained from the soft X-ray reflectivity analyses. Details of the interface reaction taking place near the substrate region of boron carbide/Cr sample are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Distribution of Neutrons in a Cold Lucite Moderator

A series of experiments was performed in 2010 to study the behavior of neutrons in a cold Lucite moderator. The diffusion of neutrons in a cold moderator is described as a two-step process. In the first step, the fast neutrons produced by a neutron source collide with the nuclei of Lucite, losing energy until they reach thermal energies. In the second step, the thermal neutrons continue to diffuse without any significant loss of energy until they are absorbed or leak out of the system. However, because the Lucite moderator is at low temperature, there is less up scattering in the neutron energy compared to a Lucite moderator at room temperature. Experiments were performed to study the behavior of neutrons and their distribution in a cold Lucite moderator and their results compared to simulations. Based on these measurements, an estimation of the absorption cross section at -78.5 °C yielded a value of 0.54 barns ± 0.02 barns.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Substorm Driven Chorus Waves: Decay Timescales and Implications for Pulsating Aurora

Abstract Energetic electron precipitation (EEP) associated with pulsating aurora can transfer greater than 30 keV electrons from the outer radiation belt region into the upper atmosphere and can deplete atmospheric ozone via collisions that produce NOx and HOx molecules. Our knowledge of exactly how EEP occurs is incomplete. Previous studies have shown that pitch angle scattering between electrons and lower‐band chorus waves can cause pulsating aurora associated with EEP and that substorms play an important role. In this work, we quantify the timescale of chorus wave decay following substorms and compare that to previously determined timescales. We find that the chorus decay e‐folding time varies based on magnetic local time (MLT), magnetic latitude, and wave frequency. The shortest timescales occur for lower‐band chorus in the 21 to 9 MLT region and compares, within uncertainty, to the energetic pulsating aurora timescale of Troyer et al. (2022, https://doi.org/10.3389/fspas.2022.1032552 ) for energetic pulsating aurora. We are able to further support this connection by modeling our findings in a quasi‐linear diffusion simulation. These results provide observations of how chorus waves behave after substorms and add additional statistical evidence linking energetic pulsating aurora to substorm driven lower‐band chorus waves.

Troyer, Riley N.↗
DOE OSTI2022

The discovery of a radio galaxy of at least 5 Mpc

Context. Giant radio galaxies (GRGs, or colloquially ‘giants’) are the Universe’s largest structures generated by individual galaxies. They comprise synchrotron-radiating active galactic nucleus ejecta and attain cosmological (megaparsec-scale) lengths. However, the main mechanisms that drive their exceptional growth remain poorly understood. Aims. To deduce the main mechanisms that drive a phenomenon, it is usually instructive to study extreme examples. If there exist host galaxy characteristics that are an important cause for GRG growth, then the hosts of the largest GRGs are likely to possess them. Similarly, if there exist particular large-scale environments that are highly conducive to GRG growth, then the largest GRGs are likely to reside in them. For these reasons, we aim to perform a case study of the largest GRG available. Methods. We reprocessed the LOFAR Two-Metre Sky Survey DR2 by subtracting compact sources and performing multi-scale CLEAN de-convolutions at 60" and 90" resolution. The resulting images constitute the most sensitive survey yet for radio galaxy lobes, whose diffuse nature and steep synchrotron spectra have allowed them to evade previous detection attempts at higher resolution and shorter wavelengths. We visually searched these images for GRGs. Results. We have discovered Alcyoneus, a low-excitation radio galaxy with a projected proper length l p = 4.99 ± 0.04 Mpc. Both its jets and lobes are detected at very high significance, and the SDSS-based identification of the host, at spectroscopic redshift z spec = 0.24674 ± 6 × 10 -5 , is unambiguous. The total luminosity density at ν = 144 MHz is L ν = 8 ± 1 × 10 25 W Hz -1 , which is below average, though near median (percentile 45 ± 3%) for GRGs. The host is an elliptical galaxy with a stellar mass M* = 2.4 ± 0.4 × 10 11 M ⊙ and a super-massive black hole mass M• = 4 ± 2 × 10 8 M ⊙ , both of which tend towards the lower end of their respective GRG distributions (percentiles 25 ± 9% and 23 ± 11%). The host resides in a filament of the Cosmic Web. Through a new Bayesian model for radio galaxy lobes in three dimensions, we estimate the pressures in the megaparsec-cubed-scale northern and southern lobes to be P min,1 = 4.8 ± 0.3 × 10 -16 Pa and P min,2 = 4.9 ± 0.6 × 10 -16 Pa, respectively. The corresponding magnetic field strengths are B min,1 = 46 ± 1 pT and B min,2 = 46 ± 3 pT. Conclusions. We have discovered what is in projection the largest known structure made by a single galaxy – a GRG with a projected proper length l p = 4.99 ± 0.04 Mpc. The true proper length is at least l min = 5.04 ± 0.05 Mpc. Beyond geometry, Alcyoneus and its host are suspiciously ordinary: the total low-frequency luminosity density, stellar mass, and super-massive black hole mass are all lower than, though similar to, those of the medial GRG. Thus, very massive galaxies or central black holes are not necessary to grow large giants, and, if the observed state is representative of the source over its lifetime, neither is high radio power. A low-density environment remains a possible explanation. The source resides in a filament of the Cosmic Web, with which it might have significant thermodynamic interaction. The pressures in the lobes are the lowest hitherto found, and Alcyoneus therefore represents the most promising radio galaxy yet to probe the warm–hot inter-galactic medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Validation of the HyRAM+ physics models for use with propane

The Hydrogen Risk Assessment Model Plus (HyRAM+) toolkit combines quantitative risk assessment with simulations of unignited dispersion, ignited turbulent diffusion flames, and indoor accumulation with delayed ignition of fuels. HyRAM+ is differentiated from HyRAM in that it includes models and leak data for other alternate fuels. The models of the physical phenomena need to be validated for each of the fuels in the toolkit. This report shows the validation for propane which is being used as a surrogate for autogas, which is a mixture of propane and butane and used in internal combustion engines in vehicles. For flame length comparisons, five previously published experiments from peer reviewed journals were used to validate our models. The validation looked at flame lengths and flame widths with respect to different leak diameters, mass flow rates, and source pressures. Most of the sources included more than one set of experimental data, which were collected using different methods (CCD cameras, IR visualization etc.). In general, HyRAM+ overpredicts the flame lengths by around 65%. For heat and radiation models, we compared the heat flux and radiation data reported from two different sources to the values calculated by HyRAM+. For higher mass flow rates, the HyRAM+ calculated flame length results gave a better estimate of what is found in the experiments (65% error), but a higher error (85%) is observed between the HyRAM+ calculated lengths and the experimental flame lengthsfor lower mass flows. Some differences can be attributed to outdoor environmental effects (i.e. wind speed) and uncertainties in jet flame shapes. The propane flame trajectory is predicted for a high Reynolds number case with Re = 12,500 and a low Reynolds number case where Re = 2,000. The Re=12,500 case which is momentum dominated matches well with the experimental flame trajectory, but the agreement for the bouancy driven low Reynolds number case is not as good. Dispersion modeling for unignited propane was also analyzed. We compared the mole fraction, mixture fraction, mean velocity, concentration half width, and inverse mass concentration over an axial distance from different credible journals to the values calculated by HyRAM+. The results display good agreement but generally, HyRAM+ predicts a wider profile for mole fraction and mixture fraction experiments. Overall, HyRAM+’s results are reasonable for predicting the flame length, heat flux, flame trajectory, and dispersion for propane and can be used in risk analyses

08 HYDROGEN↗
DOE OSTI2021

Capture efficiency and bias from the defect dynamics near grain boundaries in BCC Fe using mesoscale simulations

The capture efficiency describes the capability of a sink, such as a grain boundary (GB), dislocation, and void, to absorb point defects (PDs). The bias defines the difference in capture efficiency between the absorption of a vacancy and dumbbell at a sink. Complete kinetic information on PDs, including diffusion barriers and diffusion orientations, as well as accurate saddle points, are needed to determine the capture efficiency and bias at a sink accurately, which is computationally demanding. Here, in the present study, the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method was used to investigate the defect dynamics of PDs near different types of grain boundaries (GBs) (with both $\langle$100$\rangle$ and $\langle$110$\rangle$ families) accurately in body-centered cubic (BCC) iron (Fe). The capture efficiency, sink strength, and bias factor of different types of GBs were determined in Fe, which, different from traditional rate theory estimation, showed a distinct capture efficiency, sink strength, and bias in different GBs. The results demonstrate a strong positive correlation between the capture efficiency and the GB strain width, instead of the GB misorientation, GB energy, or GB-PD binding energy, which have been investigated previously. This work provides valuable insight into the radiation-induced microstructural evolution of GBs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Design of a high energy density experiment to measure the suppression of hydrodynamic instability in an applied magnetic field

A transverse magnetic field can suppress the hydrodynamic instability growth of an unstable plasma interface. This effect, of interest to inertial confinement fusion and astrophysics, has mostly been studied via simulation. Here, we present the design of an experiment at the National Ignition Facility (NIF) to demonstrate this effect in a laboratory. Simulations indicate that the timescale for the diffusion of the magnetic field across the mixing region should be at least comparable to the timescale of the instability growth in order to have a measurable suppression effect. This motivates the use of lower density target materials than usual high energy density (HED) hydrodynamics experiments to permit faster hydrodynamics and higher plasma conductivities (through higher temperature), for a given laser drive and magnetic field. We discuss a target design for creating a Rayleigh–Taylor unstable HED plasma interface that uses 320 mg/cc iodine-doped carbon foam as the heavy material, 20 mg/cc carbon foam as the light material, and a 6 lm amplitude, a 120 lm wavelength ripple machined at the interface, which shows a measurable suppression effect with a nominal NIF drive and a 30 T magnetic field (the present facility limit). Here, models indicate lower density foams that may display even larger suppression effects, as the Hohlraum drive also radiatively preheats the foam to permit even higher temperatures and, hence, higher conductivities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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

Mechanistic Insights into the Surface Instabilities of TiNb 2 O 7, a High‐Power Li‐Ion Anode

TiNb 2 O 7 (TNO) is a promising Li-ion battery anode for high-power applications, such as implantable medical devices and heavy-duty equipment. Hailed as being safe due to its elevated operating potential near 1.6 V, TNO has long been assumed to be highly stable in the carbonate-based electrolytes used in Li-ion batteries. Herein, all mechanisms occurring at the surface of both TNO and Nd-doped TNO are identified, and both materials in fact show significant gassing. CO 2 is even released at open circuit conditions, demonstrating the poor chemical stability of the material in the electrolyte even prior to battery operation. Such extreme instability is a critical safety concern. In addition, it was found that Ti dissolves from the surface of TNO particles at low voltage (below 1.4 V vs Li), and in fact deposits on the counter electrode. Ti further inside TNO particles then diffuses to the Ti-poor surface during discharge. Partial carbon-coating as a mitigating measure has also been tested and found to exacerbate these processes. The findings identify novel reactions occurring within TNO, and clearly highlight the need to stabilize the surfaces of TNO in order to prevent such aggressive deterioration at the surface of the particles.

25 ENERGY STORAGE↗