Multiscale Modeling of Zr Corrosion: Predicting Microstructures of Oxide Films
Abstract
A first-principles-informed multiscale modeling approach is reported for oxide formation during corrosion of zirconium. The present simulations reproduce the columnar corrosion film microstructure and grain size distribution seen in micrographs of oxidized Zr surfaces. This structural feature emerges from a hexagonal ZrO layer adjacent to the α-Zr/ZrO2−x interface, which is dynamically stabilized by continuous heterogeneous grain nucleation and grain boundary mass diffusion. This ZrO layer is identified as a significant oxygen diffusion barrier and a potential contributor to sub-parabolic pretransition corrosion kinetics. The simulations also reproduce the undulatory nature of the metal/oxide interface, leading to regions with tensile stress components capable of nucleating cracks, despite the overall state of compressive hydrostatic stress within the oxide film. First principles density functional theory (DFT) computations of fundamental ZrOx properties were combined with molecular dynamics simulations using a machine-learned potential based on a DFT training set to provide extensive input to an advanced phase field model for oxide formation and its microstructure, as implemented in the new MedeA PhaseField module. The phase field simulations reported herein incorporated computed stoichiometry-dependent thermochemical data, anisotropic elastic properties, interface energies, bulk and grain boundary mass transport processes, multiphase heterogeneous grain nucleation, and strain energy-dependent phase stabilities.
Versions
Related
- Role of defect trapping in the enhanced corrosion resistance of Ni-Mo alloys with yttrium and yttrium oxide additions
- Data-Driven Modeling of Molten Chloride Corrosion Using Static Test Data
- From Airspace Degradation to Deconfliction: Predicting Safety-Critical Traffic States through the Evolution of Multiscale Complexity Metrics
- Multiscale modeling of moisture diffusion in polymer/clay nanocomposites
- Contagion-preserving compression for multiscale epidemic modeling
BibTeX
@misc{starkey2026multiscale,
title = {{Multiscale Modeling of Zr Corrosion: Predicting Microstructures of Oxide Films}},
author = {K. Starkey and L. Kahle and M. Christensen and V. Eyert and A. El-Azab and C. B. Geller and A. P. Shivprasad and E. Wimmer},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.28493},
eprint = {2610.28493},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.28493},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Multiscale Modeling of Zr Corrosion: Predicting Microstructures of Oxide Films
AU - K. Starkey
AU - L. Kahle
AU - M. Christensen
AU - V. Eyert
AU - A. El-Azab
AU - C. B. Geller
AU - A. P. Shivprasad
AU - E. Wimmer
PY - 2026
DA - 2026/10/20
PB - arXiv
DO - 10.48550/arxiv.2610.28493
UR - https://arxiv.org/abs/2610.28493
AB - A first-principles-informed multiscale modeling approach is reported for oxide formation during corrosion of zirconium. The present simulations reproduce the columnar corrosion film microstructure and grain size distribution seen in micrographs of oxidized Zr surfaces. This structural feature emerges from a hexagonal ZrO layer adjacent to the α-Zr/ZrO_2−x interface, which is dynamically stabilized by continuous heterogeneous grain nucleation and grain boundary mass diffusion. This ZrO layer is identified as a significant oxygen diffusion barrier and a potential contributor to sub-parabolic pretransition corrosion kinetics. The simulations also reproduce the undulatory nature of the metal/oxide interface, leading to regions with tensile stress components capable of nucleating cracks, despite the overall state of compressive hydrostatic stress within the oxide film. First principles density functional theory (DFT) computations of fundamental ZrO_x properties were combined with molecular dynamics simulations using a machine-learned potential based on a DFT training set to provide extensive input to an advanced phase field model for oxide formation and its microstructure, as implemented in the new MedeA PhaseField module. The phase field simulations reported herein incorporated computed stoichiometry-dependent thermochemical data, anisotropic elastic properties, interface energies, bulk and grain boundary mass transport processes, multiphase heterogeneous grain nucleation, and strain energy-dependent phase stabilities.
ER -