“…Our procedure, based on a particle method and a splitting strategy, is an evolution of the approach introduced in [6]. The substantial difference of our problem from the one studied in [6] is that the dust particles in planetary rings, unlike gas molecules, being affected by the gravitational acceleration due to the planet and moons, satisfy the Vlasov equation ( 1), which is more complicated to deal with than the free transport equation. On the other hand, the boundary condition are, in this problem, particularly easy to handle, because the surfaces of the satellites are modelled with an absorbing boundary: the particles that collide with it disappear from the problem's domain.…”