Error analysis of SAV finite element method to phase field crystal model

TL;DR

This paper develops an energy stable SAV finite element scheme for the phase field crystal model, combining advanced numerical techniques, and verifies its accuracy and efficiency through rigorous analysis and numerical experiments.

Contribution

It introduces a novel energy stable SAV finite element scheme with adaptive strategy for the PFC model, improving computational efficiency and accuracy.

Findings

Scheme is first-order in time and second-order in space.

Numerical experiments confirm theoretical stability and accuracy.

Adaptive strategy effectively captures phase interfaces.

Abstract

In this paper, we construct and analyze an energy stable scheme by combining the latest developed scalar auxiliary variable (SAV) approach and linear finite element method (FEM) for phase field crystal (PFC) model, and show rigorously that the scheme is first-order in time and second-order in space for the L 2 and H -1 gradient flow equations. To reduce efficiently computational cost and capture accurately the phase interface, we give a simple adaptive strategy, equipped with a posteriori gradient estimator, i.e. L 2 norm of the recovered gradient. Extensive numerical experiments are presented to verify our theoretical results and to demonstrate the effectiveness and accuracy of our proposed method.