第17章：多相流基础

当你喝可乐时看到上升的气泡，当你在雨中开车时看到挡风玻璃上的水滴，或者当你使用喷雾器时观察到细小的液滴云雾——这些都是多相流动的日常例子。多相流动无处不在，从内燃机的燃油喷射到血管中的血液流动，从化工反应器中的气液混合到海浪的破碎。本章将帮助你建立对这些复杂现象的物理直觉，掌握工程实践中的经验法则。

17.1 气泡与液滴动力学

17.1.1 单个气泡的上升

想象一个气泡在静止的水中上升。它的运动速度由什么决定？这是一个看似简单却蕴含丰富物理的问题。

力的平衡

气泡受到三个主要力的作用：

  • 浮力：$F_b = \rho_l V g$（向上）
  • 重力：$F_g = \rho_g V g$（向下）
  • 阻力：$F_d = \frac{1}{2} C_D \rho_l A V_t^2$（向下）

其中 $V$ 是气泡体积，$A$ 是迎流面积，$V_t$ 是终端速度。

达到终端速度时，力平衡给出： $$V_t = \sqrt{\frac{4gd(\rho_l - \rho_g)}{3C_D\rho_l}}$$ 对于小气泡（$Re < 1$），阻力系数 $C_D = 24/Re$，得到斯托克斯上升速度： $$V_t = \frac{gd^2(\rho_l - \rho_g)}{18\mu_l}$$ 形状演化

气泡的形状随尺寸变化呈现明显的规律：

直径 < 1mm：     球形（表面张力主导）
         O

1-6mm：         椭球形（开始变形）
        ___
       (   )
        ‾‾‾

6-20mm：        球帽形（形状振荡）
       _____
      (     )
       ‾‾‾‾‾

> 20mm：        不规则/破碎
      ～～～
     ～   ～
      ～～～

经验法则：

  • 直径 < 1mm 的气泡：上升速度 ∝ d²
  • 直径 1-10mm 的气泡：上升速度 ∝ d^0.5
  • 直径 > 10mm 的气泡：上升速度基本恒定 ≈ 25 cm/s

17.1.2 液滴的运动与变形

雨滴下落是最常见的液滴运动。你可能认为雨滴是泪滴形的，但实际上：

实际雨滴形状：

小雨滴（< 2mm）：    球形
                    O

中等雨滴（2-5mm）：  扁球形
                   ___
                  (   )
                   ‾‾‾

大雨滴（> 5mm）：    汉堡形（底部凹陷）
                  _____
                 (_   _)
                  ‾‾‾‾‾

Weber数的作用

液滴变形程度由Weber数决定： $$We = \frac{\rho V^2 d}{\sigma}$$

  • $We < 1$：表面张力主导，保持球形
  • $1 < We < 10$：开始变形
  • $We > 10$：显著变形
  • $We > 100$：破碎

临界Weber数：

  • 振荡开始：$We_{crit} ≈ 4$
  • 袋状破碎：$We_{crit} ≈ 10-20$
  • 剪切破碎：$We_{crit} ≈ 100$

17.1.3 气泡/液滴的相互作用

聚并（Coalescence）

两个气泡/液滴接近时，中间液膜的排出决定是否聚并：

接近 → 液膜变薄 → 临界厚度 → 破裂 → 聚并
 O O     O|O      O:O      OO      O

聚并时间尺度： $$t_{coal} \sim \frac{\mu R}{\sigma}$$ 破碎（Breakup）

在湍流中，涡旋可以撕裂气泡/液滴：

临界条件（Kolmogorov-Hinze理论）： $$d_{max} \sim \left(\frac{\sigma}{\rho}\right)^{0.6} \varepsilon^{-0.4}$$ 其中 $\varepsilon$ 是湍流耗散率。

17.2 表面张力效应

17.2.1 毛细现象的物理本质

表面张力源于分子间作用力的不平衡。在液体内部，分子受到各向同性的吸引力；而在界面处，分子受力不平衡，产生指向液体内部的净力。

Young-Laplace方程

弯曲界面两侧的压力差： $$\Delta p = \sigma \left(\frac{1}{R_1} + \frac{1}{R_2}\right)$$ 对于球形界面：$\Delta p = \frac{2\sigma}{R}$

