Particulate Flow Simulation

Two neutrally buoyant prolate ellipsoids interact in a bounded shear flow: Passing over/under trajectories

Two neutrally buoyant prolate ellipsoids pass over/under each other in a bounded shear flow of a Newtonian fluid under creeping conditions. The mass centers return to their initial heights, respectively, after passing each other.

The dimension of computational domain = (-1.5, 1.5)x(-1,1)x(-0.5,0.5).
The semi-axes of the prolate ellipsoids are 0.2, 0.1, and 0.1.
The density of the balls = 1.
The density of fluid = 1.
The viscosity = 1.
The shear rate = 1.
Two walls move in opposite directions with the same speed.
the mass centers are (-0.5, 0, 0.05) and (0.5, 0, -0.05).
The mesh size for the velocity field = 1/48
The mesh size for the pressure = 1/24
The time step=0.001.