Scale Dependence of Contact Line Computations

The shape and velocity of a sliding droplet are computed by solving the Navier{ Stokes equation with free interface boundary conditions. The Galerkin flnite element method is implemented in a 2D computation domain discretized using an unstructured mesh with triangular elements. The mesh is reflned recursively at the corners (contact points). The stationary sliding velocity is found to be strongly dependent on grid reflnement, which is a consequence of the contact line singularity resolved through the efiective slip across the flnite elements adjacent to the contact point. For small droplets, this dependence is well approximated by a theoretical estimates obtained using multiscale expansion and matching technique in lubrication approximation, where the corner element size is used as a microscale parameter. For larger droplets, the shape is also dependent on grid reflnement. This ques- tions the validity of numerous computations of ∞ows with moving contact line where grids are invariably much more coarse than molecular scales on which the singularity is resolved. It is suggested that extrapolation to molecular scales should be used to obtain realistic results.