An Energy Stable Linear Diffusive Crank-Nicolson Scheme for the Cahn-Hilliard Gradient Flow

04/10/2020
by   Lin Wang, et al.
0

We propose and analyze a linearly stabilized semi-implicit diffusive Crank–Nicolson scheme for the Cahn–Hilliard gradient flow. In this scheme, the nonlinear bulk force is treated explicitly with two second-order stabilization terms. This treatment leads to linear elliptic system with constant coefficients and provable discrete energy dissipation. Rigorous error analysis is carried out for the fully discrete scheme. When the time step-size and the space step-size are small enough, second order accuracy in time is obtained with a prefactor controlled by some lower degree polynomial of 1/ε. Here ε is the thickness of the interface. Numerical results together with an adaptive time stepping are presented to verify the accuracy and efficiency of the proposed scheme.

READ FULL TEXT
research
06/15/2022

A linear adaptive BDF2 scheme for phase field crystal equation

In this paper, we present and analyze a linear fully discrete second ord...
research
08/24/2023

A second-order length-preserving and unconditionally energy stable rotational discrete gradient method for Oseen-Frank gradient flows

We present a second-order strictly length-preserving and unconditionally...
research
04/01/2022

An implicit–explicit second order BDF numerical scheme with variable steps for gradient flows

In this paper, we propose and analyze an efficient implicit–explicit (IM...
research
07/28/2023

An Unconditionally Energy-Stable and Orthonormality-Preserving Iterative Scheme for the Kohn-Sham Gradient Flow Based Model

We propose an unconditionally energy-stable, orthonormality-preserving, ...
research
07/18/2019

Doubly-Adaptive Artificial Compression Methods for Incompressible Flow

This report presents adaptive artificial compression methods in which th...
research
07/30/2020

An oscillation-free fully partitioned scheme for the numerical modeling of cardiac active mechanics

In silico models of cardiac electromechanics couple together mathematica...
research
09/06/2022

Second order, unconditionally stable, linear ensemble algorithms for the magnetohydrodynamics equations

We propose two unconditionally stable, linear ensemble algorithms with p...

Please sign up or login with your details

Forgot password? Click here to reset