Discontinuous Finite Volume Element Method for a Coupled Navier-Stokes-Cahn-Hilliard Phase Field Model


In this paper, we propose a discontinuous finite volume element method to solve a phase field model for two immiscible incompressible fluids. In this finite volume element scheme, discontinuous linear finite element basis functions are used to approximate the velocity, phase function, and chemical potential while piecewise constants are used to approximate the pressure. This numerical method is efficient, optimally convergent, conserving the mass, convenient to implement, flexible for mesh refinement, and easy to handle complex geometries with different types of boundary conditions. We rigorously prove the mass conservation property and the discrete energy dissipation for the proposed fully discrete discontinuous finite volume element scheme. Using numerical tests, we verify the accuracy, confirm the mass conservation and the energy law, test the influence of surface tension and small density variations, and simulate the driven cavity, the Rayleigh-Taylor instability.


Mathematics and Statistics

Research Center/Lab(s)

Center for High Performance Computing Research

Keywords and Phrases

Discontinuous finite volume element methods; Discrete energy dissipation; Navier-Stokes-Cahn-Hilliard equation; Phase field model

International Standard Serial Number (ISSN)

1019-7168; 1572-9044

Document Type

Article - Journal

Document Version


File Type





© 2020 Springer, All rights reserved.

Publication Date

01 Apr 2020