A new boundary-integral algorithm for the motion of a particle between two parallel plane walls in Poiseuille flow at low Reynolds number was developed to study the translational and rotational velocities for a broad range of particle sizes and depths in the channel. Instead of the free-space Green’s function more commonly employed in boundary-integral equations, we used the Green’s function for the domain between two infinite plane walls [Liron and Mochon, J. Eng. Math. 10, 287 (1976)]. This formulation allows us to directly incorporate the effects of the wall interactions into the stress tensor, without discretizing the bounding walls, and use well-established iterative methods. Our results are in good agreement with previous computations [Ganatos et al., J. Fluid Mech. 99, 755 (1980)] and limiting cases, over their range of application, with additional results obtained for very small particle–wall separations of less than 1% of the particle radius. In addition to the boundary-integral solution in the mobility formulation, we used the resistance formulation to derive the near-field asymptotic forms for the translational and rotational velocities, extending the results to even smaller particle–wall separations. The decrease in translational velocity from the unperturbed fluid velocity increases with particle size and proximity of the particle to one or both of the walls. The rotational velocity exhibits a maximum magnitude between the centerline and either wall, due to the competing influences of wall retardation and the greater fluid velocity gradient near the walls. The average particle velocity for a uniform distribution of particles was generally found to exceed the average fluid velocity, due in large part to exclusion of the particle centers from the region of slowest fluid near the walls. The maximum average particle velocity is 18% greater than the average fluid velocity and occurs for particle diameters that are 42% of the channel height; particles with diameters greater than 82% of the channel height have smaller average velocities than does the fluid, due to the retarding effects of the nearby walls. Additionally, the translational and rotational velocities of oblate and prolate ellipsoids were calculated using the boundary-integral algorithm. The proximity of the walls to the ellipsoids was found to have a strong effect on particle velocity, so that a prolate spheroid aligned on the channel centerline moves faster than an oblate spheroid of the same volume, because the edge of the latter spheroid is closer to the channel walls. For off-centerline locations for a prolate spheroid with its major axis at an angle to the walls, a second translational velocity component normal to the walls is present. For each lower-wall gap examined, the maximum normal translational velocity occurs for smaller gaps from the upper wall (i.e., larger angle). The direction of this velocity component changes sign for mid-range angles studied, due to increased interactions with the lower wall that prevent the ellipsoid from rotating upward and, hence, yield a negative velocity perpendicular to the walls. The rotational velocity changes direction for particles at the smallest angles studied, due to the competition between the Poiseuille fluid velocity profile, which pushes the particle clockwise, and the lubrication forces, which impede this rotation.
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.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.