Search NASA⌕ Search

SEARCH · Search NASA

Results for “interface & surface thermodynamics”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Kinetic versus Thermodynamic Stability of LLZO in Contact with Lithium Metal

Li 7 La 3 Zr 2 O 12 (LLZO) garnet-based oxides are a promising class of solid electrolytes used as the separator in all-solid-state batteries (ASSBs). While LLZO is considered to have a wide electrochemical stability window, its intrinsic stability in contact with lithium metal is not sufficiently well understood, and there is still a debate on the key question of whether LLZO does or does not form passivation layers before and during cycling. Utilizing both in situ and operando X-ray photoelectron spectroscopy techniques, in this work we reveal the presence of a kinetic barrier to the reduction of LLZO by Li metal, with the extent of oxygen-deficient interphase (ODI) formation depending sensitively on the energetics of Li metal arriving at the Li vertical bar LLZO interface. Despite the clear presence of a kinetic barrier to reduction, the electrochemical response of the Li vertical bar LLZO interface is unchanged by the presence of the ODI, indicating that ODI formation during electrochemical cycling does not hinder charge transfer across the Li vertical bar LLZO interface. Overall, these results reveal that the reactivity of LLZO with Li metal depends not only on the material properties of the adjoining phases (i.e., surface purity and active contact) and their resulting thermodynamic stability but also on the energy input at the interface and the resulting reaction kinetics. Furthermore, the presence of a kinetic barrier to reduction highlights the additional complexities governing the reactivity of solid-state interfaces in ASSBs and underscores the importance of operando characterization of interfacial stability to design more robust, high-performance protection strategies for solid electrolytes in contact with reactive electrodes.

25 ENERGY STORAGE↗
DOE OSTI2026

Remote sensing of planetary boundary layer from ground and space: Structure, thermodynamics, and boundary layer clouds

The planetary boundary layer (PBL) is a critical interface between Earth's surface and atmosphere, influencing atmospheric convection, weather patterns, and air quality. Recognized by the 2017 National Academies of Sciences, Engineering, and Medicine Earth Science Decadal Survey as an Incubation Targeted Observable, high-quality and effective observation of the PBL has become a priority. This paper reviews surface-based and satellite remote sensing techniques for characterizing PBL features, including PBL height (PBLH), boundary layer thermodynamics, turbulence, and boundary layer clouds (BLCs). These elements are treated as interconnected aspects of the PBL system, while recognizing that different instruments retrieve different physical manifestations of the PBL. The review summarizes recent advances and limitations in Micro-Pulse Lidar (MPL), Doppler lidar, Raman lidar, Differential Absorption Lidar (DIAL), ceilometers, wind profilers, GNSS Radio Occultation, radar, and hyperspectral sounders. These observations help trace dynamic processes within the PBL and link PBL structure to broader weather and climate processes. The review also highlights persistent observational gaps over oceans, remote land, and polar regions, where continuous surface-based profiling is sparse. By comparing instrument capabilities, practical limitations, and interpretation issues, this review emphasizes the need for integrated remote sensing approaches and careful definition of the retrieved PBL quantities.

Boundary layer clouds↗
DOE OSTI2020

Pulsed power accelerator surface Joule heating models

Understanding the effects of contaminant plasmas generated within the Z machine at Sandia is critical to understanding current loss mechanisms. The plasmas are generated at the accelerator electrode surfaces and include desorbed species found in the surface and substrate of the walls. These desorbed species can become ionized. The timing and location of contaminant species desorbed from the wall surface depend non-linearly on the local surface temperature. For accurate modeling, it is necessary to utilize wall heating models to estimate the amount and timing of material desorption. One of these heating mechanisms is Joule heating. In this paper, we propose several extended semi-analytic magnetic diffusion heating models for computing surface Joule heating and demonstrate their effects for several representative current histories. We quantitatively assess under what circumstances these extensions to classical formulas may provide a validatable improvement to the understanding of contaminant desorption timing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Characterizing Structure-Dependent TiS 2 /Water Interfaces Using Deep-Neural-Network-Assisted Molecular Dynamics

