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

by   Ying Wei, et al.
Institute of Software, Chinese Academy of Sciences

In this paper, we present a parallel numerical algorithm for solving the phase field crystal equation. In the algorithm, a semi-implicit finite difference scheme is derived based on the discrete variational derivative method. Theoretical analysis is provided to show that the scheme is unconditionally energy stable and can achieve second-order accuracy in both space and time. An adaptive time step strategy is adopted such that the time step size can be flexibly controlled based on the dynamical evolution of the problem. At each time step, a nonlinear algebraic system is constructed from the discretization of the phase field crystal equation and solved by a domain decomposition based, parallel Newton--Krylov--Schwarz method with improved boundary conditions for subdomain problems. Numerical experiments with several two and three dimensional test cases 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 11

page 14

page 16

page 18

page 19


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

In this paper, we study numerical methods for solving the coupled Allen-...

Energy stable schemes for gradient flows based on the DVD method

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

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

A massively parallel semi-Lagrangian solver for the six-dimensional Vlasov-Poisson equation

This paper presents an optimized and scalable semi-Lagrangian solver for...

A spatio-temporal adaptive phase-field fracture method

We present an energy-preserving mechanic formulation for dynamic quasi-b...

Parallel finite volume simulation of the spherical shell dynamo with pseudo-vacuum magnetic boundary conditions

In this paper, we study the parallel simulation of the magnetohydrodynam...

Please sign up or login with your details

Forgot password? Click here to reset