An efficient numerical method for evolving microstructures with strong elastic inhomogeneity

In this paper, we consider a fast and efficient numerical method for the modified Cahn-Hilliard equation with a logarithmic free energy for microstructure evolution. Even though it is physically more appropriate to use a logarithmic free energy, a quartic polynomial approximation is typically used for the logarithmic function due to a logarithmic singularity. In order to overcome the singularity problem, we regularize the logarithmic function and then apply an unconditionally stable scheme to the Cahn-Hilliard part in the model. We present computational results highlighting the different dynamic aspects from two different bulk free energy forms. We also demonstrate the robustness of the regularization of the logarithmic free energy, which implies the time-step restriction is based on accuracy and not stability.