Stable Phase Field Approximations of Anisotropic Solidification

Abstract

We introduce unconditionally stable finite element approximations for a phase field model for solidification, which take highly anisotropic surface energy and kinetic effects into account. We hence approximate Stefan problems with anisotropic Gibbs--Thomson law with kinetic undercooling, and quasi-static variants thereof. The phase field model is given by {align*} \vartheta\,w_t + λ\,\varrho(φ)\,φ_t & = \nabla \,.\, (b(φ)\,\nabla\, w) \,, \cPsi\,\tfrac{a}α\,\varrho(φ)\,w & = ε\,\tfracρα\,μ(\nabla\,φ)\,φ_t -ε\,\nabla \,.\, A'(\nabla\, φ) + ε^{-1}\,Ψ'(φ) {align*} subject to initial and boundary conditions for the phase variable $\phi$ and the temperature approximation $w$. Here $\epsilon >0$ is the interfacial parameter, $\Psi$ is a double well potential, $\cPsi = \int_{-1}^1 \sqrt{2\,Ψ(s)}\;{\rm d}s$, $\varrho$ is a shape function and $A(\nabla \phi )=\frac{1}{2}|\gamma (\nabla \phi ){|}^2$, where $\gamma$ is the anisotropic density function. Moreover, $\vartheta \ge 0$, $\lambda >0$, $a>0$, $\alpha >0$ and $\rho \ge 0$ are physical parameters from the Stefan problem, while $b$ and $\mu$ are coefficient functions which also relate to the sharp interface problem. On introducing the novel fully practical finite element approximations for the anisotropic phase field model, we prove their stability and demonstrate their applicability with some numerical results.

40 pages, 23 figures

References in corpus (2)

• Numerical computations of facetted pattern formation in snow crystal growth
• On the stable discretization of strongly anisotropic phase field models with applications to crystal growth

Cited by in corpus (12)

• Parametric finite element approximations of curvature driven interface evolutions
• Finite Element Approximation for the Dynamics of Fluidic Two-Phase Biomembranes
• Cahn-Hilliard-Brinkman systems for tumour growth
• A diffuse-interface approach for solid-state dewetting with anisotropic surface energies
• A structure preserving front tracking finite element method for the Mullins--Sekerka problem
• Phase Field Models versus Parametric Front Tracking Methods: Are they accurate and computationally efficient?
• Doubly Degenerate Diffuse Interface Models of Anisotropic Surface Diffusion
• A novel finite element approximation of anisotropic curve shortening flow
• The anisotropic Cahn--Hilliard equation: regularity theory and strict separation properties
• Optimal control of a quasilinear parabolic equation and its time discretization
• An unconditionally stable finite element scheme for anisotropic curve shortening flow
• Overhang penalization in additive manufacturing via phase field structural topology optimization with anisotropic energies