Search NASA⌕ Search

SEARCH · Search NASA

Results for “corrosion simulation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Simulating microgalvanic corrosion in alloys using the PRISMS phase-field framework

In this prospective paper, we first review the existing simulation tools to simulate microgalvanic corrosion during free immersion. Then, we describe a recently developed application that employs PRISMS-PF, an open-source, high-performance phase-field modeling framework. The model employed in the application accounts for the electrochemical reaction at the metal/electrolyte interface and ionic migration in the electrolyte to determine the evolution of the corrosion front. We present the implementation details for the application and discuss its features such as super-linear parallel scaling performance for a sufficiently large system. Finally, we demonstrate the capability of the application by simulating corrosion of the matrix phase of an alloy near a secondary phase particle in two and three dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Interface stability of ultrasonic additively manufactured Zircaloy-4 during hydrothermal corrosion

Simulated pressurized water reactor conditions (330 °C, 15.6 MPa, ~20 ppb oxygen) without irradiation were used to investigate the hydrothermal corrosion behavior of ultrasonic additively manufactured Zircaloy-4 up to 1000 h. X-ray computed tomography allowed for visualization of defects from processing and their progression after corrosion experiments. The specimens were found to have clear variability in the mass change data, compared to typical wrought Zircaloy-4 specimens. The variation in the mass change after exposure was attributed to weld defects connected to the specimen surface which allowed ingress of oxidant into the samples. Defects visualized by computed tomography were found via metallography and characterized. In conclusion, ultrasonic additively manufactured Zircaloy-4 was found to have comparable corrosion behavior as wrought Zircaloy-4 for specimens which did not have clear surface defects along weld interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Localised instability of titanium during its erosion-corrosion in simulated acidic hydrometallurgical slurries

Electrochemical techniques were used to investigate the erosion-corrosion of titanium in simulated acidic mineral leaching slurries. Erosion-corrosion of titanium was caused by solid particle impingement. Electrochemical noise revealed that solid particle impacts resulted in localised fracture of the passive film, and erosion-corrosion of titanium proceeded in the form of current transients. As conditions become more abrasive, erosion-corrosion is an increasing threat to titanium equipment exposed to acidic slurries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Development of a numerical model for simulating stress corrosion cracking in spent nuclear fuel canisters

Abstract Prediction and detection of the chloride-induced stress corrosion cracking (CISCC) in Type 304 stainless steel spent nuclear fuel canisters are vital for the lifetime extension of dry storage canisters. This paper conducts a critical review that focuses on the numerical modeling and simulation on the research progress of the CISCC. The numerical models emphasizing the residual stress, susceptible microstructure, and corrosive environment are summarized individually. Meanwhile, the simulation studies on the role of hydrogen-assisted cracking are reviewed. Finally, a multi-physical numerical model, which combines the different fields is proposed based on our recent investigation.

Materials Science↗
DOE OSTI2023

Phase Field Modeling of Corrosion for Design of Next-Generation Magnesium-Aluminum Vehicle Joints

In this study, WPI and its partners sought to understand galvanic corrosion in Mg-Al friction stir weld (FSW) joints. Unlike fusion welds, FSW is done in the solid state. This limits (but does not eliminate) the formation of brittle Al-Mg intermetallics, resulting in better joint strength and corrosion resistance. This study aimed to understand corrosion behavior of Mg-Al FSW joints by using a Cahn-Hilliard phase field electrochemistry model based on the work of Pongsaksawad et al. Its context challenge problem was an ultra-light door designed by Magna with 6061-Al sheet outer and ZEK100-Mg sheet inner panels, joined by an adhesive and FSW. Unlike most FSW joints, which weld through the softer material into the harder one, the Al sheet’s hem joint around Mg required the FSW tool to through the harder Al into the softer Mg, which was itself a challenge. Despite this challenge, PNNL succeeded in consistently making very strong FSW lap joints, by using a triflute tool and power control. Characterization showed mostly separate Al and Mg regions in the weld, usually with a hook protruding from Al into Mg, which likely contributed to high joint strength. Nanohardness mapping showed higher hardness in thin Mg swirls into the Al nugget, likely due to the intermetallics. The study used the industry standard SAE J2334 Cyclic Corrosion Test (CCT) to simulate corrosion conditions over the life of a vehicle, as well as linear polarization testing and ASTM G71 pitting corrosion testing. Corrosion reactions were: H 2 O + ½O 2 + 2e - → 2OH - at the Al cathode, and either Mg → Mg 2+ + 2e - or Mg + 2OH - → Mg(OH) 2 + 2e - at the Mg anode, with electron transfer through the joint. Mass loss increased with Al section length, supporting the hypothesis that galvanic corrosion in Al-Mg joints is limited by the cathodic reaction. Both pitting corrosion (first anode reaction) and hydroxide film formation (second reaction) were observed. The J2334 CCT test method was slightly modified, adding a chromium solution cleaning step after each week of corrosion testing, in order to remove Mg oxides and hydroxides and accurately measure mass loss over time. That said, FSW joints treated in this way did not did not exhibit significant reduction in lap shear strength of the joints vs. newly welded samples which had not undergone corrosion. Corrosion appeared to consist mainly of pitting in Mg in a way which did not directly affect the bond between the materials where the joints failed. A Cahn-Hilliard phase field model coupled with electrical potential described the electronically mediated galvanic corrosion electrochemical reactions described above in the four-component Al-Mg-H 2 O-(H 2 O+½O 2 ) system. This 2-D model successfully predicted Mg(OH)₂ formation at the anode in some circumstances, and Mg pitting corrosion in others, and showed the correct electric field directions in both cases. The model predicted Mg pit corrosion depth within a factor of two of measured pits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

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

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

