2022/09/16 by Jizu Huang, Chao Yang, Huang, Jizu +1
Engineering · Materials Science · #FOS: Mathematics #Fluid Dynamics and Thin Films #Magnetic Properties and Applications #Numerical Analysis (math.NA) #Solidification and crystal growth phenomena
paper · pdf · doi:10.48550/arxiv.2209.08001
openalex publication_date 2022/09/16 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
In this paper, we investigate numerical methods for solving Nickel-based phase field system related to free energy, including the elastic energy and logarithmic type functionals. To address the challenge posed by the particular free energy functional, we propose a semi-implicit scheme based on the discrete variational derivative method, which is unconditionally energy stable and maintains the energy dissipation law and the mass conservation law. Due to the good stability of the semi-implicit scheme, the adaptive time step strategy is adopted, which can flexibly control the time step according to the dynamic evolution of the problem. A domain decomposition based, parallel Newton--Krylov--Schwarz method is introduced to solve the nonlinear algebraic system constructed by the discretization at each time step. Numerical experiments show that the proposed algorithm is energy stable with large time steps, and highly scalable to six thousand processor cores.