流動層シミュレーター 戻る
化学工学

流動層シミュレーター — 最小流動化速度・圧力損失

粒子径・密度・流体粘度を調整してアルキメデス数・Umf・Ut をリアルタイム計算。Ergun 式による圧力損失曲線と粒子床アニメーションで流動化の物理を直感的に理解できます。

粒子・流体パラメータ
粒子径 dp
μm
粒子密度 ρp
kg/m³
流体密度 ρf
kg/m³
動粘度 µ 0.018 mPa·s
空隙率 ε
表面速度 U 0.10 m/s
計算結果
Ar
Umf (m/s)
Ut (m/s)
ΔP/L (Pa/m)
粒子床アニメーション
固定床 — U < Umf
圧力損失曲線(Ergun式)
圧力
理論・主要公式
Ergun式(固定床):
$\frac{\Delta P}{L}= \frac{150\mu U(1-\varepsilon)^2}{d_p^2\varepsilon^3}+ \frac{1.75\rho_f U^2(1-\varepsilon)}{d_p\varepsilon^3}$

アルキメデス数:
$Ar = \frac{d_p^3\rho_f(\rho_p-\rho_f)g}{\mu^2}$

流動層シミュレーターとは

🙋
「流動層」って何ですか?粒子が浮いてるあの状態ですよね。このシミュレーターで何がわかるんですか?
🎓
そうだね、大まかに言うと、下から流体(空気や水)を吹き上げて粒子を浮遊・混合させる状態だ。このツールでは、その状態を作り始めるのに必要な最低の流速「最小流動化速度(Umf)」や、粒子が吹き飛ばされる「終末速度(Ut)」を計算できるんだ。左側のパラメータで粒子の大きさや密度を変えてみると、アニメーションの粒子の動きやグラフの曲線がリアルタイムで変わるよ。
🙋
え、そうなんですか!グラフの縦軸が「圧力損失」で、横軸が「表面速度」って書いてありますね。このカーブは何を表しているんですか?
🎓
良いところに気づいたね。このカーブは「Ergun式」で計算される、固定床(粒子が動かない状態)の圧力損失だ。流速を上げていくと、圧力損失も上がる。でも、ある点(Umf)で粒子が動き始めると、圧力損失はほぼ一定になるんだ。シミュレーターでスライダーをゆっくり右に動かして流速Uを上げていくと、グラフ上を点が移動して、その転換点が視覚的に確認できるよ。
🙋
なるほど!で、計算結果に出てくる「アルキメデス数」って何ですか?名前は聞いたことありますが…。
🎓
これはとっても重要な無次元数なんだ。粒子にかかる「浮力」と流体の「粘性力」の比を表している。これが小さい(例えばAr<36)と粘性が支配的で粒子はゆっくり沈む。大きいと慣性が支配的でバサッと落ちる。この値によって、Umfや終末速度Utの計算式自体が切り替わるんだ。ツールで粒子径dpを大きくしてみてごらん、アルキメデス数が一気に跳ね上がるから。

よくある質問

Umfは粒子が流動し始める最低限の流体速度で、固定床から流動層へ遷移する点です。一方、Utは粒子が流体によって層外へ飛び出し始める速度です。Umfは粒子の流動化開始、Utは粒子の流出限界を表し、通常Umf < Utの関係があります。
空隙率が小さいほど粒子間の隙間が狭くなり、Ergun式の粘性・慣性損失が増大するため、同じ流速でも圧力損失が大きくなります。特に固定床領域では空隙率の影響が顕著で、流動化開始後の圧力損失はほぼ一定になります。
粒子径が大きくなると、粒子の重量が増すため、流動化に必要な流体の抗力も大きくなります。その結果、最小流動化速度(Umf)は増加します。具体的には、Umfは粒子径の約2乗に比例するため、径を2倍にするとUmfは約4倍になります。
本ツールは教育・理解目的でErgun式と簡易モデルに基づいており、実機設計には適用できません。実際の流動層では粒子形状の不均一性や壁効果、凝集性などが影響するため、実験データやより詳細なシミュレーションと併用することを推奨します。

実世界での応用

化学プラントの反応装置:触媒粒子と原料ガスを接触させる「流動床反応器」の設計に不可欠です。Umfを下回ると粒子が流動せず、Utを超えると触媒が装置外に飛び出してしまうため、適正な運転速度の決定に本シミュレーターの計算が活用されます。

石炭・バイオマスの燃焼・ガス化:微粉炭を空気で流動化させながら燃焼させる「循環流動層ボイラー」では、炉内の粒子濃度分布や滞留時間を制御するために、原料粒子のUmfとUtが重要な設計パラメータとなります。

製薬・食品の乾燥・造粒:粉末や顆粒を温風で流動させながら乾燥したり、コーティングしたりする装置(流動層造粒機)で使われます。製品の品質を均一にするため、粒子が均一に流動する速度域を把握する必要があります。

廃棄物処理・サンドブラスト:廃棄物の流動層による熱分解や、研磨材を空気で加速して行うサンドブラストでも原理は同じです。処理効率や研磨効果を最大化する流体速度の選定に、これらの基礎計算が役立ちます。

よくある誤解と注意点

このシミュレーターを使い始める際、特に現場経験が浅いエンジニアが陥りがちな落とし穴がいくつかあります。まず一つ目は、「計算されたUmfがそのまま運転速度」だと思ってしまうこと。実は、最小流動化速度はあくまで「流動が始まる」点。実際の装置では、粒子を活発に混合させたり、熱伝達を良くしたりするために、Umfの2倍から10倍程度の速度で運転することがほとんどです。例えば、粒径100μmのFCC触媒でUmfが0.02 m/sと計算されても、実際の流動床反応器では0.1 m/s以上で動かしている、といった具合です。

二つ目は「粒子径」の定義と現実です。シミュレーターでは均一な球体を想定しますが、実在する粉末は大小さまざまな粒径分布(粒度分布)を持ち、形もいびつです。このため、計算値はあくまで目安。実機設計では、代表径としてサウター平均径やメディアン径を使うなど、目的に応じて適切な値を選ぶ必要があります。また、湿潤条件下では粒子が凝集して「見かけの粒径」が大きくなるため、計算結果と実現象が大きく乖離することがあります。

三つ目は圧力損失の「一定値」への過信。流動化開始後、圧力損失がほぼ一定になる、と学びますが、これは理想的な均一流動層の場合。実際には、大型装置では気泡の発生や粒子の偏析により、圧力損失は変動します。シミュレーターの美しいカーブは「理論的な骨格」と理解し、実機では計測値とのすり合わせが不可欠です。

使い方ガイド

  1. 粒子径(dp)をμm単位で入力します。例えばシリカゲル50μm、石英砂200μmなど
  2. 粒子密度(rhoParticle)と流体密度(rhoFluid)を設定します。水のΔT=20℃時ρ=998kg/m³、空気1.2kg/m³が標準値です
  3. 動粘度(mu)をmPa・s単位で入力後、空隙率(eps)を0.4~0.5の範囲で調整すると最小流動化速度Umfが自動計算されます
  4. Ergun式とアルキメデス数(Ar)から圧力損失ΔP/Lが導出されます

具体的な計算例

石英砂流動層(dp=150μm、ρparticle=2650kg/m³、水中)の場合:流体密度998kg/m³、動粘度0.001Pa·sを入力。アルキメデス数Ar≈1850が算出され、最小流動化速度Umf≈0.028m/s、圧力損失ΔP/L≈8200Pa/mが得られます。これは高さ2mの流動層で約16.4kPaの圧力差に相当します。

実務での注意点