A decoupled, stable, and linear FEM for a phase-field model of variable density two-phase incompressible surface flow

The paper considers a thermodynamically consistent phase-field model of a two-phase flow of incompressible viscous fluids. The model allows for a non-linear dependence of fluid density on the phase-field order parameter. Driven by applications in biomembrane studies, the model is written for tangential flows of fluids constrained to a surface and consists of (surface) Navier–Stokes–Cahn–Hilliard type equations. We apply an unfitted finite element method to discretize the system and introduce a fully discrete time-stepping scheme with the following properties: (i) the scheme decouples the fluid and phase-field equation solvers at each time step, (ii) the resulting two algebraic systems are linear, and (iii) the numerical solution satisfies the same stability bound as the solution of the original system under some restrictions on the discretization parameters. Numerical examples are provided to demonstrate the stability, accuracy, and overall efficiency of the approach. Our computational study of several two-phase surface flows reveals some interesting dependencies of flow statistics on the geometry.

[1]  M. Olshanskii.  A low order Galerkin finite element method for the Navier–Stokes equations of steady incompressible flow: a stabilization issue and iterative methods , 2002 .

[2]  Alexander Zhiliakov,et al.  Recycling augmented Lagrangian preconditioner in an incompressible fluid solver , 2020, Numer. Linear Algebra Appl..

[3]  Maxim A. Olshanskii,et al.  Inf-sup stability of the trace P2-P1 Taylor-Hood elements for surface PDEs , 2019, ArXiv.

[4]  Xiaobing Feng,et al.  Fully Discrete Finite Element Approximations of the Navier-Stokes-Cahn-Hilliard Diffuse Interface Model for Two-Phase Fluid Flows , 2006, SIAM J. Numer. Anal..

[5]  Xiaoming He,et al.  A diffuse interface model and semi-implicit energy stable finite element method for two-phase magnetohydrodynamic flows , 2019, Computer Methods in Applied Mechanics and Engineering.

[6]  Jia Zhao,et al.  An energy stable algorithm for a quasi-incompressible hydrodynamic phase-field model of viscous fluid mixtures with variable densities and viscosities , 2017, Comput. Phys. Commun..

[7]  Christoph Lehrenfeld,et al.  Divergence-free tangential finite element methods for incompressible flows on surfaces , 2019, ArXiv.

[8]  Junseok Kim,et al.  A phase-field model and its efficient numerical method for two-phase flows on arbitrarily curved surfaces in 3D space , 2020 .

[9]  David Kay,et al.  Finite element approximation of a Cahn−Hilliard−Navier−Stokes system , 2008 .

[10]  H. Garcke,et al.  Existence of weak solutions for a diffuse interface model for two-phase flow with surfactants , 2017, Communications on Pure & Applied Analysis.

[11]  Xiaogang Yang,et al.  Fully Discrete Second-Order Linear Schemes for Hydrodynamic Phase Field Models of Binary Viscous Fluid Flows with Variable Densities , 2018, SIAM J. Sci. Comput..

[12]  Maxim A. Olshanskii,et al.  A Finite Element Method for the Surface Stokes Problem , 2018, SIAM J. Sci. Comput..

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

[14]  Jan Giesselmann,et al.  A quasi-incompressible diffuse interface model with phase transition , 2012 .

[15]  I. Nitschke,et al.  A finite element approach to incompressible two-phase flow on manifolds , 2012, Journal of Fluid Mechanics.

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

[17]  Maxim A. Olshanskii,et al.  Incompressible fluid problems on embedded surfaces: Modeling and variational formulations , 2017, Interfaces and Free Boundaries.

[18]  A. Quaini,et al.  A computational study of lateral phase separation in biological membranes , 2018, International journal for numerical methods in biomedical engineering.

[19]  Shuyu Sun,et al.  Fully discrete energy stable scheme for a phase-field moving contact line model with variable densities and viscosities , 2018, Applied Mathematical Modelling.

[20]  Jie Shen,et al.  A time-stepping scheme involving constant coefficient matrices for phase-field simulations of two-phase incompressible flows with large density ratios , 2012, J. Comput. Phys..

[21]  Dominic Breit,et al.  Weak solutions for a non-Newtonian diffuse interface model with different densities , 2015, 1509.05663.