As a promising layered electrode material, TiS 2 -based capacitive deionization (CDI) devices for water desalination have attracted significant attention. However, TiS 2 /H 2 O interfacial features, potentially important for device optimization, remain unidentified. Using Deep Potential Molecular Dynamics (DPMD), we characterized distinct aqueous interfaces introduced by four TiS 2 terminations expected to be present as water intercalates into TiS 2 , namely, Armchair, Zigzag, Zigzag-L, and Zigzag-R. First, we assessed important representative physical properties of the system to validate the deep potentials (DPs). DPMD simulations agree well with experiments and first-principles simulations, suggesting the DPs are accurate and reliable. Subsequent simulations of these TiS 2 /water interfaces revealed how TiS 2 surface termination influences the structure of interfacial water. This effect is most evident in the first and second water layers close to the TiS 2 surface, and more pronounced when spontaneous dissociative adsorption of water occurs. The extent of water dissociation on each surface was evaluated using enhanced sampling. Zigzag-L is the only interface where proton transfer from adsorbed water to TiS 2 surface S atoms is thermodynamically and kinetically favored. The coexistence of surface four-fold-coordinated Ti (Ti 4c ) and one-fold-coordinated S (S 1c ) is found to be essential to making proton transfer feasible on the Zigzag-L surface. Furthermore, remaining unprotonated S 1c atoms can act as good proton acceptors after water dissociation. Thus, TiS 2 with Zigzag-L termination may be a surface to avoid in CDI device construction, given that pH fluctuations adversely affect performance. Furthermore, this work provides new understanding of TiS 2 /H 2 O interfacial features that could aid future design and optimization of TiS 2 -based CDI devices for water desalination.

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

Excitonic Effects in X-ray Absorption Spectra of Fluoride Salts and Their Surfaces

Given their natural abundance and thermodynamic stability, fluoride salts may appear as evolving components of electrochemical interfaces in Li-ion batteries and emergent multivalent ion cells. This is due to the practice of employing electrolytes with fluorine-containing species (salt, solvent, or additives) that electrochemically decompose and deposit on the electrodes. Operando X-ray absorption spectroscopy (XAS) can probe the electrode–electrolyte interface with a single-digit nanometer depth resolution and offers a wealth of insights into the evolution and Coulombic efficiency or degradation of prototype cells, provided that the spectra can be reliably interpreted in terms of local oxidation state, atomic coordination, and electronic structure about the excited atoms. Here we explore fluorine K-edge XAS of mono- (Li, Na, and K) and di-valent (Mg, Ca, and Zn) fluoride salts from a theoretical standpoint and discover a surprising level of detailed electronic structure information about these materials despite the relatively predictable oxidation state and ionicity of the fluoride anion and the metal cation. Utilizing a recently developed many-body approach based on the ΔSCF method, we calculate the XAS using density functional theory and experimental spectral profiles are well reproduced despite some experimental discrepancies in energy alignment within the literature, which we can correct for in our simulations. We outline a general methodology to explain shifts in the main XAS peak energies in terms of a simple exciton model and explain line-shape differences resulting from the mixing of core-excited states with metal d character (for K and Ca specifically). Given ultimate applications to evolving interfaces, some understanding of the role of surfaces and their terminations in defining new spectral features is provided to indicate the sensitivity of such measurements to changes in interfacial chemistry.

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

Thermal interfaces in dynamic compression experiments

Temperature is a difficult thermodynamic variable to measure in dynamic compression experiments. Optical pyrometry is a general-purpose technique for measuring temperature from a radiant surface, but that surface is often the interface between distinct materials with temperatures that vary spatially along the loading direction. This leads to a fundamental problem: how does the measured interface temperature relate to this temperature profile along the compression axis? Numerical analysis of loading history and thermal diffusion at these interfaces shows that seemingly subtle changes in experiment geometry can lead to very different temperature profiles. We compare these results to laboratory temperature measurements of shock-compressed tin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

