Numerical Study of Drop Motion on a Surface with Wettability Gradient and Contact Angle Hysteresis

Abstract

In this work, the motion of a 2-D drop on a surface with given wettability gradient is studied numerically by a hybrid lattice-Boltzmann finite-difference method using the multiple-relaxation-time collision model. We incorporate the geometric wetting boundary condition that allows accurate implementation of a contact angle hysteresis model. The method is first validated through three benchmark tests, including the layered Poiseuille flow with a viscosity contrast, the motion of a liquid column in a channel with specified wettability gradient and the force balance for a static drop attached to a surface with hysteresis subject to a body force. Then, simulations of a drop on a wall with given wettability gradient are performed under different conditions. The effects of the Reynolds number, the viscosity ratio, the wettability gradient, as well as the contact angle hysteresis on the drop motion are investigated in detail. It is found that the capillary number of the drop in steady state is significantly affected by the viscosity ratio, the magnitudes of the wettability gradient and the contact angle hysteresis, whereas it only shows very weak dependence on the Reynolds number