A Precise Computation of Drag Coefficients of a Sphere

We present a computational method for drag coefficients of axisymmetric bodies. It is a kind of consistent flux method but the introduction of a proper test function enables us to establish an error estimate under some assumption. Applying the method, we obtain drag coefficients of a sphere for Reynolds numbers between 10 and 200, which are found between numerical upper and lower bounds.