[22]  Qiang Du,et al.  Finite element approximation of the Cahn–Hilliard equation on surfaces , 2011 .

[23]  Shuyu Sun,et al.  Efficient energy-stable schemes for the hydrodynamics coupled phase-field model , 2019, Applied Mathematical Modelling.

[24]  Ping Lin,et al.  Mass conservative and energy stable finite difference methods for the quasi-incompressible Navier–Stokes–Cahn–Hilliard system: Primitive variable and projection-type schemes , 2017, 1703.06606.

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

[26]  Morton E. Gurtin,et al.  A continuum theory of elastic material surfaces , 1975 .

[27]  Charles M. Elliott,et al.  Numerical simulation of dealloying by surface dissolution via the evolving surface finite element method , 2008, J. Comput. Phys..

[28]  Franck Boyer,et al.  A theoretical and numerical model for the study of incompressible mixture flows , 2002 .

[29]  Chang Shu,et al.  Diffuse interface model for incompressible two-phase flows with large density ratios , 2007, J. Comput. Phys..

[30]  Béatrice Rivière,et al.  An efficient numerical algorithm for solving viscosity contrast Cahn-Hilliard-Navier-Stokes system in porous media , 2020, J. Comput. Phys..

[31]  E. Harald van Brummelen,et al.  Diffuse-interface two-phase flow models with different densities: a new quasi-incompressible form and a linear energy-stable method , 2016, 1603.06475.

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

[33]  Maxim A. Olshanskii,et al.  A Finite Element Method for Elliptic Equations on Surfaces , 2009, SIAM J. Numer. Anal..

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

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

[36]  Günther Grün,et al.  Two-phase flow with mass density contrast: Stable schemes for a thermodynamic consistent and frame-indifferent diffuse-interface model , 2012, J. Comput. Phys..

[37]  M. Olshanskii,et al.  A Stabilized Trace Finite Element Method for Partial Differential Equations on Evolving Surfaces , 2017, SIAM J. Numer. Anal..

[38]  Arnold Reusken,et al.  Analysis of trace finite element methods for surface partial differential equations , 2015 .

[39]  Maxim A. Olshanskii,et al.  A Penalty Finite Element Method for a Fluid System Posed on Embedded Surface , 2018, Journal of Mathematical Fluid Mechanics.

[40]  H. Abels,et al.  Thermodynamically Consistent, Frame Indifferent Diffuse Interface Models for Incompressible Two-Phase Flows with Different Densities , 2011, 1104.1336.

[41]  Christoph Lehrenfeld,et al.  Analysis of a High-Order Trace Finite Element Method for PDEs on Level Set Surfaces , 2016, SIAM J. Numer. Anal..

[42]  Charles M. Elliott,et al.  Finite element methods for surface PDEs* , 2013, Acta Numerica.

[43]  Annalisa Quaini,et al.  Experimental validation of a phase-field model to predict coarsening dynamics of lipid domains in multicomponent membranes , 2020, Biochimica et biophysica acta. Biomembranes.

[44]  G. Grün,et al.  On Convergent Schemes for Diffuse Interface Models for Two-Phase Flow of Incompressible Fluids with General Mass Densities , 2013, SIAM J. Numer. Anal..

[45]  A. Bonito,et al.  A divergence-conforming finite element method for the surface Stokes equation , 2019, SIAM J. Numer. Anal..

[46]  Ricardo H. Nochetto,et al.  A diffuse interface model for two-phase ferrofluid flows , 2016, 1601.06824.

[47]  Arnold Reusken,et al.  Error analysis of higher order Trace Finite Element Methods for the surface Stokes equation , 2019, J. Num. Math..

[48]  P. Hohenberg,et al.  Theory of Dynamic Critical Phenomena , 1977 .

[49]  Heike Emmerich,et al.  The Diffuse Interface Approach in Materials Science: Thermodynamic Concepts and Applications of Phase-Field Models , 2003 .

[50]  L. Quartapelle,et al.  A projection FEM for variable density incompressible flows , 2000 .

[51]  Harald Garcke,et al.  A stable and linear time discretization for a thermodynamically consistent model for two-phase incompressible flow , 2014, 1402.6524.

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