Search NASA⌕ Search

SEARCH · Search NASA

Results for “interface structure prediction”

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

Design of lightweight BCC multi-principal element alloys with enhanced hydrogen storage using a machine learning-driven genetic algorithm

Body-centered cubic (BCC) based multi-principal element alloy (MPEA) hydrides have demonstrated significant potential for compact and efficient hydrogen storage. In this work, we first leverage machine learning (ML) models to predict the hydrogen affinity, storage capacity and phase stability of BCC MPEAs, creating a unique hydrogen-to-metal (H/M) predictor for materials with unprecedented performance. We developed a metaheuristic optimizer high-throughput framework by interfacing ML models with a genetic algorithm for the accelerated search of {Mg, Al, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Nb, Mo} based lightweight BCC MPEAs with improved hydrogen storage characteristics. We report five new MPEAs with a predicted gravimetric hydrogen storage capacity of around 3.5 wt% or more, including Cr 0.09 Mg 0.73 Ti 0.18 (4.25 wt% H) and Cr 0.21 Nb 0.11 Ti 0.35 V 0.33 (3.5 wt% H). The electronic structure of the top-performing composition, Cr 0.09 Mg 0.73 Ti 0.18 , was analyzed using density functional theory (DFT) to understand the reasons for its improved hydrogen storage properties compared to TiFe (1.90 wt% H), LaNi 5 (1.37 wt% H) or BCC MPEAs like TiVNbCr (3.70 wt% H). Temperature-dependent molecular dynamics (MD) studies were further performed on optimized BCC MPEAs to qualitatively study hydrogen mobility and analyze the effect of different elemental composition on bulk hydrogen diffusion. Our findings demonstrate how a ML assisted genetic algorithm framework can be used for efficient search of stable, lightweight and cost-effective MPEAs while minimizing the need for expensive ab initio calculations.

DFT↗
DOE OSTI2026

Development of Machine-Learned Interatomic Potentials to Predict Structure, Transport, and Reactivity in Platinum-Based Fuel Cells

Machine-learned interatomic potentials (MLIPs) have rapidly progressed in accuracy, speed, and data efficiency in recent years. However, training robust MLIPs in multicomponent systems remains a challenge. In this work, we train an MLIP to describe hydrated Nafion ionomers and platinum catalysts, which are important components of fuel cells, by constructing a diverse training set to describe the bulk polymer and interfacial catalyst–polymer interactions well. We use our trained MLIP to study the properties of the platinum–Nafion system, including polymer structure, proton mobility in a bulk Nafion polymer and near a platinum-Nafion interface, and reactions near and far from the interface, finding excellent results for structure and reactions contained within our training set. Transport seems to be well described, with both vehicular transport and Grotthuss hopping captured, although converged calculations of diffusivities were not computed because they require calculations of tens of nanoseconds that are challenging with current state-of-the-art MLIPs. The combined insights that this model provides can be leveraged to optimize fuel cell performance, and the approach can be applied to other chemical processes and devices where structure, transport, and reactivity all contribute to the overall observed performance.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Topological grain boundary segregation transitions

Engineering the structure of grain boundaries (GBs) by solute segregation is a promising strategy to tailor the properties of polycrystalline materials. Solute segregation triggering phase transitions at GBs has been suggested theoretically to offer different pathways to design interfaces, but an understanding of their intrinsic atomistic nature is missing. Here, we combined atomic resolution electron microscopy and atomistic simulations to discover that iron segregation to GBs in titanium stabilizes icosahedral units (“cages”) that form robust building blocks of distinct GB phases. Owing to their five-fold symmetry, the iron cages cluster and assemble into hierarchical GB phases characterized by a different number and arrangement of the constituent icosahedral units. Our advanced GB structure prediction algorithms and atomistic simulations validate the stability of these observed phases and the high excess of iron at the GB that is accommodated by the phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting Adhesion Energies of Metal Nanoparticles to Support Surfaces, Which Determines Metal Chemical Potential versus Particle Size and Thus Catalyst Performance

