In this paper, the process of the movement of a direct shock wave from a pure gas into a dusty medium is numerically modeled. The mathematical model took into account the viscosity, compressibility and thermal conductivity of the carrier phase. Also, the modeling technique made it possible to describe the interphase force interaction, which included the Stokes force, the dynamic force of Archimedes, the strength of the attached masses. In addition, interfacial interaction included heat transfer between the carrier and dispersed phases. The numerical solution was carried out using the explicit finite-difference method, with the subsequent application of the nonlinear correction scheme for the grid function. As a result of numerical calculations, it was revealed that with an increase in the linear particle size of the gas suspension, the velocity slip between the carrier and dispersed phases increases. Numerical modeling also showed that the absolute value of the difference between the velocities of the carrier and the dispersed phase reaches the largest value at the leading edge of the compression wave. The revealed regularities can be explained by the fact that the particles of the dispersed phase are assumed to be spherical in shape. Due to this, a multiple increase in particle size leads to a three-fold increase in their mass, a twofold increase in the area of one particle and a three-fold decrease in the number of particles. Thus, an increase in particle size leads to a decrease in the area of interfacial contact and an increase in the inertia of the particles, which in turn affects the interfacial velocity slip.