A Computational Multiscale Model for Contact Line Dynamics

Abstract: The conventional no-slip boundary condition leads to a non-integrable stress singularity at a moving contact line. This makes numerical simulations challenging, especially when capillary effects are essential for the dynamics of the flow. This paper presents a new boundary methodology, suitable for numerical simulation of flow of two immiscible and incompressible fluids in the presence of moving contact points. The methodology is based on combining a relation between the apparent contact angle and the contact point velocity with the similarity solution for Stokes flow at a planar interface. The relation between angle and velocity can be determined by theoretical arguments, or from simulations using a more detailed model. The approach here uses the phase field model in a micro domain, with physically relevant parameters for molecular diffusion and interface thickness. The methodology is used to formulate a new boundary condition for the velocity. Numerical results illustrate the usefulness.