Search NASA⌕ Search

SEARCH · Search NASA

Results for “matrix diffusion”

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 685 records · Page 38

Unraveling Shuttle Effect and Suppression Strategy in Lithium/Sulfur Cells by In Situ/Operando X-ray Absorption Spectroscopic Characterization

The polysulfides shuttle effect represents a great challenge in achieving high capacity and long lifespan of lithium/sulfur (Li/S) cells. A comprehensive understanding of the shuttle-related sulfur speciation and diffusion process is vital for addressing this issue. Herein, we employed in situ/operando X-ray absorption spectroscopy (XAS) to trace the migration of polysulfides across the Li/S cells by precisely monitoring the sulfur chemical speciation at the cathodic electrolyte-separator and electrolyte-anode interfaces, respectively, in a real-time condition. After we adopted a shuttle-suppressing strategy by introducing an electrocatalytic layer of twinborn bismuth sulfide/bismuth oxide nanoclusters in a carbon matrix (BSOC), we found the Li/S cell showed greatly improved sulfur utilization and longer life span. The operando S K-edge XAS results revealed that the BSOC modification was bi-functional: trapping polysulfides and catalyzing conversion of sulfur species simultaneously. We elucidated that the polysulfide trapping-and-catalyzing effect of the BSOC electrocatalytic layer resulted in an effective lithium anode protection. Finally, our results could offer potential stratagem for designing more advanced Li/S cells.

25 ENERGY STORAGE↗

Extensive nanoprecipitate morphology transformation in a nanostructured ferritic alloy due to extreme thermomechanical processing

Nano-oxide precipitates in a modern nanostructured ferritic alloy were investigated after extreme thermomechanical processing into a thin-walled tube geometry. It was found that the morphology of the precipitates changed from spherical to rod-shaped, with some increasing to aspect ratios of up to 9, despite the precipitate volume fraction (0.3%) and number density (> 10 23 m -3 ) of precipitates remaining unchanged. High-resolution electron microscopy showed that the precipitates likely remained coherent with the Fe-matrix, while atom probe tomography confirmed that the precipitate compositions remained unaffected by the transformation. The morphological change was attributed to the shearable nature of the (Y,Ti,O)-rich precipitates, indicating they should be considered as “soft” obstacles to dislocation motion. The elongation was most pronounced in larger (>5 nm) precipitates, which may be caused by preferential dissolution of the smallest (1–3 nm) precipitates followed by the competition between re-precipitation and solute diffusion to larger precipitates during recovery heat treatments.

36 MATERIALS SCIENCE↗

In situ TEM probing of crystallization form-dependent sodiation behavior in ZnO nanowires for sodium-ion batteries

Development of sodium-ion battery (SIB) electrode materials currently lags behind electrodes in commercial lithium-ion batteries (LIBs). However, in the long term, development of SIB components is a valuable goal. Their similar, but not identical, chemistries require careful identification of the underlying sodiation mechanism in SIBs. Here in this study, we utilize in situ transmission electron microscopy to explore quite different sodiation behaviors even in similar electrode materials through real-time visualization of microstructure and phase evolution. Upon electrochemical sodiation, single-crystalline ZnO nanowires (sc-ZNWs) are found to undergo a step-by-step electrochemical displacement reaction, forming crystalline NaZn 13 nanograins dispersed in a Na 2 O matrix. This process is characterized by a slowly propagating reaction front and the formation of heterogeneous interfaces inside the ZNWs due to non-uniform sodiation amorphization. In contrast, poly-crystalline ZNWs (pc-ZNWs) exhibited an ultrafast sodiation process, which can partly be ascribed to the availability of unobstructed ionic transport pathways among ZnO nanograins. Thus the reaction front and heterogeneous interfaces disappear. The in situ TEM results, supported by calculation of the ion diffusion coefficient, provide breakthrough insights into the dependence of ion diffusion kinetics on crystallization form. This points toward a goal of optimizing the microstructure of electrode materials in order to develop high performance SIBs.

25 ENERGY STORAGE↗

A Fast Algebraic Multigrid Solver and Accurate Discretization for Highly Anisotropic Heat Flux I: Open Field Lines

