Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

A Gas-Kinetic Method for Hyperbolic-Elliptic Equations and Its Application in Two-Phase Fluid Flow

A gas-kinetic method for the hyperbolic-elliptic equations is presented in this paper. In the mixed type system, the co-existence and the phase transition between liquid and gas are described by the van der Waals-type equation of state (EOS). Due to the unstable mechanism for a fluid in the elliptic region, interface between the liquid and gas can be kept sharp through the condensation and evaporation process to remove the "averaged" numerical fluid away from the elliptic region, and the interface thickness depends on the numerical diffusion and stiffness of the phase change. A few examples are presented in this paper for both phase transition and multifluid interface problems.

Xu, Kun↗
DOE OSTI2023

Heterostructure engineering in electrode materials for sodium-ion batteries: Recent progress and perspectives

Sodium-ion batteries (SIBs) have stepped into the spotlight as a promising alternative to lithium-ion batteries for large-scale energy storage systems. However, SIB electrode materials, in general, have inferior performance than their lithium counterparts because Na + is larger and heavier than Li + . Heterostructure engineering is a promising strategy to overcome this intrinsic limitation and achieve practical SIBs. We provide a brief review of recent progress in heterostructure engineering of electrode materials and research on how the phase interface influences Na + storage and transport properties. Efficient strategies for the design and fabrication of heterostructures (in situ methods) are discussed, with a focus on the heterostructure formation mechanism. The heterostructure's influence on Na + storage and transport properties arises primarily from local distortions of the structure and chemomechanical coupling at the phase interface, which may accelerate ion/electron diffusion, create additional active sites, and bolster structural stability. Finally, we offer our perspectives on the existing challenges, knowledge gaps, and opportunities for the advancement of heterostructure engineering as a means to develop practical, high-performance sodium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Morphology–Diffusivity Relationships in Fluorine-Free Random Terpolymers for Proton-Exchange Membranes

Here, using atomistic molecular dynamics simulations, we investigate the morphology and transport properties of a new family of fluorine-free terpolymers designed as proton-exchange membranes. Simulated random terpolymers consist of three monomers with a 5-carbon backbone with a phenylsulfonate, phenyl, or no pendant group and have ion exchange capacities (IECs) ranging from 1.06–4.14 mmol/g. At a hydration level of 9, cluster analysis reveals macrophase separation between water and terpolymers with IEC < 2.1 mmol/g and continuous, percolated hydrophilic and hydrophobic nanoscale domains at higher IECs. Channel width distribution analysis of the percolated morphologies revealed that more hydrophobic units produce less uniform channels. Decreasing the surface area per sulfonate group and increasing the fractal dimension of the hydrophilic domains correlate with increased water diffusivity, due to a more acidic interface and more isotropic water channels. Relative to the previously studied phenylsulfonate homopolymer, these terpolymers with lower IECs have only modestly lower water diffusion, and we anticipate other advantages related to processability.

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

Mechanistic origins of corrosion protection of aluminum alloys by graphene/polyetherimide nanocomposite coatings

Light-weighting vehicular components through adoption of light-metal structural alloys holds promise for reducing the fuel consumption of internal combustion engine vehicles and increasing the range of battery electric vehicles. However, the alloyed microstructure and surface precipitates of aluminum alloys render these materials susceptible to corrosion under modest excursions from neutral pH. Traditional chromium-based anodic passivation layers are subject to increasingly stringent environmental regulations, whereas options for sacrificial cathodic films are sparse for electropositive metals. While hybrid nanocomposite coatings have shown initial promise, mechanistic underpinnings remain poorly understood. Here, a fully imidized polyetherimide (PEI) resin is utilized as the continuous phase with inclusion of unfunctionalized exfoliated graphite (UFG). A comprehensive investigation of the mechanisms of corrosion protection reveals key fundamental design principles underpinning corrosion inhibition. First, strong interfacial adhesion, which for PEI is facilitated by binding of imide carbonyl moieties to Lewis acidic sites on Al surfaces. Second, the miscibility of ion-impervious nanoscopic UFG fillers and stabilization of a substantial interphase region at UFG/PEI boundaries that result in minimizing the free volume at the filler/polymer interface. Finally, extended tortuosity of ion diffusion pathways imbued by the below-percolation-threshold 2D fillers. These three design principles help govern and modulate ion transport from electrolyte/coating interfaces to the coating/metal interface and are crucial for the extended preservation of barrier properties. The results suggest an approach to systematically activate multiple modes of corrosion inhibition through rational design of hybrid nanocomposite coatings across hard-to-abate sectors where light metal alloys are likely to play an increasingly prominent role.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structural measurement of electron-phonon coupling and electronic thermal transport across a metal-semiconductor interface

