Parallel energy-stable solver for a coupled Allen-Cahn and Cahn-Hilliard system

by   Jizu Huang, et al.
Institute of Software, Chinese Academy of Sciences
Chinese Academy of Science
Peking University

In this paper, we study numerical methods for solving the coupled Allen-Cahn/Cahn-Hilliard system associated with a free energy functional of logarithmic type. To tackle the challenge posed by the special free energy functional, we propose a method to approximate the discrete variational derivatives in polynomial forms, such that the corresponding finite difference scheme is unconditionally energy stable and the energy dissipation law is maintained. To further improve the performance of the algorithm, a modified adaptive time stepping strategy is adopted such that the time step size can be flexibly controlled based on the dynamical evolution of the problem. To achieve high performance on parallel computers, we introduce a domain decomposition based, parallel Newton-Krylov-Schwarz method to solve the nonlinear algebraic system constructed from the discretization at each time step. Numerical experiments show that the proposed algorithm is second-order accurate in both space and time, energy stable with large time steps, and highly scalable to over ten thousands processor cores on the Sunway TaihuLight supercomputer.


page 12

page 13

page 14


Parallel energy-stable phase field crystal simulations based on domain decomposition methods

In this paper, we present a parallel numerical algorithm for solving the...

Parallel energy stable phase field simulations of Ni-based alloys system

In this paper, we investigate numerical methods for solving Nickel-based...

Energy stable schemes for gradient flows based on the DVD method

The existing discrete variational derivative method is only second-order...

Energy dissipation of Strang splitting for Allen–Cahn

We consider a class of second-order Strang splitting methods for Allen-C...

Parallel Multi-Stage Preconditioners with Adaptive Setup for the Black Oil Model

The black oil model is widely used to describe multiphase porous media f...

On a sinc-type MBE model

We introduce a new sinc-type molecular beam epitaxy model which is deriv...

Finite element solution of nonlocal Cahn-Hilliard equations with feedback control time step size adaptivity

In this study, we evaluate the performance of feedback control-based tim...

Please sign up or login with your details

Forgot password? Click here to reset