抗力係数（c_D）は、既知の流動条件下で抗力を測定し、その結果を無次元相関式によって適用することで決定されます。 これは次のように計算されます。
[ c_D=\frac{D}{\tfrac{1}{2}\rho V^2S} ]
ここで、(D)は抗力、(\rho)は流体密度、(V)は選択した基準速度、(S)は代表正面面積です。流体・固体相互作用を伴うパイロットプラントでは、技術者は主にレイノルズ数の関数として(c_D)を用い、粒子形状、充填層の空隙率、壁面効果、さらに沈降または流動化における重力も考慮します。
(c_D)の実用的な目的は、測定または相関によって得られた流体・固体抵抗を、スケールに依存しない設計量へ変換することです。 適切な速度、面積、レイノルズ数、および幾何形状を定義すれば、この係数を用いて圧力損失、沈降挙動、終端速度、流動化条件、スケールアップ性能を推定できます。
抗力係数が表すもの
抗力は圧力効果と表面効果の組み合わせ
固体に作用する総抗力は、次の2つの寄与による流れ方向の合力です。
- 圧力抗力または形状抗力：物体周囲の圧力が不均一になることによって生じます。
- 表面摩擦抗力：流体・固体界面における粘性せん断によって生じます。
抗力係数は、これらの効果を1つの無次元量にまとめたものです。その値は流動状態および固体または装置の幾何形状に依存します。
基準となる式は慎重に定義する必要がある
標準的な定義は次のとおりです。
[ c_D=\frac{D}{\tfrac{1}{2}\rho V^2S} ]
基準面積（S）には、通常、単独粒子の投影正面面積を用います。直径(d_p)の球形粒子の場合、
[ S=\frac{\pi d_p^2}{4} ]
速度(V)も明確に定義する必要があります。これは、単独粒子周囲の自由流速度、充填層を通過する見かけ速度、またはより詳細な多相流モデルにおける局所速度／間隙速度の場合があります。
パイロットプラントで(c_D)を決定する方法
力またはその水力学的等価量を測定する
直接実験では、パイロットプラントに既知の流体速度を与え、その結果生じる抗力を測定します。力は天びんまたはロードセルで測定できますが、多くの単位操作では、圧力損失または力の釣り合いから抗力を間接的に求めます。
充填層では、層の両端で測定される圧力損失は、単独粒子1個に作用する抗力ではなく、多数の粒子による総合的な抵抗を表します。したがって、得られたデータは単独粒子の(c_D)だけでなく、充填層の圧力損失モデルを用いて解釈されます。
関連する流体物性を確定する
実験には、流体密度(\rho)および粘度(\mu)、または動粘度(\nu=\mu/\rho)が必要です。これらの物性値は、実際の運転温度および組成において評価する必要があります。
非ニュートン流体、ガス混合物、懸濁液、または圧縮性が大きい系では、選択した物性モデルが計算されるレイノルズ数および得られる抗力相関式に大きな影響を与える可能性があります。
適切なレイノルズ数を計算する
単独粒子の場合、粒子レイノルズ数は一般に次のように表されます。
[ Re_p=\frac{\rho V d_p}{\mu} =\frac{Vd_p}{\nu} ]
ここで、(d_p)は代表粒径、(V)は流体と粒子の相対速度です。
装置規模の流れでは、技術者は塔径、水力直径、または代表的な流路寸法に基づく装置レイノルズ数を用いることもあります。適切な選択は、粒子抗力、層抵抗、または装置全体の流体力学のどれを特徴づけるかによって決まります。
(c_D)と(Re)の関係を構築する
測定された抗力を、定義式を用いて(c_D)に変換します。さまざまな速度範囲で試験を繰り返すと、次のような関係が得られます。
[ c_D=f(Re) ]
非拘束流れ中にある幾何学的に相似な単独物体では、レイノルズ数が支配的な変数となることが多いです。しかし実際のパイロットプラントでは、この関係は粒子の真球度、表面粗さ、濃度、空隙率、配向、および閉塞の影響にも依存します。
流体・固体単位操作における(c_D)の適用
充填層および充填吸収塔
充填吸収塔または充填反応器では、流体が充填要素間の空隙を通過します。その結果生じる抗力は、巨視的には層内の圧力損失として現れます。
単独粒子の抗力係数はメカニズムの説明に役立ちますが、充填層の設計では通常、粒径、見かけ速度、層の空隙率、流体物性などを含む相関式のような、層規模の圧力損失式を用います。パイロットプラントでは、実際の充填物および運転範囲に対してその相関式を検証するために必要なデータを取得します。
したがって、係数は間接的に適用されます。すなわち、粒子スケールの抵抗をモデル内に表現し、充填塔全体の圧力損失を予測します。
流動層
流動層では、上向きの流体抗力が粒子重量および粒子間力に対抗します。速度が増加すると層は膨張し、最終的に最小流動化状態に達することがあります。
流動化開始時には、層全体の圧力損失は、単位断面積当たりの固体の有効重量とおおむね関係します。抗力相関式は、この状態に近づくために必要な速度の決定や、流動化前後における圧力損失の変化予測に役立ちます。
流動層では、(c_D)だけでは不十分です。粒度分布、固体密度、層の空隙率、粒子形状、および固体循環が観測される挙動に影響します。
粒子沈降塔
沈降塔では、粒子は抗力が有効重量と釣り合うまで加速します。静止流体中の単独粒子が終端沈降速度に達したとき、
[ D=(\rho_p-\rho_f)V_p g ]
となります。ここで、(\rho_p)および(\rho_f)は粒子密度と流体密度、(V_p)は粒子体積、(g)は重力加速度です。
この力の釣り合いを抗力の定義式に代入することで、技術者は測定された終端速度から(c_D)を決定できます。逆に、既知の、または相関式で求めた(c_D)を用いて沈降速度を推定できます。
濃厚懸濁液では、妨害沈降および粒子間相互作用を補正する必要があります。希薄でない懸濁液に、単独粒子の係数を修正せずに適用してはなりません。
空気輸送およびスラリー輸送
輸送システムでは、抗力が粒子の気体または液体速度の変化への応答速度を決定します。抗力は、粒子の加速、すべり速度、懸濁状態、および固体輸送に必要な圧力に影響します。
適切な速度は、必ずしも管内のバルク速度ではなく、流体と粒子の相対速度です。管内では、流れ場、粒子濃度、衝突、および乱流により、実効抗力が単独粒子の場合とは大きく異なる可能性があります。
パイロットプラントのスケールアップにおける相似則の利用
支配的な無次元群を一致させる
パイロットプラントが有用なのは、工業規模の寸法を完全に再現しなくても、支配的な無次元関係を再現できるためです。抗力が支配的な系では、通常、レイノルズ数が最初に検討すべき相似パラメータです。
パイロットと実機の幾何形状が相似で、レイノルズ数も同程度であれば、流動パターンと抗力挙動も類似する可能性があります。この原理により、パイロットでの測定結果を、実機の圧力損失および固体運動の予測モデルの構築または検証に利用できます。
レイノルズ数だけでは不十分な場合を認識する
流体・固体操作は、重力、浮力、粒子慣性、および相間相互作用にも支配される場合があります。関連する無次元群には、次のものがあります。
- 粒子レイノルズ数：粒子周囲の粘性効果および慣性効果を表します。
- アーチメデス数：粒子と流体の密度差、重力、粒径、および粘度の複合効果を表します。
- フルード数：自由表面または粒子運動を重力と慣性の効果が支配する場合に用います。
- ペクレ数：対流輸送と分散が競合する場合に用います。
- シャーウッド数、ビオ数、およびスタントン数：固体・流体界面での物質移動または熱移動が重要な場合に用います。
例えば、(Re_p)を一致させることで局所的な抗力挙動は再現できても、重力に関する相似性が大きく異なる場合、沈降または流動化を再現できない可能性があります。
パイロットデータは単なるスケール換算ではなく、検証に利用する
パイロットプラントの最も有効な活用方法は、盲目的な幾何学的スケールアップではありません。測定した圧力損失、速度分布、層膨張、沈降速度、および固体ホールドアップを、物理的に適切なモデルと比較することです。
モデルが想定する運転範囲全体でパイロットの挙動を再現できれば、明示した仮定と不確かさのもとで実機性能を推定できます。
パイロットプラントデータの測定と解釈
速度を一貫して定義する
よくある誤差の原因は、抗力の式に誤った速度を使用することです。充填層および流動層では、報告上は見かけ速度が便利ですが、空隙を通過する実際の流体速度はそれより高くなります。
選択した速度は相関式と一致していなければなりません。見かけ速度、間隙速度、局所速度、相対速度を混同すると、測定自体が正しくても、係数が矛盾しているように見える可能性があります。
粒子スケール量と層スケール量を区別する
単独粒子の(c_D)は、指定された条件下で1つの物体に作用する抵抗を表します。充填層の圧力損失係数は、多孔質集合体全体の総合抵抗を表します。
これらの量には関係がありますが、互換性はありません。層の空隙率、粒子配列、接触点、チャネリング、および壁面効果によって、個々の粒子抗力と層全体の圧力損失との関係が変化します。
信頼性の高い面積と粒子寸法を使用する
非球形粒子では、投影面積は配向によって変化します。技術者は等価直径、真球度補正、または測定した投影面積を使用できますが、その基準を明示する必要があります。
粒度分布がある場合、1つの代表径だけでは固体の抗力挙動を十分に表せないことがあります。再現性を確保するには、粒径の基準および形状特性を報告することが不可欠です。
トレードオフを理解する
完全な動的相似は実現困難な場合がある
パイロットプラントと工業装置の間ですべての無次元群を一致させることは、多くの場合不可能です。装置寸法、利用可能な流体、粒径、運転圧力などにより、レイノルズ数、フルード数、重力パラメータを同時に一致させられないことがあります。
実務的には、支配的なメカニズムを特定し、最も重要な無次元群を一致させ、一致しない無次元群の影響を定量化します。
単独粒子相関式の適用範囲には限界がある
単一の滑らかな球について導出された相関式は基準値として有用ですが、粗い粒子、多孔質粒子、不規則粒子、または凝集する粒子には不正確な場合があります。また、濃厚な多相流を自動的に表現できるわけでもありません。
したがって、粒子形状、濃度、または相互作用が重要な場合、パイロット試験では実際の固体および流体系を反映させる必要があります。
圧力損失データだけでは不均一流を見逃す可能性がある
単一の圧力損失測定から得られるのは全体的な抵抗であり、チャネリング、流量分布不良、バイパス、または局所的な流動化を明らかにできない場合があります。2つのシステムで平均圧力損失が類似していても、内部の流動構造は大きく異なる可能性があります。
測定された抗力を正しく解釈するには、速度分布、層膨張、固体ホールドアップ、および目視またはトレーサーによる診断が必要になる場合があります。
スケールアップによって流動状態が変化する可能性がある
パイロットプラントと実機では、異なるレイノルズ数または粒径と装置寸法の比で運転される場合があります。これにより、系が層流、遷移流、乱流の間で移行し、適用される抗力関係が変化する可能性があります。
したがって、スケールアップ計算では、どこでも1つの一定の(c_D)が適用できると仮定せず、想定される流動状態を確認する必要があります。
プロジェクトへの適用方法
次のワークフローに従うことで、抗力係数を物理的に意味のある、スケールアップに有用なものとして扱えます。
- 主な関心が充填層の圧力損失である場合：さまざまな見かけ速度で圧力損失を測定し、粒径と空隙率を特性評価したうえで、層を1つの単独粒子として扱うのではなく、充填層相関式を使用します。
- 主な関心が流動化である場合：粒子レイノルズ数と最小流動化挙動を決定し、粒子密度、空隙率、粒度分布、および層膨張を解析に含めます。
- 主な関心が粒子沈降である場合：終端速度を測定し、粒子重量の釣り合いから抗力を計算します。懸濁液が希薄でない場合は、濃度または壁面効果を補正します。
- 主な関心がスケールアップである場合：適切なレイノルズ数を特に重視して、支配的な無次元群を一致させ、得られたモデルをパイロットでの圧力、速度、および固体移動データに対して検証します。
- 主な関心が熱移動または物質移動である場合：(c_D)を用いて流体力学的挙動を確立しますが、抗力だけでは移動性能を決定できないため、ペクレ数、シャーウッド数、ビオ数、または関連する輸送無次元群も評価します。
適切に定義された(c_D)は、単にパイロットプラントで測定された数値ではありません。それは、流体・固体力学、運転条件、および信頼性の高い単位操作のスケールアップを結び付ける、体系的なリンクです。
まとめ表：
| 用途 | cDの決定方法または使用方法 | 主な要因 |
|---|---|---|
| 充填層 | 圧力損失を粒子および層の抵抗と関連付ける | 見かけ速度、空隙率、粒径、壁面効果 |
| 流動層 | 抗力、最小流動化、および層膨張を推定する | 粒子密度、粒度分布、形状、固体循環 |
| 沈降塔 | 終端沈降速度からcDを計算する | 密度差、重力、濃度、壁面効果 |
| 空気輸送またはスラリー輸送 | 粒子加速、すべり速度、および輸送圧力を予測する | 相対速度、乱流、濃度、粒子衝突 |
| パイロットプラントのスケールアップ | 圧力、速度、および固体移動データを用いてモデルを検証する | レイノルズ数、アーチメデス数、フルード数、および輸送相似性 |
| 測定の基準 | 定義した速度と面積を用いてcD = D/(0.5rhoV^2S)を使用する | 流体物性、基準面積、粒子形状、流動状態 |
流体・固体プロセスまたはスケールアップ試験を計画していますか？ LABPARKにお問い合わせください。LABPARKは、化学工学、バイオプロセス・バイオテクノロジー、環境・水処理分野向けに、教育・職業訓練用の単位操作パイロットプラントを提供しています。当社のシステムは、大学、研究機関、企業が、管理された条件下で圧力損失、沈降、流動化、輸送挙動、および関連する熱移動・物質移動性能を測定できるよう支援します。お客様に適したパイロットプラント構成の選定、意味のある運転変数と無次元群の定義、流体力学モデルの検証、さらに教育、研究、プロセス開発、産業スケールアップに活用できる信頼性の高いデータの取得をサポートします。