Numerical Simulation for Moving Contact Line with Continuous Finite Element Schemes

In this paper, we compute a phase field (diffuse interface) model of Cahn-Hilliard type for moving contact line problems governing the motion of isothermal multiphase incompressible fluids. The generalized Navier boundary condition proposed by Qian et al. [1] is adopted here. We discretize model equations using a continuous finite element method in space and a modified midpoint scheme in time. We apply a penalty formulation to the continuity equation which may increase the stability in the pressure variable. Two kinds of immiscible fluids in a pipe and droplet displacement with a moving contact line under the effect of pressure driven shear flow are studied using a relatively coarse grid. We also derive the discrete energy law for the droplet displacement case, which is slightly different due to the boundary conditions. The accuracy and stability of the scheme are validated by examples, results and estimate order.

[1]  Hang Ding,et al.  Numerical Simulations of Flows with Moving Contact Lines , 2014 .

[2]  David Jacqmin,et al.  Contact-line dynamics of a diffuse fluid interface , 2000, Journal of Fluid Mechanics.

[3]  Jie Shen,et al.  Numerical simulations of jet pinching-off and drop formation using an energetic variational phase-field method , 2006, J. Comput. Phys..

[4]  Ping Lin,et al.  Continuous finite element schemes for a phase field model in two-layer fluid Bénard-Marangoni convection computations , 2014, Comput. Phys. Commun..

[5]  Roger T. Bonnecaze,et al.  Displacement of a two-dimensional immiscible droplet adhering to a wall in shear and pressure-driven flows , 1999, Journal of Fluid Mechanics.

[6]  Brian G. Thomas,et al.  Inclusion removal by bubble flotation in a continuous casting mold , 2004 .

[7]  J. Li,et al.  Numerical simulation of moving contact line problems using a volume-of-fluid method , 2001 .

[8]  P. Spelt,et al.  Validation and modification of asymptotic analysis of slow and rapid droplet spreading by numerical simulation , 2013, Journal of Fluid Mechanics.

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

[10]  Yi Sui,et al.  An efficient computational model for macroscale simulations of moving contact lines , 2013, J. Comput. Phys..

[11]  R. Cerro,et al.  Surface forces and inertial effects on moving contact lines , 2007 .

[12]  S. Gavrilyuk,et al.  Modeling of the moving deformed triple contact line: influence of the fluid inertia. , 2006, Journal of colloid and interface science.

[13]  Rüdiger Weiner,et al.  Half-explicit Runge-Kutta methods for semi-explicit differential-algebraic equations of index 1 , 1993 .

[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]  Rainer Mertz,et al.  Numerical simulation of bubble formation on orifice plates with a moving contact line , 2009 .

[16]  Robert Altmann,et al.  Finite Element Decomposition and Minimal Extension for Flow Equations , 2015 .

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

[18]  Yi Shi,et al.  A finite element method for the numerical solution of the coupled Cahn-Hilliard and Navier-Stokes system for moving contact line problems , 2012, J. Comput. Phys..

[19]  M. Miksis,et al.  Gas bubble with a moving contact line rising in an inclined channel at finite Reynolds number , 2004 .

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

[21]  Pierre Seppecher,et al.  Moving contact lines in the Cahn-Hilliard theory , 1996 .

[22]  Qi Wang,et al.  Energy law preserving C0 finite element schemes for phase field models in two-phase flow computations , 2011, J. Comput. Phys..

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

[24]  Ping Lin.  A Sequential Regularization Method for Time-Dependent Incompressible Navier--Stokes Equations , 1997 .

[25]  Chunfeng Zhou,et al.  Sharp-interface limit of the Cahn–Hilliard model for moving contact lines , 2010, Journal of Fluid Mechanics.

[26]  Hongkai Zhao,et al.  An augmented method for free boundary problems with moving contact lines , 2010 .

[27]  Jie Shen,et al.  An efficient moving mesh spectral method for the phase-field model of two-phase flows , 2009, J. Comput. Phys..

[28]  Chun Liu,et al.  Simulations of singularity dynamics in liquid crystal flows: A C0 finite element approach , 2006, J. Comput. Phys..

[29]  James J. Feng,et al.  Wall energy relaxation in the Cahn–Hilliard model for moving contact lines , 2009 .

[30]  Gunilla Kreiss,et al.  A conservative level set method for contact line dynamics , 2009, J. Comput. Phys..

[31]  Hui Zhang,et al.  An energy law preserving C0 finite element scheme for simulating the kinematic effects in liquid crystal dynamics , 2007, J. Comput. Phys..

[32]  M. Arnold.  Half-explicit Runge-Kutta methods with explicit stages for differential-algebraic systems of index 2 , 1998 .

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

[34]  N. Hadjiconstantinou.  Regular Article: Hybrid Atomistic–Continuum Formulations and the Moving Contact-Line Problem , 1999 .

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