圧縮性多相流とは?#

流体力学において**多相流(multiphase flow)**とは、2つ以上の相(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$$

問題は数値的に界面付近で2つのEOSが**混合(mixing)**された時に発生します。 保存変数ベースで単純に平均を取ると、非物理的な圧力振動(pressure oscillation)が生じます。

なぜ難しいか:中核的な3つの難題#

1. 界面における圧力振動#

最も有名な問題です。1994年にAbgrallが指摘して以来、 この問題を解決することが多相流数値手法の出発点となりました。

“

「圧縮性多相流数値解析の歴史は、すなわち圧力振動との戦いの歴史である。」

単純な保存型有限体積法は、界面で $\gamma$ が中間値を取ることでspurious oscillationを生み出します。

2. 状態方程式の多様性#

実際の問題では理想気体だけでは不十分です。 水はTait EOSまたはstiffened gas EOS、 高圧状態ではJones-Wilkins-Lee (JWL) EOS、 反応性流では化学反応と結合したEOSが必要です。

各EOSに対してRiemann solverを一般化する必要があるため、実装の複雑さが大きく上がります。

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

界面をどう扱うかは大きく2つのアプローチがあります:

なぜ重要か#

圧縮性多相流シミュレーションは多様な工学分野で核心的です:

• 水中爆発(UNDEX): 艦艇の生存性評価
• 超音速燃焼: スクラムジェットエンジン内部の燃料噴射/混合
• 医療分野: 衝撃波結石破砕(lithotripsy)における気泡挙動
• 宇宙推進: 極低温推進剤のスロッシング(sloshing)

次の記事では#

この分野の核心的ツールであるRiemann問題とGodunov型スキームを扱います。 単相Euler方程式のRiemann問題をまず理解した上で、 これが多相流にどう拡張されるかを見ていきます。

α≈0.5 付近にスライダを動かすと Wood 音速が単相の 1/100 まで下がる。