A Second Order in Time, Decoupled, Unconditionally Stable Numerical Scheme for the Cahn-Hilliard-Darcy System
Abstract
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.
Recommended Citation
D. Han and X. Wang, "A Second Order in Time, Decoupled, Unconditionally Stable Numerical Scheme for the Cahn-Hilliard-Darcy System," Journal of Scientific Computing, vol. 77, no. 2, pp. 1210 - 1233, Springer Verlag, Nov 2018.
The definitive version is available at https://doi.org/10.1007/s10915-018-0748-0
Department(s)
Mathematics and Statistics
Research Center/Lab(s)
Center for High Performance Computing Research
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
Citation
File Type
text
Language(s)
English
Rights
© 2018 Springer Verlag, All rights reserved.
Publication Date
01 Nov 2018
Comments
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.