A decoupled, stable, and linear FEM for a phase-field model of variable density two-phase incompressible surface flow
arXiv:2104.08996 · doi:10.1016/j.cma.2021.114167
Abstract
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.
arXiv comment22 pages, 5 figures, 1 table
References in corpus (9)
- Incompressible fluid problems on embedded surfaces: Modeling and variational formulations
- Mass Conservative and Energy Stable Finite Difference Methods for the Quasi-incompressible Navier-Stokes-Cahn-Hilliard system: Primitive Variable and Projection-Type Schemes
- A diffuse interface model for two-phase ferrofluid flows
- A fully discrete energy stable scheme for a phase-field moving contact line model with variable densities and viscosities
- Divergence-free tangential finite element methods for incompressible flows on surfaces
- Weak Solutions for a Non-Newtonian Diffuse Interface Model with Different Densities
- Inf-sup stability of the trace P2-P1 Taylor-Hood elements for surface PDEs
- A penalty finite element method for a fluid system posed on embedded surface
- Recycling augmented Lagrangian preconditioner in an incompressible fluid solver