Recently, the kinetic equation for two-particle distribution function under the Boltzmann approximations was proposed. Presenting scattering operator in divergence form enables to simulate the system of molecules by an equivalent system of quasiparticles. Effect of collisions is accounted by rotating quasiparticle pairs with the angular velocity that depends on the distribution function. In the paper we present new expressions for angular velocity of quasiparticle pair rotation, which are convenient for numerical calculations. We perform test simulations in the differential approximation for scattering operator. The approach expected to be useful in the range between continuum gas dynamics and Boltzmann equation descriptions.