Snapshot of the simulation (random initial positions) for different area fractions ϕ at dimensionless time t/t^=106 . Note that the image shows only the center 44a × 44a region of the simulation domain. The black lines on the particle indicate the instantaneous velocity direction.