An Isotropic Discretization with Semi-implicit Approach for Phase Field Model of Alloy Solidification

TL;DR

This paper introduces a hexagonal mesh discretization strategy combined with a semi-implicit algorithm for phase field models, significantly reducing anisotropy and improving computational efficiency in alloy solidification simulations.

Contribution

It presents a novel isotropic discretization method using a hexagonal mesh and semi-implicit approach, enhancing accuracy and efficiency over existing methods.

Findings

Reduces lattice-induced anisotropy in simulations

Increases accuracy compared to previous discretization methods

Enhances computational efficiency with semi-implicit algorithm

Abstract

Quantitative phase field models have been extensively used to study the solidification behavior of alloys under different conditions. However, a longstanding challenge of phase field models is the directional bias caused by the discretization-induced lattice effects. In particular, widely used discretization methods may introduce significant spurious anisotropy for simulations of polycrystalline solidification. In this paper, we demonstrate a feasible 2D discretization strategy utilizing a hexagonal mesh to reduce the lattice-induced anisotropy of the phase field model. The leading differential terms of the 2D discretization methods are analyzed by using known methods in Fourier space. Using Taylor expansion of discrete Fourier Transform up to sixth order, we found that the proposed discretization strategy is more accurate and isotropic than other methods, including the isotropic discretization recently proposed by Ji et al.[1]. Additionally, the proposed 2D discretization method can be easily incorporated into a semi-implicit algorithm to solve phase field equations, thereby greatly reducing time step constraints and improving computational efficiency compared to explicit approaches. To prove the accuracy and efficiency of the proposed isotropic discretization with semi-implicit algorithm, 2D simulations of alloy solidification with different discretization schemes were performed and compared. We show that the proposed discretization using a hexagonal mesh can drastically reduce grid-induced anisotropy compared to conventional methods.