Diffusion-induced stress concentrations in diffusional creep

The lattice continuum formulation for diffusional creep is implemented into the computational framework based on the finite difference method for space and the predictor–corrector algorithm for time discretization, to solve the coupled elasticity-diffusion problem with moving boundaries. The numerical scheme is implemented and tested considering 3D periodic structures without grain boundary sliding. Of primary interest in this paper are the stress non-uniformities resulting from nonuniform composition eigenstrains.