A Second Order in Time, Decoupled, Unconditionally Stable Numerical Scheme for the Cahn-Hilliard-Darcy System


We propose a novel second order in time, fully decoupled and unconditionally stable numerical scheme for solving the Cahn-Hilliard-Darcy system which models two-phase flow in porous medium or in a Hele-Shaw cell. The scheme is based on the ideas of second order convex-splitting for the Cahn-Hilliard equation and pressure-correction for the Darcy equation. The computation of order parameter, pressure and velocity is completely decoupled in our scheme. We show that the scheme is uniquely solvable, unconditionally energy stable and mass-conservative. Ample numerical results are presented to gauge the efficiency and robustness of our scheme.


Mathematics and Statistics

Research Center/Lab(s)

Center for High Performance Computing Research


The work of X. Wang is supported in part by DMS 1715504 and grants from Fudan University. The work of D. Han is supported by a seed fund from the Material Research Center at Missouri University of Science and Technology. The authors wish to thank Wenbin Chen, Mike Jolly and Jie Shen for helpful discussions.

Keywords and Phrases

Porous materials; Decoupling; Energy law; Numerical results; Numerical scheme; Order parameter; Pressure correction; Unconditional stability; Unconditionally stable; Two phase flow; Cahn-Hilliard-Darcy; Pressure-correction

International Standard Serial Number (ISSN)

0885-7474; 1573-7691

Document Type

Article - Journal

Document Version


File Type





© 2018 Springer Verlag, All rights reserved.

Publication Date

01 Nov 2018