English

Transition in a numerical model of contact line dynamics and forced dewetting

Fluid Dynamics 2018-10-30 v4 Computational Physics

Abstract

We investigate the transition to a Landau-Levich-Derjaguin film in forced dewetting using a quadtree adaptive solution to the Navier-Stokes equations with surface tension. We use a discretization of the capillary forces near the receding contact line that yields an equilibrium for a specified contact angle θΔ\theta_\Delta called the numerical contact angle. Despite the well-known contact line singularity, dynamic simulations can proceed without any explicit additional numerical procedure. We investigate angles from 1515^\circ to 110110^\circ and capillary numbers from 0.000850.00085 to 0.20.2 where the mesh size Δ\Delta is varied in the range of 0.00350.0035 to 0.060.06 of the capillary length lcl_c. To interpret the results, we use Cox's theory which involves a microscopic distance rmr_m and a microscopic angle θe\theta_e. In the numerical case, the equivalent of θe\theta_e is the angle θΔ\theta_\Delta and we find that Cox's theory also applies. We introduce the scaling factor or gauge function ϕ\phi so that rm=Δ/ϕr_m = \Delta/\phi and estimate this gauge function by comparing our numerics to Cox's theory. The comparison provides a direct assessment of the agreement of the numerics with Cox's theory and reveals a critical feature of the numerical treatment of contact line dynamics: agreement is poor at small angles while it is better at large angles. This scaling factor is shown to depend only on θΔ\theta_\Delta and the viscosity ratio qq. In the case of small θe\theta_e, we use the prediction by Eggers [Phys. Rev. Lett., vol. 93, pp 094502, 2004] of the critical capillary number for the Landau-Levich-Derjaguin forced dewetting transition. We generalize this prediction to large θe\theta_e and arbitrary qq and express the critical capillary number as a function of θe\theta_e and rmr_m. An analogy can be drawn between rmr_m and the numerical slip length.

Keywords

Cite

@article{arxiv.1703.07038,
  title  = {Transition in a numerical model of contact line dynamics and forced dewetting},
  author = {S. Afkhami and J. Buongiorno and A. Guion and S. Popinet and Y. Saade and R. Scardovelli and S. Zaleski},
  journal= {arXiv preprint arXiv:1703.07038},
  year   = {2018}
}

Comments

This version of the paper includes the corrections indicated in Ref. [1]