Engineering Papers⌕ Search

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.

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Memoirs of Mass Accretion: Probing the Edges of Intracluster Light in Simulated Galaxy Clusters

The diffuse starlight extending throughout massive galaxy clusters, known as intracluster light (ICL), has the potential to be read as a memoir of mass accretion: informative, individual, and yet imperfect. Here, we combine dark-matter-only zoom-in simulations from the Symphony suite with the Nimbus “star-tagging” model of the stellar halo to assess how much information about the mass assembly of an individual galaxy cluster can be gleaned from idealized measurements of ICL outskirts. We show that the edges of a cluster’s stellar profile—the primary (R sp⋆,1 ) and secondary (R sp⋆,2 ) stellar “splashback” radii—are sensitive to both continuous mass accretion histories (MAHs) and discrete merger events, making them potentially powerful probes of a cluster’s past. We find that R sp⋆,1 strongly correlates with the cluster’s mass ∼1 dynamical time ago, while R sp⋆,2 traces more recent MAH to a slightly lesser degree. In combination, these features can further distinguish between clusters that have and have not undergone a major merger within the past dynamical time. We use both to predict realistic cluster MAHs with the MultiCAM framework. These outer ICL features are significantly more sensitive to mass accretion and merger histories than the stellar mass gap and halo concentration, and perform comparably to the commonly used X-ray-based tracer of relaxedness, x off . While our analysis is idealized, the relevant ICL features are potentially detectable in next-generation deep imaging of nearby clusters. This work highlights the promise of ICL measurements and lays the groundwork for more detailed forecasts of their power.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

High-energy neutrino constraints on primordial black holes as dark matter

Primordial black holes (PBHs) are one of the most appealing dark matter candidates over a wide range of masses and abundances. This broad parameter space has been constrained by a variety of observational probes. In this work, for the first time, we use data from high-energy neutrino telescopes, like IceCube and ANTARES, to constrain sub-asteroid mass ($\lesssim 10^{18}\,\mathrm{g}$) Schwarzschild PBHs with extended mass functions. We derive limits from the diffuse high-energy neutrino flux produced by the direct evaporation of PBHs, as well as from the transient signatures associated with PBHs passing in the vicinity of the Earth. While our bounds are slightly weaker than existing constraints from gamma-ray observations, they provide an independent and complementary probe based on observational high-energy neutrino data. We further show that future detectors such as IceCube-Gen2 and KM3NeT can significantly improve these constraints, potentially excluding PBHs with masses up to $\sim \mathrm{few} \times 10^{18}\,\mathrm{g}$ composing the entirety of dark matter.

Mukhopadhyay, Mainak [Fermilab; Chicago U., KICP] ↗
DOE OSTI2023

A Novel Experimental Approach to Understand the Transport of Nanodrugs

Nanoparticle-based drugs offer attractive advantages like targeted delivery to the diseased site and size and shape-controlled properties. Therefore, understanding the particulate flow of the nanodrugs is important for effective delivery, accurate prediction of required dosage, and developing efficient drug delivery platforms for nanodrugs. In this study, the transport of nanodrugs including flow velocity and deposition is investigated using three model metal oxide nanodrugs of different sizes including iron oxide, zinc oxide, and combined Cu-Zn-Fe oxide synthesized via a modified polyol approach. The hydrodynamic size, size, morphology, chemical composition, crystal phase, and surface functional groups of the water-soluble nanodrugs were characterized via dynamic light scattering, transmission electron microscopy, scanning electron microscopy-energy dispersive X-ray, X-ray diffraction, and fourier transform infrared spectroscopy, respectively. Two different biomimetic flow channels with customized surfaces are developed via 3D printing to experimentally monitor the velocity and deposition of the different nanodrugs. A diffusion dominated mechanism of flow is seen in size ranges 92 nm to 110 nm of the nanodrugs, from the experimental velocity and mass loss profiles. The flow velocity analysis also shows that the transport of nanodrugs is controlled by sedimentation processes in the larger size ranges of 110–302 nm. However, the combined overview from experimental mass loss and velocity trends indicates presence of both diffusive and sedimentation forces in the 110–302 nm size ranges. It is also discovered that the nanodrugs with higher positive surface charges are transported faster through the two test channels, which also leads to lower deposition of these nanodrugs on the walls of the flow channels. The results from this study will be valuable in realizing reliable and cost-effective in vitro experimental approaches that can support in vivo methods to predict the flow of new nanodrugs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of water droplet growth dynamics on electrode current in fuel-cell catalyst layers