Scattering of energetic charge carriers and their coupling to lattice vibrations (phonons) in dielectric materials and semiconductors are crucial processes that determine the functional limits of optoelectronics, photovoltaics, and photocatalysts. The strength of these energy exchanges is often described by the electron-phonon coupling coefficient, which is difficult to measure due to the microscopic time- and length-scales involved. In the present study, we propose an alternate means to quantify the coupling parameter along with thermal boundary resistance and electron conductivity by performing a high angular-resolution time-resolved X-ray diffraction measurement of propagating lattice deformation following laser excitation of a nanoscale, polycrystalline metal film on a semiconductor substrate. Our data present direct experimental evidence for identifying the ballistic and diffusive transport components occurring at the interface, where only the latter participates in thermal diffusion. This approach provides a robust measurement that can be applied to investigate microscopic energy transport in various solid-state materials.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Numerical Study of Mixing of Two Fluids Under Low Gravity

The mixing characteristics of two fluids inside a cavity due to buoyancy driven flow fields for low gravity conditions is investigated via numerical experiments. The buoyancy driven flow, depending on the parametric region, stretches and deforms the material interface into a wave morphological pattern. The morphological pattern affects the resulting stratification thickness of the mixed region. Three basic mixing regimes occur: convective, diffusive, and chaotic. In the convective regime, an overturning motion occurs which gives rise to a stable wave formation. This wave oscillates and its decay leads to a stable stratification. Whereas, in the diffusive regime, the length of the interface remains constant while mixing occurs. This limiting behavior is very important to materials processing in space, and it admits a closed form solution corresponding to vanishing convective terms which agrees with computational results. Finally, in the chaotic regime, the material interface continuously stretches and folds on itself similar to a horseshoe map. The length of stretch of the interface increases exponentially. Internal wavebreaking occurs for this case. This wavebreaking generates local turbulence, and provides an effective mechanism for mixing.

Duval, Walter M. B.↗
NASA NTRS1999

Surface Tension and Fingering of Miscible Interfaces

Experiments on miscible, buoyantly unstable reaction-diffusion fronts and non-reacting displacement fronts in Hele-Shaw cells show a fingering-type instability whose wavelengths (lambda*) are consistent with an interfacial tension (sigma) at the front caused by the change in chemical composition, even though the solutions are miscible in all proportions. In conjunction with the Saffman-Taylor model, the relation sigma = K/tau, where tau is the interface thickness and K approximately equal 4 +/- 2 x 10(exp -6) dyne, enables prediction of our measured values of lambda* as well as results from prior experiments on miscible interfaces. These results indicate that even for miscible fluids, surface tension is generally a more significant factor than diffusion in interfacial stability and flow characteristics.

Abib, Mohammed↗
DOE OSTI2025

Micro-tensile testing of neutron-irradiated Al/Zr and Zr/U-Mo diffusion bonds

This study focuses on micro-tensile testing of neutron-irradiated Al/Zr and Zr/U-Mo diffusion bonds, which are integral to the structure and performance of nuclear fuel plates used in the fuel being developed for U.S. high-power research reactors. During fabrication, two distinct diffusion bonds occur—one at the Zr/U-10Mo interface via roll bonding, and the Zr/Al interface via hot isostatic pressing at 833 K. This work utilizes micro-tensile testing specimens approximately 7 × 7 × 18 μm 3 , which are sufficiently large to encompass the identified diffusion zones to evaluate the failure strength and identify the failure location. For the Zr/Al interface, all failures occurred in the bulk aluminum. The average yield strength, ultimate strength, and strain at failure for the two parent samples were: 125 ± 10 MPa, 139 ± 24 MPa, 156 ± 21 MPa, 185 ± 14 MPa, 37.9 ± 5.6 %, and 29.7 ± 3.4 %, respectively. For specimens across the Zr/U-Mo interface, failures occurred either in the bulk Zr or U-Mo, with no failures observed in the diffusion region. For specimens that failed in the bulk Zr the failure strength and strain for the two parent samples were 570 MPa, 710 ± 131 MPa, and 43.3 %, 27.3 ± 6.6 % respectively. The specimens that failed in the bulk U-Mo had failure strengths and strains of 596 ± 32 MPa, 548 ± 124 MPa, and 2.2 ± 1.3 %, 1.5 ± 0.6 % respectively. In conclusion, these findings support the development of accurate thermo-mechanical models for irradiated fuel performance by identifying the mechanical limits and failure locations within the bonded regions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Phosphorus diffusion and deactivation during SiGe oxidation

