A gradient stable scheme for a phase field model for the moving contact line problem

Min Gao, Xiao-Ping Wang

Research output: Contribution to journalArticlepeer-review

75 Scopus citations


In this paper, an efficient numerical scheme is designed for a phase field model for the moving contact line problem, which consists of a coupled system of the Cahn-Hilliard and Navier-Stokes equations with the generalized Navier boundary condition [1,2,4]. The nonlinear version of the scheme is semi-implicit in time and is based on a convex splitting of the Cahn-Hilliard free energy (including the boundary energy) together with a projection method for the Navier-Stokes equations. We show, under certain conditions, the scheme has the total energy decaying property and is unconditionally stable. The linearized scheme is easy to implement and introduces only mild CFL time constraint. Numerical tests are carried out to verify the accuracy and stability of the scheme. The behavior of the solution near the contact line is examined. It is verified that, when the interface intersects with the boundary, the consistent splitting scheme [21,22] for the Navier Stokes equations has the better accuracy for pressure. © 2011 Elsevier Inc.
Original languageEnglish (US)
Pages (from-to)1372-1386
Number of pages15
JournalJournal of Computational Physics
Issue number4
StatePublished - Feb 2012
Externally publishedYes

Bibliographical note

KAUST Repository Item: Exported on 2020-10-01
Acknowledged KAUST grant number(s): SA-C0040/UK-C0016
Acknowledgements: This publication was based on work supported in part by Award No. SA-C0040/UK-C0016, made by King Abdullah University of Science and Technology (KAUST), Hong Kong RGC-GRF Grants 603107, 604209 and the National Basic Research Program Project of China under project 2009CB623200. Min Gao is a Ph.D. student under HKUST-Shanghai Jiaotong University collaboration program.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.


Dive into the research topics of 'A gradient stable scheme for a phase field model for the moving contact line problem'. Together they form a unique fingerprint.

Cite this