This study is devoted to investigating a flow around a stationary or moving sphere by using direct numerical simulation with immersed boundary method (IBM) for the three-dimensional compressible Navier-Stokes equations. A hybrid scheme developed to solve both shocks and turbulent flows is employed to solve the flow around a sphere in the equally spaced Cartesian mesh. Drag coefficients of the spheres are compared with reliable values obtained from highly accurate boundary-fitted coordinate (BFC) flow solver to clarify the applicability of the present method. As a result, good agreement was obtained between the present results and those from the BFC flow solver. Moreover, the effectiveness of the hybrid scheme was demonstrated to capture the wake structure of a sphere. Both advantages and disadvantages of the simple IBM were investigated in detail.