TY - JOUR
T1 - A novel diffuse-interface model and a fully-discrete maximum-principle-preserving energy-stable method for two-phase flow with surface tension and non-matching densities
AU - ten Eikelder, M. F.P.
AU - Akkerman, I.
PY - 2021
Y1 - 2021
N2 - Two well-established classes of the interface capturing models are the level-set and phase-field models. Level-set formulations satisfy the maximum principle for the density but are not energy-stable. On the other hand, phase-field models do satisfy the second law of thermodynamics but lack the maximum principle for the density. In this paper we derive a novel model for incompressible immiscible two-phase flow with non-matching densities and surface tension that is both energetically-stable and satisfies the maximum principle for the density. The model finds its place at the intersection of level-set and phase-field models. Its derivation is based on a diffusification of the incompressible two-phase Navier–Stokes equations with non-matching densities and surface tension and involves functional entropy variables. Additionally, we present an associated fully-discrete energy-stable method. Isogeometric analysis is used for the spatial discretization and the temporal-integration is performed with a new time-integration scheme that is a perturbation of the second-order midpoint scheme. The fully-discrete scheme is unconditionally energy-dissipative, pointwise divergence-free and satisfies the maximum principle for the density. Numerical examples in two and three dimensions verify the energetic-stability of the methodology.
AB - Two well-established classes of the interface capturing models are the level-set and phase-field models. Level-set formulations satisfy the maximum principle for the density but are not energy-stable. On the other hand, phase-field models do satisfy the second law of thermodynamics but lack the maximum principle for the density. In this paper we derive a novel model for incompressible immiscible two-phase flow with non-matching densities and surface tension that is both energetically-stable and satisfies the maximum principle for the density. The model finds its place at the intersection of level-set and phase-field models. Its derivation is based on a diffusification of the incompressible two-phase Navier–Stokes equations with non-matching densities and surface tension and involves functional entropy variables. Additionally, we present an associated fully-discrete energy-stable method. Isogeometric analysis is used for the spatial discretization and the temporal-integration is performed with a new time-integration scheme that is a perturbation of the second-order midpoint scheme. The fully-discrete scheme is unconditionally energy-dissipative, pointwise divergence-free and satisfies the maximum principle for the density. Numerical examples in two and three dimensions verify the energetic-stability of the methodology.
KW - Energy dissipation
KW - Incompressible two-phase flow
KW - Isogeometric analysis
KW - Level-set formulations
KW - Phase-field models
KW - Surface tension
UR - http://www.scopus.com/inward/record.url?scp=85102509657&partnerID=8YFLogxK
U2 - 10.1016/j.cma.2021.113751
DO - 10.1016/j.cma.2021.113751
M3 - Article
AN - SCOPUS:85102509657
SN - 0045-7825
VL - 379
JO - Computer Methods in Applied Mechanics and Engineering
JF - Computer Methods in Applied Mechanics and Engineering
M1 - 113751
ER -