Drag Coefficients of Viscous Spheres at Intermediate and High Reynolds Numbers

A finite-difference scheme is used to solve the Navier-Stokes equations for the steady flow inside and outside viscous spheres in a fluid of different properties. Hence, the hydrodynamic force and the steady-state drag coefficient of the spheres are obtained. The Reynolds numbers of the computations range between 0.5 and 1000 and the viscosity ratio ranges between 0 (inviscid bubble) and infinity (solid particle). Unlike the numerical schemes previously implemented in similar studies (uniform grid in a stretched coordinate system) the present method introduces a two-layer concept for the computational domain outside the sphere. The first layer is a very thin one $[ORe−1/2]$ and is positioned at the interface of the sphere. The second layer is based on an exponential function and covers the rest of the domain. The need for such a double-layered domain arises from the observation that at intermediate and large Reynolds numbers a very thin boundary layer appears at the fluid-fluid interface. The computations yield the friction and the form drag of the sphere. It is found that with the present scheme, one is able to obtain results for the drag coefficient up to 1000 with relatively low computational power. It is also observed that both the Reynolds number and the viscosity ratio play a major role on the value of the hydrodynamic force and the drag coefficient. The results show that, if all other conditions are the same, there is a negligible effect of the density ratio on the drag coefficient of viscous spheres.