This paper is devoted to the development and quantitative evaluation of the properties of the numerical algorithm of the Cartesian grid method for three-dimensional (3D) simulation of shock-wave propagation in areas of varying shape. The detailed description of the algorithm is presented. The algorithm is relatively simple to implement and does not require solving the problem of determination of the shape of the body’s boundary intersection with regular computational cell. The accuracy of the algorithm is demonstrated by comparing the simulated and experimental data in the problems of the interaction of a shock wave (SW) with a nonmoving sphere and a moving particle.