Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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

Material properties of aqueous potassium acetate solutions for use in liquid desiccant air conditioning

We present thermophysical properties of aqueous potassium acetate solutions relevant to their application as a liquid desiccant in air-conditioning systems. Liquid desiccant air conditioners could provide up to 80% energy savings compared to high-efficiency vapor compression air conditioners. However, the commonly used liquid desiccants, particularly chloride salt solutions, are highly corrosive. This precludes the use of metallic components, necessitating specialized plastics and thereby driving up cost, weight, and limiting operational temperature and pressure ranges. Potassium acetate has been identified as a less corrosive alternative, potentially enabling broader material compatibility in LDAC components. In this study we systematically measured and modelled key thermophysical properties of KAc solutions density, vapor pressure, viscosity, specific heat capacity, phase behavior, and diffusion coefficient across a range of concentrations (10 - 70wt%) and temperatures. For instance, KAc exhibits a vapor pressure of 3.1 kPa at 50wt% and 40degrees C, demonstrating dehumidification capability comparable to chloride salts. By providing an extended dataset and associated models, this work advances the understanding needed to design energy-efficient, corrosion mitigating LDAC system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Statistical Learning for Nonlinear Model Reduction from Local Simulations of Stochastic and Particle- and Agent-Based Systems

Stochastic physical systems across the sciences that have very high-dimensional state spaces, with a large number of fast degrees of freedom that force direct simulators to proceed by integration steps that are orders of magnitude smaller than events of interests (e.g., particle collisions). Examples range from molecular motion to dynamics of large populations of cells. A grand challenge in the simulation and understanding of such systems is the systematic construction of accurate, interpretable, reduced models, enabling faster simulations, revealing fundamental properties of the dynamics, and predicting phenomena of interest that the original simulator could not reached with sufficient accuracy or within a given computational budget. In this projected we developed novel statistical estimation/machine learning techniques for analyzing and building empirical reduced models for important families of high-dimensional stochastic systems, in particular: - we developed techniques for estimating interaction kernels in interacting particle- and agent-based systems, which are ubiquitous in Physics, Biology and many other sciences, given observed trajectories of the system; - we developed techniques for nonlinear model reduction for high-dimensional stochastic systems that have a small number of unknown, nonlinear slow variables, and a large number of fast modes, that are possibly of large magnitude, given observed short trajectories of the system in the form of bursts of trajectories from different initial conditions; - we developed novel techniques for estimating linear dynamical systems on graphs when both the dynamics and the underlying graph are unknown, and we have a sparse set of space-time observations; - we considered the problem of estimating an unknown nonlinear observation function of a standard process (e.g. Brownian motion), so that we can recognized if an observed dynamics is "just" a nonlinear version of a known dynamics; we also developed benchmarks for learning algorithms aimed at learning and classifying diffusion processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Modeling and Validation of Sodium Plugging for Heat Exchangers in Sodium-cooled Fast Reactor Systems