毛细上升

细管中的液面上升高度： $$h = \frac{2\sigma \cos\theta}{\rho g r}$$ 其中 $\theta$ 是接触角，$r$ 是管径。

日常应用：

  • 毛巾吸水：纤维间的毛细作用
  • 植物输水：导管中的毛细上升
  • 墨水书写：笔尖的毛细供墨

17.2.2 表面张力的尺度效应

Bond数（或Eötvös数）

重力与表面张力的比值： $$Bo = \frac{\rho g L^2}{\sigma}$$

  • $Bo << 1$：表面张力主导（水滴保持球形）
  • $Bo >> 1$：重力主导（大水体表面平坦）

临界尺度

毛细长度：$l_c = \sqrt{\frac{\sigma}{\rho g}} ≈ 2.7mm$（水）

这解释了为什么：

  • 小于3mm的水滴是球形的
  • 昆虫能在水面行走
  • 微重力下液体行为完全改变

17.2.3 马兰戈尼效应

温度或浓度梯度引起的表面张力梯度驱动流动：

马兰戈尼数： $$Ma = \frac{\Delta\sigma L}{\mu \alpha}$$ 其中 $\alpha$ 是热扩散系数。

日常例子：

  • 酒的"挂杯"现象（酒精浓度梯度）
  • 咖啡环效应（蒸发引起的流动）
  • 焊接熔池的流动模式

17.3 空化现象

17.3.1 空化的物理机理

空化是液体中局部压力降低到饱和蒸汽压以下时产生气泡的现象。这不是沸腾（温度驱动），而是压力驱动的相变。

空化的发生条件

当局部压力 $p_{local} < p_v$（饱和蒸汽压）时发生空化。

利用伯努利方程： $$p_{local} = p_∞ - \frac{1}{2}\rho V^2$$ 空化开始的临界速度： $$V_{crit} = \sqrt{\frac{2(p_∞ - p_v)}{\rho}}$$ 空化数

无量纲参数描述空化倾向： $$\sigma = \frac{p_∞ - p_v}{\frac{1}{2}\rho V^2}$$

  • $\sigma > 1$：不会空化
  • $\sigma ≈ 1$：初生空化
  • $\sigma < 0.5$：严重空化

17.3.2 空化的类型

游移空化（Travelling Cavitation）

单个气泡随流动移动
 o → o → o → ×（溃灭）

片状空化（Sheet Cavitation）

翼型吸力面的连续气膜
  ___________
 /═══════════＼

云状空化（Cloud Cavitation）

周期性脱落的气泡云
 ≈≈≈≈ → ≈≈≈≈ → ≈≈≈≈

超空化（Supercavitation）

完全包裹物体的大空泡
  ▶════════════⟩

17.3.3 空化的工程影响

破坏作用

空泡溃灭产生的局部高压（可达GPa量级）：

  • 材料剥蚀
  • 噪声振动
  • 性能下降

溃灭压力估算（Rayleigh-Plesset）： $$p_{max} \sim \rho c \sqrt{\frac{p_∞ - p_v}{\rho}}$$ 其中 $c$ 是声速。

工程实例：

  1. 船舶螺旋桨：叶梢空化限制转速
  2. 水轮机：尾水管空化影响效率
  3. 阀门节流：空化噪声和振动
  4. 柴油喷嘴：空化改善雾化

防空化设计原则：

  • 增加系统压力（如泵的灌注高度）
  • 降低流速
  • 优化型线减少局部低压
  • 使用抗空蚀材料

17.3.4 空化的有益应用

超声清洗

  • 频率：20-40 kHz
  • 原理：微射流冲刷污垢

医疗碎石

  • 聚焦超声产生空化
  • 气泡溃灭粉碎结石

空化水射流

  • 切割加工
  • 表面处理

17.4 喷雾与雾化

17.4.1 雾化机理

将连续液体分解成细小液滴的过程涉及多种不稳定性：

Rayleigh不稳定性

圆柱液流的轴对称扰动增长：

━━━━━━━ → ～～～～ → ●●●●●
稳定液流    波动增长    断裂成滴