In-situ thermodynamics measurements at metal oxides-solution interfaces using Flow Adsorption Microcalorimetry.

Mineral-fluid interfaces are the principal sites of geochemical processes near Earth’s surface, hosting chemical reactions that play a fundamental role in (bio)geochemical cycles and in the fate and transport of anthropogenic contaminants, and that effectively control the compositions of soil and water environments. These complex interfaces are critical for our energy and environmental future. The mineral-fluid interface has been studied in an unprecedented level of detail with both experimental and computational approaches, separately and in combination. However, conspicuously missing from studies of the mineral-fluid interface are direct measurements of the energies of ion sorption and exchange. The literature on energetics and enthalpies of exchange, adsorption, dissolution, precipitation, and surface protonation reactions, especially those directly supported by experimental data, remains scarce despite their fundamental nature. The overarching goal of this project is to complete a systematic study of the thermodynamics properties of interfacial reactions at four MO surfaces (Rutile (α-TiO2), Quartz (SiO2), boehmite (γ-AlOOH) and goethite (α-FeOOH)) through the application and construction of novel flow adsorption microcalorimetry techniques and instrumentations. These unique and specialized microcalorimeters will operate at various temperatures and solution chemical compositions allowing for in-situ measurements across metal oxides and ligands of various characteristics. The overall research goal will be accomplished by completing the following three specific objectives (O): O1) Determine the energetics of surface protonation and deprotonation, ion exchange and ligand sorption reactions; O2) Investigate the surface charge thermodynamic properties under a range of temperature and solution chemical compositions; and O3) Develop predictive trends of the interplay between MO structure, surface coverage and surface reactivity. In addition to key thermodynamics parameters, calorimetric measurements provide a wealth of mechanistic information about reactions energetics and kinetics, surface charge characteristics, and structure-reactivity or selectivity relationships, all obtained in-situ and in real-time. This report includes science highlights from various projects completed over the performance period of the project. Also listed are dissemination opportunities, people supported on the grant and the impact on available physical resources and the discipline as a whole.

58 GEOSCIENCES↗
DOE OSTI2025

Multicationic interactions mitigating lattice strain in sodium layered cathodes

Transition-metal (TM) layered oxides have emerged as the primary cathode choice for sodium-ion batteries (SIBs) due to their high energy density and sustainable chemistry using non-critical elements. However, their anisotropic lattice strain and stress accumulation during (de)sodiation lead to severe structural degradation, yet an intrinsic strain-depressant approach remains elusive. Herein, we propose entropy regulation with zero Li/Co usage to mitigate harmful lattice displacements and enhance the electrochemical performance of sodium layered cathodes. Our findings demonstrate that high entropy design effectively inhibits TMO 6 octahedra distortions upon cycling, as evidenced by hard X-ray absorption spectroscopy, greatly reducing near-surface structural deconstruction and interface side reactions. Furthermore, multicationic interactions driven by configurational entropy thermodynamically mitigate the formation of oxygen defects and strengthen ligand-to-metal coordination. The complementarity inherent in charge compensation within complex systems is unveiled and the restrained lattice parameters deviations without interior volume residuals are successfully achieved. As a result, the multicationic cathode exhibits improved cycling stability and Na + diffusion kinetics in both half and full cells. The cathode chemistries outlined here broaden the prospects for lattice engineering to alleviate bulk fatigue and open up the possibility to develop an economically viable layered oxides with long durability.