Dopant profiles near the semiconductor–oxide interface are critical for microelectronic device performance. As the incorporation of Si 1−x Ge x into transistors continues to increase, it is necessary to understand the behavior of dopants in Si 1−x Ge x . In this paper, the diffusion and electrical activation of phosphorus within a strained, single-crystal Si 0.7 Ge 0.3 layer on Si during oxidation are reported. Both layers were uniformly doped, in situ, with an average phosphorus concentration of 4 × 10 19 atoms/cm 3 . After high-temperature oxidation, secondary ion mass spectrometry measurements revealed that the bulk of the phosphorus diffuses out of only the SiGe layer and segregates at the oxidizing SiGe–SiO 2 interface. Hall effect measurements corroborate the observed phosphorus loss and show that the phosphorus diffusing to the oxidizing interface is electrically inactive. Through density functional theory (DFT) calculations, it is shown that phosphorus interstitials prefer sites near the SiGe–SiO 2 interface. Finally, based on a combination of experimental data and DFT calculations, we propose that the phosphorus atoms are displaced from their lattice sites by Ge interstitials that are generated during SiGe oxidation. The phosphorus atoms then migrate toward the SiGe–SiO 2 interface through a novel mechanism of hopping between Ge sites as P–Ge split interstitials. Once they reach the interface, they are electrically inactive, potentially in the form of interstitial clusters or as part of the reconstructed interface or oxide.

Thornton, Chappel S.↗
NASA NTRS2005

Structures Self-Assembled Through Directional Solidification

Nanotechnology has created a demand for new fabrication methods with an emphasis on simple, low-cost techniques. Directional solidification of eutectics (DSE) is an unconventional approach in comparison to low-temperature biomimetic approaches. A technical challenge for DSE is producing microstructural architectures on the nanometer scale. In both processes, the driving force is the minimization of Gibb's free energy. Selfassembly by biomimetic approaches depends on weak interaction forces between organic molecules to define the architectural structure. The architectural structure for solidification depends on strong chemical bonding between atoms. Constituents partition into atomic-level arrangements at the liquid-solid interface to form polyphase structures, and this atomic-level arrangement at the liquid-solid interface is controlled by atomic diffusion and total undercooling due to composition (diffusion), kinetics, and curvature of the boundary phases. Judicious selection of the materials system and control of the total undercooling are the keys to producing structures on the nanometer scale. The silicon-titanium silicide (Si-TiSi2) eutectic forms a rod structure under isothermal cooling conditions. At the NASA Glenn Research Center, directional solidification was employed along with a thermal gradient to promote uniform rods oriented with the thermal gradient. The preceding photomicrograph shows the typical transverse microstructure of a solidified Si-TiSi2 eutectic composition. The dark and light gray regions are Si and TiSi2, respectively. Preferred rod orientation along the thermal gradient was poor. The ordered TiSi2 rods have a narrow distribution in diameter of 2 to 3 m, as shown. The rod diameter showed a weak dependence on process conditions. Anisotropic etch behavior between different phases provides the opportunity to fabricate structures with high aspect ratios. The photomicrographs show the resulting microstructure after a wet chemical etch and a dry plasma etch. The wet chemical etches the silicon away, exposing the TiSi2 rods, whereas plasma etching preferentially etches the Si-TiSi2 interface to form a crater. The porous architectures are applicable to fabricating microdevices or creating templates for part fabrication. The porous rod structure can serve as a platform for fabricating microplasma devices for propulsion or microheat exchangers and for fabricating microfilters for miniatured chemical reactors. Although more work is required, self-assembly from DSE can have a role in microdevice fabrication.

