Unconditionally energy stable DG schemes for the Swift-Hohenberg equation

10/01/2019
by   Hailiang Liu, et al.
0

The Swift-Hohenberg equation as a central nonlinear model in modern physics has a gradient flow structure. Here we introduce fully discrete discontinuous Galerkin (DG) schemes for a class of fourth order gradient flow problems, including the nonlinear Swift-Hohenberg equation, to produce free-energy-decaying discrete solutions, irrespective of the time step and the mesh size. We exploit and extend the mixed DG method introduced in [H. Liu and P. Yin, J. Sci. Comput., 77: 467–501, 2018] for the spatial discretization, and the "Invariant Energy Quadratization" method for the time discretization. The resulting IEQ-DG algorithms are linear, thus they can be efficiently solved without resorting to any iteration method. We actually prove that these schemes are unconditionally energy stable. We present several numerical examples that support our theoretical results and illustrate the efficiency, accuracy and energy stability of our new algorithm. The numerical results on two dimensional pattern formation problems indicate that the method is able to deliver comparable patterns of high accuracy.

READ FULL TEXT

page 17

page 18

page 19

research
12/21/2019

Unconditionally energy stable discontinuous Galerkin schemes for the Cahn-Hilliard equation

In this paper, we introduce novel discontinuous Galerkin (DG) schemes fo...
research
11/04/2022

A numerical study of vortex nucleation in 2D rotating Bose-Einstein condensates

This article introduces a new numerical method for the minimization unde...
research
06/25/2021

Numerical analysis of the LDG method for large deformations of prestrained plates

A local discontinuous Galerkin (LDG) method for approximating large defo...
research
08/27/2020

On the SAV-DG method for a class of fourth order gradient flows

For a class of fourth order gradient flow problems, integration of the s...
research
03/21/2023

On the Stability of IMEX Upwind gSBP Schemes for Linear Advection-Diffusion Equations

A fully discrete energy stability analysis is carried out for linear adv...
research
10/21/2021

Q-tensor gradient flow with quasi-entropy and discretizations preserving physical constraints

We propose and analyze numerical schemes for the gradient flow of Q-tens...

Please sign up or login with your details

Forgot password? Click here to reset