A finite element method for the numerical solution of the coupled Cahn-Hilliard and Navier-Stokes system for moving contact line problems

In this paper, a semi-implicit finite element method is presented for the coupled Cahn-Hilliard and Navier-Stokes equations with the generalized Navier boundary condition for the moving contact line problems. In our method, the system is solved in a decoupled way. For the Cahn-Hilliard equations, a convex splitting scheme is used along with a P1-P1 finite element discretization. The scheme is unconditionally stable. A linearized semi-implicit P2-P0 mixed finite element method is employed to solve the Navier-Stokes equations. With our method, the generalized Navier boundary condition is extended to handle the moving contact line problems with complex boundary in a very natural way. The efficiency and capacity of the present method are well demonstrated with several numerical examples.

[1]  Ping Sheng,et al.  Molecular scale contact line hydrodynamics of immiscible flows. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  Qinjun Kang,et al.  Displacement of a two-dimensional immiscible droplet in a channel , 2002 .

[3]  Ping Sheng,et al.  Moving contact line over undulating surfaces , 2006 .

[4]  Ping Sheng,et al.  Moving contact line on chemically patterned surfaces , 2008, Journal of Fluid Mechanics.

[5]  Andreas Prohl,et al.  Error analysis of a mixed finite element method for the Cahn-Hilliard equation , 2004, Numerische Mathematik.

[6]  S. H. Davis,et al.  On the motion of a fluid-fluid interface along a solid surface , 1974, Journal of Fluid Mechanics.

[7]  Chunfeng Zhou,et al.  Spontaneous shrinkage of drops and mass conservation in phase-field simulations , 2007, J. Comput. Phys..

[8]  P. Sheng,et al.  Molecular Hydrodynamics of the Moving Contact Line in Two-Phase Immiscible Flows , 2005, cond-mat/0510403.

[9]  Roland Glowinski,et al.  A least-squares/finite element method for the numerical solution of the Navier-Stokes-Cahn-Hilliard system modeling the motion of the contact line , 2011, J. Comput. Phys..

[10]  Per-Olof Persson,et al.  A Simple Mesh Generator in MATLAB , 2004, SIAM Rev..

[11]  Vivette Girault,et al.  Finite Element Methods for Navier-Stokes Equations - Theory and Algorithms , 1986, Springer Series in Computational Mathematics.

[12]  Xiaobing Feng,et al.  Fully Discrete Finite Element Approximations of the Navier-Stokes-Cahn-Hilliard Diffuse Interface Model for Two-Phase Fluid Flows , 2006, SIAM J. Numer. Anal..

[13]  J.-F. Gerbeau,et al.  Generalized Navier boundary condition and geometric conservation law for surface tension , 2008, 0804.1563.

[14]  Min Gao,et al.  A gradient stable scheme for a phase field model for the moving contact line problem , 2012, J. Comput. Phys..

[15]  Timothy A. Davis,et al.  A column pre-ordering strategy for the unsymmetric-pattern multifrontal method , 2004, TOMS.

[16]  P. Spelt.  Shear flow past two-dimensional droplets pinned or moving on an adhering channel wall at moderate Reynolds numbers: a numerical study , 2006, Journal of Fluid Mechanics.