| Popis: |
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. |