Dynys, Frederick W.↗
DOE OSTI2023

Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale (Final Report)

This project titled “Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale” represents a collaborative effort between the University of Michigan (Jason Demers, PI) and Oak Ridge National Laboratory (Scott Brooks, co-I). Much has been learned about mercury (Hg) cycling in stream ecosystems, and East Fork Poplar Creek (EFPC) in particular, through decades of previous research. Nevertheless, some of the most fundamental questions regarding the sources of bioavailable Hg and its transformation to toxic methylmercury (MeHg) have remained unanswered. These fundamental questions include: (1) what are the sources and biogeochemical processes that lead to the input of dissolved Hg to stream water across critical subsurface interfaces within stream ecosystems, and EFPC in particular? and (2) what are the sources and biogeochemical processes that control the production and fate of bioaccumulative MeHg within stream ecosystems, and in EFPC in particular? To address these fundamental questions, our project aimed to couple laboratory experiments and field observations, both utilizing natural abundance Hg stable isotope techniques, to identify the processes responsible for generating mobile, bioavailable dissolved Hg from recalcitrant legacy sources within critical subsurface zones (e.g., streambed hyporheic zone, riparian floodplain subsurface). We used the isotopic signature of this bioavailable dissolved Hg to track its mobilization across these critical interfaces in order to link diffuse subsurface sources of dissolved Hg with increases in surface water dissolved Hg flux measured at the watershed scale. Additionally, our research aimed to determine the isotopic composition of MeHg within these same critical subsurface zones. We directly assessed the isotopic composition of MeHg within biota in order to gain insight into which subsurface sources of inorganic Hg and toxic MeHg are available for bioaccumulation within the EFPC ecosystem. Net fluxes of dissolved Hg along the flow path of EFPC were shown to vary spatially and temporally. In the Upper EFPC, within the Y12 boundary, stream water flux of dissolved Hg consistently decreased between the outfall and the downstream boundary of Y12 (57% ± 29%, 1SD). Within the Upper EFPC, an assessment of Hg isotopic composition suggested that losses were strongly reaction-driven, although isotopic diagnostics did not conform to any known processes. Downstream of Y12, in the upper reach of the Lower EFPC, dissolved Hg fluxes tended to increase during the dormant season (net gain of 11-120%), and decrease during the growing season (net loss of 23% +/- 18%, 1SD). In the downstream-most reach of Lower EFPC, dissolved Hg fluxes increased by 12-108% in 9 out of 10 monthly assessments. Overall, diffuse fluxes from the non-Y12 watershed accounted for 34% (+/- 17%, 1SD) of all dissolved Hg exported during base flow. Within Lower EFPC, an assessment of Hg isotopic composition was consistent with the contribution of diffuse Hg inputs from high-concentration hotspots within riparian floodplains and streambed hyporheic pore water. To investigate remobilization of recalcitrant Hg from legacy sediment sources, we developed procedures that coupled isotopic analysis with sequential extractions of streambed sediment. We found that the proportion of weakly-bound Hg within EFPC streambed sediment was relatively small, but could still account for a large proportion of the annual flux of dissolved Hg from EFPC. These sequential extractions also showed that this weakly-bound Hg fraction could be replenished from the much larger fraction of recalcitrant Hg in sediment. The isotopic composition of these weakly-bound and remobilized recalcitrant Hg fractions within the sediment was consistent with high-concentration dissolved Hg hotspots within hyporheic pore water. Thus, this research provided novel evidence that legacy mercury sources within streambed sediment could provide an ongoing contribution of dissolved Hg to surface waters. Finally, we developed new methods for the direct determination of the MeHg isotopic composition of organisms, which allowed a more direct evaluation of inorganic Hg and MeHg sources accumulating in the food web. We found that fish and aquatic invertebrates in both EFPC and a regional background site obtained inorganic Hg and MeHg from multiple isotopically distinct sources, including sediment, suspended particulates, and periphyton. Photodemethylation was found to be an important reaction influencing MeHg dynamics at both sites. However, the balance of microbial methylation and demethylation processes differed between the two streams, with fractionation resulting from methylation and demethylation processes being relatively in balance within the regional background site, whereas microbial methylation appeared to be dominant over microbial demethylation within the EFPC ecosystem. Broadly, the application of Hg isotopic analysis in this study led to numerous novel insights regarding the biogeochemical cycling of Hg in stream ecosystems. This research project promoted the development of two new approaches, including the coupling of sequential extractions with Hg isotopic analysis to assess remobilization of recalcitrant Hg within sediments, and a new method for the isotopic analysis of MeHg isolated from environmental samples. Both of these efforts represent advances in capacity for the field of mercury isotopic analysis and environmental assessment. Overall, this study demonstrates that the application of Hg stable isotope techniques continues to provide new insights into the biogeochemical cycling of Hg in complex aquatic environments, both within the EFPC and beyond.

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

