A flow containing multiple particles and the shock wave is investigated by the direct numerical simulation with immersed boundary method. The shock Mach number and the Reynolds numbers of particle behind the shock wave are set to be 1.5 to 2.0 and 300 to 600, respectively. The comparison of the present results with one-dimensional simulation results, shows good agreement. From the results, we clarified characteristic flow structure at different shock Mach and Reynolds number. The turbulence kinetic energy was enhanced from the vortex structure in the wake of particles for the high Reynolds number case. The drag coefficient from the present simulation and the previous prediction models shows almost the same values at Mach number 1.5. At Mach number 2.0, however, discrepancy is obtained for the drag coefficient between the present flow simulation and the previous prediction models.