In this paper, we consider the numerical approximation for a phase field model of the coupled two-phase free flow and two-phase porous media flow. This model consists of Cahn—Hilliard—Navier—Stokes equations in the free flow region and Cahn—Hilliard—Darcy equations in the porous media region that are coupled by seven interface conditions. The coupled system is decoupled based on the interface conditions and the solution values on the interface from the previous time step. A fully discretized scheme with finite elements for the spatial discretization is developed to solve the decoupled system. In order to deal with the difficulties arising from the interface conditions, the decoupled scheme needs to be constructed appropriately for the interface terms, and a modified discrete energy is introduced with an interface component. Furthermore, the scheme is linearized and energy stable. Hence, at each time step one need only solve a linear elliptic system for each of the two decoupled equations. Stability of the model and the proposed method is rigorously proved. Numerical experiments are presented to illustrate the features of the proposed numerical method and verify the theoretical conclusions. © 2018 Society for Industrial and Applied Mathematics.
Y. Gao et al., "Decoupled, Linear, and Energy Stable Finite Element Method for the Cahn-Hilliard-Navier-Stokes-Darcy Phase Field Model," SIAM Journal on Scientific Computing, vol. 40, no. 1, pp. B110 - B137, Society for Industrial and Applied Mathematics (SIAM), Jan 2018.
The definitive version is available at https://doi.org/10.1137/16M1100885
Mathematics and Statistics
Center for High Performance Computing Research
Keywords and Phrases
Cahn-Hilliard-Navier-Stokes-Darcy; Diffuse interface; Energy stability; Finite element method
International Standard Serial Number (ISSN)
Article - Journal
© 2018 Society for Industrial and Applied Mathematics (SIAM), All rights reserved.
01 Jan 2018