A gradient stable scheme for a phase field model for the moving contact line problem
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 CahnHilliard and NavierStokes equations with the generalized Navier boundary condition [1], [2]and[4]. The nonlinear version of the scheme is semi-implicit in time and is based on a convex splitting of the CahnHilliard free energy (including the boundary energy) together with a projection method for the NavierStokes 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]and[22] for the Navier Stokes equations has the better accuracy for pressure.
phase field moving contact line navierstokes equations convex splitting cahnhilliard equation