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.

[1]  John Lowengrub,et al.  Microstructural Evolution in Inhomogeneous Elastic Media , 1997 .

[2]  Harald Garcke,et al.  The Cahn-Hilliard equation with elasticity-finite element approximation and qualitative studies , 2001 .

[3]  Qiang Du,et al.  An iterative-perturbation scheme for treating inhomogeneous elasticity in phase-field models , 2005 .

[4]  Jie Shen,et al.  Morphological evolution during phase separation and coarsening with strong inhomogeneous elasticity , 2001 .

[5]  Yunxian Liu,et al.  A class of stable spectral methods for the Cahn-Hilliard equation , 2009, J. Comput. Phys..

[6]  Dietmar Gross,et al.  The equilibrium shape of an elastically inhomogeneous inclusion , 1997 .

[7]  Long-Qing Chen,et al.  Computer simulation of structural transformations during precipitation of an ordered intermetallic phase , 1991 .

[8]  Shenyang Y. Hu,et al.  A phase-field model for evolving microstructures with strong elastic inhomogeneity , 2001 .

[9]  P. K. Chan.  Effect of concentration gradient on the thermal-induced phase separation phenomenon in polymer solutions , 2006 .

[10]  Jeffrey W. Bullard,et al.  Computational and mathematical models of microstructural evolution , 1998 .

[11]  Dietmar Gross,et al.  The effect of elastic inhomogeneity on equilibrium and stability of a two particle morphology , 1998 .

[12]  J. Warren,et al.  Controlling the accuracy of unconditionally stable algorithms in the Cahn-Hilliard equation. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  J. K. Lee.  Elastic Stress and Microstructural Evolution , 1998 .

[14]  Jong K. Lee,et al.  A study on coherency strain and precipitate morphologyvia a discrete atom method , 1996 .

[15]  Steven M. Wise,et al.  EFFICIENT PHASE-FIELD SIMULATION OF QUANTUM DOT FORMATION IN A STRAINED HETEROEPITAXIAL FILM , 2004 .

[16]  Long-Qing Chen.  Phase-Field Models for Microstructure Evolution , 2002 .

[17]  Qiang Du,et al.  Coarsening Kinetics of a Two Phase Mixture with Highly Disparate Diffusion Mobility , 2010 .

[18]  John W. Barrett,et al.  Finite element approximation of the Cahn-Hilliard equation with concentration dependent mobility , 1999, Math. Comput..

[19]  A. G. Khachaturi︠a︡n.  Theory of structural transformations in solids , 1983 .

[20]  Perry H Leo,et al.  A diffuse interface model for microstructural evolution in elastically stressed solids , 1998 .

[21]  E. Favvas,et al.  What is spinodal decomposition , 2008 .

[22]  Khachaturyan,et al.  Elastic strain energy of inhomogeneous solids. , 1995, Physical review. B, Condensed matter.

[23]  Yinnian He,et al.  On large time-stepping methods for the Cahn--Hilliard equation , 2007 .

[24]  E. Mello,et al.  Numerical study of the Cahn–Hilliard equation in one, two and three dimensions , 2004, cond-mat/0410772.

[25]  Richard Welford,et al.  A multigrid finite element solver for the Cahn-Hilliard equation , 2006, J. Comput. Phys..

[26]  Yunzhi Wang,et al.  Kinetics of strain-induced morphological transformation in cubic alloys with a miscibility gap , 1993 .

[27]  H. Garcke,et al.  Spinodal Decomposition in the Presence of Elastic Interactions , 2003 .

[28]  Shenyang Hu.  Phase-field Models of Microstructure Evolution in a System with Elastic Inhomogeneity and Defects , 2004 .

[29]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[30]  P. Fratzl,et al.  A Possible Criterion for Slowing Down of Precipitate Coarsening due to Elastic Misfit Interactions , 1995 .

[31]  John W. Cahn,et al.  On Spinodal Decomposition , 1961 .