$ C_D= \begin{cases} \frac{24}{Re}(1+0.15Re^{0.687}&Re<1000 \\ 0.44&Re>1000 \end{cases} $ $ Re=\frac{\rho_q|\boldsymbol{U}_p-\boldsymbol{U_q}|d_pj}{\mu_q} $