We present a novel solver technique for the anisotropic heat flux equation, aimed at the high level of anisotropy seen in magnetic confinement fusion plasmas. Such problems pose two major challenges: (i) discretization accuracy and (ii) efficient implicit linear solvers. We simultaneously address each of these challenges by constructing a new finite element discretization with excellent accuracy properties, tailored to a novel solver approach based on algebraic multigrid (AMG) methods designed for advective operators. We pose the problem in a mixed formulation, introducing the directional temperature gradient as an auxiliary variable. The temperature and auxiliary fields are discretized in a scalar discontinuous Galerkin space with upwinding principles used for discretizations of advection. We demonstrate the proposed discretization’s superior accuracy over other discretizations of anisotropic heat flux, achieving error 1000x smaller for anisotropy ratio of 10 9 , for closed field lines. The block matrix system is reordered and solved in an approach where the two advection operators are inverted using AMG solvers based on approximate ideal restriction, which is particularly efficient for upwind discontinuous Galerkin discretizations of advection. To ensure that the advection operators are nonsingular, in this paper we restrict ourselves to considering open (acyclic) magnetic field lines for the linear solvers. We demonstrate fast convergence of the proposed iterative solver in highly anisotropic regimes where other diffusion-based AMG methods fail.

97 MATHEMATICS AND COMPUTING↗

Kinetics and transport of hydrogen in graphite at high temperature and the effects of oxidation, irradiation and isotopics

The kinetics of uptake and desorption impact the performance of graphite as a vector for tritium in high-temperature fission reactors and in the blanket of fusion reactors. Graphite components in these reactors are exposed to temperatures > 500 °C and H 2 partial pressures of few Pa and desorption temperatures are limited to < 1600 °C; limited data is available at these conditions. Here we review the mechanisms for uptake in, transport and desorption of hydrogen from graphite at high temperature, compiling data on uptake rates, diffusion coefficients and activation energies and providing a discussion of the impact of irradiation, pre-oxidation and isotope. At FHR conditions, trapping impacts uptake rates, leading to a reduction in apparent diffusivity by 35 to 80% compared to higher partial-pressure uptake. Timelines for desorption are not clearly defined; extrapolating from available data, at 1150 °C desorbing 80% of tritium uptaken at FHR conditions may take from 100 to 10,000 h.

36 MATERIALS SCIENCE↗

Photocured epoxy/graphene nanocomposites with enhanced water vapor barrier properties

A transparent, water vapor barrier film made of an epoxy resin and graphene oxide (GO) was synthesized by photopolymerization process. The epoxy/GO film with just 0.05 wt% GO gives a 93% WVTR reduction with respect to the pristine polymer, reaching barrier properties better than other polymer composites containing higher amounts of graphene. The excellent water vapor barrier is attributed to the good dispersion of GO in the polymer matrix. Moreover, GO significantly enhances the toughness and the damping capacity of the epoxy resins. The hybrid film can have potential applications in anticorrosive coatings, electronic devices, pharmaceuticals and food packaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A bi-layer barrier design for 122-type iron-based superconducting wires and tapes

Iron-based superconducting wires and tapes hold great promise for high-field magnet applications. A promising design for 122-type wires and tapes based on the powder-in-tube method is using silver and copper double-layer sheaths. For this design a heat treatment temperature below ~ 779 °C is required to prevent Ag-Cu liquid formation. However, this may be below the optimal heat treatment temperature for the critical current density, and still cannot prevent Ag-Cu interdiffusion occurring in the solid state. In this work we propose adding a niobium or tantalum or vanadium (or their alloys) barrier layer between the Ag and Cu to solve the Ag-Cu interdiffusion issue, given that the group-VB metals (vanadium, niobium, tantalum) are relatively inert to both Ag and Cu. To investigate the effectiveness of this design, BaFe 1.84 Co 0.16 As 2 wires and tapes with Ag/Cu and Ag/Ta/Cu sheaths, as well as Ba 0.6 K 0.4 Fe 2 As 2 wires and tapes with Ag/Cu and Ag/Nb/Cu sheaths, were fabricated. It was found that both the Ta and Nb layers kept integral after wire drawing, but after a large flat-rolling reduction the Ta layer broke while the Nb layer kept integral. In the tapes with Ag/Cu sheaths (without the Ta or Nb layer) Cu diffused through the Ag layer and into the powder cores during 740 °C heat treatment, while in the tapes with Ag/Nb/Cu sheaths the Nb layer effectively blocked Ag-Cu interdiffusion even at 900 °C. Further, this work demonstrates that Ta is a suitable barrier material for 122-type wires, while Nb is suitable for both wires and tapes. In this design using Ag/Nb (or Ta)/Cu sheaths, we can regard the outer Cu as the conductor matrix while the Ag and Nb (or Ta) serve as two layers of barriers that suppress reactions between the components. Thus, we call this design a “bi-layer barrier” design for 122-type wires and tapes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pore-scale prediction of transport properties in reconstructed nanostructures of organic matter in shales