Fuel cells are a promising next-generation energy-conversion technology designed to replace internal combustion engines in transportation applications. However, much work remains to optimize them. Operation at high humidities causes liquid water droplet formation on Pt catalyst particles during oxygen reduction, potentially impeding reactant arrival to the reactive electrode. In this work, four different cases of water droplet growth in fuel-cell catalyst layers are considered: pinned or advancing droplets on a bare Pt surface, advancing droplets on a Nafion film, and water-layer growth in carbon nanopores. Transient drop growth is captured with a combination of mass, species mass, and momentum balances, and the subsequent limiting current is determined via oxygen diffusion and Tafel kinetics. Further, water droplets are found not to be mass-transfer limiting due to the relatively large liquid-gas area compared to the Pt nanoparticle. Mass-transfer-limited behavior is calculated in carbon nanopores.

25 ENERGY STORAGE↗
DOE OSTI2021

Abell 1430: A merging cluster with exceptional diffuse radio emission

Diffuse radio emission has been found in many galaxy clusters, predominantly in massive systems which are in the state of merging. The radio emission can usually be classified as relic or halo emission, which are believed to be related to merger shocks or volume-filling turbulence, respectively. Recent observations have revealed radio bridges for some pairs of very close galaxy clusters. The mechanisms that may allow one to explain the high specific density of relativistic electrons, which are necessary to explain the radio luminosity of these bridge regions, have been poorly explored until now. When inspecting the first data release of the LOFAR Two-Metre Sky Survey (LoTSS), we discovered diffuse radio emission in the galaxy cluster Abell 1430. Here, for this work, we aim to determine the dynamical state of the cluster and characterise the diffuse radio emission. We analysed the LoTSS data in detail and complemented them with recent Karl G. Jansky Very Large Array observations in the L-band. To study the dynamical state of the cluster, we analysed XMM-Newton data, Chandra data, and Sloan Digital Sky Survey data. Moreover, we compared our results to clusters extracted from THE THREE HUNDRED PROJECT cosmological simulation. We find that Abell 1430 consists of two components, namely A1430-A and A1430-B, with a mass ratio of about 2:1. The massive component shows diffuse radio emission which can be classified as radio halo which shows a low radio power at 1.4 GHz with respect to the mass of the cluster. Most interestingly, there is extended diffuse radio emission in the following dubbed as the ‘Pillow’ according to its morphology, which is apparently related to A1430-B and which is neither typical halo nor typical relic emission. The origin of this emission is puzzling. We speculate that the two components of Abell 1430 undergo an off-axis merger. In this scenario, A1430-B is moving towards the main cluster component and may have compressed and stirred the medium in the filament between the two cluster components. We have discovered evidence for diffuse radio emission related to the low-density intracluster or intergalactic medium in Abell 1430. To date, only a few examples of emission originating from such regions are known. These discoveries are crucial to constrain possible acceleration mechanisms which may allow us to explain the presence of relativistic electrons in these regions. In particular, our results indicate a spectral index of α 144 MHz 1.5 GHz = -1.4±0.5 for the Pillow. If upcoming observations confirm a slope as flat as -1.4 or even flatter, this would pose a challenge for the electron acceleration scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Lattice B -field correlators for heavy quarks

We analyze the color-magnetic (or “ B ”) field two-point function that encodes the finite-mass correction to the heavy-quark momentum-diffusion coefficient. The simulations are done on fine isotropic lattices in the quenched approximation at 1.5 T c , using a range of gradient flow times for noise suppression and operator renormalization. The continuum extrapolation is performed at fixed flow time followed by a second extrapolation to zero flow time. Perturbative calculations to next-to-leading order of this correlation function, matching gradient-flowed correlators to MS ¯ , are used to resolve nontrivial renormalization issues. We perform a spectral reconstruction based on perturbative model fits to estimate the coefficient κ B of the finite-mass correction to the heavy-quark momentum-diffusion coefficient. The approach we present here yields high-precision data for the correlator with all renormalization issues incorporated at next-to-leading order and is also applicable for actions with dynamical fermions. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

