Moneytizerによる広告を掲載しています。粒子が大きくなると、流体中の粒子は慣性の効果により流体と同じ動きができなくなります。粒子を含む流体中にフィルター等の障害物を設置すると、流体は障害物を避けて流れますが、粒子は障害物を避けきれずに衝突し捕集されます。このような慣性力を利用した捕集機構を慣性捕集機構といいます。慣性捕集は数ある捕集機構の中でも主要なものの1つです。よる粒子の捕集は、基本的に粒子の運動方程式を解くことで求められます。粒子レイノルズ数Rep≦2のストークス則が成り立つ領域では、$$F_{d}=3πμd_{p}(u-v)・・・(1)$$Fd:流体の抵抗力[N]、μ:流体粘度[Pa・s]、dp:粒子径[m]u:流体速度[m/s]、v:粒子速度[m/s](1)式が成立します。(u-v)は粒子の相対速度を表わします。粒子が大きくなるとストークス則で計算される抵抗力にずれが生じますから、通常はカニンガムの補正係数で(1)式を補正します。$$F_{d}=\frac{3πμd_{p}(u-v)}{C_{c}}・・・(2)$$Cc:カニンガムの補正係数[-](2)式を利用して、粒子の運動方程式を立てると、$$m_{p}\frac{dv}{dt}=\frac{3πμd_{p}(u-v)}{C_{c}}+F_{e}・・・(3)$$mp:粒子重量[kg]、Fe:粒子の受ける外力[N](3)式となります。次に(3)式を(4)~(6)式を使用して無次元化します。$$v_{n}=\frac{v}{u_{0}}・・・(4)$$$$u_{n}=\frac{u}{u_{0}}・・・(5)$$$$t_{n}=\frac{u_{0}t}{d_{f}}・・・(6)$$vn:無次元粒子速度[-]、un:無次元流体速度[-]、tn:無次元時間[-]u0:流れの平均速度[m/s]、d_{f}:繊維径[m]また、粒子重量mpは、$$m_{p}=ρ_{p}・\frac{4}{3}π(\frac{d_{p}}{2})^3=\frac{π}{6}ρ_{p}{d_{p}}^{3}・・・(7)$$ρp:粒子密度[kg/m3](7)式で表されます。(4)~(7)式を(3)式に代入すると、$$\frac{π}{6}ρ_{p}{d_{p}}^{3}\frac{d(u_{0}v_{n})}{d(d_{f}t_{n}/u_{0})}=\frac{3πμd_{p}(u_{0}u_{n}-u_{0}v_{n})}{C_{c}}+F_{e}$$$$\frac{C_{c}ρ_{p}{d_{p}}^{2}u_{0}}{18μd_{f}}\frac{dv_{n}}{dt_{n}}=(u_{n}-v_{n})+\frac{F_{e}}{u_{0}}\frac{C_{c}}{3πμd_{p}}$$$$ψ\frac{dv_{n}}{dt_{n}}=(u_{n}-v_{n})+\frac{F_{e}}{u_{0}}\frac{C_{c}}{3πμd_{p}}・・・(8)$$$$ψ=\frac{C_{c}ρ_{p}{d_{p}}^{2}u_{0}}{18μd_{f}}・・・(9)$$ψ:慣性パラメーター[-](8)式となります。(8)式中のψは慣性パラメーターと呼ばれており、$$Stk=2×ψ・・・(10)$$慣性パラメーターの2倍がストークス数Stkとなります。粒子の捕集効率は、理論的には(8)式の粒子の運動方程式を数値的に解くことで求められます。(8)式を解くのは大変です。ここでは、慣性捕集効率ηI[-]のいくつか有名な推算式を紹介します。Stechkinaらの式$$η_{I}=\frac{Stk}{4{K_{u}}^{2}}(29.6-28α^{0.62})R^2-27.5R^{2.8}・・・(11)$$$$K_{u}=-0.5{\rm{ln}}α+α-0.25α^2-0.75・・・(12)$$$$R=\frac{d_{p}}{d_{f}}・・・(13)$$$$Stk=\frac{C_{c}ρ_{p}{d_{p}}^{2}u_{0}}{9μd_{f}}・・・(14)$$Ku:水力学因子[-]、α:フィルターの充填率[-]、R:さえぎりパラメーター[-]Stechkinaらは慣性効果が比較的小さい領域で成り立つ式として、(11)式を発表しています。0.0035<α<0.111、0.01<R<0.4が適用範囲とされています。α=0.01、Rを3×10-5<R<3.0まで変化させたときのグラフを以下に示します。適用範囲外の領域で計算すると、慣性捕集効率ηIが負の値になったり、100%を超えたりするので気を付けましょう。Hermannの式$$η_{I}=\frac{Stk^3}{Stk^3+1.54Stk^2+1.76}・・・(15)$$$$Stk=\frac{C_{c}ρ_{p}{d_{p}}^{2}u_{0}}{9μd_{f}}・・・(16)$$$$Re_{f}=\frac{ρd_{f}u}{μ}・・・(17)$$LandahlとHermannは繊維レイノルズ数Ref=10で、(15)式を発表しています。Ref=10のときのグラフを以下に示します。式の形として、ηIが0~100%の間に収まる点は良いですが、Ref≒10でどれほどの推算精度を保てるかはわかりません。高橋の式$$η_{I}={\rm{exp}}(-0.462α^{-0.5}Stk^{-1.0})・・・(18)$$$$Stk=\frac{C_{c}ρ_{p}{d_{p}}^{2}u_{0}}{9μd_{f}}・・・(19)$$α:フィルターの充填率[-]高橋は0.005≦α≦0.1の範囲で慣性捕集効率ηIの線図をまとめており、定式化したものが(18)式です。α=0.01、Ref=10のときのグラフを以下に示します。Hermannの式と似たような値を示します。まとめ慣性捕集機構について解説しました。この記事に示したもの以外にも推算式はたくさんあります。適用範囲をよく確認して使用するようにしましょう。コメントをキャンセルメールアドレスが公開されることはありません。※サイト次回のコメントで使用するためブラウザーに自分の名前、メールアドレス、サイトを保存する。date*日本語が含まれない投稿は無視されますのでご注意ください。定圧変化における系のエネルギーや熱量を表わす指標をエンタルピーといいます。定圧変化においては、状態量ではない熱量Qを状態量であるエンタルピーの変化ΔHから算出することができます。棚段塔に使用される液を保持するための板をトレイ(棚段)といいます。トレイ上に保持した液と塔底から登ってくるガスを接触させるために、トレイにはガスの通り穴が多数空いています。流体が壁面との摩擦や流体同士の摩擦によってエネルギーを損失することを圧力損失といいます。この記事では実務で計算する頻度の多い配管流路の圧力損失について計算方法を紹介します。【Fickの法則】を解説:物質移動・拡散の基礎式Fickの法則は物質拡散に関する法則です。ガス吸収におけるモデルの1つである二重境膜説でも、ベースとなっている法則なので覚えておきましょう。撹拌槽内の流動は撹拌翼のフローパターンによって変化するため、槽壁の境膜伝熱係数もそれに応じて変化します。先人達は境膜伝熱係数を比較的簡単な相関式で求められるよう研究してきました。以下ではその相関式について紹介します。【ニュートンの粘性法則】を解説:ニュートン流体・非ニュートン流体の違いとは?現実の流体は粘性を持っており、流体が動くときには相互運動を妨げる方向に力が働きます。【さえぎり捕集機構】について解説:障害物による粒子のさえぎりプロフィールこのブログを運営しているルートです。化学メーカーで働くケミカルエンジニアです。