Time-dependent Ginzburg-Laudau (TDGL) equation is the simplest kinetic equation for the temporal evolution of a continuum field, which assumes that the rate of evolution of the field is linearly proportional to the thermodynamical driving force. The computation model based on this equation is also called phase field model. Phase field simulation can predict quite beautiful patterns of microstructures of material. It has been widely applied to simulating the evolution of microstructure by choosing different field variables. For example, using the single conserved field (concentration field), continuum phase field models has been employed to describe the pattern formation in phase-separating alloys (Nishimori and Onuki, 1990 Phys. Rev. B, 42,980) and the nanoscale pattern formation of an epitaxial monolayer (Lu and Suo, 2001 J. Mech. Phys. Solids, 49,1937). On the other hand, using the nonconserved field (polarization field), the phase field model has been utilized to simulating the formation of domain structure in ferroelectrics (Li et al. 2002 Acta Mater, 50,395). The thermodynamical driving force is usually nonlinear with respect to the field variable. In the case of nonlinearity, the solution to TDGL equation may not be unique. Different grid density, length of iteration step, initial state and random term (introduced to describe the nucleation process) may induce different results in the simulation. Does anyone investigate the effect of these factors on the final pattern? I wonder whether we can prove the solution is unique or not.
Recent comments