最不稳定波长：$\lambda = 4.51d_{jet}$ 液滴直径：$d_{drop} ≈ 1.89d_{jet}$

Kelvin-Helmholtz不稳定性

高速气流剪切作用：

气 →→→→→
━━━━━━━━ → ～～～～ → 撕裂
液

临界Weber数：$We_{gas} = \frac{\rho_g V_{rel}^2 d}{\sigma} > 12$

17.4.2 雾化方式分类

压力雾化

     ╱╲
    ╱  ╲    锥形液膜
   ╱    ╲
  │  ○  │   旋流室
  └──┬──┘
     │      供液

特点：

  • SMD ∝ (压力)^(-0.25)
  • 喷雾锥角：60°-120°
  • 应用：柴油机、燃气轮机

气动雾化（双流体）

气 →→ ╱╲ ←← 气
    ╱液╲
   ╱╱╲╲
  ╱╱  ╲╲

特点：

  • 更细的液滴（SMD < 50μm）
  • 气液比：0.1-10
  • 应用：涂装、医疗雾化器

旋转雾化

    ╱─────╲
   │   ●   │  高速旋转
    ╲─────╱
     ╲ ╱
    ╲ · ╱    离心甩出
   ╲ · · ╱

特点：

  • 均匀液滴尺寸
  • SMD ∝ (转速)^(-0.8)
  • 应用：喷雾干燥、农药喷洒

17.4.3 液滴尺寸分布

Sauter平均直径（SMD或D32）

体积-表面积平均直径： $$D_{32} = \frac{\sum n_i d_i^3}{\sum n_i d_i^2}$$ 物理意义：具有相同体积-表面积比的均匀液滴直径。

经验关联式

压力雾化（Lefebvre）： $$SMD = 2.25\sigma^{0.25}\mu^{0.25}\dot{m}^{0.25}\Delta p^{-0.5}\rho_a^{-0.25}$$ 气动雾化： $$SMD = 0.95\left(\frac{\sigma}{\rho_a V_{rel}^2}\right)^{0.5}d_{liquid} + 0.13\left(\frac{\mu^2}{\sigma\rho_l}\right)^{0.5}d_{liquid}$$ 典型SMD值：

  • 农药喷洒：200-400 μm
  • 汽车喷漆：20-50 μm
  • 柴油喷射：10-100 μm
  • 医疗雾化：1-5 μm

17.4.4 喷雾特性参数

喷雾锥角

影响因素：

  • 旋流强度（压力雾化）
  • 气液速度比（气动雾化）
  • 环境压力

估算公式（压力雾化）： $$\theta = 2\arctan\left(K\sqrt{\frac{\Delta p}{\rho}}\right)$$ 其中K是喷嘴几何常数。

穿透深度

高压喷射进入静止环境： $$L_{penetration} = C\left(\frac{\rho_l}{\rho_g}\right)^{0.25}d_{nozzle}\sqrt{\frac{\Delta p t^2}{\rho_l}}$$ 蒸发时间

单个液滴的D²定律： $$d^2(t) = d_0^2 - Kt$$ 蒸发常数：$K = \frac{8D_v\rho_v}{\rho_l}\ln(1+B)$

其中B是传质驱动势。

历史人物：约翰·威廉·斯特拉特，第三代瑞利男爵（1842-1919）

瑞利男爵不仅因发现氩气获得1904年诺贝尔物理学奖，他在流体不稳定性理论方面的贡献同样深远。1878年，他发表了关于液体射流不稳定性的开创性论文，奠定了现代雾化理论的基础。

主要贡献：

  1. 瑞利不稳定性（1878） - 证明了圆柱形液流本质上是不稳定的 - 预测了最不稳定波长：λ = 4.51d - 解释了为什么水龙头的水流会断裂成液滴

  2. 瑞利-泰勒不稳定性（1883） - 研究了密度分层流体的界面不稳定性 - 应用：超新星爆炸、惯性约束聚变

  3. 空泡动力学（1917） - 推导了空泡崩溃的经典解 - Rayleigh崩溃时间：$t_c = 0.915R_0\sqrt{\frac{\rho}{p_∞}}$

经典实验：瑞利用小提琴弦的振动激发水射流的受控扰动，首次定量验证了理论预测。这个简单而优雅的实验至今仍在流体力学教学中使用。

