TY - JOUR
T1 - On phase-field modeling with a highly anisotropic interfacial energy
AU - Fleck, M.
AU - Mushonger, L.
AU - Pilipenko, D.
AU - Ankit, K.
AU - Emmerich, H.
PY - 2011/10
Y1 - 2011/10
N2 - We report on phase-field approaches that allow for anisotropies sufficiently high so that the interface develops sharp corners due to missing crystallographic orientations. This implies the necessity of a regularization that enforces local equilibrium at the corners, and we use the method of Eggelston et al. (Physica D 150, 91 (2001)), generalized to arbitrary crystal symmetries and rotations of the crystalline axes. Two different anisotropic phase-field formulations are presented and discussed: The classical model that allows the interface to vary with orientation, and another more recent formulation that has a constant interface width. We develop an explicit finite-difference scheme that combines a two-step differentiation with a stagnation grid formulation. The presented numerical implementation is stable and accurate enough to account for odd crystal symmetries and high angle rotations of the initial crystalline orientation. Even in the case of highly anisotropic interfacial energies, both formulations show excellent agreement with the well-known Wulff construction of the equilibrium shape of a particle embedded in a matrix.
AB - We report on phase-field approaches that allow for anisotropies sufficiently high so that the interface develops sharp corners due to missing crystallographic orientations. This implies the necessity of a regularization that enforces local equilibrium at the corners, and we use the method of Eggelston et al. (Physica D 150, 91 (2001)), generalized to arbitrary crystal symmetries and rotations of the crystalline axes. Two different anisotropic phase-field formulations are presented and discussed: The classical model that allows the interface to vary with orientation, and another more recent formulation that has a constant interface width. We develop an explicit finite-difference scheme that combines a two-step differentiation with a stagnation grid formulation. The presented numerical implementation is stable and accurate enough to account for odd crystal symmetries and high angle rotations of the initial crystalline orientation. Even in the case of highly anisotropic interfacial energies, both formulations show excellent agreement with the well-known Wulff construction of the equilibrium shape of a particle embedded in a matrix.
UR - http://www.scopus.com/inward/record.url?scp=84860820936&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84860820936&partnerID=8YFLogxK
U2 - 10.1140/epjp/i2011-11095-7
DO - 10.1140/epjp/i2011-11095-7
M3 - Article
AN - SCOPUS:84860820936
VL - 126
SP - 1
EP - 11
JO - European Physical Journal Plus
JF - European Physical Journal Plus
SN - 2190-5444
IS - 10
M1 - 95
ER -