Stokes flows in a sessile hemispherical drop due to evaporation and surface tension gradient

本文针对具有钉扎接触线的蒸发半球形液滴，通过解析方法研究了无滑移与全滑移边界条件下由蒸发和表面张力梯度（如非均匀冷却或光敏表面活性剂引起）驱动的斯托克斯流动，揭示了无滑移条件下蒸发率与表面张力梯度的刚性关系，并重新阐释了从毛细流向马兰戈尼对流转变的临界马兰戈尼数。

要点

这篇文章探讨了一个非常有趣且贴近生活的物理现象：为什么一滴水在蒸发时，里面的液体会像有生命一样流动？

想象一下，你在一块玻璃上滴了一滴墨水或咖啡。随着水分蒸发，你会看到边缘出现了一圈深色的痕迹，这就是著名的“咖啡环效应”。这篇文章就是深入研究了这滴“小水珠”内部到底发生了什么，以及为什么有时候水流会突然变得不一样。

作者把这个问题拆解成了几个关键部分，我们可以用一些生动的比喻来理解：

1. 舞台设定：半球形的小水滴

想象这滴水是一个完美的半球形（像半个西瓜皮扣在桌子上）。

• 蒸发：就像水分子从表面“逃跑”到空气中。
• 接触线（Contact Line）：水滴边缘和桌子接触的那一圈线。在这篇文章里，作者假设这圈线是被钉死的（Pinned），也就是说，水滴在蒸发变小时，边缘不会收缩，而是像被胶水粘住一样，只能变扁。

2. 两种主要的“水流舞蹈”

水滴内部的液体流动主要有两种模式，就像两种不同的舞蹈：

A. 德根流（Deegan Flow）：补位舞者

• 起因：水滴边缘蒸发得最快（因为边缘薄，空气接触多）。为了补充边缘流失的水分，底部的液体就会像补位一样，从中心流向边缘。
• 结果：这就是造成“咖啡环”的罪魁祸首。液体带着溶质（比如咖啡粉）冲向边缘，水分蒸发后，粉末就留在了那里。
• 比喻：就像一场接力赛，边缘的选手跑得太快把水耗尽了，后方的选手（中心液体）必须拼命跑过去填补空缺。

B. 马兰戈尼流（Marangoni Flow）：表面张力舞者

• 起因：如果水滴表面温度不均匀（比如一边冷一边热），或者表面张力不一样，液体就会从“张力低”的地方流向“张力高”的地方。
• 结果：这会产生一种旋转的涡流，像是一个漩涡，把液体从边缘拉回中心，或者形成复杂的循环。
• 比喻：想象你在浴缸里滴入一滴肥皂水，肥皂水周围的表面张力变小，水就会把肥皂水推开。这就是马兰戈尼效应。

3. 核心发现：地板是“粘”还是“滑”？

这篇文章最精彩的地方在于它发现：水滴底部的地板（基底）是“粘人”的还是“打滑”的，决定了这两种舞蹈是混在一起跳，还是分开跳。

情况一：地板很“粘”（无滑移条件，No-slip）

• 设定：液体分子紧紧贴在地板上，完全不动。
• 现象：在这种条件下，“补位舞”（德根流）和“漩涡舞”（马兰戈尼流）是锁死的，分不开！
• 比喻：就像你穿着粘满胶水的鞋子在地板上走。你想往边缘跑（补位），但地板的摩擦力强迫你必须同时产生一种旋转（马兰戈尼效应）。
• 结论：作者发现，如果地板很粘，那么蒸发速度和表面张力梯度之间有一个死板的数学关系。这意味着，只要水滴在蒸发，它内部就必然伴随着某种程度的马兰戈尼对流。你无法只让“补位舞”单独存在。这就像是你无法只让一个人既在跑步又在原地转圈，除非地板允许他打滑。

情况二：地板很“滑”（滑移条件，Slip）