$S^5$: Tidal Disruption in Crater 2 and Formation of Diffuse Dwarf Galaxies in the Local Group

We present results of a spectroscopic campaign around the diffuse dwarf galaxy Crater 2 (Cra2) and its tidal tails as part of the Southern Stellar Stream Spectroscopic Survey ($S^5$). Cra2 is a Milky Way dwarf spheroidal satellite with extremely cold kinematics, but a huge size similar to the Small Magellanic Cloud, which may be difficult to explain within collisionless cold dark matter. We identify 143 Cra2 members, of which 114 belong to the galaxy's main body and 29 are deemed part of its stellar stream. We confirm that Cra2 is dynamically cold (central velocity dispersion $2.51^{+0.33}_{-0.30}\,{\rm km\,s^{-1}}$) and also discover a $\approx$7$σ$ velocity gradient consistent with its tidal debris track. We separately estimate the stream velocity dispersion to be $5.74^{+0.98}_{-0.83}\,{\rm km\,s^{-1}}$. We develop a suite of $N$-body simulations with both cuspy and cored density profiles on a realistic Cra2 orbit to compare with $S^5$ observations. We find that the velocity dispersion ratio between Cra2 stream and galaxy ($2.30^{+0.41}_{-0.35}$) is difficult to reconcile with a cuspy halo with fiducial concentration and an initial mass predicted by standard stellar mass$-$halo mass relationships. Instead, either a cored halo with relatively small core radius or a low-concentration cuspy model can reproduce this ratio. Despite tidal mass loss, Cra2 is metal-poor ($\langle \rm[Fe/H]\rangle=-2.16\pm0.04$) compared to the stellar mass$-$metallicity relation for its luminosity. Other diffuse dwarf galaxies similar to Cra2 in the Local Group (Antlia 2 and Andromeda 19) also challenge galaxy formation models. Finally, we discuss possible formation scenarios for Cra2, including ram-pressure stripping of a gas-rich progenitor combined with tides.

Limberg, Guilherme [Chicago U., KICP; Chicago U.] ↗
DOE OSTI2024

Compression Acceleration of Protons and Heavier Ions at the Heliospheric Current Sheet

Abstract Recent observations by the Parker Solar Probe (PSP) suggest that protons and heavier ions are accelerated to high energies by magnetic reconnection at the heliospheric current sheet (HCS). By solving the energetic particle transport equation in large-scale MHD simulations, we study the compression acceleration of protons and heavier ions in the reconnecting HCS. We find that the acceleration of multispecies ions results in nonthermal power-law distributions with a spectral index consistent with the PSP observations. Our study shows that the high-energy cutoff of protons can reach E max ∼ 0.1 –1 MeV depending on the particle diffusion coefficients. We also study how the high-energy cutoff of different ion species scales with the charge-to-mass ratio E max ∝ ( Q / M ) α . When determining the diffusion coefficients from the quasi-linear theory with a Kolmogorov magnetic power spectrum, we find that α ∼ 0.4, which is somewhat smaller than α ∼ 0.7 observed by PSP.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Understanding Mass Transport in Copper Electrolyte-Based Dye-Sensitized Solar Cells

