Dynamic simulation of spheroid motion between two parallel plane walls in low-Reynolds-number Poiseuille flow

A novel boundary-integral algorithm is used to study the general, three-dimensional motion of neutrally buoyant prolate and oblate spheroids in a low-Reynolds-number Poiseuille flow between parallel plates. Adaptive meshing of the spheroid surface assists in obtaining accurate numerical results for particle–wall gaps as small as 1.3% of the spheroid's major axis. The resistance formulation and lubrication asymptotic forms are then used to obtain results for arbitrarily small particle–wall separations. Spheroids with their major axes shorter than the channel spacing experience oscillating motion when the spheroid's centre is initially located in or near the midplane of the channel. For both two-dimensional and three-dimensional oscillations, the period length decreases with an increase in the initial inclination of the spheroid's major axis with respect to the lower wall. These spheroids experience tumbling motions for centre locations further from the midplane of the channel, with a period length that decreases as the spheroid is located closer to a wall. The transition from two-dimensional oscillating motion to two-dimensional tumbling motion occurs for an initial centre location closer to a wall as the initial inclination of the major axis is increased. For these spheroids, the average translational velocity along the channel length for two-dimensional oscillating motion decreases for an increase in the initial inclination of the major axis, and the average translational velocity for two-dimensional tumbling motion decreases for a decrease in the initial centre location. A prolate spheroid with its major axis 50% longer than the channel spacing and confined to the ($x_2$, $x_3$)-plane (where $x_2$ is the primary flow direction and $x_3$ is normal to the walls) cannot experience two-dimensional tumbling; instead, the spheroid becomes wedged between the walls for initial centre locations near the midplane of the channel when the initial inclination of the large spheroid's major axis is steep, and experiences two-dimensional oscillations for initial centre locations near a wall. When this spheroid's major axis is not confined to the ($x_2$, $x_3$)-plane, it experiences three-dimensional oscillations for initial centre locations in or near the midplane of the channel, and three-dimensional tumbling for initial centre locations near a wall.