3D printed transparent ceramic YAG laser rods: Matching the core-clad refractive index

Yttrium Aluminum Garnet (YAG) solid state laser gain media rods with an active Neodymium-doped core and an optically-clear cladding region were additively manufactured via direct-ink-writing (DIW), followed by sintering and hot isostatic pressing to form fully dense optical ceramics. Lutetium and Gadolinium were chosen as optically-inert co-doping ions in the clad to match the increase in refractive index caused by the Neodymium in the core. Additionally, either 11.6% Lutetium or 3.8% Gadolinium can be used to match the index change from 2% neodymium; however, differences in diffusion distances across the core-clad interface lead to large fluctuations in index in that region. These index fluctuations can be minimized either by matching dopants with similar diffusion distances, or by implementing a gradual gradient in the doping profile, possible through DIW, rather than a sharp compositional interface attainable via more standard fabrication methods. This improvement in index homogeneity resulted in a 40% improvement in lasing performance compared with that of a core-clad rod fabricated with a sharp interface between the doped and undoped regions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational Methodology to Simulate Pyrometallurgical Processes in a Secondary Lead Furnace

Pyrometallurgical recovery of nonferrous metals involves a combination of thermally intensive transformations during exothermic gas-phase reactions, endothermic decomposition of solid charge, and melting of simpler solids. In the recovery of secondary lead, simultaneous thermal effects in a reverberatory-style furnace cause a melt pool to accumulate at the bottom, with lighter solids (slag) floating above and gaseous products from decomposition of the charge diffusing through the gas–slag interface. Species from oxy-fuel combustion of natural gas, species profiles from smelting reactions, and the formation of a melt pool consisting primarily of lead are simulated via a time-averaged formulation. Predictions of outflow are compared with preset inflow profiles to ensure conservation of mass. Thermal profiles for solid, liquid, and gas phases are presented by species. A novel method is implemented to model the latent heat of fusion using a heterogeneous chemical reaction. The simulation is conducted in Simcenter STAR-CCM+ v. 16.02.009-R8.

Rao, Vivek↗
DOE OSTI2023

A comparative study of solute trapping in Fe-(33–45 at%) Cu alloys manufactured by laser directed energy deposition

The high cooling rates in laser directed energy deposition (DED-LB) of alloys lead to substantial amounts of solute trapping as solute atoms cannot diffuse away from the solid/liquid interface before it advances. In some concentrated alloys, this results in supersaturated solid phases that form nanoscale hierarchical microstructures when the solute atoms precipitate out during reheating from subsequent laser passes. We choose the Iron-Copper (Fe-Cu) binary alloy as model system as it is chemically homogeneous in the liquid phase and has negligible solid solubility at room temperature. Two alloys with nominal compositions in atomic (at.) %, Fe67Cu33 and Fe55Cu45, were manufactured using DED-LB. Scanning transmission electron microscopy (STEM), energy dispersive spectroscopy (EDS) and wavelength dispersive spectroscopy (WDS) were used to characterize the nanostructures and heterogeneous chemical compositions. A non-equilibrium solute partitioning model was used to compute the supersaturated chemistries of the constituent phases and validated with experimentally measured compositions. The measured phase compositions of the two alloys were very similar, at roughly 12 at.% Cu and 4 at.% Fe in the α(bcc)-Fe and ε(fcc)-Cu phases respectively, despite having different processing parameters and mechanical behavior. This indicates that the total thermal history, that depends on both the processing parameters and laser scan pattern, plays a stronger role on the final microstructure evolution than just the initial quantity of trapped solute. In addition, we find that current non-equilibrium solute partitioning models applied on the continuum scale fall short of predicting accurate quantitative phase compositions in concentrated alloys, although the qualitative trends are captured correctly.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Understanding formation mechanisms of intermetallic compounds in dissimilar Al/steel joint processed by resistance spot welding