Size, morphology and distributions of pores in organic matter of shale matrix are discussed based on high resolution images from experiments in the literature. 250 nanoscale structures of the organic matter are then reconstructed by randomly placing pore spheres with different diameters and overlap tolerances. Effects of porosity, the mean diameter and the overlap tolerance on void space connectivity and pore size distribution are studied. Furthermore, a pore-scale model based on the lattice Boltzmann method developed in a previous study is used to predict the Knudsen diffusivity and permeability of the reconstructed organic matter. The simulation results show that the mean pore diameter and overlap tolerance significantly affect the transport properties. The predicted Knudsen effective diffusivity is compared with Bruggeman equation and it is found that this equation underestimates the tortuosity. A modified Bruggeman equation is proposed based on the simulation results. The predicted intrinsic permeability is in acceptable agreement with Kozeny–Carman (KC) equation. In addition, the apparent permeability is determined based on Knudsen diffusivity and intrinsic permeability predicted. The apparent permeability is compared with that obtained with various correlations in the literature. Lastly, Knudsen’s correlations match best with our numerical results and are recommended for calculating the apparent permeability.

04 OIL SHALES AND TAR SANDS↗

Accelerating network layouts using graph neural networks

Graph layout algorithms used in network visualization represent the first and the most widely used tool to unveil the inner structure and the behavior of complex networks. Current network visualization software relies on the force-directed layout (FDL) algorithm, whose high computational complexity makes the visualization of large real networks computationally prohibitive and traps large graphs into high energy configurations, resulting in hard-to-interpret “hairball” layouts. Here we use Graph Neural Networks (GNN) to accelerate FDL, showing that deep learning can address both limitations of FDL: it offers a 10 to 100 fold improvement in speed while also yielding layouts which are more informative. We analytically derive the speedup offered by GNN, relating it to the number of outliers in the eigenspectrum of the adjacency matrix, predicting that GNNs are particularly effective for networks with communities and local regularities. Finally, we use GNN to generate a three-dimensional layout of the Internet, and introduce additional measures to assess the layout quality and its interpretability, exploring the algorithm’s ability to separate communities and the link-length distribution. The novel use of deep neural networks can help accelerate other network-based optimization problems as well, with applications from reaction-diffusion systems to epidemics.

97 MATHEMATICS AND COMPUTING↗

Electrochemical Modeling and Experimental Verification of Lithiation Gradients in Oxide Cathodes of Lithium-Ion Cells

Lithiated nickel-cobalt-manganese oxides, such as NCM523, are used in the positive electrode (cathode) of Li-ion cells. Using operando X-ray diffraction profilometry, lithiation gradients in the cathode matrix can be observed and quantified by expansion into Legendre polynomials with time-dependent weights. These weights (referred to as gradients) increase in magnitude when electric current flows through the cell, decrease during potentiostatic hold and finally relax to zero when the current is interrupted during open circuit rest. Both physics-based electrochemical models and operando X-ray experiments suggest that the time constants for gradient growth and abatement are primarily determined by ionic diffusion in the oxide particles, which in turn depends on their lithium content. In contrast, the magnitude of gradients depends mainly on the applied current. The X-ray profilometry provides a way of directly probing the formation and disappearance of Li gradients across the cathode during fast cycling, which can help to diagnose the effects of material degradation in the cells.

25 ENERGY STORAGE↗

Calculation of Higher Eigenpairs of the Transport Equation Using IRAM Based on Domain Decomposition

The higher eigenpairs of the neutron transport (or diffusion) equation are useful in core flux distribution reconstruction, transient analysis, stability analysis, and loading pattern optimization. There are various methods for the calculation of higher eigenpairs, such as the power iteration method with filtering technique, the Generalized Modified Power Method (GMPM), the power iteration method with innovative boundary conditions for symmetric cores, and Krylov subspace methods. In regard to Krylov subspace methods, the Implicitly Restarted Arnoldi Method (IRAM) and ORTHOMIN were implemented to calculate higher eigenpairs for both transport and diffusion equations. In Warsa's work, however, there was a significant drawback that problems with energy upscatter could not be treated efficiently. In Modak's work, the mode to which his method was likely to converge depended on the shape of initial guess flux, and only one mode could be computed each time, making it impractical for real use. Furthermore, the computational burden for higher eigenpairs is generally much larger than for only the first eigenpair, especially for neutron transport equations. Therefore, the implementation of parallel computation for higher eigenpairs is essential to improve the efficiency, and no such work has been reported in literature so far. In this summary, we propose an innovative method to calculate higher eigenpairs of the neutron transport equation based on high performance computation. In this method, 2D Matrix Method of Characteristic (MMOC) is used to discretize the neutron transport equation; IRAM is used to evaluate higher eigenpairs; spatial domain decomposition is used to realize parallel computation. The highlight of this method is using the MGMD PGMRES (Multi-Group and Multi-Domain coupled Parallel Generalized Minimal Residual Method) algorithm to simultaneously accelerate the upscatter iterations and the inner interface angular fluxes iterations, avoiding the drawback in Warsa's method and the efficiency problem of conventional Jacobi iteration. To the best of our knowledge, our work represents the first application of IRAM based on domain decomposition for higher eigenpairs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗

