Numerical study on the stick-slip motion of contact line moving on heterogeneous surfaces

Abstract

We present a numerical study of a moving contact line (CL) crossing the intersecting region of hydrophilic and hydrophobic patterns on a solid wall using lattice Boltzmann methods (LBMs). To capture the interface between the two phases properly, we applied a phase field model coupled with the LBM. The evolutions of the CL velocity, dynamic contact angle, and apparent contact angle are analyzed for the so-called “stick” and “slip” processes. In the two processes, the evolution of the quantities follows different rules shortly after the initial quick transition, which is probably caused by finite interfacial thickness or non-equilibrium effects. For the stick process, the CL is almost fixed and energy is extracted from the main flow to rebuild the meniscus’ profile. The evolution of the meniscus is mainly governed by mass conservation. The CL is depinned after the apparent contact angle surpasses the dynamic one, which implies that the interfacial segment in the vicinity of contact line is bended. For the slip process, the quantities evolve with features of relaxation. In the microscopic scale, the velocity of the CL depends on the balance between unbalanced Young’s capillary force and viscous drag. To predict the apparent contact angle evolution, a model following the dynamics of an overdamped spring-mass system is proposed. Our results also show that the capillary flows in a channel with heterogeneous wall can be described generally with the Poiseuille flow superimposed by the above transient one.