工程影响：瑞利的工作直接影响了：

  • 喷墨打印技术（控制液滴大小）
  • 燃料喷射系统设计
  • 农药喷洒技术
  • 3D打印中的液滴沉积

瑞利的格言："理论指导实验，实验验证理论"，完美体现了他的研究哲学。

高级话题：界面追踪方法（VOF/Level-Set）

在数值模拟多相流时，准确追踪相界面是关键挑战。两种主流方法各有优劣：

Volume of Fluid (VOF) 方法

基本思想： 用体积分数 $\alpha$ 表示网格单元中液相的占比：

  • $\alpha = 0$：纯气相
  • $\alpha = 1$：纯液相
  • $0 < \alpha < 1$：界面单元

界面重构：

PLIC (分段线性)：
┌─────┐     ┌─────┐
│ 0.7 │ →   │  ╱  │
└─────┘     └─────┘

优点：

  • 严格守恒
  • 处理拓扑变化（破碎、聚并）自然
  • 计算效率高

缺点：

  • 界面法向和曲率计算困难
  • 数值扩散导致界面模糊

Level-Set 方法

基本思想： 用符号距离函数 $\phi$ 表示到界面的距离：

  • $\phi < 0$：液相内部
  • $\phi = 0$：界面
  • $\phi > 0$：气相内部

演化方程： $$\frac{\partial\phi}{\partial t} + \vec{u} \cdot \nabla\phi = 0$$ 优点：

  • 界面几何信息准确（法向、曲率）
  • 处理表面张力自然
  • 高阶精度容易实现

缺点：

  • 质量不守恒
  • 需要重新初始化

混合方法 CLSVOF

结合两种方法的优点：

  1. VOF保证质量守恒
  2. Level-Set提供几何信息
  3. 相互校正提高精度

工程选择指南：

  • 破碎/雾化问题：VOF
  • 表面张力主导：Level-Set
  • 大变形+精确界面：CLSVOF
  • 工业CFD：通常使用VOF（稳健性）

本章小结

多相流动涉及气液（或液液、气固）界面的复杂相互作用。关键概念包括：

  1. 尺度效应：表面张力在小尺度占主导（毛细长度 ~3mm）
  2. 无量纲数： - Weber数：惯性力/表面张力 - Bond数：重力/表面张力 - 空化数：压力余量/动压
  3. 不稳定性机理： - Rayleigh：表面张力驱动 - Kelvin-Helmholtz：剪切驱动
  4. 工程经验： - 气泡上升速度：小泡∝d²，大泡≈常数 - 空化判据：σ < 1 - 雾化液滴：SMD ∝ We^(-0.5)

记住：多相流的复杂性来自界面的动态特性，而界面行为由局部力平衡决定。

练习题

基础题

17.1 气泡上升速度估算 一个直径2mm的空气泡在20°C的水中上升，估算其终端速度。水的性质：ρ = 1000 kg/m³，μ = 0.001 Pa·s，σ = 0.073 N/m。

Hint：先判断流动区域（Re数），选择合适的阻力系数关联式。

答案

首先估算Re数量级。假设速度~0.1 m/s： Re = ρVd/μ = 1000×0.1×0.002/0.001 = 200

这是过渡区，使用经验关联式。对于2mm气泡，形状接近球形。 根据经验法则，1-10mm气泡：V ∝ d^0.5 典型值：V ≈ 0.2-0.25 m/s

精确计算：使用Mendelson公式 $$V_t = \sqrt{\frac{2\sigma}{ρd} + \frac{gd}{2}}$$ $$V_t = \sqrt{\frac{2×0.073}{1000×0.002} + \frac{9.81×0.002}{2}} = 0.27 m/s$$

17.2 毛细上升高度 玻璃毛细管（内径0.5mm）插入水中，计算毛细上升高度。假设完全润湿（θ = 0°）。

Hint：使用Young-Laplace方程和静压平衡。

答案

毛细上升公式： $$h = \frac{2\sigma\cos\theta}{\rho g r}$$ 代入数值（θ = 0°，cos θ = 1）： $$h = \frac{2×0.073×1}{1000×9.81×0.00025} = 0.0595 m ≈ 60 mm$$ 这解释了为什么细管中液面明显高于容器液面。

