中文

接触线动力学与强制去湿数值模型中的转变

流体动力学 2018-10-30 v4 计算物理

摘要

我们利用带有表面张力的 Navier-Stokes 方程的四叉树自适应解,研究了强制去湿中向 Landau-Levich-Derjaguin 薄膜的转变。我们对后退接触线附近的毛细力进行离散化,该离散化针对指定接触角 θΔ\theta_\Delta(称为数值接触角)产生平衡。尽管存在众所周知的接触线奇异性,动态模拟仍可在无任何显式额外数值过程的情况下进行。我们研究了从 1515^\circ110110^\circ 的角度以及从 0.000850.000850.20.2 的毛细数,其中网格尺寸 Δ\Delta 在毛细长度 lcl_c0.00350.00350.060.06 范围内变化。为解释结果,我们使用涉及微观距离 rmr_m 和微观角 θe\theta_e 的 Cox 理论。在数值情形下,θe\theta_e 的等价量是角 θΔ\theta_\Delta,并且我们发现 Cox 理论同样适用。我们引入缩放因子或规范函数 ϕ\phi 使得 rm=Δ/ϕr_m = \Delta/\phi,并通过将我们的数值结果与 Cox 理论比较来估计该规范函数。该比较提供了数值结果与 Cox 理论一致性的直接评估,并揭示了接触线动力学数值处理的一个关键特征:在小角度下一致性较差,而在大角度下较好。该缩放因子被证明仅依赖于 θΔ\theta_\Delta 和黏度比 qq。在小 θe\theta_e 情形下,我们使用 Eggers [Phys. Rev. Lett., vol. 93, pp 094502, 2004] 对 Landau-Levich-Derjaguin 强制去湿转变临界毛细数的预测。我们将该预测推广到大 θe\theta_e 和任意 qq,并将临界毛细数表示为 θe\theta_ermr_m 的函数。可以在 rmr_m 与数值滑移长度之间进行类比。

关键词

引用

@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}
}

备注

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