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

Abstract 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.

[1]  Pao-Hsiung Chiu,et al.  A conservative phase field method for solving incompressible two-phase flows , 2011, J. Comput. Phys..

[2]  Anna C Balazs,et al.  Convection-driven pattern formation in phase-separating binary fluids. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[4]  L. Antanovskii.  A phase field model of capillarity , 1995 .

[5]  Chang Shu,et al.  Diffuse interface model for incompressible two-phase flows with large density ratios , 2007, J. Comput. Phys..

[6]  Jie Shen,et al.  A Phase-Field Model and Its Numerical Approximation for Two-Phase Incompressible Flows with Different Densities and Viscosities , 2010, SIAM J. Sci. Comput..

[7]  J. Lowengrub,et al.  Quasi–incompressible Cahn–Hilliard fluids and topological transitions , 1998, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[8]  Abner J. Salgado,et al.  A splitting method for incompressible flows with variable density based on a pressure Poisson equation , 2009, J. Comput. Phys..

[9]  L. Quartapelle,et al.  A projection FEM for variable density incompressible flows , 2000 .

[10]  Yuying Yan,et al.  A lattice Boltzmann method for incompressible two-phase flows on partial wetting surface with large density ratio , 2007, J. Comput. Phys..

[11]  P. Sheng,et al.  A variational approach to moving contact line hydrodynamics , 2006, Journal of Fluid Mechanics.

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

[13]  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..

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

[15]  Jie Shen,et al.  Gauge-Uzawa methods for incompressible flows with variable density , 2007, J. Comput. Phys..

[16]  Xiaofeng Yang,et al.  Mass and Volume Conservation in Phase Field Models for Binary Fluids , 2013 .

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

[18]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .