High-temperature oxidation is an important property to evaluate thermal protection materials. However, since oxidation is a complex process involving microstructure evolution, its quantitative analysis has always been a challenge. In this work, a phase field method （PFM） based on the thermodynamics theory is developed to simulate the oxidation behavior and oxidation induced growth stress. It involves microstructure evolution and solves the problem of quantitatively computational analysis for the oxidation behavior and growth stress. Employing this method, the diffusion process, oxidation performance, and stress evolution axe predicted for Fe-Cr-A1-Y alloys. The numerical results agree well with the experimental data. The linear relationship between the maximum growth stress and the environment oxygen concentration is found. PFM provides a powerful tool to investigate high-temperature oxidation in complex environments.