The project “Modeling and Validation of Sodium Plugging for Heat Exchangers in Sodium-cooled Fast Reactor Systems” was conducted jointly by Westinghouse Electric Company (Westinghouse) and Argonne National Laboratory (ANL), over the period October 1, 2013- March 31, 2016. The project’s motivation was the need to provide designers of Sodium Fast Reactors (SFRs) with a validated, state-of-the-art computational tool for the prediction of sodium oxide (Na 2 O) deposition in small-diameter sodium heat exchanger (HX) channels, such as those in the diffusion bonded HXs proposed for SFRs coupled with a supercritical CO 2 (sCO 2 ) Brayton cycle power conversion system. In SFRs, Na 2 O deposition can potentially occur following accidental air ingress in the intermediate heat transport system (IHTS) sodium and simultaneous failure of the IHTS sodium cold trap. In this scenario, oxygen can travel through the IHTS loop and reach the coldest regions, represented by the cold end of the sodium channels of the HXs, where Na 2 O precipitation may initiate and continue. In addition to deteriorating HX heat transfer and pressure drop performance, Na 2 O deposition can lead to channel plugging especially when the size of the sodium channels is small, which is the case for diffusion bonded HXs whose sodium channel hydraulic diameter is generally below 5 mm. Sodium oxide melts at a high temperature well above the sodium melting temperature such that removal of a solid plug such as through dissolution by pure sodium could take a lengthy time. The Sodium Plugging Phenomena Loop (SPPL) was developed at ANL, prior to this project, for investigating Na 2 O deposition phenomena within sodium channels that are prototypical of the diffusion bonded HX channels envisioned for SFR-sCO 2 systems. In this project, a Computational Fluid Dynamic (CFD) model capable of simulating the thermal-hydraulics of the SPPL test section and provided with Na 2 O deposition prediction capabilities, was developed. This state-of-the-art computational tool incorporates a first-principles Na 2 O deposition model developed by ANL, and combines it with predictive capabilities for the spatial and temporal variation of temperature, velocity, dissolved oxygen concentration, and wall temperature under flowing sodium conditions. The CFD model was validated under no-deposition conditions using experimental data collected with the SPPL, demonstrating the model’s capability to predict the thermal-hydraulics of the SPPL test section within the measurement uncertainty characterizing the SPPL instrumentation. The model’s deposition prediction capability was not, however, validated as the SPPL could not be operated under plugging conditions during the project, resulting in the lack of deposition data with adequate pedigree for a CFD model validation. Two novel diagnostic techniques to detect and characterize Na 2 O deposits, i.e. Ultrasonic Time Domain Reflectometry (UTDR) and Potential Drop (PD) techniques, were developed to ultimately assist in the validation effort under plugging conditions, which can be performed once the SPPL becomes operational. This development effort consisted first in demonstrating, analytically and/or computationally, the capability of these techniques to diagnose Na 2 O deposits inside of small channels (particularly the deposit’s thickness), and subsequently in the fabrication and testing of prototypical UTDR and PD instrumentation. The testing, performed on mockups of the SPPL test section, demonstrated the capability of these techniques to detect and characterize material discontinuities like those induced by sodium oxide deposition on stainless steel channel walls. Because of the mentioned impossibility to run the SPPL in a plugging mode, the developed instrumentation could not be tested in-situ, i.e. at the SPPL while deposits are being formed inside of the SPPL test section. Recommended future work includes a possible enhancement in the CFD modeling technique and installation of the developed UTDR and PD instrumentation on the test section, followed by plugging tests to be conducted with the SPPL. The installation of the UTDR and PD diagnostic instrumentation on the SPPL test section will allow collection of Na 2 O deposition data after the onset of deposition to nearly complete channel plugging, which can ultimately be used for the validation of the CFD model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗
DOE OSTI2024

Understanding particle size effect on fast-charging behavior of graphite anode using ultra-thin-layer electrodes

Extreme fast charging (≤ 15 min) of lithium-ion batteries is highly desirable to accelerate mass-market adoption of electric vehicles. However, significant capacity fading, as well as safety issues due to the lithium plating caused by the fast charging rate, limit its implementation. In this study, we investigated the fast-charging capability of graphite materials with various particle sizes. To eliminate the Li + ion concentration gradient effect across the thickness of the electrode, ultra-thin-layer graphite electrodes were developed to investigate the "real" fast-charging capability of graphite at the particle level. Electrochemical assessments as well as microscopic characterizations revealed that smaller particles exhibited superior fast-charging performance, featuring enhanced capacity reversibility, faster charging rate, and less lithium plating under the same fast-charging conditions. It is shown that small-particle graphite (mean radius of 3.3 μm) could withstand a 4C charge (to 80 % state-of-charge) without plating, with minimal plating occurring at 6C. Thicker particles exhibited plating at lower C-rates. Since the experimental data could not directly explain whether intra-particle diffusion limitations or interfacial reaction limitations dominated the plating mechanism, the pseudo-2-dimensional model was used to evaluate the most likely plating mechanism. The model suggested that particle-level diffusion is the dominant mechanism contributing to plating at high rates. Finally, this work provides comprehensive insights into the particle size effects on fast-charging capability, offering a better understanding of fast-charging behavior and valuable guidance for designing optimal electrode architecture for high-rate lithium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2019

