An efficient scheme for a phase field model for the moving contact line problem with variable density and viscosity

Min Gao, Xiao-Ping Wang

Research output: Contribution to journalArticlepeer-review

50 Scopus citations


In this paper, we develop an efficient numerical method for the two phase moving contact line problem with variable density, viscosity, and slip length. The physical model is based on a phase field approach, which consists of a coupled system of the Cahn-Hilliard and Navier-Stokes equations with the generalized Navier boundary condition [1,2,5]. To overcome the difficulties due to large density and viscosity ratio, the Navier-Stokes equations are solved by a splitting method based on a pressure Poisson equation [11], while the Cahn-Hilliard equation is solved by a convex splitting method. We show that the method is stable under certain conditions. The linearized schemes are easy to implement and introduce only mild CFL time constraint. Numerical tests are carried out to verify the accuracy, stability and efficiency of the schemes. The method allows us to simulate the interface problems with extremely small interface thickness. Three dimensional simulations are included to validate the efficiency of the method. © 2014 Elsevier Inc.
Original languageEnglish (US)
Pages (from-to)704-718
Number of pages15
JournalJournal of Computational Physics
StatePublished - Sep 2014
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), the Hong Kong RGC-GRF Grants 605311, 605513 and NNSF of China Grant 91230102.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.


Dive into the research topics of 'An efficient scheme for a phase field model for the moving contact line problem with variable density and viscosity'. Together they form a unique fingerprint.

Cite this