• 设定：液体在地板上可以像溜冰一样滑动。
• 现象：这时，两种舞蹈彻底分家了！
• 比喻：就像你穿上了溜冰鞋。你可以选择只往边缘跑（纯粹的补位流，形成咖啡环），也可以只转圈圈（纯粹的马兰戈尼流）。它们互不干扰，可以独立存在。
• 结论：只有在允许滑动的情况下，我们才能在理论上把“咖啡环效应”和“马兰戈尼对流”完全分开研究。

4. 为什么这很重要？（临界点）

文章提出了一个全新的视角来理解临界马兰戈尼数（一个判断对流是否发生的指标）：

• 传统观点：我们通常认为，只有当温度差大到一定程度，马兰戈尼流才会突然“爆发”。
• 本文观点：其实，只要地板是“粘”的，马兰戈尼流一直存在，只是它被蒸发流“伪装”了，混在一起了。
• 转折点：当蒸发非常剧烈，或者液体变得很热，导致靠近地板的液体受到的剪切力（摩擦力）太大时，原本“粘”在地板上的液体可能会突然**“打滑”**。
  • 一旦打滑，原本混在一起的两种流就会分离。
  • 这时候，系统可能会从一种混乱的混合状态，突然转变为清晰的、独立的马兰戈尼对流状态。

总结

这篇文章就像是在研究一滴水里的“交通系统”：

1. 蒸发是交通的驱动力。
2. 地板的摩擦力（粘还是滑）决定了交通规则。
3. 如果地板太粘，所有的车（液体）都被迫走一条混合路线，你分不清哪辆车是为了补位，哪辆车是为了旋转。
4. 如果地板变滑了，车流就会分道扬镳，有的专门补位，有的专门旋转。

对实验人员的启示：
如果你在做实验，发现水滴里的流动模式突然变了，不要只盯着温度看。也许是因为蒸发太快，导致液体和桌子的摩擦力变了（从“粘”变成了“滑”），从而彻底改变了内部的流动模式。

这就解释了为什么有时候微小的环境变化（比如温度升高一点点），会导致水滴内部发生剧烈的、非线性的流动变化。

技术摘要

这是一份关于论文《Stokes flows in a sessile hemispherical drop due to evaporation and surface tension gradient》（蒸发和表面张力梯度引起的半球形液滴中的斯托克斯流动）的详细技术总结。

1. 研究问题 (Problem Statement)

本文旨在研究缓慢蒸发的半球形驻留液滴（sessile droplet）内部的粘性流体动力学流动。核心问题在于阐明蒸发引起的毛细补偿流（即“咖啡环效应”中的 Deegan 流）与表面张力梯度引起的马兰戈尼流（Marangoni flow）之间的相互作用机制。

研究特别关注液滴与基底之间的边界条件对流动结构的影响：

• 无滑移条件 (No-slip)：液滴在基底处速度为零。
• 全滑移条件 (Full-slip)：允许液滴在基底上滑动。

作者试图解决一个理论难题：在无滑移条件下，蒸发流和热毛细流是否本质上是分离的？还是存在某种刚性耦合？

2. 方法论 (Methodology)

• 物理模型：

  • 假设液滴为半球形（接触角 $\pi/2$），接触线固定（Pinned contact line）。
  • 流体为不可压缩粘性流体，流动处于低雷诺数状态（Stokes 流/蠕动流）。
  • 考虑蒸发导致的体积损失以及由此产生的表面温度梯度和表面张力梯度。

• 数学工具：

  • 使用球坐标系下的轴对称斯托克斯方程（Stokes equations）。
  • 利用勒让德多项式 (Legendre polynomials) 展开速度场、流函数和压力场，以获得解析解。
  • 分别处理两种边界条件情况：
    1. 无滑移情况：基底处径向和切向速度均为零。
    2. 滑移情况：基底处仅满足不可渗透条件，允许切向滑动。

• 热传递分析：

  • 分析了低佩克莱特数（$Pe \ll 1$，分子热传导主导）和高佩克莱特数（$Pe \gg 1$，对流主导）两种热传递机制。
  • 推导了与特定流动模式相对应的特征温度分布和表面张力梯度。

3. 主要贡献与关键发现 (Key Contributions & Results)