Improved catalysts and electrocatalysts composed of transition metal nanoparticles dispersed on high-area supports are essential for energy and environmental technologies. The chemical potential of the metal atoms in these supported nanoparticles is an important descriptor that correlates with both their catalytic activity and deactivation rate. This descriptor (μ M ) is predictably determined by the particle size and the adhesion energy per unit area at the metal/support interface (E adh ). We show here that the adhesion energies for different metals on a given support scale linearly with a simple property of the metal: for oxides, it is proportional to the metal oxophilicity, and for the carbon support, it increases linearly with metal carbophilicity (both divided by the area per metal atom). Furthermore, these relationships allow predicting E adh for other metal/support combinations, thus allowing estimation of μ M versus particle size and thereby better structure-based predictions of catalysts’ performance, which can aid in designing improved catalysts.

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

OpenFacadeControl: enabling integration of automated facades with other building systems

Automated facades are, for the most part, still considered as separate from other building systems throughout the design, installation, commissioning, operation, and maintenance cycle. This takes place despite the fact that their energy and comfort performance are deeply interlinked with the operation of lighting and HVAC systems. Over the last two decades, research has shown that there are significant advantages from operating facades as an integrated system with the rest of the building. Nevertheless, significant barriers prevent this type of integration becoming more common. One of them is the lack of a platform that is inexpensive to implement and that easily allows the practical implementation of integrated control algorithms across fenestration and other building systems, using a variety of communications protocols. This is particularly challenging when automated facades are installed in existing buildings, where interaction with legacy building systems that were installed over the past lifetime of the building can require a high degree of interoperability. OpenFacadeControl (OFC) is an open-source controls framework aimed at unified control of facades and other building systems, including the sharing of third-party sensor information. Through leveraging the Volttron controls platform, it allows the integration of systems and sensors that are manufactured by different companies and that use different communications protocols into an ensemble that functions as a single system. OFC is designed to enable integrated control algorithms of varying degrees of complexity, ranging from simple, heuristic controls to more sophisticated approaches like model-predictive control. Use of a research version to test advanced lighting and shading strategies in a full-scale experimental testbed has demonstrated the ease of deploying advanced control solutions using OpenFacadeControl. This paper presents the structure of OpenFacadeControl and a demonstration case showing the use of OFC in laboratory tests of advanced lighting and fenestration controls that coordinated motorized shades communicating via the BACnet building communications standard and lights communicating via internet-protocol-based application programming interface (API), based on the readings of a shared light level sensor communicating via a different API.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

ReaxFF Parameter Set for Boron Clusters and Icosahedral Boron Crystals: Comparison with Density Functional Theory and Machine-Learning Potentials

Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters and similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B 80 clusters, including a low-energy core–shell structure containing a B 12 icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B 80 clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.

boron↗
DOE OSTI2024

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN↗
DOE OSTI2025

Phase-field modeling of diffusion bonding in 316H stainless steel: Impact of processing conditions on grain morphology and bonding quality

A novel multi-phase, multi-component phase‐field model is presented to study the diffusion bonding of 316H stainless steel. Combined with targeted experimental investigations, this model simulates the bond-growth process and predicts the bonding quality. Unlike previous models, our approach captures the simultaneous evolution of voids and grain structures, while quantifying bonding quality using defined bonding ratio. A comprehensive analysis of bond process control is performed by changing temperature, pressure and surface roughness observing the resulting bond structure, which is consistent with experimental observations and analytical predictions. Temperature is determined to be the dominant factor, with the transition from a flat to a robust bond occurring between 1000 °C and 1050 °C. At the ideal bonding temperature of 1050 °C, a surface roughness exceeding 0.6 μm or an applied stress below 4 MPa results in poor bonding quality. Beyond this, higher pressures and smoother surfaces reduce void size, accelerate void shrinkage, and lead to improved bond integrity. This diffuse-interface model can be extended to other material systems if supplied with appropriate thermodynamic and kinetic data. In conclusion, this makes it an effective modeling platform for optimizing high-temperature diffusion bonding and developing reliable bonded components such as compact heat exchangers.

Diffusion bonding↗
DOE OSTI2024

Microstructure development during rapid alloy solidification

