Abstract
The classical thermoelastic diffusion coupling equations in equilibrium states are extended to the nonequilibrium status introducing the interface energy combining with phase field method to simulate a full coupling thermoelastic diffusion in nonisothermal process without phase transition. Thermodynamic and kinetic coefficients are chosen through the parametric analysis of the asymptotic phase interface to meet the requirements of sharp interface model. This method could evaluate the interaction influence among temperature, deformation, and concentration in the process of microstructural evolution. As an example, a four-sphere diffusion model under heating and applied pressure is utilized to simulate microstructural evolution.