17.3 空化判断 水泵入口处压力80 kPa（绝对），流速5 m/s，水温20°C（饱和蒸汽压2.34 kPa）。判断是否会发生空化。

Hint：计算空化数并与临界值比较。

答案

空化数： $$\sigma = \frac{p_∞ - p_v}{0.5\rho V^2} = \frac{80000 - 2340}{0.5×1000×25} = 6.21$$ σ = 6.21 >> 1，不会发生空化。

安全余量充足。若流速增加到11 m/s，σ ≈ 1.28，接近空化临界。

17.4 喷雾液滴尺寸 压力雾化喷嘴，喷射压差5 MPa，燃油流量0.1 kg/s，估算Sauter平均直径。燃油：σ = 0.025 N/m，μ = 0.003 Pa·s，空气密度1.2 kg/m³。

Hint：使用Lefebvre关联式。

答案

Lefebvre公式： $$SMD = 2.25\sigma^{0.25}\mu^{0.25}\dot{m}^{0.25}\Delta p^{-0.5}\rho_a^{-0.25}$$ 代入数值： $$SMD = 2.25×(0.025)^{0.25}×(0.003)^{0.25}×(0.1)^{0.25}×(5×10^6)^{-0.5}×(1.2)^{-0.25}$$ $$SMD ≈ 45 μm$$ 这是典型的柴油喷射液滴尺寸。

挑战题

17.5 雨滴的最大尺寸 为什么很少看到直径超过6mm的雨滴？从流体力学角度解释雨滴尺寸的上限。

Hint：考虑Weber数和气动破碎。

答案

雨滴下落达到终端速度时，气动力与重力平衡。大雨滴的终端速度更高，导致Weber数增加：

对于6mm雨滴：

  • 终端速度 V ≈ 9 m/s
  • Weber数：We = ρ_air×V²×d/σ = 1.2×81×0.006/0.073 ≈ 8

当We > 10-12时，气动力克服表面张力，雨滴变形并破碎。这就是雨滴尺寸的自然上限。

更大的水滴会经历"袋状破碎"：

  1. 底部压力使中心凹陷
  2. 形成袋状结构
  3. 袋壁变薄破裂
  4. 分裂成多个小滴

这种机制确保地面降雨的雨滴尺寸分布相对均匀。

17.6 啤酒泡沫的稳定性 为什么啤酒泡沫比肥皂泡持续时间长？分析液膜排水和泡沫稳定机制。

Hint：考虑表面活性剂、粘度和马兰戈尼效应。

答案

啤酒泡沫稳定性来自多个因素：

  1. 蛋白质稳定：啤酒中的蛋白质是天然表面活性剂，降低表面张力并形成粘弹性膜。

  2. 高粘度：CO₂和溶解物增加液体粘度，减缓液膜排水： 排水时间 ∝ μ/ρg

  3. 马兰戈尼效应：液膜变薄处表面张力增大，驱动液体回流，自我修复。

  4. CO₂补充：持续的CO₂释放维持泡内压力，延缓破裂。

相比之下，肥皂泡：

  • 液膜排水快（低粘度）
  • 蒸发导致变薄
  • 无内部气体补充

典型寿命：肥皂泡~秒级，啤酒泡沫~分钟级。

17.7 超疏水表面的液滴弹跳 液滴撞击超疏水表面（接触角>150°）可以完全弹跳。分析弹跳的能量转换和临界条件。

Hint：考虑动能、表面能和粘性耗散的平衡。

答案

液滴弹跳涉及复杂的能量转换：

能量分析：

  1. 初始动能：$E_k = \frac{1}{2}mV^2$
  2. 最大铺展时转化为表面能：$E_s = \sigma\Delta A$
  3. 回缩时表面能转回动能
  4. 粘性耗散：$E_d \sim \mu V \Delta A t_{contact}$

弹跳条件： Weber数适中：5 < We < 50

  • We太小：动能不足，gentle deposition
  • We太大：破碎或粘附

接触时间（Rayleigh时间尺度）： $$t_{contact} \sim \sqrt{\frac{\rho R^3}{\sigma}} ≈ 2.6\sqrt{\frac{\rho R^3}{\sigma}}$$