Abstract Solidification processing of structural alloys can take place over an extremely wide range of solid–liquid interface velocities spanning six orders of magnitude, from the low-velocity constitutional supercooling limit of microns/s to the high-velocity absolute stability limit of m/s. In between these two limits, the solid–liquid interface is morphologically unstable and typically forms cellular-dendritic microstructures, but also other microstructures that remain elusive. Rapid developments in additive manufacturing have renewed the interest in modeling the high-velocity range, where approximate analytical theories provide limited predictions. In this article, we discuss recent advances in phase-field modeling of rapid solidification of metallic alloys, including a brief description of state-of-the-art experiments used for model validation. We describe how phase-field models can cope with the dual challenge of carrying out simulations on experimentally relevant length- and time scales and incorporating nonequilibrium effects at the solid–liquid interface that become dominant at rapid rates. We present selected results, illustrating how phase-field simulations have yielded unprecedented insights into high-velocity interface dynamics, shedding new light on both the absolute stability limit and the formation of banded microstructures that are a hallmark of rapid alloy solidification near this limit. We also discuss state-of-the-art experiments used to validate those insights. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2024

INSPIRED: Inelastic neutron scattering prediction for instantaneous results and experimental design

Inelastic neutron scattering (INS) has unique advantages in probing how atoms vibrate and how the vibrations propagate and interact. Such dynamic information is crucial in understanding various material properties, from heat capacity, thermal conductivity, phase transitions, and chemical reactions to more exotic quantum behavior. The analysis and interpretation of the INS spectra often start from a model structure of the sample, followed by a series of calculations to obtain the simulated spectra to compare with experiments. The conventional way to perform such calculations usually requires significant time, computing resources, and specialized expertise. Here, we present a new program named INSPIRED (Inelastic Neutron Scattering Prediction for Instantaneous Results and Experimental Design), which enables users to perform rapid INS simulations in several different ways on their personal computers in just a few clicks, with the crystal structure as the only input file. Specifically, the users can choose a pre-trained symmetry-aware neural network (coupled with an autoencoder) to predict the phonon density of states (DOS), 1D S(E) and 2D S(|Q|,E) spectra for any given structure. One can also choose an existing density functional theory (DFT) calculation from a database (containing over 12,000 crystals), and quickly obtain the simulated INS spectra for single crystals and powders. It is also possible to use pre-trained universal machine learning force fields to relax a given crystal structure, calculate the phonon dispersion and DOS, and, subsequently, the INS spectra. All these functions are implemented with a PyQt graphic user interface. Finally, we expect these new tools will benefit broad user communities and significantly improve the efficiency of experiment design, execution, and data analysis for INS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Surface dependence of electronic growth of Cu(111) on MoS 2

Scanning tunneling microscopy shows that copper deposited at room temperature onto a freshly exfoliated MoS 2 surface forms Cu(111) clusters with periodic preferred heights of 5, 8, and 11 atomic layers. These height intervals correlate with Fermi nesting regions along the necks of the bulk Cu Fermi surface, indicating a connection between physical and electronic structures. Density functional theory calculations of freestanding Cu(111) films support this as well, predicting a lower density of states at the Fermi level for these preferred heights. This is consistent with other noble metals deposited on MoS 2 that exhibit electronic growth, in which the metal films self-assemble as nanostructures minimizing quantum electronic energies. Here, we have discovered that it is critical for the metal deposition to begin on a clean MoS 2 surface. If copper is deposited onto an already Cu coated surface, even if the original film displays electronic growth, the resulting Cu film lacks quantization. Instead, the preferred heights of the Cu clusters simply increase linearly with the amount of Cu deposited upon the surface. We believe this is due to different bonding conditions during the initial stages of growth. Newly deposited copper would bond strongly to the already present copper clusters, rather than the weak bonding, which exists to the van der Waals terminated surface of MoS 2 . The stronger bonding with previously deposited clusters hinders additional Cu atoms from reaching their lowest quantum energy state. The interface characteristics of the van der Waals surface enable surface engineering of self-assembled structures to achieve different applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

On the importance of species immiscibility in mixing-layer dynamics at supercritical pressures