Ramteke, Rajat Rajat Durgesh Ramteke [University o↗
DOE OSTI2020

In situ synchrotron X-ray tomography of 304 stainless steels undergoing chlorine-induced stress corrosion cracking

To understand the mechanisms behind chlorine-induced stress corrosion cracking (CISCC) in 304 stainless steel, synchrotron high energy x-rays were used to perform in-situ x-ray microtomography on a pre-cracked sample under load in a corrosive simulated marine environment. The tomography scans showed the crack morphology evolving into a branching crack around the surface of the sample. Finite element analysis and stress intensity analysis are conducted to explain the observed crack branching.

304 Stainless Steel↗
DOE OSTI2022

MOSCATO Development and Integration in Fiscal Year 2022

During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. We built upon the work for FY22 and improved the code from several different perspectives. First, we improved the user interface by adding a new component to the official Nek5000 input file (.par). This new part contains documents parameters like, salt properties (density, viscosity, Cp, thermal conductivity), diffusion coefficients, standard potential, etc. Second, we built a conversion script to extract salt properties from the MSTDB-TP salt database and write to MOSCATO input file. Third, we migrated the code to NekRS, which is the GPU branch of Nek5000 and suitable for next generation supercomputers. Verification and Validation (V&V) work was also continued in FY22. Two tasks were performed. The first V&V task involved the validation of the Poisson-Nernst-Planck equation solver and Butler-Volmer electrode kinetics, by comparing with numerical and experimental data about thermoelectric cells. The second task involved the comparisons to corrosion results from a thermal convection loop run during the MSRE era. Satisfactory agreement was obtained from both tasks.

Yuan, Haomin↗
DOE OSTI2021

MOSCATO Solver Development and Integration Plan

During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. Although capable of large-scale simulations, MOSCATO is not well-suited to parametric studies of complete reactor geometries. These types of simulations are instead better handled by reduced-order modeling codes such as ORNL’s Mole code. Reduced-order simulation tools like Mole, however, are dependent on high fidelity correlations to account for complex, coupled three-dimensional phenomena that they do not directly simulate. Tools such as MOSCATO must therefore be used to create these correlations, as suitable empirical relationships are not available for most molten salt systems. Toward that end, we used the Nek-derived version of MOSCATO to create new mass transfer correlations for three relevant cases including tubular, tube bank, and subchannel geometries. These new correlations are more accurate than any existing ones and can be readily integrated into any reduced-order modeling tools that are targeting full-scale MSR simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Controlled thermal pre-treatment of ZMG232G10® for corrosion mitigation under simulated SOFC interconnect exposure conditions

The IT-SOFC candidate interconnect material ZMG232G10® was thermally pre-treated in air and in H 2 – 3% H 2 O to mitigate dual atmosphere corrosion. Bare and pre-treated steel samples were exposed to a dual atmosphere exposure condition (dry air/metal/4% H 2 – N 2 (bal.) + 3% H 2 O) for 500h at 600°C. Corrosion products and oxide morphology developed on the air-exposed surfaces were analyzed. Pre-treatments of the metal in oxidizing (air) and reducing (H 2 – 3% H 2 O) gas atmospheres provide improved passivation against dual atmosphere corrosion compared to bare steel. Further, the uniform and iron-free oxide scale developed during the pre-treatment of alloy in the low pO 2 H 2 – 3% H 2 O gas atmosphere provides an effective diffusion barrier against the outward transport of cations and oxidation of iron observed under dual atmosphere conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Iron surface corrosion in supercritical CO2 at atomic scale investigated by molecular dynamics simulations

Understanding the corrosion behavior of steels in supercritical carbon dioxide (S-CO2) is essential for ensuring the safe application of S-CO2 as a heat-transfer fluid in high-temperature energy systems, including advanced nuclear reactors. In this work, molecular dynamics (MD) simulations using ReaxFF potential are performed to explore the atomic-scale corrosion mechanisms of body-centered cubic iron (BCC-Fe) in S-CO2. The results show that CO2 molecules in S-CO2 decompose at the Fe surface, generating free C and O atoms that form Fe-C and Fe-O bonds and subsequently produce oxides and carbides. Concurrently, Fe atoms dissolve from the surface and diffuse into the S-CO2 region, resulting in interdiffusion of Fe, C and O atoms at the interface. The corrosion-layer thickness calculations show that high pressure and temperature induced by S-CO2 have stronger effects than surface orientation on the corrosion process. In addition, surface Fe atoms undergo substantial displacement under S-CO2 exposure, further accelerating corrosion. When a radiation-induced void is introduced near the Fe surface, the corrosion is enhanced. The void-matrix interface expands the reaction surface area and simultaneously induces corrosion reactions inside the bulk, resulting in a deeper penetration of C and O and thicker corrosion layers. All these results indicate that high-temperature, high-pressure and radiation-induced voids can seriously affect the corrosion of Fe in S-CO2, and must be considered to better use S-CO2 in nuclear facilities.

Li, Wenhua↗
DOE OSTI2024

Critical role of slags in pitting corrosion of additively manufactured stainless steel in simulated seawater

Abstract Pitting corrosion in seawater is one of the most difficult forms of corrosion to identify and control. A workhorse material for marine applications, 316L stainless steel (316L SS) is known to balance resistance to pitting with good mechanical properties. The advent of additive manufacturing (AM), particularly laser powder bed fusion (LPBF), has prompted numerous microstructural and mechanical investigations of LPBF 316L SS; however, the origins of pitting corrosion on as-built surfaces is unknown, despite their utmost importance for certification of LPBF 316L SS prior to fielding. Here, we show that Mn-rich silicate slags are responsible for pitting of the as-built LPBF material in sodium chloride due to their introduction of deleterious defects such as cracks or surface oxide heterogeneities. In addition, we explain how slags are formed in the liquid metal and deposited at the as-built surfaces using high-fidelity melt pool simulations. Our work uncovers how LPBF changes surface oxides due to rapid solidification and high-temperature oxidation, leading to fundamentally different pitting corrosion mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accelerating charge estimation in molecular dynamics simulations using physics-informed neural networks: corrosion applications

Molecular Dynamics (MD) simulations are used to understand the effects of corrosion on metallic materials in salt brine. Reactive force fields in classical MD enable accurate modeling of bond formation and breakage in the aqueous medium and at the metal-electrolyte interface, while also facilitating dynamic partial charge equilibration. However, MD simulations are computationally intensive and unsuitable for modeling the long time scales characteristic of corrosive phenomena. To address this, we develop reduced-order machine learning models that provide accurate and efficient predictions of charge density in corrosive environments. Specifically, we use Long Short-Term Memory (LSTM) networks to forecast charge density evolution based on atomic environments represented by Smooth Overlap of Atomic Positions (SOAP) descriptors. A physics-informed loss function enforces charge neutrality and electronegativity equivalence. The atomic charges predicted by the deep learning model trained on this work were obtained two orders of magnitude faster than those from molecular dynamics (MD) simulations, with an error of less than 3% compared to the MD-obtained charges, even in extrapolative scenarios, while adhering to physical constraints. This demonstrates the excellent accuracy, computational efficiency, and validity of the developed model. Lastly, even though developed for corrosion, these protocols are formulated in a phenomenon-agnostic manner, allowing application to various variable-charge interatomic potentials and related fields.

Atomistic models↗
DOE OSTI2020

Drying Wet Stored and Corroded Magnox Fuel for Interim Dry Storage - 20187

To investigate the feasibility of drying corroded Magnesium alloy clad nuclear fuel, the corrosion products of inactive and unirradiated Magnox simulant dry corrosion product and dried corroded Magnox sludge were characterised by X-ray diffraction and thermogravimetric analysis. XRD identified brucite in the dried corrosion product and both brucite and hydromagnesite in the sludge. TGA identified mass loss on dehydration of ∼9.1% up to 250 deg. C and 27.9% from dehydroxylation up to 440 deg. C, with total mass loss of 41.9% up to 800 deg. C for the dried corrosion product. The sludge TGA showed 8.3% mass loss up to 250 deg. C from dehydration, 7.5% up to 330 deg. C from dehydroxylation and 31.4% mass loss up to 47.2% for decarbonation, with total mass loss of 49.2% up to 800 deg. C. Water removal up to 1.8 g was performed by cold vacuum drying (40 deg. C-120 deg. C) on a ≅14.8 g (wet) corroded Magnox sample. The process was monitored by observing pressure, dew point, temperature and gas flow changes supported by measuring the sample mass loss as water is removed. From these tests, it was observed that the dried corroded Magnox displayed some hygroscopic capacity which has implications for water retention following vacuum drying and the length of exposure required to achieve the desired dryness. Whilst at higher temperatures the water was removed faster, the increased temperature also risks the fuel cladding undergoing unwanted chemical reactions with the residual water during the process. However, at lower temperatures the achievable levels of dryness were reduced. Therefore, from these observations this work has identified that there is a temperature balance that may be necessary to optimise the drying process with respect to allowing the greatest level of dryness whilst restricting unwanted chemical reactions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