神奇的是，接触时间几乎与撞击速度无关！这在自清洁表面设计中很重要。

应用：防冰涂层、自清洁玻璃、防雾表面。

17.8 喷墨打印的精度极限 现代喷墨打印机能产生约20μm的液滴。分析更小液滴面临的物理限制。

Hint：考虑表面张力、蒸发和静电效应。

答案

小液滴面临多重物理限制：

  1. Rayleigh不稳定性： 喷嘴直径d_n → 液滴直径d_d ≈ 1.89d_n 要产生20μm液滴，需要~10μm喷嘴，制造困难。

  2. 表面张力压力： ΔP = 2σ/R 对于10μm液滴：ΔP = 2×0.073/5×10^-6 = 29 kPa 需要更高驱动压力。

  3. 快速蒸发： 蒸发时间 ∝ d² 20μm水滴在空气中：t_evap ~ 毫秒级 可能在到达纸面前蒸发。

  4. 气动偏移： Stokes数：St = ρ_d×d²×V/(18μ_g×L) 小液滴St小，易受气流干扰，定位精度降低。

  5. 静电效应： 小液滴电荷密度高，库仑斥力可能超过惯性，导致轨迹不稳定。

解决方案：

  • 使用低挥发性墨水
  • 控制环境湿度
  • 静电辅助定位
  • 减小喷射距离

物理极限：约5-10μm（实验室条件），商业应用：20-50μm。

常见陷阱与错误

概念误区

  1. "雨滴是泪滴形的" - 错误：受空气动力影响，大雨滴实际是汉堡形 - 正确：只有非常小的雨滴（<1mm）保持球形

  2. "气泡上升速度与大小成正比" - 错误：这只对很小的气泡（<1mm）成立 - 正确：大气泡（>1cm）上升速度几乎恒定

  3. "表面张力总是使液滴球形" - 错误：当其他力（重力、惯性）超过表面张力时失效 - 正确：需要检查Bond数和Weber数

计算错误

  1. 混淆不同平均直径 - D10（算术平均）≠ D32（Sauter平均） - SMD用于传热传质计算，D10用于数量统计

  2. 忽略尺度效应 - 微米级：表面张力主导 - 毫米级：重力和表面张力竞争 - 厘米级：重力主导

  3. 空化数计算使用表压 - 必须使用绝对压力 - p_v是绝对压力

工程陷阱

  1. 忽视多相流的时间尺度 - 气泡聚并：毫秒级 - 液滴蒸发：秒级 - 泡沫排水：分钟级

  2. 简单外推单相流关联式 - 两相流压降≠单相×倍数 - 需要专门的两相流关联式

最佳实践检查清单

设计阶段

☐ 识别流型

  • 确定是分散流还是分离流
  • 估算体积分数范围
  • 判断主导作用力

☐ 无量纲数分析

  • 计算We、Bo、Ca数
  • 确定流动区域
  • 选择合适的关联式

☐ 尺度评估

  • 毛细长度vs系统尺寸
  • Kolmogorov尺度vs气泡/液滴尺寸
  • 时间尺度匹配

系统设计

☐ 防空化措施

  • 计算NPSH裕量
  • 避免急剧压降
  • 考虑温度影响

☐ 雾化系统

  • 匹配液滴尺寸需求
  • 考虑蒸发损失
  • 评估分布均匀性

☐ 气液分离

  • 重力分离器尺寸
  • 旋流器切割粒径
  • 除雾器效率

运行优化

☐ 监测参数

  • 压力（空化预警）
  • 流型转变
  • 液滴尺寸分布

☐ 故障诊断

  • 异常噪声→空化
  • 效率下降→流型改变
  • 振动→两相流不稳定

☐ 维护要点

  • 定期检查空蚀
  • 清理气液分离器
  • 校准两相流量计

安全考虑

☐ 压力波动

  • 段塞流冲击载荷
  • 水锤效应放大
  • 安全阀选型

☐ 相变风险

  • 闪蒸可能性
  • 冷凝锤
  • 热分层

记住：多相流比单相流复杂得多，经验关联式的不确定性通常在±30%，设计时需要充足的安全裕量。