Here is a gif animated of a simulation with 1 million particles. The initial condition was an isolated, uniform,  prolate ellipsoid with a some rigid rotation. The system is evolved for 50 dynamical times. This is a projection on the XY plane where rotation is around the Z axis. The color code is proportional to the log of the density.

More: Long-lived transient structure in collisionless self-gravitating systems David Benhaiem, Francesco Sylos Labini, and Michael Joyce, Phys. Rev. E 99, 022125, 2019