The scalar auxiliary variable (SAV) approach for gradient flows

Abstract We propose a new approach, which we term as scalar auxiliary variable (SAV) approach, to construct efficient and accurate time discretization schemes for a large class of gradient flows. The SAV approach is built upon the recently introduced IEQ approach. It enjoys all advantages of the IEQ approach but overcomes most of its shortcomings. In particular, the SAV approach leads to numerical schemes that are unconditionally energy stable and extremely efficient in the sense that only decoupled equations with constant coefficients need to be solved at each time step. The scheme is not restricted to specific forms of the nonlinear part of the free energy, so it applies to a large class of gradient flows. Numerical results are presented to show that the accuracy and effectiveness of the SAV approach over the existing methods.

[1]  Santiago Badia,et al.  Finite element approximation of nematic liquid crystal flows using a saddle-point structure , 2011, J. Comput. Phys..

[2]  Daozhi Han,et al.  Linearly first- and second-order, unconditionally energy stable schemes for the phase field crystal model , 2017, J. Comput. Phys..

[3]  James J. Feng,et al.  A diffuse-interface method for simulating two-phase flows of complex fluids , 2004, Journal of Fluid Mechanics.

[4]  Xiaofeng Yang,et al.  Numerical approximations for a phase-field moving contact line model with variable densities and viscosities , 2017, J. Comput. Phys..

[5]  Tao Tang,et al.  Long Time Numerical Simulations for Phase-Field Problems Using p-Adaptive Spectral Deferred Correction Methods , 2015, SIAM J. Sci. Comput..

[6]  Lili Ju,et al.  Efficient linear schemes with unconditional energy stability for the phase field elastic bending energy model , 2017 .

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

[8]  L. Greengard,et al.  Spectral Deferred Correction Methods for Ordinary Differential Equations , 2000 .

[9]  Junseok Kim.  Phase-Field Models for Multi-Component Fluid Flows , 2012 .

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

[11]  S. Cox,et al.  Exponential Time Differencing for Stiff Systems , 2002 .

[12]  Xiaofeng Yang,et al.  Numerical approximations for a three-component Cahn–Hilliard phase-field model based on the invariant energy quadratization method , 2017 .

[13]  Lili Ju,et al.  Linear and unconditionally energy stable schemes for the binary fluid–surfactant phase field model , 2017, 1701.07446.

[14]  Martin Grant,et al.  Modeling elasticity in crystal growth. , 2001, Physical review letters.

[15]  Jie Shen,et al.  Numerical simulations of jet pinching-off and drop formation using an energetic variational phase-field method , 2006, J. Comput. Phys..

[16]  S. Edwards,et al.  The Theory of Polymer Dynamics , 1986 .

[17]  Bath Ba.  THE GLOBAL DYNAMICS OF DISCRETE SEMILINEAR PARABOLIC EQUATIONS , 1993 .

[18]  M. Gurtin,et al.  TWO-PHASE BINARY FLUIDS AND IMMISCIBLE FLUIDS DESCRIBED BY AN ORDER PARAMETER , 1995, patt-sol/9506001.

[19]  F. M. Leslie.  Theory of Flow Phenomena in Liquid Crystals , 1979 .

[20]  Xiaofeng Yang,et al.  Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends , 2016, J. Comput. Phys..

[21]  Jun Li,et al.  Energy Stable Numerical Schemes for a Hydrodynamic Model of Nematic Liquid Crystals , 2016, SIAM J. Sci. Comput..

[22]  Jie Shen,et al.  Modeling and simulations of drop pinch-off from liquid crystal filaments and the leaky liquid crystal faucet immersed in viscous fluids , 2013, J. Comput. Phys..

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

[24]  Xiaofeng Yang,et al.  A novel linear second order unconditionally energy stable scheme for a hydrodynamic Q-tensor model of liquid crystals , 2017 .

[25]  Xiaofeng Yang,et al.  Error analysis of stabilized semi-implicit method of Allen-Cahnequation , 2009 .

[26]  Dong Li,et al.  On Second Order Semi-implicit Fourier Spectral Methods for 2D Cahn–Hilliard Equations , 2017, J. Sci. Comput..

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

[28]  Tao Tang,et al.  Stability Analysis of Large Time-Stepping Methods for Epitaxial Growth Models , 2006, SIAM J. Numer. Anal..

[29]  D. M. Anderson,et al.  DIFFUSE-INTERFACE METHODS IN FLUID MECHANICS , 1997 .

[30]  Jie Shen,et al.  A decoupled energy stable scheme for a hydrodynamic phase-field model of mixtures of nematic liquid crystals and viscous fluids , 2016, J. Comput. Phys..

[31]  Jie Shen,et al.  Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[32]  J. Cahn,et al.  A microscopic theory for antiphase boundary motion and its application to antiphase domain coasening , 1979 .

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

[34]  Bo Li,et al.  Epitaxial Growth Without Slope Selection: Energetics, Coarsening, and Dynamic Scaling , 2004, J. Nonlinear Sci..