Plasma kinetic effects on interfacial mix and burn rates in multispatial dimensions

The physics of mixing in plasmas is of great significance to inertial confinement fusion and high energy density laboratory experiments. Two- and three-dimensional (2D and 3D) particle-in-cell simulations with a binary collision model are used to explore kinetic effects arising during the mixing of plasma media. The applicability of the one-dimensional (1D) ambipolarity condition is evaluated in 2D and 3D simulations of a plasma interface with a sinusoidal perturbation. The 1D ambipolarity condition is found to remain valid in 2D and 3D, as electrons and ions flow together required for $J$ = 0. Simulations of perturbed interfaces show that diffusion-induced total pressure imbalance and hydroflows flatten fine interface structures and drive rapid atomic mix. The atomic mix rate from a structured interface is faster than the ~$\sqrt{\bar{t}}$ scaling obtained from 1D theory in the small-Knudsen-number limit. Plasma kinetic effects inhibit the growth of the Rayleigh-Taylor instability at small wavelengths and result in a nonmonotonic growth rate scaling with wavenumber k with a maximum at a low k value, much different from $\sqrt{Agk}$ (where A is the Atwood number and g is the gravitational constant) as expected in the absence of plasma kinetic effects. Simulations under plasma conditions relevant to MARBLE separated-reactant experiments on Omega and the NIF show kinetic modification of DT fusion reaction rates. With non-Maxwellian distributions and relative drifts between D and T ions, DT reactivity is higher than that inferred from rates using stationary Maxwellian distributions. Reactivity is additionally revealed to be reduced in the presence of finite-Knudsen-layer losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Transport Processes in a Li-ion Cell during an Internal Short-Circuit

Internal short-circuit in a lithium-ion cell causes an abrupt increase in cell temperature and triggers subsequent thermal runaway. In this work, we present a detailed electrochemical-thermal model to investigate the physical behavior during an internal short-circuit. Simulations at wide range of heat transfer coefficients and short-circuit resistances are conducted to illustrate electrochemical and thermal behavior under a wide range of conditions. The Joule heating at the shorted region promotes electrochemical reactions nearby, causing in-plane non-uniformity of electrolyte and active material transport. Furthermore, it is found that diffusion in solid active materials plays a significant role at very high shorting currents (~20 C), as electrochemical reactions rate are being controlled progressively more by availability of Li + at the interface, due to limitations in diffusion through the active material with increasing discharge rates. This diffusion limitation causes a drop in available energy, and subsequently a decrease in cell equilibrium potential and the heat generation rate at the location of the short. On the other hand, rapid depletion of lithium concentration in the electrolyte and accumulation on the electrode surface results in highly non-uniform transport properties resulting in higher heat generation rates. Hence, the heating regime shifts from "local heating" to "global heating". Based on the findings, important design parameters for battery safety are discussed.

25 ENERGY STORAGE↗
DOE OSTI2022

Macrovoid resolved simulations of transport through HPRO relevant membrane geometries

Modeling the transport properties such as diffusivity and permeability of high pressure reverse osmosis (HPRO) membranes is critical for the selection and manufacture of membranes suitable for operation under high pressure. These properties can be significantly affected by the changes in heterogeneous pore structures due to compaction. Here, the modeling platform presented in this work resolves the two scale porosity in HPRO relevant membranes. A synthetic membrane geometry is constructed based on available experimental visualizations of pore structures. The simulations directly capture the flow channeling that results from a combination of material properties and geometric features of the macrovoids. A parametric study is presented to account for the transition of flow characteristics from a material governed regime to a macrovoid governed regime. Permeability of the membrane is evaluated using the simulation data and compared with an existing model that scales with the square of tortuosity over a range of material properties. The model is found to perform well under a narrow range of tortuosity while deviating from the calculated permeabilities at extreme conditions. The effective permeability of the membrane is found to vary by at least two orders of magnitude between the two flow regimes. It is also observed that tortuosity is a bounded property with its upper limit determined by the macrovoid geometry. Consequently, the tortuosity based correlations fail near a flow regime that is mainly governed by the macrovoids. The modeled permeability can be more than an order of magnitude smaller than the simulation result. A new model based on flux-weighted porosity of a membrane is introduced and its correlation with tortuosity is studied. The model agrees with the simulated data as, in addition to tortuosity, it also accounts for the flux partition within and outside the flow channels. Such correlations enable extending the existing understanding of flow characteristics to enhance predictability of porous media models.

