Boundary element simulations of dynamic wetting with a mesoscale contact line model

It is known that numerical simulations of moving contact lines are challenging owing to the fact that multiple scales are inherently involved. In this paper, we propose an efficient boundary element method for numerical simulations of dynamic wetting/dewetting. The flow domain is truncated in a mesoscopic scale, where boundary conditions resulted from a wedge flow and the asymptotic theory of the intermediate region are imposed. This procedure avoids the high resolution near the contact line in full-scale simulations and hence significantly reduces the computational cost. Numerical tests for dip coating problem show that the meniscus profiles and slopes produced by the proposed method agree well with high-resolution full-scale simulations as well as the local asymptotic theory.