終端沈降速度シミュレーター 戻る
流体力学シミュレーター

終端沈降速度シミュレーター — 球形粒子の沈降と抗力係数

球形粒子が流体中を重力沈降する終端速度を Reynolds 域別に反復計算。粒子径・密度・粘度を変えて、Stokes 域・中間域・Newton 域での沈降挙動と抗力係数の遷移を直感的に学べます。

パラメータ設定
粒子径 D
mm
粒子密度 ρ_p
kg/m³
流体密度 ρ_f
kg/m³
流体粘度 μ
Pa·s

重力加速度 g = 9.81 m/s²。反復は最大 20 回、|ΔU/U|<1e−4 で収束判定。

計算結果
終端沈降速度 U_t
終端 Reynolds 数 Re_t
抗力係数 C_D
流れ域
Re ― U_t 関係(等粒径理論線)

横軸 = Reynolds 数(対数)/縦軸 = 終端速度 U_t(対数)/灰線 = D=0.01,0.1,1,10mm の理論線/黄点 = 現在の (Re, U_t)/縦点線 = 域境界 Re=0.1, 1000

理論・主要公式

球形粒子の重力沈降では、実効重量と流体抗力がつり合うとき終端速度に達します。

力のつり合い(重力 − 浮力 = 抗力)から得られる一般式:

$$U_t = \sqrt{\dfrac{4\,g\,D\,(\rho_p-\rho_f)}{3\,\rho_f\,C_D}}, \qquad Re=\dfrac{\rho_f\,U_t\,D}{\mu}$$

抗力係数 C_D の Reynolds 域別表現:

$$C_D = \begin{cases} 24/Re & (Re<0.1)\\ \dfrac{24}{Re}\,(1+0.15\,Re^{0.687}) & (0.1\le Re<1000)\\ 0.44 & (1000\le Re<2\times10^5) \end{cases}$$

Stokes 域では解析解が直接得られます:

$$U_t^{\text{Stokes}} = \dfrac{g\,D^2(\rho_p-\rho_f)}{18\,\mu}$$

中間域では C_D が U_t に依存するため反復計算(Re→C_D→U_t を更新)で解きます。

終端沈降速度シミュレーターとは

🙋
水の中に砂を落とすと、最初は速くなって、途中から一定の速さで沈んでいきますよね?あれって何が起きているんですか?
🎓
ざっくり言うと、重力で加速していくうちに、流体抗力がだんだん大きくなって、最後は重力(から浮力を引いたもの)と抗力がぴたっとつり合うんだ。そのつり合った瞬間からは加速がゼロ、つまり一定速度——これが「終端沈降速度 U_t」だよ。上のシミュレーターで「砂/水」プリセットを押すと、D=1mm の砂粒が水中で約 0.14 m/s で沈むのが見える。
🙋
じゃあ、粒子が小さい時と大きい時で、式も変わるんですか?右のグラフが折れ曲がっているのが気になって…
🎓
いいところに気づいたね。粒子周りの流れの状態は Reynolds 数 Re=ρUD/μ で決まる。Re<0.1 の遅い流れ(Stokes 域)では、粘性が支配的で C_D=24/Re。式は $U_t=gD^2\Delta\rho/(18\mu)$ になり、なんと径の2乗に比例するんだ。逆に Re>1000 の Newton 域では C_D≈0.44 で一定になり、$U_t \propto \sqrt{D}$ で平方根でしか効かなくなる。
🙋
でも、さっきの砂粒は Re=138 でしたよね。中間域だと、どうやって計算するんですか?
🎓
そこが面白いところで、中間域では C_D が U_t に依存し、U_t を求めるのに C_D が必要、という「ニワトリと卵」状態なんだ。だから反復で解く。Stokes 解で初期推定 → Re を計算 → C_D を計算 → U_t を更新、を繰り返す。本ツールでは |ΔU/U|<1e−4 を満たすまで最大 20 回回している。試しに「鋼球/油」プリセットを押すと、5mm の鋼球が粘い油の中で何 m/s で沈むか出るよ。粘度を上げると Re が落ちて Stokes 域に近づくのも面白い。

よくある質問