42 ENGINEERING↗
DOE OSTI2023

Impact of phase state and non-ideal mixing on equilibration timescales of secondary organic aerosol partitioning

Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Gas Flow and Ion Transfer in Heated ESI Capillary Interfaces

Transfer capillaries are the preferred means to transport ions, generated by electrospray ionization, from ambient conditions to vacuum. During the transfer of ions through the narrow, long tubes into vacuum, substantial losses are typical. However, recently it was demonstrated that these losses can be avoided altogether. To understand the experimental observation and provide a general model for the ion transport, here, we investigate the ion transport through capillaries by numerical simulation of interacting ions. The simulation encompasses all relevant factors, such as space charge, diffusion, gas flow, and heating. Special attention is paid to the influence of the gas flow on the transmission and especially the change imposed by heating. The gas flow is modeled by a one-dimensional gas dynamics description. A large number of ions are treated as point particles in this gas flow. This allows to investigate the influence of the capillary heating on the gas flow and by this on the ion transport. The results are compared with experimental findings. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Enhanced water use efficiency in global terrestrial ecosystems under increasing aerosol loadings

Aerosols play a crucial role in the climate system, affecting incoming radiation and cloud formation. Based on a modelling framework that couples ecosystem processes with the atmospheric transfer of radiation, we analyze the effect of aerosols on surface incoming radiation, gross primary productivity (GPP), water losses from ecosystems through evapotranspiration (ET) and ecosystem water use efficiency (WUE, defined as GPP/ET) for 2003–2010 and validate them at global FLUXNET sites. The total diffuse radiation increases under relatively low or intermediate aerosol loadings, but decreases under more polluted conditions. We find that aerosol-induced changes in GPP depend on leaf area index, aerosol loading and cloudiness. Specifically, low and moderate aerosol loadings cause increases in GPP for all plant types, while heavy aerosol loadings result in enhancement (decrease) in GPP for dense (sparse) vegetation. On the other hand, ET is mainly negatively affected by aerosol loadings due to the reduction in total incoming radiation. Finally, WUE shows a consistent rise in all plant types under increasing aerosol loadings. Overall, the simulated daily WUE compares well with observations at 43 eddy-covariance tower sites (R 2 =0.84 and RMSE=0.01gC (kg H 2 O) -1 ) with better performance at forest sites. In addition to the increasing portions of diffuse light, the rise in WUE is also favored by the reduction in radiation- and heat-stress caused by the aerosols, especially for wet and hot climates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Dopant-controlled single-electron pumping through a metallic island

We investigate a hybrid metallic island/single dopant electron pump based on fully depleted silicon-on-insulator technology. Electron transfer between the central metallic island and the leads is controlled by resonant tunneling through single phosphorus dopants in the barriers. Top gates above the barriers are used to control the resonance conditions. Applying radio frequency signals to the gates, non-adiabatic quantized electron pumping is achieved. A simple deterministic model is presented and confirmed by comparing measurements with simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Hydra-TH Extensions for Multispecies and Thermosolutal Convection

This report describes the Consortium for Advanced Simulation of Light Water Reactors (CASL) work conducted for completion of the Thermal Hydraulics Methods (THM) Level 3 Milestone THM.CFD.P11.02: Hydra-TH Extensions for Multispecies and Thermosolutal Convection. A critical requirement for modeling reactor thermal hydraulics is to account for species transport within the fluid. In particular, this capability is needed for modeling transport and diffusion of boric acid within water for emergency, reactivity-control scenarios. To support this need, a species transport capability has been implemented in Hydra-TH for binary systems (for example, solute within a solvent). A species transport equation is solved for the species (solute) mass fraction, and both thermal and solutal buoyancy effects are handled with specification of a Boussinesq body force. Species boundary conditions can be specified with a Dirichlet condition on mass fraction or a Neumann condition on diffusion flux. To enable enhanced species/fluid mixing in turbulent flow, the molecular diffusivity for the binary system is augmented with a turbulent diffusivity in the species transport calculation. The new capabilities are demonstrated by comparison of Hydra-TH calculations to the analytic solution for a thermosolutal convection problem, and excellent agreement is obtained.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

