二相流シミュレーター 戻る
流体工学・化学工学

二相流シミュレーター(気液流動)

気液二相流のフローパターン(気泡流・スラグ流・環状流等)をベーカー図で判定。ロックハート-マルチネッリ相関で圧力損失を計算。

流体パラメータ

計算結果

計算結果

フローパターン-
ボイド率 α-
クオリティ x-
スリップ比 S-
M-Mパラメータ X_tt-
二相流乗数 φ²_L-
圧力損失/m [Pa/m]-
ベーカー図(フローパターンマップ)
圧力損失 vs クオリティ x
理論・主要公式

$X_{tt}= \left(\frac{1-x}{x}\right)^{0.9}\left(\frac{\rho_G}{\rho_L}\right)^{0.5}\left(\frac{\mu_L}{\mu_G}\right)^{0.1}$
$\phi_L^2 = 1 + \frac{C}{X_{tt}}+ \frac{1}{X_{tt}^2}$

二相流シミュレーター(気液流動)とは

🙋
二相流って何ですか?水と空気が一緒に流れるようなことですか?
🎓
その通り!気体と液体が一緒に管の中を流れる現象だね。例えば、エアコンの冷媒配管や石油プラントの流送管では日常的に起こっているよ。このシミュレーターでは、流量や物性を上のスライダーで変えると、流れの様子(フローパターン)がどう変わるか、ベーカー図上で即座に確認できるんだ。
🙋
フローパターンが変わるとなぜ重要なのですか?
🎓
非常に重要!流れ方で圧力損失や熱の伝わり方が何十倍も変わってしまうんだ。例えば「スラグ流」だと大きな気泡が断続的に流れるから、配管がガタガタ振動する原因になる。シミュレーターで「液相質量速度」を大きくして確認してみて。「気泡流」から「環状流」に変わるのが分かるよ。設計ではこの予測が命なんだ。
🙋
圧力損失の計算で出てくる「ロックハート-マルチネッリ相関」って、どうやって使うんですか?
🎓
大まかに言うと、気相だけ、液相だけの圧力損失をまず計算して、それらを結びつける「二相流乗数」を求める便利な経験式だ。シミュレーター右側のグラフで、蒸気乾き度「x」のスライダーを0.1から0.9に動かしてみて。パラメータ$X_{tt}$が小さくなって、圧力損失比$\phi_L^2$が急激に大きくなるのが見えるだろう?これが二相流の難しいところだね。

よくある質問

気相と液相の質量流量[kg/s]、各相の密度[kg/m³]、大気密度、管断面積[m²]が必要です。これらの値から修正マスフラックスG_GとG_Lを計算し、ベーカー図上の領域を参照して流動様式(気泡流・スラグ流・環状流など)を判定します。
二相流乗数は、単相の圧力損失に対して二相流による増加率を表す無次元数です。この乗数と単相圧力損失を掛け合わせることで、気液二相流の実際の圧力損失を簡便に推定できます。入力条件から計算されるパラメータXを用いて求めます。
ベーカー図は水平管・低圧の実験データに基づくため、垂直管や高圧条件では誤差が大きくなります。また、ロックハート-マルチネッリ相関は環状流や高粘度流体で精度が低下するため、適用範囲を確認してご使用ください。
化学プラントの配管設計、原子炉の冷却系解析、石油・ガスの輸送ライン設計、熱交換器内の二相流評価などに利用できます。特に、圧力損失予測と流動様式の判定が重要なプロセス設計の初期検討に適しています。

実世界での応用

化学プラント・石油精製:反応器からの生成物や加熱炉内の流体の流送管設計に不可欠です。フローパターンを誤ると、配管の腐食・磨耗が特定部位で集中したり、流量計測が不正確になったりするため、シミュレーターによる事前検討が行われます。

空調・冷凍サイクル:エアコンや冷蔵庫の冷媒配管では、蒸発器や凝縮器内で気液二相流が発生します。環状流は熱伝達に優れますが、圧力損失も大きいため、効率とポンプ動力のトレードオフをこのようなツールで最適化します。

地熱・原子力発電:地熱発電の蒸気・熱水混合流体や、原子炉冷却材の沸騰流の挙動を予測するために使用されます。特にスラグ流は流量や圧力の脈動を引き起こすため、安全設計上その発生領域を避ける必要があります。

宇宙機・航空機の燃料供給:無重力状態や高加速度条件下での燃料(液体)とヘリウムガス(加圧用)の二相流挙動は複雑です。地球上での実験データと相関式を組み合わせ、信頼性の高い供給システム設計に応用されています。

よくある誤解と注意点

このシミュレーターを使い始めるとき、いくつかつまずきやすいポイントがあるから気をつけて。まず「物性値の入力はデフォルトのままでいいや」という考え。これは危険だ。例えば、高温高圧の蒸気-水の二相流を常温常圧の空気-水の物性で計算すると、密度比や粘度比が大きく異なるから、フローパターン判定も圧損計算も大きく外れる。実プロセスに近い条件、できれば物性データベースで調べた値を入れよう。

次に「ベーカー図の境界線は絶対的なもの」という誤解。あの図はあくまで代表的な実験データに基づく「目安」だ。実際の配管の傾きや内面粗さ、入口形状で遷移は変わる。シミュレーターで「スラグ流」と出ても、配管レイアウトを少し変えるだけで「気泡流」に落ち着くかもしれない。ツールの出力は「可能性の一つ」として受け止め、特に危険領域と判定された時は安全側の設計を心がけよう。

最後に、圧力損失計算の「C定数」の選択。ロックハート-マルチネッリ相関の $C$ 値は流動様式で変わるけど、実際の流れは遷移領域で曖昧だ。例えば、計算上は「環状流」と「波状流」の境界付近だと、両方の $C$ 値で計算して悪い方(圧損が大きい方)を採用するのが実務の鉄則だ。ツールは自動判定してくれるが、その背後にある理屈を知っておくと、判断に自信が持てるよ。

使い方ガイド

  1. 液相質量流量(kg/s)と気相質量流量(kg/s)を入力し、液体密度(kg/m³)と気体密度(kg/m³)を設定します
  2. 配管内径(mm)を指定してから、液相速度と気相速度を計算開始ボタンで自動算出します
  3. ベーカー図にプロットされたフローパターン(層流・波状流・栓流・環状流)を確認し、ロックハート-マルチネッリ相関式で圧力損失(Pa/m)を取得します

具体的な計算例

水-空気系で配管内径25.4mm、液相流量0.5kg/s(ρ=998kg/m³)、気相流量0.08kg/s(ρ=1.2kg/m³)の場合、液相速度2.0m/s、気相速度17.6m/sとなります。ベーカー線図上でこれらの速度パラメータはフルード数で無次元化されフローパターンが環状流と判定され、相関式により圧力損失は約340Pa/mと算出されます。

実務での注意点