Energy stable arbitrary order ETD-MS method for gradient flows with Lipschitz nonlinearity

02/22/2021
by   Wenbin Chen, et al.
0

We present a methodology to construct efficient high-order in time accurate numerical schemes for a class of gradient flows with appropriate Lipschitz continuous nonlinearity. There are several ingredients to the strategy: the exponential time differencing (ETD), the multi-step (MS) methods, the idea of stabilization, and the technique of interpolation. They are synthesized to develop a generic k^th order in time efficient linear numerical scheme with the help of an artificial regularization term of the form Aτ^k∂/∂ tℒ^p(k)u where ℒ is the positive definite linear part of the flow, τ is the uniform time step-size. The exponent p(k) is determined explicitly by the strength of the Lipschitz nonlinear term in relation to ℒ together with the desired temporal order of accuracy k. To validate our theoretical analysis, the thin film epitaxial growth without slope selection model is examined with a fourth-order ETD-MS discretization in time and Fourier pseudo-spectral in space discretization. Our numerical results on convergence and energy stability are in accordance with our theoretical results.

READ FULL TEXT

page 17

page 18

page 22

page 23

page 24

research
10/16/2019

Arbitrarily High-order Linear Schemes for Gradient Flow Models

We present a paradigm for developing arbitrarily high order, linear, unc...
research
03/04/2022

A general class of linear unconditionally energy stable schemes for the gradient flows

This paper studies a class of linear unconditionally energy stable schem...
research
10/31/2020

A third order BDF energy stable linear scheme for the no-slope-selection thin film model

In this paper we propose and analyze a (temporally) third order accurate...
research
11/09/2021

An F-modulated stability framework for multistep methods

We introduce a new 𝐅-modulated energy stability framework for general li...
research
03/01/2021

Stability and conservation properties of Hermite-based approximations of the Vlasov-Poisson system

Spectral approximation based on Hermite-Fourier expansion of the Vlasov-...
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...

Please sign up or login with your details

Forgot password? Click here to reset