A new Lagrange Multiplier approach for gradient flows

Abstract We propose a new Lagrange multiplier approach to design unconditional energy stable schemes for gradient flows. The new approach leads to unconditionally energy stable schemes that are as accurate and efficient as the recently proposed SAV approach (Shen, Xu, and Yang 2018), but enjoys two additional advantages: (i) schemes based on the new approach dissipate the original energy, as opposed to a modified energy in the recently proposed SAV approach (Shen, Xu, and Yang 2018);and (ii) they do not require the nonlinear part of the free energy to be bounded from below as is required in the SAV approach. The price we pay for these advantages is that a nonlinear algebraic equation has to be solved to determine the Lagrange multiplier. We present ample numerical results to validate the new approach, and, as a particular example of challenging applications, we consider a block copolymer (BCP)/coupled Cahn–Hilliard model, and carry out new and nontrivial simulations which are consistent with experiment results.

[1]  Q. Du,et al.  The phase field method for geometric moving interfaces and their numerical approximations , 2019, Geometric Partial Differential Equations - Part I.

[2]  H. Frieboes,et al.  Three-dimensional multispecies nonlinear tumor growth--I Model and numerical method. , 2008, Journal of theoretical biology.

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

[4]  Jie Shen,et al.  Highly Efficient and Accurate Numerical Schemes for the Epitaxial Thin Film Growth Models by Using the SAV Approach , 2019, J. Sci. Comput..

[5]  Y. Nishiura,et al.  Frustrated phases under three-dimensional confinement simulated by a set of coupled Cahn-Hilliard equations. , 2016, Soft matter.

[6]  A. Karma,et al.  Phase-Field Simulation of Solidification , 2002 .

[7]  J. Villain.  Continuum models of crystal growth from atomic beams with and without desorption , 1991 .

[8]  A. A. Wheeler,et al.  Thermodynamically-consistent phase-field models for solidification , 1992 .

[9]  Qing Cheng,et al.  Multiple Scalar Auxiliary Variable (MSAV) Approach and its Application to the Phase-Field Vesicle Membrane Model , 2018, SIAM J. Sci. Comput..

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

[11]  Yu U. Wang,et al.  Phase field microelasticity modeling of dislocation dynamics near free surface and in heteroepitaxial thin films , 2003 .

[12]  Christian Miehe,et al.  Thermodynamically consistent phase‐field models of fracture: Variational principles and multi‐field FE implementations , 2010 .

[13]  Cv Clemens Verhoosel,et al.  A phase-field description of dynamic brittle fracture , 2012 .

[14]  J. Lahann,et al.  Surface‐Reactive Patchy Nanoparticles and Nanodiscs Prepared by Tandem Nanoprecipitation and Internal Phase Separation , 2018, Advanced Functional Materials.

[15]  Bruce T. Murray,et al.  Adaptive phase-field computations of dendritic crystal growth , 1997 .

[16]  A. Karma,et al.  Phase-field simulations of dendritic crystal growth in a forced flow. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[18]  F. Schmid,et al.  Liquid-vapor phase behavior of a symmetrical binary fluid mixture , 1998, cond-mat/9801265.

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

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

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

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

[23]  R. Kobayashi.  Modeling and numerical simulations of dendritic crystal growth , 1993 .

[24]  A. Karma,et al.  Phase-field method for computationally efficient modeling of solidification with arbitrary interface kinetics. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[25]  J. Tinsley Oden,et al.  GENERAL DIFFUSE-INTERFACE THEORIES AND AN APPROACH TO PREDICTIVE TUMOR GROWTH MODELING , 2010 .

[26]  Alain Karma,et al.  Spiral Surface Growth without Desorption , 1998, cond-mat/9809358.

[27]  Phase-field simulations for evaporation with convection in liquid-vapor systems , 2005 .

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

[29]  Y. Nishiura,et al.  Transformation of Block Copolymer Nanoparticles from Ellipsoids with Striped Lamellae into Onionlike Spheres and Dynamical Control via Coupled Cahn–Hilliard Equations , 2018, ACS omega.