Numerical approximations for a three-component Cahn–Hilliard phase-field model based on the invariant energy quadratization method

How to develop efficient numerical schemes while preserving energy stability at the discrete level is challenging for the three-component Cahn–Hilliard phase-field model. In this paper, we develop a set of first- and second-order temporal approximation schemes based on a novel “Invariant Energy Quadratization” approach, where all nonlinear terms are treated semi-explicitly. Consequently, the resulting numerical schemes lead to well-posed linear systems with a linear symmetric, positive definite at each time step. We prove that the developed schemes are unconditionally energy stable and present various 2D and 3D numerical simulations to demonstrate the stability and the accuracy of the schemes.

[1]  P. Gennes.  Dynamics of fluctuations and spinodal decomposition in polymer blends , 1980 .

[2]  John S. Rowlinson,et al.  Molecular Theory of Capillarity , 1983 .

[3]  K. Binder.  Collective diffusion, nucleation, and spinodal decomposition in polymer mixtures , 1983 .

[4]  Charles M. Elliott,et al.  `A generalised diffusion equation for phase separation of a multi-component mixture with interfacial free energy' , 1991 .

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

[6]  Charles M. Elliott,et al.  Numerical analysis of a model for phase separation of a multi- component alloy , 1996 .

[7]  Harald Garcke,et al.  Diffusional phase transitions in multicomponent systems with a concentration dependent mobility matrix , 1997 .

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

[9]  J. Lowengrub,et al.  Quasi–incompressible Cahn–Hilliard fluids and topological transitions , 1998, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

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

[11]  Gunduz Caginalp,et al.  Convergence of the phase field model to its sharp interface limits , 1998, European Journal of Applied Mathematics.

[12]  Harald Garcke,et al.  A MultiPhase Field Concept: Numerical Simulations of Moving Phase Boundaries and Multiple Junctions , 1999, SIAM J. Appl. Math..

[13]  D. Jacqmin.  Regular Article: Calculation of Two-Phase Navier–Stokes Flows Using Phase-Field Modeling , 1999 .

[14]  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.

[15]  Harald Garcke,et al.  On Fully Practical Finite Element Approximations of Degenerate Cahn-Hilliard Systems , 2001 .

[16]  K R Elder,et al.  Sharp interface limits of phase-field models. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  John W. Barrett,et al.  An improved error bound for a finite element approximation of a model for phase separation of a multi-component alloy with a concentration dependent mobility matrix , 2001, Numerische Mathematik.

[18]  Andreas Prohl,et al.  Numerical analysis of the Allen-Cahn equation and approximation for mean curvature flows , 2003, Numerische Mathematik.

[19]  Jie Shen,et al.  A phase field model for the mixture of two incompressible fluids and its approximation by a Fourier-spectral method , 2003 .

[20]  Junseok Kim,et al.  CONSERVATIVE MULTIGRID METHODS FOR TERNARY CAHN-HILLIARD SYSTEMS ∗ , 2004 .

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

[22]  B. Stinner,et al.  Multicomponent alloy solidification: phase-field modeling and simulations. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[23]  Junseok Kim,et al.  Phase field modeling and simulation of three-phase flows , 2005 .

[24]  James J. Feng,et al.  Diffuse-interface simulations of drop coalescence and retraction in viscoelastic fluids , 2005 .

[25]  Harald Garcke,et al.  Second order phase field asymptotics for multi-component systems , 2006 .

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

[27]  Franck Boyer,et al.  Study of a three component Cahn-Hilliard flow model , 2006 .

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

[29]  Yinnian He,et al.  Analysis of finite element approximations of a phase field model for two-phase fluids , 2006, Math. Comput..

[30]  Junseok Kim,et al.  A second-order accurate non-linear difference scheme for the N -component Cahn–Hilliard system , 2008 .

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

[32]  Jie Shen,et al.  An efficient moving mesh spectral method for the phase-field model of two-phase flows , 2009, J. Comput. Phys..

[33]  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.

[34]  Xiaofeng Yang,et al.  Energy stable schemes for Cahn-Hilliard phase-field model of two-phase incompressible flows , 2010 .

[35]  Christian Miehe,et al.  A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits , 2010 .

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

[37]  Jie Shen,et al.  A Phase-Field Model and Its Numerical Approximation for Two-Phase Incompressible Flows with Different Densities and Viscosities , 2010, SIAM J. Sci. Comput..

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

[39]  Xiaofeng Yang,et al.  Journal of Non-newtonian Fluid Mechanics Shear Cell Rupture of Nematic Liquid Crystal Droplets in Viscous Fluids , 2011 .

[40]  V Mironov,et al.  Engineering a 3D, biological construct: representative research in the South Carolina Project for Organ Biofabrication , 2011, Biofabrication.

[41]  Steven M. Wise,et al.  An Energy Stable and Convergent Finite-Difference Scheme for the Modified Phase Field Crystal Equation , 2011, SIAM J. Numer. Anal..

[42]  Franck Boyer,et al.  Numerical schemes for a three component Cahn-Hilliard model , 2011 .

[43]  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..

[44]  V. Mironov,et al.  Modeling fusion of cellular aggregates in biofabrication using phase field theories. , 2012, Journal of theoretical biology.

[45]  Xiaofeng Yang,et al.  LCP droplet dispersions: a two-phase, diffuse-interface kinetic theory and global droplet defect predictions , 2012 .

[46]  Cheng Wang,et al.  A Linear Energy Stable Scheme for a Thin Film Model Without Slope Selection , 2012, J. Sci. Comput..

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

[48]  Qi Wang,et al.  A phase field approach for multicellular aggregate fusion in biofabrication. , 2013, Journal of biomechanical engineering.

[49]  Xiaofeng Yang,et al.  Mass and Volume Conservation in Phase Field Models for Binary Fluids , 2013 .

[50]  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..

[51]  Jia Zhao,et al.  A 3D Hydrodynamic Model for Cytokinesis of Eukaryotic Cells , 2014 .

[52]  Jie Shen,et al.  Decoupled Energy Stable Schemes for Phase-Field Models of Two-Phase Complex Fluids , 2014, SIAM J. Sci. Comput..

[53]  Keith Promislow,et al.  High accuracy solutions to energy gradient flows from material science models , 2014, J. Comput. Phys..

[54]  Jie Shen,et al.  Efficient energy stable numerical schemes for a phase field moving contact line model , 2015, J. Comput. Phys..

[55]  Jie Shen,et al.  Decoupled, Energy Stable Schemes for Phase-Field Models of Two-Phase Incompressible Flows , 2015, SIAM J. Numer. Anal..

[56]  Jie Shen,et al.  Decoupled Energy Stable Schemes for a Phase-Field Model of Two-Phase Incompressible Flows with Variable Density , 2014, Journal of Scientific Computing.

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

[58]  Xiaoming Wang,et al.  A second order in time, uniquely solvable, unconditionally stable numerical scheme for Cahn-Hilliard-Navier-Stokes equation , 2014, J. Comput. Phys..

[59]  Jia Zhao,et al.  A 3D Multi-Phase Hydrodynamic Model for Cytokinesis of Eukaryotic Cells , 2016 .

[60]  Jia Zhao,et al.  Modeling the Excess Cell Surface Stored in a Complex Morphology of Bleb-Like Protrusions , 2016, PLoS Comput. Biol..

[61]  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..

[62]  Ya Shen,et al.  A 3D numerical study of antimicrobial persistence in heterogeneous multi-species biofilms. , 2016, Journal of theoretical biology.

[63]  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..

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

[65]  Jia Zhao,et al.  Numerical approximations to a new phase field model for two phase flows of complex fluids , 2016 .

[66]  Xiaofeng Yang,et al.  Numerical Approximations for Allen-Cahn Type Phase Field Model of Two-Phase Incompressible Fluids with Moving Contact Lines , 2017 .

[67]  Xiaofeng Yang,et al.  Numerical approximations for a phase field dendritic crystal growth model based on the invariant energy quadratization approach , 2017 .

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

[69]  Daozhi Han,et al.  Numerical Analysis of Second Order, Fully Discrete Energy Stable Schemes for Phase Field Models of Two-Phase Incompressible Flows , 2017, J. Sci. Comput..

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

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

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

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

[74]  Jie Shen,et al.  Efficient and accurate numerical schemes for a hydro-dynamically coupled phase field diblock copolymer model , 2017, J. Comput. Phys..

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

[76]  Jia Zhao,et al.  Decoupled Energy Stable Schemes for a Phase Field Model of Three-Phase Incompressible Viscous Fluid Flow , 2016, Journal of Scientific Computing.