Moving contact line dynamics: from diffuse to sharp interfaces

• Published in: Journal of fluid mechanics Vol. 788; pp. 209 - 227
• Main Authors: Kusumaatmaja, H., Hemingway, E. J., Fielding, S. M.
• Format: Journal Article
• Language: English
• Published: Cambridge, UK Cambridge University Press 10.02.2016
• Subjects: Fluid mechanics; Interfaces; Interfacial bonding; multiphase flow; contact lines; interfacial flows (free surface)
• ISSN: 0022-1120; 1469-7645
• DOI: 10.1017/jfm.2015.697

We reconcile two scaling laws that have been proposed in the literature for the slip length associated with a moving contact line in diffuse interface models, by demonstrating each to apply in a different regime of the ratio of the microscopic interfacial width $l$ and the macroscopic diffusive length $l_{D}=(M{\it\eta})^{1/2}$ , where ${\it\eta}$ is the fluid viscosity and $M$ the mobility governing intermolecular diffusion. For small $l_{D}/l$ we find a diffuse interface regime in which the slip length scales as ${\it\xi}\sim (l_{D}l)^{1/2}$ . For larger $l_{D}/l>1$ we find a sharp interface regime in which the slip length depends only on the diffusive length, ${\it\xi}\sim l_{D}\sim (M{\it\eta})^{1/2}$ , and therefore only on the macroscopic variables ${\it\eta}$ and $M$ , independent of the microscopic interfacial width $l$ . We also give evidence that modifying the microscopic interfacial terms in the model’s free energy functional appears to affect the value of the slip length only in the diffuse interface regime, consistent with the slip length depending only on macroscopic variables in the sharp interface regime. Finally, we demonstrate the dependence of the dynamic contact angle on the capillary number to be in excellent agreement with the theoretical prediction of Cox (J. Fluid Mech., vol. 168, 1986, p. 169), provided we allow the slip length to be rescaled by a dimensionless prefactor. This prefactor appears to converge to unity in the sharp interface limit, but is smaller in the diffuse interface limit. The excellent agreement of results obtained using three independent numerical methods, across several decades of the relevant dimensionless variables, demonstrates our findings to be free of numerical artefacts.