当你喝可乐时看到上升的气泡，当你在雨中开车时看到挡风玻璃上的水滴，或者当你使用喷雾器时观察到细小的液滴云雾——这些都是多相流动的日常例子。多相流动无处不在，从内燃机的燃油喷射到血管中的血液流动，从化工反应器中的气液混合到海浪的破碎。本章将帮助你建立对这些复杂现象的物理直觉，掌握工程实践中的经验法则。
想象一个气泡在静止的水中上升。它的运动速度由什么决定？这是一个看似简单却蕴含丰富物理的问题。
力的平衡
气泡受到三个主要力的作用：
其中 $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： 不规则/破碎
～～～
～ ～
～～～
经验法则：
雨滴下落是最常见的液滴运动。你可能认为雨滴是泪滴形的，但实际上：
实际雨滴形状：
小雨滴（< 2mm）： 球形
O
中等雨滴（2-5mm）： 扁球形
___
( )
‾‾‾
大雨滴（> 5mm）： 汉堡形（底部凹陷）
_____
(_ _)
‾‾‾‾‾
Weber数的作用
液滴变形程度由Weber数决定： \(We = \frac{\rho V^2 d}{\sigma}\)
临界Weber数：
聚并（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$ 是湍流耗散率。
表面张力源于分子间作用力的不平衡。在液体内部，分子受到各向同性的吸引力；而在界面处，分子受力不平衡，产生指向液体内部的净力。
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$ 是管径。
日常应用：
Bond数（或Eötvös数）
重力与表面张力的比值： \(Bo = \frac{\rho g L^2}{\sigma}\)
临界尺度
毛细长度：$l_c = \sqrt{\frac{\sigma}{\rho g}} ≈ 2.7mm$（水）
这解释了为什么：
温度或浓度梯度引起的表面张力梯度驱动流动：
马兰戈尼数： \(Ma = \frac{\Delta\sigma L}{\mu \alpha}\)
其中 $\alpha$ 是热扩散系数。
日常例子：
空化是液体中局部压力降低到饱和蒸汽压以下时产生气泡的现象。这不是沸腾（温度驱动），而是压力驱动的相变。
空化的发生条件
当局部压力 $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}\)
游移空化（Travelling Cavitation）
单个气泡随流动移动
o → o → o → ×（溃灭）
片状空化（Sheet Cavitation）
翼型吸力面的连续气膜
___________
/═══════════＼
云状空化（Cloud Cavitation）
周期性脱落的气泡云
≈≈≈≈ → ≈≈≈≈ → ≈≈≈≈
超空化（Supercavitation）
完全包裹物体的大空泡
▶════════════⟩
破坏作用
空泡溃灭产生的局部高压（可达GPa量级）：
溃灭压力估算（Rayleigh-Plesset）： \(p_{max} \sim \rho c \sqrt{\frac{p_∞ - p_v}{\rho}}\)
其中 $c$ 是声速。
工程实例：
防空化设计原则：
超声清洗
医疗碎石
空化水射流
将连续液体分解成细小液滴的过程涉及多种不稳定性：
Rayleigh不稳定性
圆柱液流的轴对称扰动增长：
━━━━━━━ → ～～～～ → ●●●●●
稳定液流 波动增长 断裂成滴
最不稳定波长：$\lambda = 4.51d_{jet}$ 液滴直径：$d_{drop} ≈ 1.89d_{jet}$
Kelvin-Helmholtz不稳定性
高速气流剪切作用：
气 →→→→→
━━━━━━━━ → ～～～～ → 撕裂
液
临界Weber数：$We_{gas} = \frac{\rho_g V_{rel}^2 d}{\sigma} > 12$
压力雾化
╱╲
╱ ╲ 锥形液膜
╱ ╲
│ ○ │ 旋流室
└──┬──┘
│ 供液
特点：
气动雾化（双流体）
气 →→ ╱╲ ←← 气
╱液╲
╱╱╲╲
╱╱ ╲╲
特点：
旋转雾化
╱─────╲
│ ● │ 高速旋转
╲─────╱
╲ ╱
╲ · ╱ 离心甩出
╲ · · ╱
特点：
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值：
喷雾锥角
影响因素：
估算公式（压力雾化）： \(\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是传质驱动势。
瑞利男爵不仅因发现氩气获得1904年诺贝尔物理学奖，他在流体不稳定性理论方面的贡献同样深远。1878年，他发表了关于液体射流不稳定性的开创性论文，奠定了现代雾化理论的基础。
主要贡献：
经典实验：瑞利用小提琴弦的振动激发水射流的受控扰动，首次定量验证了理论预测。这个简单而优雅的实验至今仍在流体力学教学中使用。
工程影响：瑞利的工作直接影响了：
瑞利的格言：”理论指导实验，实验验证理论”，完美体现了他的研究哲学。
在数值模拟多相流时，准确追踪相界面是关键挑战。两种主流方法各有优劣：
基本思想： 用体积分数 $\alpha$ 表示网格单元中液相的占比：
界面重构：
PLIC (分段线性)：
┌─────┐ ┌─────┐
│ 0.7 │ → │ ╱ │
└─────┘ └─────┘
优点：
缺点：
基本思想： 用符号距离函数 $\phi$ 表示到界面的距离：
演化方程： \(\frac{\partial\phi}{\partial t} + \vec{u} \cdot \nabla\phi = 0\)
优点：
缺点：
结合两种方法的优点：
工程选择指南：
多相流动涉及气液（或液液、气固）界面的复杂相互作用。关键概念包括：
记住：多相流的复杂性来自界面的动态特性，而界面行为由局部力平衡决定。
17.1 气泡上升速度估算 一个直径2mm的空气泡在20°C的水中上升，估算其终端速度。水的性质：ρ = 1000 kg/m³，μ = 0.001 Pa·s，σ = 0.073 N/m。
Hint：先判断流动区域（Re数），选择合适的阻力系数关联式。
17.2 毛细上升高度 玻璃毛细管（内径0.5mm）插入水中，计算毛细上升高度。假设完全润湿（θ = 0°）。
Hint：使用Young-Laplace方程和静压平衡。
17.3 空化判断 水泵入口处压力80 kPa（绝对），流速5 m/s，水温20°C（饱和蒸汽压2.34 kPa）。判断是否会发生空化。
Hint：计算空化数并与临界值比较。
17.4 喷雾液滴尺寸 压力雾化喷嘴，喷射压差5 MPa，燃油流量0.1 kg/s，估算Sauter平均直径。燃油：σ = 0.025 N/m，μ = 0.003 Pa·s，空气密度1.2 kg/m³。
Hint：使用Lefebvre关联式。
17.5 雨滴的最大尺寸 为什么很少看到直径超过6mm的雨滴？从流体力学角度解释雨滴尺寸的上限。
Hint：考虑Weber数和气动破碎。
17.6 啤酒泡沫的稳定性 为什么啤酒泡沫比肥皂泡持续时间长？分析液膜排水和泡沫稳定机制。
Hint：考虑表面活性剂、粘度和马兰戈尼效应。
17.7 超疏水表面的液滴弹跳 液滴撞击超疏水表面（接触角>150°）可以完全弹跳。分析弹跳的能量转换和临界条件。
Hint：考虑动能、表面能和粘性耗散的平衡。
17.8 喷墨打印的精度极限 现代喷墨打印机能产生约20μm的液滴。分析更小液滴面临的物理限制。
Hint：考虑表面张力、蒸发和静电效应。
☐ 识别流型
☐ 无量纲数分析
☐ 尺度评估
☐ 防空化措施
☐ 雾化系统
☐ 气液分离
☐ 监测参数
☐ 故障诊断
☐ 维护要点
☐ 压力波动
☐ 相变风险
记住：多相流比单相流复杂得多，经验关联式的不确定性通常在±30%，设计时需要充足的安全裕量。