Terminal-velocity calculations assume a sphere density ρ_p = 2700 kg/m³ (aluminum) and gravitational acceleration g = 9.81 m/s².
When an aluminum sphere (ρ_p = 2700 kg/m³) is released from rest, it accelerates and approaches the terminal velocity v_t. Left: settling motion and force vectors; right: operating point on the C_D–Re curve.
Horizontal axis = Re (log); vertical axis = C_D (log); solid blue = composite C_D(Re) curve; red circle = current operating point; dashed lines = regime boundaries
Flow around a sphere changes substantially with Reynolds number Re, and drag coefficient C_D is represented in four regimes.
Reynolds number. ρ is fluid density, U is flow speed, D is sphere diameter, and μ is dynamic viscosity:
$$Re = \frac{\rho\,U\,D}{\mu}$$Composite drag-coefficient equation (simplified model by regime):
$$C_D = \begin{cases} 24/Re & (Re \lt 0.1) \\ \dfrac{24}{Re}\bigl(1 + 0.15\,Re^{0.687}\bigr) & (0.1 \le Re \lt 10^3) \\ 0.44 & (10^3 \le Re \lt 2\times 10^5) \\ 0.10 & (Re \ge 2\times 10^5) \end{cases}$$Drag F_D. A = πD²/4 is the projected area:
$$F_D = \tfrac{1}{2}\,C_D\,\rho\,U^2\,A$$Terminal velocity U_t under gravity (particle density ρ_p; g is gravitational acceleration):
$$U_t = \sqrt{\tfrac{4}{3}\,\frac{(\rho_p-\rho)\,g\,D}{\rho\,C_D}}$$Because C_D depends on U_t itself, U_t is found iteratively (this tool uses five iterations). Near Re ≈ 2×10⁵, the boundary layer becomes turbulent and C_D drops sharply in the “drag crisis.”