A positivity preserving and conservative variational scheme for phase-field modeling of two-phase flows

Abstract

We present a positivity preserving variational scheme for the phase-field modeling of incompressible two-phase flows with high density ratio and using meshes of arbitrary topology. The variational finite element technique relies on the Allen-Cahn phase-field equation for capturing the phase interface on a fixed mesh with a mass conservative and energy-stable discretization. Mass is conserved by enforcing a Lagrange multiplier which has both temporal and spatial dependence on the solution of the phase-field equation. The spatial part of the Lagrange multiplier is written as a mid-point approximation to make the scheme energy-stable. This enables us to form a conservative, energy-stable and positivity preserving scheme. The proposed variational technique reduces spurious and unphysical oscillations in the solution while maintaining second-order spatial accuracy. To model a generic two-phase free-surface flow, we couple the Allen-Cahn phase-field equation with the Navier-Stokes equations. Comparison of results between standard linear stabilized finite element method and the present variational formulation shows a remarkable reduction of oscillations in the solution while retaining the boundedness of the phase-indicator field. We perform a standalone test to verify the accuracy and stability of the Allen-Cahn two-phase solver. Standard two-phase flow benchmarks such as Laplace-Young law and sloshing tank problem are carried out to assess the convergence and accuracy of the coupled Navier-Stokes and Allen-Cahn solver. Two- and three-dimensional dam break problem are then solved to assess the scheme for the problem with topological changes of the air-water interface on unstructured meshes. Finally, we demonstrate the phase-field solver for a practical problem of wave-structure interaction in offshore engineering using general three-dimensional unstructured meshes.

35 pages, 22 figures, 5 tables

Cited by in corpus (13)

• A conservative diffuse interface method for two-phase flows with provable boundedness properties
• An adaptive variational procedure for the conservative and positivity preserving Allen-Cahn phase-field model
• A variational interface-preserving and conservative phase-field method for the surface tension effect in two-phase flows
• Transverse flow-induced vibrations of a sphere in the proximity of a free surface: A numerical study
• A hybrid variational Allen-Cahn/ALE scheme for the coupled analysis of two-phase fluid-structure interaction
• A consistent and conservative volume distribution algorithm and its applications to multiphase flows using Phase-Field models
• An interface and geometry preserving phase-field method for fully Eulerian fluid-structure interaction
• Eulerian framework for contact between solids represented as phase fields
• Tusas: A fully implicit parallel approach for coupled phase-field equations
• Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows
• Low Order Finite Element Methods for the Navier-Stokes-Cahn-Hilliard Equations
• A robust and accurate finite element framework for cavitating flows with fluid-structure interaction
• A fractional step lattice Boltzmann model for two phase flows with large density differences