本ツールは球形を仮定しています。実際の粉体や砂の粒は不規則な形状で、抗力係数は同じ Re でも球より大きくなります。実務では「形状係数」や「等積球径」「等速球径」などの補正を施すか、Haider-Levenspiel 相関のような非球形相関式を使います。設計では本ツールの値を上限の目安とし、安全側に補正してください。
本ツールは無限希釈(粒子が単独で沈む)を仮定しています。粒子濃度が高くなると上昇する流体の戻り流れや、粒子間相互作用で見かけの沈降速度が低下します。これを「干渉沈降」と呼び、Richardson-Zaki の式 U_int=U_t·(1−φ)^n(φ=粒子体積分率、n は Re 依存)で補正します。沈殿池や濃縮槽の設計では干渉項が支配的です。
ρ_p < ρ_f のときは粒子は浮力で上昇します。本ツールは速度の絶対値で表示し、Reynolds 域の判定だけ行います(油中の気泡やプラスチック片の上昇速度など)。なお、気泡は表面張力や変形で球形を保てない場合が多く、Re が大きい領域では Hadamard-Rybczynski や Mendelson の波動理論など別のモデルが必要です。
Re が約 1000 を超えると、球の後流に大規模な乱れた剥離渦が発生し、その圧力差で抗力の大半が決まるようになります。粘性応力の寄与は相対的に小さくなり、流れのパターンが Re にあまり依存しなくなるため C_D≈0.44 でほぼ一定になります。さらに Re≈2×10^5 を超えると境界層が乱流に遷移し、剥離点が後方に移動して C_D が急激に低下する「ドラッグクライシス」が起こります。

実世界での応用

上下水道の沈殿池設計:浄水処理や下水処理の沈殿池では、原水中の浮遊物質(SS)の沈降速度が表面負荷率(流量/池表面積)の上限を決めます。本ツールの U_t を Stokes 域で計算し、目標除去粒径から池サイズを設計します。実際は凝集沈殿で粒子径を大きくし、より速い沈降を狙います。

サイクロン分離器・遠心分級:気流中の粉塵をサイクロンで分離する際、遠心場での終端速度(重力 g を遠心加速度 ω²r に置き換えたもの)の差で粒径を分級します。粉体プロセスや空気清浄機で広く使われ、設計には粒子の Cd-Re 関係が不可欠です。

気象・降雨物理:雨滴の終端速度は粒径とともに増加し、直径 1mm の小雨で約 4 m/s、5mm の大粒では約 9 m/s に達します(変形効果のため理論より低くなる)。降雨強度の推定や、レーダー反射エコーから粒径分布を逆推定するアルゴリズムの基礎になります。

化学プロセス・流動層:固体触媒の流動層リアクタでは、粒子の終端速度より速いガス流速で粒子が舞い上がり「飛散」します。安定運転には U_mf(流動化開始速度)と U_t の間にガス速度を保つ必要があり、本ツールの U_t は飛散開始の上限速度として使われます。

よくある誤解と注意点

最も多い誤解は、「沈降速度は粒径に比例する」と単純に考えてしまうことです。実際には Stokes 域では D の2乗、Newton 域では D の平方根、中間域はその中間の指数で増加します。シミュレーターで粒径スライダーを 0.001mm から 100mm まで動かしながら U_t と Re を見てください。径を 10 倍にしたとき U_t が何倍になるかが域によって全く違うことが体感できます。設計で「径を倍にすれば速度も倍」と決めつけると、特に砂粒以上のサイズでは大きく外します。

次に多いのが、Stokes 解だけで全粒径をカバーしようとするミスです。Stokes 式は Re<0.1 でしか使えません。シミュレーターの初期値(砂 1mm/水)でも Stokes 推定の Re は 818 と完全に範囲外で、もし Stokes 式で押し通すと U_t≈0.82 m/s と現実の約 6 倍を出してしまいます。中間域・Newton 域では必ず C_D を Re で更新する反復解を使ってください。本ツールは自動的に反復しているため、Re カードが域境界を越えたら式が切り替わっていることを意識すると良いでしょう。

最後に、「終端速度は瞬時に達成される」と思い込む点に注意が必要です。実際には粒子が静止状態から U_t の 99% に達するまでに有限の時間と距離が必要で、Stokes 域の小さな粒子(μm オーダー)では数ミリ秒、Newton 域の大粒子(cm オーダー)では数秒かかります。沈殿池や落下塔の設計では、粒子が U_t に到達する加速距離が無視できない場合があり、特に短い容器では助走区間の補正が必要になります。