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.
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
Mathematics and Statistics
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)
Article - Journal
© 2018 Springer Verlag, All rights reserved.
01 Nov 2018