Shock and pulsating bubble, as two processes for structure damage in underwater explosion, are often studied separately due to their different time scales. In this study, we would pay particular attention to the interaction between them in underwater explosion nearby a wall by finite volume method. In order to capture the shock precisely, it is an issue at priority how to properly select grid density and artificial viscosity. Then a “two-step” strategy is adopted to overcome the bottleneck of high computing requirements in CPU time and memory. The shock reflection, bubble oscillation and the interaction between them along with resultant load are presented and discussed.