A conservative diffuse interface method for two-phase flows with provable boundedness properties

Abstract

Central finite difference schemes have long been avoided in the context of two-phase flows for the advection of the phase indicator function due to numerical overshoots and undershoots associated with their dispersion errors. We will show however, for an incompressible flow, in the context of a specific diffuse interface model, one can maintain the boundedness of the phase field while also taking advantage of the low cost and ease of implementation of central differences to construct a non-dissipative discretization scheme for the advective terms. This is made possible by combining the advection and reinitialization steps of a conservative level set scheme introduced by Olsson and Kreiss [J. Comput. Phys., 210, 225 (2005)] to form a phase field equation similar to that of Chiu and Lin [J. Comput. Phys., 230, 185 (2011)]. Instead of resorting to specialized upwind methods as in these articles, we prove that the boundedness of the phase field is guaranteed for certain choices of the free parameters ($\epsilon$ and $\gamma$) for a specific central difference scheme that we propose. The proposed discretely conservative and bounded phase field equation, which is free of any reinitialization or mass redistribution, possesses desirable properties that can be leveraged in the coupled finite difference discretization of the two-phase momentum equation. Additionally, as compared to the state-of-the-art conservative and bounded two-phase flow methods, the proposed method boasts competitive accuracy-vs-cost trade-off, small memory requirements, ease of implementation and excellent parallelizability, providing a viable alternative for realistic two-phase flow calculations.

Under review for publication in Journal of Computational Physics

References in corpus (1)

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

Cited by in corpus (17)

• Comparison between the diffuse interface and volume of fluid methods for simulating two-phase flows
• A conservative diffuse-interface method for compressible two-phase flows
• Accurate conservative phase-field method for simulation of two-phase flows
• An energy-preserving level set method for multiphase flows
• Assessment of an energy-based surface tension model for simulation of two-phase flows using second-order phase field methods
• A consistent and conservative volume distribution algorithm and its applications to multiphase flows using Phase-Field models
• An accurate lattice Boltzmann method for interface capturing
• Modeling of non-equilibrium effects in intermittency region between two phases
• Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows
• High Fidelity Simulations of Micro-Bubble Shedding from Retracting Thin Gas Films in the Context of Liquid-Liquid Impact
• Discrete exterior calculus discretization of two-phase incompressible Navier-Stokes equations with a conservative phase field method
• A general pressure equation based method for incompressible two-phase flows
• On the Energy Analysis of Two-phase Flows Simulated with the Diffuse Interface Method
• An Energy Stable Nonlinear Incompressible Multi-Phase Flow Formulation
• Enhanced numerical approaches for modeling insoluble surfactants in two-phase flows with the diffuse-interface method
• Grid-agnostic volume of fluid approach with interface sharpening and surface tension for compressible multiphase flows
• Topological changes due to non-equilibrium effects by means of the statistical model of two-phase flow