STEADY STATE DUST DISTRIBUTIONS IN DISK VORTICES: OBSERVATIONAL PREDICTIONS AND APPLICATIONS TO TRANSITIONAL DISKS

The Atacama Large Millimeter Array has returned images of transitional disks in which large asymmetries are seen in the distribution of millimeter sized dust in the outer disk. The explanation in vogue borrows from the vortex literature and suggests that these asymmetries are the result of dust trapping in giant vortices, excited via Rossby wave instabilities at planetary gap edges. Due to the drag force, dust trapped in vortices will accumulate in the center and diffusion is needed to maintain a steady state over the lifetime of the disk. While previous work derived semi-analytical models of the process, in this paper we provide analytical steady-steady solutions. Exact solutions exist for certain vortex models. The solution is determined by the vortex rotation profile, the gas scale height, the vortex aspect ratio, and the ratio of dust diffusion to gas-dust friction. In principle, all of these quantities can be derived from observations, which would validate the model and also provide constrains on the strength of the turbulence inside the vortex core. Based on our solution, we derive quantities such as the gas-dust contrast, the trapped dust mass, and the dust contrast at the same orbital location. We apply our model to the recently imaged Oph IRS 48 system, finding values within the range of the observational uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Driven electrostatic phase space vortices in a 1D weakly dissipative Vlasov–Poisson system

In this paper, the effect of collisions on driven electrostatic phase space vortices is analyzed by means of Eulerian simulation for two different collision models. It was demonstrated recently that in the absence of collisions, at late times, steady state phase space vortices manifest to form a plateau in the resonant region of the particle velocity distribution function, due to trapping of particles supporting multiextrema giant phase space vortices (PSVs). In the presence of collisions, over long time, this multiextrema plateau are found to smooth out, since collisions drive the velocity distribution toward Maxwellian, irrespective of how weak the collisions are as long as they are non-zero. In these conditions, kinetic processes and collisionality are found to be in competition, and the evolution of the plasma is found, therefore, to be a result of nontrivial combination of these two effects. An attempt has been made by means of numerical simulations to study the effect of weak collisionality on the electrostatic driven phase space vortices with two types of collision operators: (1) Bhatnagar–Gross–Krook (Krook) collision operator, where the colliding particles can be treated as isolated pairs and, (2) Fokker–Planck (FP) type collision operator (Zakharov–Karpman) in one dimension, where many weak collisions lead to particle diffusion in velocity space. It is shown that depending on the collision model used, the nature of smoothing in velocity space of giant PSVs results in qualitatively very different phase space structures. However, irrespective of the collision model used, excess density fractions over 10% are retained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Q2 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q2: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $\circ$ Run CESOL with BOUT++/Hermes-3 and immersed boundary condition to directly map to wall: • Run BOUT++/Hermes-3 through the IPS workflow to find radial particle and energy diffusivities to match either the Eich or the physics-based scaling of the SOL heat flux width, and • Expand source of first wall heat flux to include charged particles, neutrals, and radiation from the core+edge. 2. Generate medium fidelity parametrized CAD: $\circ$ Develop the TRACER tool to read an existing CAD, regenerate the geometry based on vertex location and connectivity information, define vertex translation and parameters needed for scaling the CAD, and $\circ$ Utilize the FreeGS code to determine CAT PF coil placement, including minimizing the number of coils, coil current, and electromechanical stresses. 3. Utilize plasma loading for engineering analysis: $\circ$ Couple the plasma loading to input for OpenFOAM and demonstrate initial test of thermal analysis of CAT first wall loading with typical DCLL blanket component cooling boundary conditions. 4. Demonstrate nuclear analysis: $\circ$ Apply initial analysis of tritium transport in DCLL blanket by evaluating spatially resolved tritium generation rates, tritium diffusion and convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