Electrospun ZnO Nanoparticles Doped Core–Sheath Nanofibers: Characterization and Antimicrobial Properties

Coaxial electrospinning technique was used to fabricate the core–sheath composite nanofibers of ZnO nanoparticle (Nps) (10%, 20% w/w) doped polymethyl methacrylate (PMMA) (as sheath) and polyvinyl alcohol (PVA) (as core). Fourier transform infrared (FT-IR) spectra were confirmed the weak forces arise between ZnO Nps, PMMA and PVA matrixes. The hexagonal (wurtzite) structure of ZnO Nps with ~ 30.8 nm of diameter was confirmed from the X-ray diffraction pattern. The morphology and microstructure of core–sheath composite nanofibers were confirmed from the scanning electron microscopy (SEM) and transmission electron microscopy (TEM). It is clearly seen from the TEM images that the PMMA encapsulate the PVA core. Core–sheath composite nanofibers were assessed against Escherichia coli (E. coli) and Bacillus subtilis (B. subtilis) bacteria through quantitative, disk diffusion and viable cell count methods. It was found that ZnO Nps doped core–sheath nanofibers were effectively inhibit the growth of gram positive bacteria, B. subtilis.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Versatile soil gas concentration and isotope monitoring: optimization and integration of novel soil gas probes with online trace gas detection

Abstract. Gas concentrations and isotopic signatures can unveil microbial metabolisms and their responses to environmental changes in soil. Currently, few methods measure in situ soil trace gases such as the products of nitrogen and carbon cycling or volatile organic compounds (VOCs) that constrain microbial biochemical processes like nitrification, methanogenesis, respiration, and microbial communication. Versatile trace gas sampling systems that integrate soil probes with sensitive trace gas analyzers could fill this gap with in situ soil gas measurements that resolve spatial (centimeters) and temporal (minutes) patterns. We developed a system that integrates new porous and hydrophobic sintered polytetrafluoroethylene (sPTFE) diffusive soil gas probes that non-disruptively collect soil gas samples with a transfer system to direct gas from multiple probes to one or more central gas analyzer(s) such as laser and mass spectrometers. Here, we demonstrate the feasibility and versatility of this automated multiprobe system for soil gas measurements of isotopic ratios of nitrous oxide (δ18O, δ15N, and the 15N site preference of N2O), methane, carbon dioxide (δ13C), and VOCs. First, we used an inert silica matrix to challenge probe measurements under controlled gas conditions. By changing and controlling system flow parameters, including the probe flow rate, we optimized recovery of representative soil gas samples while reducing sampling artifacts on subsurface concentrations. Second, we used this system to provide a real-time window into the impact of environmental manipulation of irrigation and soil redox conditions on in situ N2O and VOC concentrations. Moreover, to reveal the dynamics in the stable isotope ratios of N2O (i.e., 14N14N16O, 14N15N16O, 15N14N16O, and 14N14N18O), we developed a new high-precision laser spectrometer with a reduced sample volume demand. Our integrated system – a tunable infrared laser direct absorption spectrometry (TILDAS) in parallel with Vocus proton transfer reaction mass spectrometry (PTR-MS), in line with sPTFE soil gas probes – successfully quantified isotopic signatures for N2O, CO2, and VOCs in real time as responses to changes in the dry–wetting cycle and redox conditions. Broadening the collection of trace gases that can be monitored in the subsurface is critical for monitoring biogeochemical cycles, ecosystem health, and management practices at scales relevant to the soil system.

54 ENVIRONMENTAL SCIENCES↗

Room temperature neutron diffraction and magnetic studies of multiferroic Pb{sub 0.9}Bi{sub 0.1}Fe{sub 0.55}Nb{sub 0.45}O{sub 3} solid solution

