重力加速度 g = 9.81 m/s²。反復は最大 20 回、|ΔU/U|<1e−4 で収束判定。
横軸 = 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 を更新)で解きます。