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

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

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

流体パラメータ

一時停止中はスライダーを動かすと結果が即座に更新されます。

流動様式リアルタイム・アニメーション(気液二相流)
流動様式
ボイド率 α
クオリティ x
圧力勾配 [Pa/m]
気相流量 Q_G を上げるとボイド率 α が増え、気泡流→スラグ流→チャーン流→環状流へ遷移します。左パネルのスライダーで変化を確認してください。
計算結果

計算結果

フローパターン-
ボイド率 α-
クオリティ 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は液相流量Q_Lと気相流量Q_Gから自動計算されるので、Q_G/Q_Lの比を変えてみて。パラメータ$X_{tt}$や圧力損失比$\phi_L^2$が変わる様子を確認できるよ。

よくある質問

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

実世界での応用

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

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

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

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

よくある誤解と注意点

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

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

最後に、圧力損失計算の「C定数」の選択。ロックハート-マルチネッリ相関の $C$ 値は、各相が層流か乱流かで変わります。本ツールではRe_LとRe_Gから、層流-層流は5、乱流-層流は10、層流-乱流は12、乱流-乱流は20を用いる簡易判定にしています。遷移域では値が不確かになるため、実務では複数相関や安全側の圧損も確認してください。

使い方ガイド

  1. 液相質量流量(kg/s)と気相質量流量(kg/s)を入力し、液体密度(kg/m³)と気体密度(kg/m³)を設定します
  2. 配管内径(mm)、液相流量Q_L(L/min)、気相流量Q_G(L/min)を指定すると、速度・クオリティ・ボイド率・圧力損失が自動算出されます
  3. ベーカー図にプロットされたフローパターン(気泡流・スラグ流・層状流・波状流・環状流)を確認し、ロックハート-マルチネッリ相関式で圧力損失(Pa/m)を取得します

具体的な計算例

デフォルトの水-空気系で配管内径25mm、液相流量10L/min、気相流量5L/min、ρ_L=1000kg/m³、ρ_G=1.2kg/m³の場合、液相速度≈0.340m/s、気相速度≈0.170m/s、クオリティx≈0.060%、ボイド率α≈0.0505です。補正後のベーカー図ではスラグ流、ロックハート-マルチネッリ相関ではX_tt≈22.04、φ²_L≈1.46、圧力損失≈111Pa/mとなります。

実務での注意点