A. 无滑移条件下的刚性耦合 (Rigid Coupling under No-Slip)

• 发现：在无滑移边界条件下，蒸发引起的 Deegan 流（补偿流）与马兰戈尼流不可分离。
• 机制：为了满足基底处的无滑移条件（速度为零），蒸发引起的径向流动必须伴随着特定的切向速度分布，这反过来要求存在一个特定的表面张力梯度。
• 结果：
  • 即使没有外部施加的温度梯度，蒸发本身也会诱导出一个等效的“特征温度差”和表面张力梯度。
  • 计算得出的有效马兰戈尼数（$Ma$）虽然数值很小（对于水约为 0.0017），但在数学上是非零的。
  • 结论：在无滑移条件下，不存在纯粹的“咖啡环流”；任何蒸发流都必然伴随着马兰戈尼效应。两者是刚性耦合的，无法通过简单的参数调整将它们分开。

B. 滑移条件下的流动分离 (Flow Separation under Slip)

• 发现：当允许液滴在基底上滑动（滑移条件）时，蒸发流和马兰戈尼流可以完全分离。
• 机制：滑移条件消除了基底处对切向速度的强制约束，使得蒸发通量（$J$）和表面张力梯度（$\beta$）成为两个独立的参数。
• 结果：
  • 流动可以分解为两部分：
    1. 纯补偿流 (Deegan flow)：仅由蒸发驱动，对应于零表面张力梯度（$Ma=0$）。
    2. 纯马兰戈尼流：仅由表面张力梯度驱动。
  • 流场可视化显示，滑移条件下的流线比无滑移条件（无滑移条件下流线呈现分形/束状结构）更加简单清晰。

C. 对临界马兰戈尼数的新视角 (New Perspective on Critical Marangoni Number)

• 传统观点认为马兰戈尼对流有一个临界阈值，超过该阈值流动才会发生。
• 本文提出：在无滑移条件下，由于蒸发流本身已经包含了非零的马兰戈尼效应，因此“从毛细流到马兰戈尼对流”的过渡并非简单的阈值跨越，而是边界条件性质的改变。
• 物理图景：随着蒸发速率增加，热负荷和粘性剪切应力增大，可能导致基底处的无滑移条件失效（转变为部分或全滑移）。这种转变使得原本耦合的流动解耦，从而出现独立的、强烈的马兰戈尼对流。

4. 结果可视化与流场特征

• 无滑移流场 (Fig. 2)：流线呈现复杂的束状和分形特征，表明流动结构高度依赖于蒸发与表面张力梯度的耦合。
• 滑移流场 (Fig. 3)：
  • 图 3(a) 展示了纯补偿流，流线平滑。
  • 图 3(b) 展示了纯马兰戈尼流（环流）。
  • 两者在滑移条件下可以独立存在和叠加。

5. 研究意义 (Significance)

1. 理论突破：挑战了将 Deegan 流和马兰戈尼流视为独立现象的传统分类。证明了在严格的无滑移边界条件下，两者是数学上不可分割的。
2. 实验指导：为实验物理学家提供了新的研究方向。在研究蒸发液滴时，必须高度关注液 - 固界面的边界条件（是粘附还是滑移）。
3. 相变机制解释：解释了为什么在某些系统中，随着蒸发加剧，流动模式会发生突变。这种突变可能不是由于马兰戈尼数超过了某个绝对临界值，而是由于高剪切应力导致基底边界条件从“无滑移”转变为“滑移”，从而释放了马兰戈尼对流。
4. 应用价值：对微流体、喷墨打印、自组装技术、蒸发光刻及薄膜涂层等领域具有指导意义，特别是在需要精确控制液滴内部混合和溶质沉积（如咖啡环效应）的场景中。

总结

该论文通过严格的解析推导，揭示了蒸发液滴中流动机制的深层物理本质：边界条件决定了流动模式的耦合程度。无滑移导致蒸发流与马兰戈尼流强制耦合，而滑移条件则允许两者解耦。这一发现为理解蒸发液滴中的复杂对流现象提供了全新的理论框架。