Energy-stable Runge-Kutta schemes for gradient flow models using the energy quadratization approach

Abstract In this letter, we present a novel class of arbitrarily high-order and unconditionally energy-stable algorithms for gradient flow models by combining the energy quadratization (EQ) technique and a specific class of Runge–Kutta (RK) methods, which is named the EQRK schemes. First of all, we introduce auxiliary variables to transform the original model into an equivalent system, with the transformed free energy a quadratic functional with respect to the new variables and the modified energy dissipative law is conserved. Then a special class of RK methods is employed for the reformulated system to arrive at structure-preserving time-discrete schemes. Along with rigorous proofs, numerical experiments are presented to demonstrate the accuracy and unconditionally energy-stability of the EQRK schemes.

[1]  Jie Shen,et al.  A New Class of Efficient and Robust Energy Stable Schemes for Gradient Flows , 2017, SIAM Rev..

[2]  Bo Li,et al.  Center for Scientific Computation And Mathematical Modeling , 2003 .

[3]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy and Free Energy of a Nonuniform System. III. Nucleation in a Two‐Component Incompressible Fluid , 2013 .

[4]  Xiaofeng Yang,et al.  Regularized linear schemes for the molecular beam epitaxy model with slope selection , 2018 .

[5]  Xiao Li,et al.  Energy stability and error estimates of exponential time differencing schemes for the epitaxial growth model without slope selection , 2018, Math. Comput..

[6]  Amanda E. Diegel,et al.  Stability and Convergence of a Second Order Mixed Finite Element Method for the Cahn-Hilliard Equation , 2014, 1411.5248.

[7]  Qiang Du,et al.  Analysis of Fully Discrete Approximations for Dissipative Systems and Application to Time-Dependent Nonlocal Diffusion Problems , 2018, J. Sci. Comput..

[8]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[9]  Jaemin Shin,et al.  Unconditionally stable methods for gradient flow using Convex Splitting Runge-Kutta scheme , 2017, J. Comput. Phys..

[10]  D. J. Eyre.  Unconditionally Gradient Stable Time Marching the Cahn-Hilliard Equation , 1998 .

[11]  Francisco Guillén-González,et al.  On linear schemes for a Cahn-Hilliard diffuse interface model , 2013, J. Comput. Phys..

[12]  J. Butcher.  Numerical methods for ordinary differential equations , 2003 .

[13]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[14]  Jun Li,et al.  Hydrodynamic Theories for Flows of Active Liquid Crystals and the Generalized Onsager Principle , 2016, Entropy.

[15]  Jia Zhao,et al.  A GENERAL STRATEGY FOR NUMERICAL APPROXIMATIONS OF NON-EQUILIBRIUM MODELS – PART I : THERMODYNAMICAL SYSTEMS , 2018 .

[16]  Jia Zhao,et al.  Linear second order in time energy stable schemes for hydrodynamic models of binary mixtures based on a spatially pseudospectral approximation , 2018, Adv. Comput. Math..

[17]  Xiaofeng Yang,et al.  Numerical approximations of Allen-Cahn and Cahn-Hilliard equations , 2010 .

[18]  Jia Zhao,et al.  Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method , 2017, J. Comput. Phys..

[19]  L. Onsager.  Reciprocal Relations in Irreversible Processes. II. , 1931 .

[20]  Steven M. Wise,et al.  Unconditionally stable schemes for equations of thin film epitaxy , 2010 .

[21]  Ignacio Romero,et al.  Thermodynamically consistent time‐stepping algorithms for non‐linear thermomechanical systems , 2009 .

[22]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .