$\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を下回ると粒子が流動せず、Utを超えると触媒が装置外に飛び出してしまうため、適正な運転速度の決定に本シミュレーターの計算が活用されます。
石炭・バイオマスの燃焼・ガス化:微粉炭を空気で流動化させながら燃焼させる「循環流動層ボイラー」では、炉内の粒子濃度分布や滞留時間を制御するために、原料粒子のUmfとUtが重要な設計パラメータとなります。
製薬・食品の乾燥・造粒:粉末や顆粒を温風で流動させながら乾燥したり、コーティングしたりする装置(流動層造粒機)で使われます。製品の品質を均一にするため、粒子が均一に流動する速度域を把握する必要があります。
廃棄物処理・サンドブラスト:廃棄物の流動層による熱分解や、研磨材を空気で加速して行うサンドブラストでも原理は同じです。処理効率や研磨効果を最大化する流体速度の選定に、これらの基礎計算が役立ちます。
このシミュレーターを使い始める際、特に現場経験が浅いエンジニアが陥りがちな落とし穴がいくつかあります。まず一つ目は、「計算されたUmfがそのまま運転速度」だと思ってしまうこと。実は、最小流動化速度はあくまで「流動が始まる」点。実際の装置では、粒子を活発に混合させたり、熱伝達を良くしたりするために、Umfの2倍から10倍程度の速度で運転することがほとんどです。例えば、粒径100μmのFCC触媒でUmfが0.02 m/sと計算されても、実際の流動床反応器では0.1 m/s以上で動かしている、といった具合です。
二つ目は「粒子径」の定義と現実です。シミュレーターでは均一な球体を想定しますが、実在する粉末は大小さまざまな粒径分布(粒度分布)を持ち、形もいびつです。このため、計算値はあくまで目安。実機設計では、代表径としてサウター平均径やメディアン径を使うなど、目的に応じて適切な値を選ぶ必要があります。また、湿潤条件下では粒子が凝集して「見かけの粒径」が大きくなるため、計算結果と実現象が大きく乖離することがあります。
三つ目は圧力損失の「一定値」への過信。流動化開始後、圧力損失がほぼ一定になる、と学びますが、これは理想的な均一流動層の場合。実際には、大型装置では気泡の発生や粒子の偏析により、圧力損失は変動します。シミュレーターの美しいカーブは「理論的な骨格」と理解し、実機では計測値とのすり合わせが不可欠です。
石英砂流動層(dp=150μm、ρparticle=2650kg/m³、水中)の場合:流体密度998kg/m³、動粘度0.001Pa·sを入力。アルキメデス数Ar≈1850が算出され、最小流動化速度Umf≈0.028m/s、圧力損失ΔP/L≈8200Pa/mが得られます。これは高さ2mの流動層で約16.4kPaの圧力差に相当します。