The effect of contact interaction between particles is decreased by the Riemann solver.Efficiency of the method is demonstrated by numerical tests of a 1D shock tube.