Numerical Analysis of Second Order, Fully Discrete Energy Stable Schemes for Phase Field Models of Two-Phase Incompressible Flows


In this paper, we propose several second order in time, fully discrete, linear and nonlinear numerical schemes for solving the phase field model of two-phase incompressible flows, in the framework of finite element method. The schemes are based on the second order Crank-Nicolson method for time discretization, projection method for Navier-Stokes equations, as well as several implicit-explicit treatments for phase field equations. The energy stability and unique solvability of the proposed schemes are proved. Ample numerical experiments are performed to validate the accuracy and efficiency of the proposed schemes.


Mathematics and Statistics


D. Han is partially supported by NSF DMS-1312701. Alex. Brylev is partially supported by NSF DMS-1418898. X. Yang is partially supported by NSF-DMS-1200487, NSF-DMS-1418898, AFOSR-FA9550-12-1-0178, NSFC-11471372, and NSFC-11571385. Z. Tan is partially supported by the NSFC-11571385, the special project "High performance computing" of National Key Research and Development Program (No. 2016YFB0200604), the Fundamental Research Funds for the Central Universities (15lgjc17), and Guangdong Province Key Laboratory of Computational Science at the Sun Yat-sen University.

Keywords and Phrases

Convergence of numerical methods; Incompressible flow; Navier Stokes equations; Numerical methods; Phase interfaces; Two phase flow; Energy stability; Explicit treatments; Numerical experiments; Phase field equation; Phase field models; Phase fields; Projection method; Time discretization; Finite element method; Navier-Stokes; Stability

International Standard Serial Number (ISSN)

0885-7474; 1573-7691

Document Type

Article - Journal

Document Version


File Type





© 2017 Springer Verlag, All rights reserved.

Publication Date

01 Mar 2017