Here, this paper confirmed the formation mechanism of intermetallic compounds (IMCs) in Al/steel resistance spot welds with transmission electron microscopy, electron backscatter diffraction, nanoindentation and thermodynamic calculations. In particular, the formation of AlFe with BCC_B2 structure, which is not commonly seen in welds, was identified. The formation mechanism of IMCs at the high welding energy region is described as follows. Al 13 Fe 4 first nucleates from the Al side, followed by Al 5 Fe 2 growth with Fe atoms accumulating at the Al 13 Fe 4 grain boundaries. Then, Al 5 Fe 2 grains grow continuously to coarse columnar grains, and AlFe forms at the interface of Al 5 Fe 2 grains and the ferrite phase. Lastly, the needle-like Al 13 Fe 4 forms in the cooling process. At the middle welding energy region, only equiaxed Al 5 Fe 2 and small Al 13 Fe 4 grains are formed at the interface because of the lower diffusion rates of Al and Fe, hence postponing the growth of Al 5 Fe 2 and Al 13 Fe 4 . At the low welding energy region, only sporadic Al 5 Fe 2 and Al 13 Fe 4 grains are formed surrounded by the Al phase.

42 ENGINEERING↗
DOE OSTI2024

Partial Solvation of Lithium Ions Enhances Conductivity in a Nanophase-Separated Polymer Electrolyte

In this study, we demonstrate that a multiblock lithium-ion-conducting polymer can be swollen with ethylene carbonate solvent to increase the conductivity relative to the dry polymer material by nearly 4 orders of magnitude. This increase is due to the partial solvation of lithium ions by ethylene carbonate, which leads to Li + diffusion along the solvent–polymer interface. This differs from the vehicular transport mechanism for lithium ions in pure solvent. We use a combination of broadband dielectric spectroscopy, X-ray scattering, and all-atom molecular dynamics simulations to probe the effect of the solvent on the polymer morphology and to elucidate the mechanism of lithium ion transport.

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

Grain boundary effects in high-temperature liquid-metal dealloying: a multi-phase field study

Abstract A multi-phase field model is employed to study the microstructural evolution of an alloy undergoing liquid dealloying, specifically considering the role of grain boundaries. A semi-implicit time-stepping algorithm using spectral methods is implemented, which enables simulating large 2D and 3D domains over long time scales while still maintaining a realistic interfacial thickness. Simulations reveal a mechanism of coupled grain–boundary migration to maintain equilibrium contact angles with the topologically complex solid–liquid interface, which locally accelerates diffusion-coupled growth of a liquid channel into the precursor. This mechanism asymmetrically disrupts the ligament connectivity of the dealloyed structure in qualitative agreement with published experimental observations. The grain boundary migration-assisted corrosion channels form even for precursors with small amounts of the dissolving alloy species, below the parting limit . The activation of this grain boundary dealloying mechanism depends strongly on grain boundary mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multilayer interface tracking model of pure tungsten oxidation

Here, we present a numerical model to predict oxide scale growth on tungsten surfaces under exposure to oxygen at high temperatures. The model captures the formation of four thermodynamically-compatible oxide sublayers, WO 2 , WO 2.72 , WO 2.9 , and WO 3 , on top of the metal substrate. Oxide layer growth is simulated by tracking the oxide/oxide and oxide/metal interfaces using a sharp-interface Stefan model coupled to diffusion kinetics. The model is parameterized using selected experimental measurements and electronic structure calculations of the diffusivities of all the oxide subphases involved. We simulate oxide growth at temperatures of 600°C and above, extracting the power law growth exponents in each case, which we find to deviate from classical parabolic growth in several cases. We conduct a comparison of the model predictions with an extensive experimental data set, with reasonable agreement at most temperatures. While many gaps in our understanding still exist, this work is a first attempt at embedding the thermodynamic and kinetic complexity of tungsten oxide growth into a comprehensive mesoscale kinetic model that attempts to capture the essential features of tungsten oxidation to fill existing knowledge gaps and guide and enhance future tungsten oxidation models.

36 MATERIALS SCIENCE↗