Particulate Flow Simulation

The motion of a neutrally buoyant prolate in a tube Poiseuille flow

The motion of a neutrally buoyant prolate in a pressure driven tube Poiseuille flow for 0 < t < 15:
The long axis is initially parallel to the central axix of the cylinder.

The flow is periodic in the horizontal direction.
The length of the cylinder is 2 and its radius is R=0.5.
The dimension of computational domain = (0,1+4h)x(0,2)X(0,1+4h).
The initial velocity field is a fully developed tube Poiseuille flow with maximal speed of 10.
Three semi-axes of the prolate are 0.065, 0.065, 0.195.
The density of tha prolate is 1.
The density of fluid is 1.
The viscosity is 0.1
The pressure drop is 16 per unit length.
The mesh size for the velocity field is h=1/80.
The mesh size for the pressure is 1/40.
The time step=0.001.
The distance from the mass center to the central axis is 0.2728 (0.5456R).
The averaged translation speed is 6.732 for 215 < t < 225.
The averaged particle Reynolds number 26.2548 for 215 < t < 225.