The mixing dynamics of injected propellants is a key factor in determining the ignition performance and combustion-instability response of rocket engines, internal combustion engines, and gas turbines. A key source of uncertainty in the prediction of phase-exchange dynamics is the immiscibility of the injected propellants at supercritical pressure conditions, which induces liquid/vapor phase separation and surface-tension dynamics. While experimental observations indicate the presence of liquid/vapor interfacial structures, this interfacial dynamics is typically neglected in numerical analyses. To address this issue, the objective of the present study is to systematically evaluate the importance of species immiscibility on the phase-exchange dynamics of cryogenic LOX/GH2 mixing layers at typical rocket engine injection conditions. This is accomplished by comparing simulations of (i) a recently developed interface-capturing Regularized-Interface Method (RIM) formulation, and (ii) the commonly employed Diffuse-Interface Method (DIM) formulation. Analysis shows that the interfacial dynamics significantly impact the atomization and mixing of the propellants in the near-injector region, which is not captured by the DIM formulation. The findings of this study extend to other immiscible injection systems, such as LOX/kerosene and hydrocarbon-fuel/air, thereby highlighting the importance of resolving species immiscibility in simulating high-pressure combustion engines.

diffused-interface method↗
DOE OSTI2024

Electric field sensitivity of molecular color centers

Molecular color centers with S=1 ground states are promising candidates for quantum sensing of electric fields. These molecules have an electronic structure similar to solid state color centers, but they allow for processing modalities that permit direct interfacing with an analyte. Currently, it is unknown how sensitive these molecules are to electric fields and what molecular properties affect their sensitivity. We perform density functional theory calculations to understand the impact of electric fields on the electronic structure of five nominally tetrahedral molecular color centers exhibiting variable transition metal chemistry and ligand densities. We then extract the Stark parameters from each of these molecules and compare them to molecular properties such as the dipole moment and inner shell stiffness and find that the dipole moment of the molecule largely governs sensitivity. We predict that polar heteroleptic molecules may have electric field sensitivities comparable to solid state color centers such as nitrogen-vacancy centers in diamond.

Physics↗
DOE OSTI2025

Feature engineering descriptors, transforms, and machine learning for grain boundaries and variable-sized atom clusters

Abstract Obtaining microscopic structure-property relationships for grain boundaries is challenging due to their complex atomic structures. Recent efforts use machine learning to derive these relationships, but the way the atomic grain boundary structure is represented can have a significant impact on the predictions. Key steps for property prediction common to grain boundaries and other variable-sized atom clustered structures include: (1) describing the atomic structure as a feature matrix, (2) transforming the variable-sized feature matrix to a fixed length common to all structures, and (3) applying a machine learning algorithm to predict properties from the transformed matrices. We examine how these steps and different combinations of engineered features impact the accuracy of grain boundary energy predictions using a database of over 7000 grain boundaries. Additionally, we assess how different engineered features support interpretability, offering insights into the physics of the structure-property relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

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

Magic Diamond: Covalent Bond Formation of Melamine and Other Amines on Nanodiamond Surfaces

High-temperature, high-pressure (HPHT) nanodiamond (ND) hosts nitrogen-vacancy (NV) centers, solid-state qubits that enable room-temperature quantum sensing by all-optical magnetometry, electrometry, and thermometry. However, the covalent surface functionalization of nanoscale diamond remains largely limited to carboxylate-based chemistries. Amine termination is particularly attractive because theoretical studies predict suppression of midgap states and extended electron-spin coherence times. Recently, chemical activation of alcohol-terminated NDs to alkyl bromides (ND-Br) using SOBr2 has enabled nucleophilic substitution through a carbocation intermediate, allowing formation of simple amine terminations. Here, we evaluate whether sterically demanding amines can form covalent diamond−nitrogen bonds on ND-Br surfaces. ND-Br was reacted with branched, linear, and cyclic amines, including polyethylenimine, diethylenetriamine, and melamine. X-ray spectroscopies were used to confirm successful and to probe the resulting electronic structure at the diamond−amine interface. These results expand the chemical toolbox for tuning diamond surface dipoles and electron affinity, providing new pathways for engineering nanodiamond surfaces for quantum sensing and photocatalysis applications.

Amines↗
DOE OSTI2024

Ion-Depleted Microenvironments During Lithium Deposition Revealed by Operando Freezing Cryogenic Electron Microscopy