Wang, Haoji [Central South University, Changsha (C↗
DOE OSTI2025

Facet-dependent structure and dissociation of water at pristine IrO 2 /water interfaces

Understanding the microscopic structure of water at metal oxide interfaces is crucial for advancing electrocatalysis. IrO 2 , specifically, has shown exceptional activity for electrochemical water oxidation, but we currently lack a fundamental understanding of how the surface structure of IrO 2 impacts water reactivity. In this work, we developed a machine learning potential trained to first-principles accuracy for modeling IrO 2 /water interfaces across different facets: (110), (100), (101), and (001). Using extensive machine learning molecular dynamics simulations, we investigated the spontaneous dissociation of water molecules at these interfaces. Our results reveal a distinct dissociation probability trend: (110) > (100) ≈ (101) > (001), which we attribute primarily to the reaction thermodynamics of surface water dissociation. A strong correlation is observed between the surface Ir–O bond distances and the dissociation probabilities, highlighting the role of surface geometry in modulating reactivity. As a consequence, the interfacial solvation structures and hydrogen bonding environments are dynamically tuned by the varying water dissociation capabilities across facets. This work elucidates how water dissociation energetics depend on surface orientation and interfacial structure, offering atomistic insights into manipulating reaction chemistry at electrocatalytic interfaces.

organic↗
DOE OSTI2020

Weak scaling of the contact distance between two fluctuating interfaces with system size

A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ* depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ* scales weakly with system size, i.e., the dependence of ℓ* on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim L→∞ (ℓ*/L) = 0. In particular, for (2+1)-dimensional interfaces ℓ* is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. Overall, this weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Thermodynamic Control of Interface Directs MnO 2 Nucleation Chemistry for Dense and Conformal Electrodeposition

Manganese dioxide (MnO 2 ) is widely recognized as a promising material for high-energy-density energy storage systems due to its broad applicability and facile electrodeposition. However, achieving uniform, thin, and high-mass-loading MnO 2 coatings on high-surface-area electrodes remains a significant challenge. Conventional electrodeposition methods typically yield nonuniform, thick layers with poor conductivity and limited material utilization, restricting their practical use. Here, we uncover a thermodynamically engineered vanadyl/pervanadyl (VO 2+ /VO 2 + ) interface that fundamentally reshapes MnO 2 electrodeposition chemistry, enabling highly uniform and dense coatings. Here, combining in situ AFM measurement, Classical Nucleation Theory, and Johnson–Mehl–Avrami–Kolmogorov modeling, we show that this interface reduces early-stage detectable MnO 2 island size by 35-fold and shifts the MnO 2 growth from diffusion-limited to reaction-limited progressive nucleation. This thermodynamically controlled interface yields highly dense and conformal MnO 2 films with record-high mass loading of 241 mg cm –2 (1607 mg cm –3 ) on 3D-printed graphene aerogels, without compromising porosity or inducing thickness gradient. As a prototype demonstration, the resulting MnO 2 electrodes deliver record-setting volumetric performance in both capacitors (106 F cm –3 ) and Zn//MnO 2 pouch cells (162 mAh cm –3 ). Beyond energy storage, our findings demonstrate the significance of thermodynamic interface control in MnO 2 nucleation chemistry for achieving dense and uniform coatings on various substrates, with implications for electrocatalysis, semiconductor processing, and advanced materials manufacturing.

Batteries↗
DOE OSTI2026

Amphiphilic Baskets for Supramolecular Nanoarchitectures at Interfaces: Inverted Monolayer Formation on Water

Interfacial chemistry of molecular baskets remains poorly understood despite their promise for supramolecular applications of detection and sequestration of toxic molecules including those of illicit drugs, organophosphorus compounds, and anticancer agents. We present a fundamental investigation of the interfacial behavior of three amphiphilic supramolecular baskets (ASB 4, 8, and 12), having increasingly longer yet linear alkyl chains at the top of their bowl-shaped cavity. The studies were completed at the air−water interface to elucidate surface activity, interfacial stability, self-assembly, and monolayer organization that drive inverted monolayer formation, in which the molecular arms orient toward the aqueous phase in a configuration opposite to that typically observed for lipids. Herein, surface pressure−area isotherms of ASB 4, 8, 12, deposited on a water surface, were performed in tandem with nonequilibrium relaxation experiments to quantify surface activity, thermodynamic stability, and monolayer compressibility of the baskets’ monolayer assembly. Brewster angle microscopy enabled direct visualization of morphological evolution, aggregation, and packing at the interface. We show that systematic extension of the hydrocarbon arms, from four to 12 methylene groups, progressively modifies intermolecular packing, drives distinct two-dimensional aggregation pathways, and increases number densities at the air−water interface. Atomistic molecular dynamics simulations corroborate many of these experimentally observed trends and provide mechanistic detail on the cooperative roles of basket topology and interfacial concentration in regulating the hydration structure and dynamics within the cavities generated by surface-adsorbing baskets, consistent with observed variations in surface activity and packing. Our results establish how the topology of these unique supramolecules and their concentration govern interfacial organization and offer a rational framework for designing amphiphiles with predictable behavior at soft interfaces.

Hydrocarbons↗
DOE OSTI2020

Combined Cr and S poisoning behaviors of La 1-x Sr x MnO 3±δ and La 1-x Sr x Co 1-y FeyO 3-δ cathodes in solid oxide fuel cells

Although the individual effects of airborne Cr and S contaminants on SOFC cathode performance degradation have been extensively studied, the combined effects of Cr and S contaminants remain largely unexplored. Under the real SOFC operating condition where the Cr and S species coexist, their effects may compete, affecting the poisoning behavior. Here, our investigation reveals that the combined Cr and S poisoning behavior of LSCF remains different from those of individual Cr and S poisonings, while the combined poisoning mechanism of LSM is equivalent to the sum of those of individual Cr and S effects. For LSCF electrode, gaseous Cr species are deposited mainly at the LSCF/GDC interface by electrochemical reduction, rather than forming SrCrO 4 on LSCF surfaces (as in the case of Cr-only poisoning), indicating no reaction between Cr vapors and SrO on the LSCF surface. Thermodynamic analysis demonstrates that the SrO on LSCF surface absorbs SO 2 (g) and thereby loses the Cr-gettering effect, allowing Cr vapors to flow through LSCF and to reach the LSCF/GDC interface where the Cr deposition occurs. Unlike LSCF, LSM electrode shows cumulative effects of Cr and S, as Cr accumulation occurs at the triple-phase boundary and S absorption takes place at localized Sr-rich regions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Revealing the Origin and Nature of the Buried Metal‐Substrate Interface Layer in Ta/Sapphire Superconducting Films

Abstract Despite constituting a smaller fraction of the qubit's electromagnetic mode, surfaces and interfaces can exert significant influence as sources of high‐loss tangents, which brings forward the need to reveal properties of these extended defects and identify routes to their control. Here, we examine the structure and composition of the metal‐substrate interfacial layer that exists in Ta/sapphire‐based superconducting films. Synchrotron‐based X‐ray reflectivity measurements of Ta films, commonly used in these qubits, reveal an unexplored interface layer at the metal‐substrate interface. Scanning transmission electron microscopy and core‐level electron energy loss spectroscopy identified an intermixing layer (≈0.65 ± 0.05 nm) at the metal‐substrate interface containing Al, O, and Ta atoms. Density functional theory modeling reveals that the structure and properties of the Ta/sapphire heterojunctions are determined by the oxygen content on the sapphire surface prior to Ta deposition for two atomic terminations of sapphire. Using a multimodal approach, we gained deeper insights into the interface layer between the metal and substrate, which suggests that the orientation of deposited Ta films depend on the surface termination of sapphire. The observed elemental intermixing at the metal‐substrate interface influences the thermodynamic stability and electronic behavior of the film, which may also affect qubit performance.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling Pb(II) Adsorption on Mineral Surfaces: Bridging Density Functional Theory and Experiment with Thermodynamic Insights

Despite decades of work on aqueous lead (Pb) adsorption on a-Fe2O3 (hematite) and a-Al2O3 (alumina), gaps between measurements and modeling obscure molecular-level understanding. Achieving well-matched geometries between theory and experimental for mineral-water interfaces is a hurdle, as surface functional group type and distribution must be accounted for in determining mechanisms. Additionally, computational methods that can describe the substrate are often not appropriate to capture aqueous effects. Progress requires focusing on well-studied and relevant systems, such as key facets (001),(012), and (110) of hematite and alumina, and ubiquitous contaminants such as aqueous Pb. In the past, bulk-parametrized bond-valence principles were used to rationalize Pb(II) adsorption trends. These approaches can break down at surfaces, where flexible bonding environments and adsorption-induced surface relaxations play a critical role. Here, we adapt and apply a density functional theory (DFT) and thermodynamics framework, integrating DFT-calculated energies with experimental data and electrochemical principles, to predict Pb(II) adsorption. Our model results capture trends across the full set of surfaces and predict that inner-sphere Pb(II) sorption on (001) alumina varies from unfavorable to weakly favorable across a range of pH conditions. This aligns with experiment insights that Pb(II) interacts at that surface through outer-sphere interactions. Extending to Fe(II) adsorption, we demonstrate a coverage-dependent site preference, potentially explaining disorder in overlayers grown by the oxidative adsorption of Fe(II) on hematite (001).

lead contamination↗
DOE OSTI2022

Nonequilibrium statistical thermodynamics of multicomponent interfaces

Nonequilibrium interfacial thermodynamics has important implications for crucial biological, physical, and industrial-scale transport processes. Here, we discuss a theory of local equilibrium for multiphase multicomponent interfaces that builds upon the “sharp” interface concept first introduced by Gibbs, allowing for a description of nonequilibrium interfacial processes such as those arising in evaporation, condensation, adsorption, etc. By requiring that the thermodynamics be insensitive to the precise location of the dividing surface, one can identify conditions for local equilibrium and develop methods for measuring the values of intensive variables at the interface. We then use extensive, high-precision nonequilibrium molecular dynamics (NEMD) simulations to verify the theory and establish the validity of the local equilibrium hypothesis. In particular, we demonstrate that equilibrium equations of state are also valid out of equilibrium, and can be used to determine interfacial temperature and chemical potential(s) that are consistent with nonequilibrium generalizations of the Clapeyron and Gibbs adsorption equations. We also show, for example, that, far from equilibrium, temperature or chemical potential differences need not be uniform across an interface and may instead exhibit pronounced discontinuities. However, even in these circumstances, we demonstrate that the local equilibrium hypothesis and its implications remain valid. These results provide a thermodynamic foundation and computational tools for studying or revisiting a wide variety of interfacial transport phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Self-healing Growth of LaNiO 3 on a Mixed-Terminated Perovskite Surface

Developing atomic-scale synthesis control is a prerequisite for understanding and engineering the exotic physics inherent to transition-metal oxide heterostructures. Thus, far, however, the number of materials systems explored has been extremely limited, particularly with regard to the crystalline substrate, which is routinely SrTiO 3 . Here, we investigate the growth of a rare-earth nickelate-LaNiO 3 -on (LaAlO 3 )(Sr 2 AlTaO 6 ) (LSAT) (001) by oxide molecular beam epitaxy (MBE). Whereas the LSAT substrates are smooth, they do not exhibit the single surface termination usually assumed necessary for control over the interface structure. Performing both nonresonant and resonant anomalous in situ synchrotron surface X-ray scattering during MBE growth, we show that reproducible heterostructures can be achieved regardless of both the mixed surface termination and the layer-by-layer deposition sequence. The rearrangement of the layers occurs dynamically during growth, resulting in the fabrication of high-quality LaNiO 3 /LSAT heterostructures with a sharp and consistent interfacial structure. This is due to the thermodynamics of the deposition window as well as the nature of the chemical species at interfaces-here, the flexible charge state of nickel at the oxide surface. This has important implications regarding the use of a wider variety of substrates for fundamental studies on complex oxide synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