Kinetic modeling of high temperature oxidation of Ni-base alloys

The rates of isothermal and cyclic oxidation and the elemental concentration profiles as a function of time of oxidation for a few Ni-base superalloys were determined through a modified Wagner’s oxidation model and the solution of coupled elemental diffusion equations. Thermodynamically calculated interfacial elemental concentrations and oxygen partial pressures for the multi-component Ni-base alloys were used as boundary conditions for the solution of Wagner’s equation and the elemental coupled diffusion equations (for Cr, Al and O). The multiple elemental diffusion and mass conservation equations were solved using a numerical procedure. The dependence of self/tracer-diffusivities of Cr, Al and O in the corundum phase on the oxygen partial pressures was deduced using a genetic algorithm based optimization procedure incorporating the experimental parabolic rate constants for several Ni-base alloys. Rates of cyclic oxidation were then deduced from the deterministic interfacial cyclic oxidation spalling model (DICOSM) developed by Smialek [1]. The calculated oxidation rates were in reasonable agreement with the experimental values for a range of multi-component Ni-base alloys.