Copper redox shuttles, particularly [Cu(tmby) 2 ] 2+/1+ (tmby = bis(4,4',6,6'-tetramethyl-2,2'-bipyridine)), proved to be among the best electrolytes for dye-sensitized solar cells (DSCs), realizing higher power conversion efficiencies both under full sun and indoor illumination than conventional iodide/triiodide and cobalt electrolytes. Even though [Cu(tmby) 2 ] 2+/1+ renders a relatively higher performance, this metal complex is bulky and is limited by mass transport. Since the regeneration of the dye ground state by Cu I and the reaction of Cu II at the counter electrode are comparatively faster processes, the efficiency of DSC involving Cu I /Cu II electrolytes under relatively high light intensities is largely governed by the diffusion of Cu I /Cu II species. Understanding mass transport in these solar cells will enable further improvements in the performance of such copper-based DSCs. Here in the present study, the role of illumination intensity on the photogenerated current and its relationship to mass transport is evaluated using the best cosensitized dye (D35:XY1) and copper electrolyte ([Cu(tmby) 2 ] 2+/1+ ) combination.

14 SOLAR ENERGY↗
DOE OSTI2023

Description of the multinucleon transfer mechanism for Ca 48 + Pu 244 and Kr 86 + Pt 198 reactions in a quantal transport approach

Multinucleon transfer (MNT) reactions involving heavy projectile and target combinations stand as a promising method for synthesizing new neutron-rich exotic nuclei, which may not be possible using hot or cold fusion reactions or fragmentation. Exploring the mechanisms behind MNT reactions is essential and it requires a comprehensive theoretical framework that can explain the physical observables in these reactions. This work aims to show that the quantal diffusion approach based on the stochastic mean-field (SMF) theory is capable of explaining the reaction dynamics observed in MNT reactions. Primary product mass distributions in 48 Ca + 244 Pu reaction at E c.m. = 203.2 MeV and 86 Kr + 198 Pt reaction at E c.m. = 324.2 MeV are calculated and compared with the available experimental data. In this work, we utilize the time-dependent Hartree-Fock (TDHF) calculations to analyze the mean-field reaction dynamics computationally in the reactions 48 Ca + 244 Pu and 86 Kr + 198 Pt for a broad range of initial angular momenta. Quantal transport description based on the SMF approach is used to calculate quantal diffusion coefficients and mass variances in 48 Ca + 244 Pu and 86 Kr + 198 Pt systems. The primary products arising from quasifission reactions are described by joint probability distribution in the SMF approach and those arising from fusion-fission are estimated by using the statistical deexcitation code gemini + +. Mean values of charge and mass numbers, scattering angles of the primary reaction products, and the total kinetic energies after the collision are calculated within the TDHF framework for a broad range of initial angular momenta. Throughout all the collisions, drift toward the mass symmetry and large mass dispersion associated with this drift are observed. Here, the calculated primary fragment and mass distributions using the SMF approach successfully explain experimental observations for the 48 Ca + 244 Pu and 86 Kr + 198 Pt systems. The primary mass distributions, mean values of binary products, and mass dispersions are determined and results are compared with the available experimental data. The observed agreement between the experimental data and SMF results highlights the effectiveness of the quantal diffusion mechanism based on the SMF approach, which does not include any adjustable parameters other than standard parameters of Skyrme energy density functional.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Thermodynamics and kinetics of Pb intercalation under graphene on SiC(0001)

SiC-supported graphene intercalated by a two-dimensional Pb monolayer can provide an appealing platform for spintronic applications. Such a monolayer structure is thermodynamically ultrastable, as observed in recent experiments. However, important fundamentals such as the structure of intercalated phases, locations of intercalated atoms, thermodynamic preference for intercalation, and intercalation pathways for this system have not yet been understood conclusively. Here, in this work, extensive density functional theory calculations are performed to assess Pb intercalation thermodynamics and kinetics under graphene on SiC(0001). We find that intercalation of isolated Pb atoms is strongly disfavored over adsorption on top of graphene. However, intercalation of complete Pb layers in the gallery between SiC and graphene buffer layer is strongly favored over supported Pb monolayers and moderately favored over formation of supported large three-dimensional Pb islands. We also find that initiation of intercalation either by individual Pb atoms directly penetrating graphene or by hopping under a static graphene step edge is energetically prohibitive at experimental temperatures. Consequently, more complex intercalation pathways are operative and further analyzed. We demonstrated that once an intercalated Pb monolayer forms around a graphene step edge, a facile Pb mass transport by Pb vacancy-mediated diffusion can be triggered for continued growth of the intercalated monolayer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Machine learning-based ethylene and carbon monoxide estimation, real-time optimization, and multivariable feedback control of an experimental electrochemical reactor

Electrochemical reduction of CO 2 gas is a novel CO 2 utilization technique that has the potential to mitigate the global climate crisis caused by anthropogenic CO 2 emissions, and enable the large-scale storage of energy generated from renewable sources in the form of carbon-based chemicals and fuels. However, due to the complexity of the electrochemical reactions, the explicit first-principles models for CO2 reduction are not available yet, and there has been a limited effort to develop process modeling, optimization and control of CO 2 electrochemical reactors. To this end, a rotating cylinder electrode (RCE) reactor has been constructed at UCLA to understand the mass transfer and reaction kinetics effects separately on the productivity. In the RCE reactor, the applied potential strongly influences the reaction energetics and the electrode rotation speed affects the hydrodynamic boundary layer and modifies the film mass transfer coefficient, which involves convective and diffusive transport. Further, the present work aims to develop a multi-input multi-output (MIMO) control scheme for the RCE reactor that integrates techniques from artificial and recurrent neural network modeling, nonlinear optimization, and process controller design. Specifically, production rates of two products from the experimental reactor, ethylene and carbon monoxide, are controlled by manipulating two inputs, applied potential and catalyst rotation speed. Process dynamics and controllability are analyzed, a feedback control strategy is designed and the controllers are tuned accordingly. The experimental electrochemical cell is employed to gather data for process modeling and implement the multivariable control system. Finally, the experimental results are presented which demonstrate excellent closed-loop performance by the control system and regulation of the outputs at three different set-points including an economically-optimal set-point.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Superior CO 2 /N 2 separation performance of highly branched Poly(1,3 dioxolane) plasticized by polyethylene glycol

With its high content of CO 2 -philic ether oxygen groups, poly(1,3-dioxolane) (PDXLA) has emerged as an attractive platform to achieve excellent CO 2 /N 2 separation properties for post-combustion carbon capture. Herein we demonstrate that the separation properties of PDXLA can be further enhanced by plasticizing with miscible polyethylene glycol (PEG)-based additives using an integrated experimentation and modeling approach. The effects of the chain end groups and loading level of the additives on the physical properties of the blends are thoroughly investigated, including glass transition temperature (T g ), fractional free volume, and gas transport properties, and the effects can be satisfactorily described using models available for homogeneous blends. Notably, a T g -integrated free volume model is adapted to successfully interpret the unified effect of the blend composition and temperature on gas diffusivity and permeability. A sample containing 45 mass% PEG dimethyl ether (PEGDME with a molecular mass of 240 g/mol) displays stable mixed-gas CO 2 permeability of 1540 Barrer and CO 2 /N 2 selectivity of 40 when challenged with a model flue gas at 60 °C, outperforming Robeson's 2008 upper bound. Elucidating how small plasticizers impact gas transport in homogeneous blends may unravel a facile way to design high-performance membranes for gas separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Bridging the Gap in Carbon Free Iron Making: How Hydrogen Affects the Reduction of Iron Ore between 900 and 1590 °C

Hydrogen-based reduction of iron ore for iron and steel production has emerged as a promising alternative to coal and natural gas. Unlike other hydrogen-based iron ore reduction studies, this research focuses on a wide temperature range across 900–1590 °C, encompassing reduction in solid, mixed, and liquid (slag) phases. For a 20 min exposure to hydrogen, the reduction degree increased monotonically from ∼35% at 900 °C to >90% at 1550 °C, except between 1100 °C and 1400 °C, where it stagnated ∼60%. This experimental work challenges the widely accepted notion that higher temperatures enhance the reduction process. Instead, it reveals an overlooked kinetic bottleneck, suggesting complex thermodynamic and mass transfer limitations influenced by phase transformations, diffusion barriers, and microstructural changes. Density functional theory-based molecular dynamics simulations indicate that oxygen diffusivity in BCC iron is 3.88 × 10 –5 cm 2 /s which is ∼5–10 times higher than that in FCC iron. This study reports that this stagnant reduction degree in the mixed solid–liquid phase is due to competition of multiple mechanisms, such as surface- and bulk-diffusion, pore collapse mechanisms, and crystallographic transitions.

08 HYDROGEN↗
DOE OSTI2021

Does inhomogeneous big bang nucleosynthesis produce an inhomogeneous element distribution today?

Inhomogeneous big bang nucleosynthesis (BBN) produces a spatially inhomogeneous distribution of element abundances at $T \sim 10^9$ K , but subsequent element diffusion will tend to erase these inhomogeneities. In this work, we calculate the cosmological comoving diffusion length for the BBN elements. This diffusion length is limited by atomic scattering and is therefore dominated by diffusion when the atoms are neutral, between the redshifts of recombination and reionization. We find that the comoving diffusion length today is $d_{\text{com}} ≈ 70$ pc for all of the elements of interest except 7 Li , for which $d_{\text{com}}$ is an order of magnitude smaller because 7 Li remains ionized throughout the relevant epoch. This comoving diffusion length corresponds to a substellar baryonic mass scale and is roughly equal to the horizon scale at BBN. These results lend support to the possibility that inhomogeneities on scales larger than the horizon at BBN could lead to a spatially inhomogeneous distribution of elements today, while purely subhorizon fluctuations at BBN can result only in a homogeneous element distribution at present.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Recovery of Alkaline-Earths into Liquid Bi in Ternary LiCl-KCl-SrCl 2 /BaCl 2 Electrolytes at 500 °C

Electrochemical reduction of Sr 2+ and Ba 2+ into liquid Bi was investigated in dilute concentrations of SrCl 2 /BaCl 2 (0-5 mol%) in LiCl-KCl electrolytes at 500 °C to ascertain the limit of liquid Bi electrodes for alkaline-earth recovery. Analysis of the electrodes after constant current electrolysis to the specific charge of 270 C g -1 showed Sr 2+ ions consuming 29% of charge at 5 mol% before dropping to 8%-10% of the total charge at 0.45-0.72 mol% SrCl 2 . Ba 2+ ions consumed 54% at 5 mol% BaCl 2 before decreasing to 22%-24% at 0.42-0.89 mol% BaCl 2 ; substantial co-deposition of Li was observed in all chemistries, consuming up to 53% of charge. Considering only 1% of the total charge was consumed for depositing Ba 2+ and Sr 2+ ions in ~0.1 mol% SrCl 2 /BaCl 2 electrolyte, the lower recovery limit of Bi for alkaline-earth elements is suggested to be at ~0.4 mol% SrCl 2 /BaCl 2 to achieve appreciable deposition of alkaline-earths (>1.5 mol% Ba/Sr in liquid Bi). The overpotentials of liquid Bi at 5 mol% of SrCl 2 /BaCl 2 were evaluated by electrochemical impedance spectroscopy. The co-deposition of Sr and Li exhibited the largest increase in charge transfer resistances implying sluggish charge transfer kinetics whereas the co-deposition of Ba and Li exhibited a large increase in mass transport resistances due to the slow diffusion of Ba 2+ ions in the electrolyte.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Correlation Between Microscopic Current Fluctuations Observed at Ultra-Microelectrodes and Macroscopic Bulk Electrolysis Performance in Redox-Active Microemulsions

Microemulsions (μEs) have been proposed as redox flow battery (RFB) electrolytes that maximize ionic conductivity and charge capacity by synergizing two immiscible phases. However, charge transfer during electrolysis in μEs is poorly understood. Here, we show that ultramicroelectrode electrolysis of ferrocene-loaded μEs –20%, 60%, and 90% water - reveals stochastic current fluctuations. These are differentiated in the scanning electrochemical microscopy (SECM) geometry, where power spectral density analysis showed distinct changes in the frequency contributions. SECM in the substrate generation-tip collection mode showed that fluctuations arise under mass-transfer control. Significant differences in the diffusion coefficient of ferrocene species were deducted from SECM approach curves, suggesting phase transfer behavior. Using bulk electrolysis, we calculated the charge accessibility and cycling behavior in the μEs. A decrease in the stochastic behavior of the μEs seems to correlate to a higher accessibility and cycling performance, with the 90% water μE displaying the best reversibility and the 60% the lowest. Altogether, these results suggest that Marangoni-type convection driven by concentration gradients and/or μE restructuring during charge transfer play a role in the electrochemical performance of μEs. This presents opportunities for screening and diagnosing the performance of these emerging RFB electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2025

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