Multiple Scalar Auxiliary Variable (MSAV) Approach and its Application to the Phase-Field Vesicle Membrane Model

We consider in this paper gradient flows with disparate terms in the free energy that cannot be efficiently handled with the scalar auxiliary variable (SAV) approach, and we develop the multiple scalar auxiliary variable (MSAV) approach to deal with these cases. We apply the MSAV approach to the phase-field vesicle membrane (PF-VMEM) model which, in addition to some usual nonlinear terms in the free energy, has two additional penalty terms to enforce the volume and surface area. The MSAV approach enjoys the same computational advantages as the SAV approach but can handle free energies with multiple disparate terms such as the volume and surface area constraints in the PF-VMEM model. The MSAV schemes are unconditional energy stable and second-order accurate in time and lead to decoupled elliptic equations with constant coefficients to solve at each time step. Hence, these schemes are easy to implement and extremely efficient when coupled with an adaptive time stepping. Ample numerical results are presented...

[1]  Jian Zhang,et al.  A phase field model for vesicle-substrate adhesion , 2009, J. Comput. Phys..

[2]  Thomas J. R. Hughes,et al.  Provably unconditionally stable, second-order time-accurate, mixed variational methods for phase-field models , 2011, J. Comput. Phys..

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

[4]  Jie Shen,et al.  Convergence and Error Analysis for the Scalar Auxiliary Variable (SAV) Schemes to Gradient Flows , 2018, SIAM J. Numer. Anal..

[5]  P. Canham.  The minimum energy of bending as a possible explanation of the biconcave shape of the human red blood cell. , 1970, Journal of theoretical biology.

[6]  W. Helfrich.  Elastic Properties of Lipid Bilayers: Theory and Possible Experiments , 1973, Zeitschrift fur Naturforschung. Teil C: Biochemie, Biophysik, Biologie, Virologie.

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

[8]  Q. Du,et al.  Modelling and simulations of multi-component lipid membranes and open membranes via diffuse interface approaches , 2006, Journal of Mathematical Biology.

[9]  Max Gunzburger,et al.  A two phase field model for tracking vesicle–vesicle adhesion , 2016, Journal of mathematical biology.

[10]  Cheng Wang,et al.  Stable and efficient finite-difference nonlinear-multigrid schemes for the phase field crystal equation , 2009, J. Comput. Phys..

[11]  Q. Du,et al.  A phase field approach in the numerical study of the elastic bending energy for vesicle membranes , 2004 .

[12]  Tao Tang,et al.  An Adaptive Time-Stepping Strategy for the Molecular Beam Epitaxy Models , 2011, SIAM J. Sci. Comput..

[13]  W. Marsden.  I and J , 2012 .

[14]  Jie Shen,et al.  On the maximum principle preserving schemes for the generalized Allen–Cahn equation , 2016 .

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

[16]  Francisco Guillén-González,et al.  Unconditionally energy stable numerical schemes for phase-field vesicle membrane model , 2018, J. Comput. Phys..

[17]  Jiang Yang,et al.  The scalar auxiliary variable (SAV) approach for gradient flows , 2018, J. Comput. Phys..

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

[19]  Qiang Du,et al.  Simulating the deformation of vesicle membranes under elastic bending energy in three dimensions , 2006, J. Comput. Phys..

[20]  Axel Voigt,et al.  Phase-field modeling of the dynamics of multicomponent vesicles: Spinodal decomposition, coarsening, budding, and fission. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  Qiang Du,et al.  Efficient and stable exponential time differencing Runge-Kutta methods for phase field elastic bending energy models , 2016, J. Comput. Phys..

[22]  Richard S. Chadwick,et al.  Axisymmetric Indentation of a Thin Incompressible Elastic Layer , 2002, SIAM J. Appl. Math..

[23]  Qiang Du,et al.  ANALYSIS OF A PHASE FIELD NAVIER-STOKES VESICLE-FLUID INTERACTION MODEL , 2007 .

[24]  Jie Shen,et al.  Second-order Convex Splitting Schemes for Gradient Flows with Ehrlich-Schwoebel Type Energy: Application to Thin Film Epitaxy , 2012, SIAM J. Numer. Anal..

[25]  Xiaofeng Yang,et al.  Decoupled energy stable schemes for phase-field vesicle membrane model , 2015, J. Comput. Phys..