Local structures and chemistry at active electrochemical interfaces are critical to determining safety, lifetime, and energy density in lithium metal batteries and other devices, but they are challenging to characterize at the nanoscale. We address this issue by developing operando freezing cryogenic electron microscopy (cryo-EM) to preserve battery interfaces in an active state for subsequent high-resolution characterization. We find that ion-depleted microenvironments form locally in the electrolyte adjacent to the lithium deposition interface in lithium metal batteries and are linked to heterogenous growth morphologies. These depleted environments arise locally even under conditions for which ion depletion is not predicted at steady state; this provides a mechanistic explanation for why dangerous lithium morphologies can still propagate in such systems and lead to thermal runaway. Operando freezing cryo-EM thus provides a method to directly visualize nanoscale heterogeneities that arise locally at electrochemical interfaces and play key roles in device failure.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,E↗
DOE OSTI2026

Cathodic Protection Modeling for Hanford Underground Double-Shell Tank Farms

Hanford stores millions of gallons of radioactive and chemically hazardous waste from the production of weapon materials in tank farms consisting of underground carbon-steel storage tanks surrounded by reinforced concrete. Six of these Hanford tank farms use double-shell storage tanks (DSTs). The DST farms were constructed from 1968 to 1986 with a planned 40–50 year design life, so some are already operating beyond their initial life expectancy. Ultrasonic testing (UT) has indicated significant thinning on the bottom of the secondary (outer) liner of these tanks, believed to arise from groundwater intrusion driving concrete side corrosion. There is no direct access to the steel/concrete interface between the tank and the concrete pad, making it difficult to apply a chemical-based mitigation strategy or to conduct repairs, but cathodic protection (CP) is a possible method to inhibit further concrete-side corrosion. Hanford already uses CP to protect below grade steel piping within the tank farms and connected to the tanks, but this system was not designed to protect the tank bottoms. CP design must account for the structures surrounding the DSTs, including the steel reinforcing bars (rebar) within the concrete pad and vault, various process lines, and the existing CP system. In this study, finite element analysis (FEA) modeling was carried out to simulate CP protection of 1) a single tank and CP anode to develop options for modeling the rebar and to compare to a simpler circuit model and 2) the entire Hanford AN tank farm as a representative example consisting of seven tanks, associated piping, and both existing and new CP anodes. Both circuit and FEA models predict that significant protective current could be delivered to the bottoms of the tanks with the addition of tank-protection anodes below the depth of the tanks. Simulations with only the existing pipe-protection anodes active confirmed that only a very small current to the tank bottoms is predicted under present conditions. Multiple simplified representations of the dome and wall rebar were tested to reduce the computational complexity of the tank-farm simulations, resulting in modeling the rebar as edge elements with a prescribed effective circumference that matches the real rebar surface area. The geometry of the rebar is also simplified into horizontal hoops around the tank walls and radial rebar over the dome with increased effective circumference to retain the target surface area. This simplification was found to greatly reduce the complexity and solution time of the models without large changes in current distributions, especially to the tank bottom. A range of values were tested for model parameters such as soil and concrete resistivities and polarization resistance to investigate their impact on the current and electric potential distributions. Depending on the parameters used, FEA simulations predict some risk of overprotection, particularly on the piping system; since overprotection can also lead to surface damage associated with hydrogen gas generation at the interface (e.g. hydrogen embrittlement or damage to coatings), this needs to be considered when refining the design of the new CP system. Comparison between the FEA models and the circuit model representation demonstrated that the circuit model could not match the predicted FEA current distribution, even when using the exact same surface areas. This discrepancy appeared to be at least partly attributable to the impact of the relative positions of the tank components and anodes to each other and to the ground surface. The FEA model accounts for the relative positions since it solves the governing equations in three dimensions, but the circuit model cannot account for the positioning. In particular, the circuit model underpredicts the current to the tank bottom and overpredicts the current to the dome compared to FEA for the baseline geometry. The FEA models omitted the electrically isolated rebar in the bottom concrete slab. However, a circuit based stray current model estimated that only 2.1% of the total current through the slab would stray into the rebar, corresponding to ~0.21 A for a target current density of 2 mA/ft2 to the tank bottom. The estimated corrosion driven by this amount of stray current is predicted to yield a lifetime of >400 years for the minimum rebar diameter, assuming an acceptable cross-section area loss of 10%.

d'Entremont, Anna [Savannah River National Laborat↗