什么是可压缩多相流?#

在流体力学中,**多相流(multiphase flow)**指两个或更多相(phase)共存的流动。 水与空气、燃料液滴与高温气体、水下爆炸中的气泡等都是典型例子。

加上可压缩性(compressibility),问题就急剧复杂化。 各相拥有不同的状态方程(EOS),并以界面(interface)为分界 产生冲击波(shock wave),发生反射、透射、折射。

与单相可压缩流的区别#

单相 Euler 方程如下:

$$\frac{\partial \mathbf{U}}{\partial t} + \frac{\partial \mathbf{F}}{\partial x} = 0$$

其中守恒变量 $\mathbf{U} = [\rho, \rho u, E]^T$, 单一 EOS(例如理想气体 $p = (\gamma - 1)\rho e$)即可闭合系统。

多相流中界面两侧 EOS 不同。 设气相区域使用 $\gamma = 1.4$ 的理想气体, 液相区域使用 stiffened gas EOS:

$$p = (\gamma - 1)\rho e - \gamma p_\infty$$

数值上的麻烦出现在界面附近两种 EOS **混合(mixing)**时。 基于守恒变量直接平均会产生非物理的压力振荡(pressure oscillation)。

为什么难:三大核心难题#

1. 界面处的压力振荡#

最著名的问题。1994 年 Abgrall 指出后, 解决该问题成为多相流数值方法的起点。

“

"可压缩多相流数值方法的历史,就是与压力振荡斗争的历史。"

普通的守恒型有限体积法,在界面处由于 $\gamma$ 取中间值会产生 spurious oscillation。

2. 状态方程的多样性#

实际问题中仅靠理想气体不够。 水使用 Tait EOS 或 stiffened gas EOS, 高压状态使用 Jones-Wilkins-Lee (JWL) EOS, 反应性流动需要与化学反应耦合的 EOS。

需要对每种 EOS 推广 Riemann solver,实现复杂度大幅上升。

3. 界面追踪 vs 界面捕捉#

如何处理界面大致分两种方法:

为什么重要#

可压缩多相流仿真在多个工程领域处于核心地位:

• 水下爆炸(UNDEX): 舰艇生存性评估
• 超音速燃烧: 冲压发动机内部燃料喷射/混合
• 医疗领域: 冲击波碎石(lithotripsy)中的气泡行为
• 航天推进: 低温推进剂的晃动(sloshing)

下篇预告#

下一篇介绍该领域的核心工具:Riemann 问题与 Godunov 型格式。 首先理解单相 Euler 方程的 Riemann 问题, 再看它如何扩展到多相流。

将探针滑到 α≈0.5 附近，Wood 声速会降到单相的 1/100。