A new class of efficient and robust energy stable schemes for gradient flows

Abstract

We propose a new numerical technique to deal with nonlinear terms in gradient flows. By introducing a scalar auxiliary variable (SAV), we construct efficient and robust energy stable schemes for a large class of gradient flows. The SAV approach is not restricted to specific forms of the nonlinear part of the free energy, and only requires to solve {\it decoupled} linear equations with {\it constant coefficients}. We use this technique to deal with several challenging applications which can not be easily handled by existing approaches, and present convincing numerical results to show that our schemes are not only much more efficient and easy to implement, but can also better capture the physical properties in these models. Based on this SAV approach, we can construct unconditionally second-order energy stable schemes; and we can easily construct even third or fourth order BDF schemes, although not unconditionally stable, which are very robust in practice. In particular, when coupled with an adaptive time stepping strategy, the SAV approach can be extremely efficient and accurate.

Cited by in corpus (12)

• A new Lagrange Multiplier approach for gradient flows
• Structure-preserving algorithms for the two-dimensional sine-Gordon equation with Neumann boundary conditions
• Arbitrarily High-order Unconditionally Energy Stable Schemes for Thermodynamically Consistent Gradient Flow Models
• An Accurate and Efficient Algorithm for The Time-fractional Molecular Beam Epitaxy Model with Slope Selection
• Second-order Decoupled Energy-stable Schemes for Cahn-Hilliard-Navier-Stokes equations
• A fully discrete energy stable scheme for a phase-field moving contact line model with variable densities and viscosities
• A structure-preserving algorithm for the fractional nonlinear Schrödinger equation based on the SAV approach
• On a SAV-MAC scheme for the Cahn-Hilliard-Navier-Stokes Phase Field Model
• A linearly implicit and local energy-preserving scheme for the sine-Gordon equation based on the invariant energy quadratization approach
• A linearly implicit structure-preserving scheme for the fractional sine-Gordon equation based on the IEQ approach
• Energy Stability and Convergence of SAV Block-centered Finite Difference Method for Gradient Flows
• The fast scalar auxiliary variable approach with unconditional energy stability for nonlocal Cahn-Hilliard equation