終端速度の計算では球密度 ρ_p = 2700 kg/m³(アルミニウム)、重力加速度 g = 9.81 m/s² を仮定しています。
横軸=Re(対数)/縦軸=C_D(対数)/青実線=C_D(Re) 合成曲線/赤丸=現在の動作点/破線=各域の境界
球周りの流れは Reynolds 数 Re によって性質が大きく変わり、抗力係数 C_D は 4 つの領域に分けて表されます。
Reynolds 数。ρ は流体密度、U は流速、D は球直径、μ は動的粘度:
$$Re = \frac{\rho\,U\,D}{\mu}$$抗力係数の合成式(領域別の簡易モデル):
$$C_D = \begin{cases} 24/Re & (Re < 0.1) \\ \dfrac{24}{Re}\bigl(1 + 0.15\,Re^{0.687}\bigr) & (0.1 \le Re < 10^3) \\ 0.44 & (10^3 \le Re < 2\times 10^5) \\ 0.10 & (Re \ge 2\times 10^5) \end{cases}$$抗力 F_D。A = πD²/4 は投影面積:
$$F_D = \tfrac{1}{2}\,C_D\,\rho\,U^2\,A$$重力下落の終端速度 U_t(粒子密度 ρ_p、g は重力加速度):
$$U_t = \sqrt{\tfrac{4}{3}\,\frac{(\rho_p-\rho)\,g\,D}{\rho\,C_D}}$$U_t は C_D が U_t 自身に依存するため、反復計算で求めます(本ツールは 5 回反復)。Re ≈ 2×10⁵ 付近では境界層が乱流化して C_D が急減する「ドラッグクライシス」が起こります。