The Pb{sub 0.9}Bi{sub 0.1}Fe{sub 0.55}Nb{sub 0.45}O{sub 3} (PBFNO) solid solution was synthesized by single step solid state reaction method and the optimized parameters are 700°C for 2hr (calcination) and 800°C for 3hr (sintering). The formation of desired material was confirmed using X-Ray Diffraction (XRD) and Neutron Diffraction (ND) studies. The structural and magnetic properties of the sintered pellets were investigated at room temperature (RT) through XRD, ND and Magnetic (M-H) studies. The structural analysis was carried out by Rietveld refinement through the Full Prof program. Rietveld refined XRD and ND patterns confirms the monoclinic structure with Cm space group and obtained cell parameters from the ND data are a = 5.651(2) Å, b = 5.658(2) Å, c = 4.005(1) Å and α = 90°, β = 89.98(1)°, γ = 90°. RT M-H curve studies have been carried out. It shows the clear opening of hysteresis (M–H) loop, is evidenced as the existence of weak ferromagnetism at RT. The M-H data shows existence weak ferromagnetism embedded in an antiferromagnetic matrix structure. The dielectric constant with frequency shows the formation of barrier layers at the grain and grain boundary interfaces gives rise to interfacial space charge polarization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Melting and solidification dynamics during laser melting of reaction-based metal matrix composites uncovered by in-situ synchrotron X-ray diffraction

Laser additive manufacturing (AM) of reaction-based metal matrix composites (MMCs) involves highly complex and non-equilibrium material transformation behavior, including melting, dissolution, precipitation, and solidification. Yet, the dynamics and interplay of these phase transformation processes remain poorly understood, posing substantial challenges in identifying the microstructure formation mechanism, and predicting and controlling the microstructure in the printed parts. Here, in this study, we performed the in-situ X-ray diffraction experiment to characterize the phase evolution dynamics of the 316L + 10 vol.%TiC system during laser melting, which provides direct and quantitative insights of the complex phase reaction and evolution dynamics under rapid heating and cooling conditions relevant to additive manufacturing of reaction-based MMCs. Further in-depth thermodynamic and kinetic calculations revealed that most of the phase evolution behavior observed in the in-situ X-ray diffraction experiment cannot be solely explained by widely used equilibrium thermodynamic models, and diffusion-controlled nonequilibrium dissolution and precipitation kinetics must be considered to elucidate the complex phase evolution behavior, including incomplete TiC dissolution, and three-step TiC precipitation. The three distinct types of precipitates generate unique hierarchical TiC micro- and nanostructures, which enhances the yield strength from 513 MPa to 877 MPa by 71 %, tensile strength from 628 MPa to 1054 MPa by 68 %, and Young's modulus from 193 GPa to 221 GPa by 14 %. The findings of our research provide the knowledge foundation for the design of unique microstructures and advanced MMC materials through laser AM.

36 MATERIALS SCIENCE↗

Novel Gamma Attenuation Methods Utilized for Waste Characterization of UO{sub 2}F{sub 2}-Bearing Materials at the Paducah GDP - 17327

The Paducah Gaseous Diffusion Plant (PGDP) is currently under deactivation with demolition in the decades to follow. UO{sub 2}F{sub 2}-bearing waste is generated at PGDP by UF{sub 6} hydrolysing to materials such as pipes, valves, PPE, etc. when UF{sub 6} is exposed to humid air. In addition, Al{sub 2}O{sub 3} and NaF trap materials contain significant amounts of Uranium. The UO{sub 2}F{sub 2} deposits within the generated waste can be of sub-gram levels to several hundreds of localized grams and are much denser than the surrounding interfering matrix materials and may be of very complex geometric and chemical forms such as with the Al{sub 2}O{sub 3} and NaF trap materials. The main objective is to characterize the fissile (U-235) content of the UO{sub 2}F{sub 2}-bearing waste so as to meet the Data Quality Objectives (DQO) for proper disposal. If a waste container cannot meet the DQO then it must be subjected to costly remediation, repackaging and re-characterization processes. We have developed, at PGDP, several novel model-based methods and techniques that vastly improve the standard gamma attenuation correction methods and do not require Working Reference Materials (WRM) for calibration. These new methods involve rotating and non-rotating 18.9 liter (5 U.S. gallon) waste containers with and without a transmission modality. A new model-based variant of differential attenuation has been developed that provides a high-quality self-absorption correction that does not require knowledge of the Uranium enrichment or isotopic ratios. For more challenging waste containers, a method using multi-view (non-rotating) 18.9 liter waste containers has been considered using a recently developed Advanced ISOCS{sup C} Uncertainty Estimator (AIUE) package. The new methods suffice to significantly improve accuracy and precision and to lower the Total Measurement Uncertainty (TMU) leading to satisfactory DQO's that were not previously achievable using the standard gamma attenuation correction methods. In this paper we present real 18.9 liter waste container and Working Reference Material (WRM) data subjected to the new methods and compared to standard methods depicting the significant improvement in the DQO that is realized. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