Display options
Cite
Dvinsky AS, Popel AS. Motion of a Rigid Cylinder Between Parallel Plates in Stokes Flow: Part 1: Motion in A Quiescent Fluid and Sedimentation. Comput Fluids. 1987;15(4):391-404doi: 10.1016/0045-7930(87)90031-4.
Dvinsky, A. S., & Popel, A. S. (1987). Motion of a Rigid Cylinder Between Parallel Plates in Stokes Flow: Part 1: Motion in A Quiescent Fluid and Sedimentation. Computers & fluids, 15(4), 391-404. https://doi.org/10.1016/0045-7930(87)90031-4
Dvinsky, A S, and Popel, A S. "Motion of a Rigid Cylinder Between Parallel Plates in Stokes Flow: Part 1: Motion in A Quiescent Fluid and Sedimentation." Computers & fluids vol. 15,4 (1987): 391-404. doi: https://doi.org/10.1016/0045-7930(87)90031-4
Dvinsky AS, Popel AS. Motion of a Rigid Cylinder Between Parallel Plates in Stokes Flow: Part 1: Motion in A Quiescent Fluid and Sedimentation. Comput Fluids. 1987;15(4):391-404. doi: 10.1016/0045-7930(87)90031-4. PMID: 28943671; PMCID: PMC5609700.
Copy Download .nbib Format: NLM
Share it on

Comput Fluids. 1987;15(4):391-404. doi: 10.1016/0045-7930(87)90031-4.

Motion of a Rigid Cylinder Between Parallel Plates in Stokes Flow: Part 1: Motion in A Quiescent Fluid and Sedimentation.

Computers & fluids

A S Dvinsky, A S Popel

Affiliations

  1. Creare Inc., P.O. Box 71, Hanover, NH 03755, U.S.A.
  2. Department of Biomedical Engineering, School of Medicine, Johns Hopkins University, Baltimore, MD 21205, U.S.A.

PMID: 28943671 PMCID: PMC5609700 DOI: 10.1016/0045-7930(87)90031-4

Abstract

A general numerical scheme for solution of two-dimensional Stokes equations in a multiconnected domain of arbitrary shape [1,2] is applied to the motion of a rigid circular cylinder between plane parallel boundaries. Numerically generated boundary-conforming coordinates are used to transform the flow domain into a domain with rectilinear boundaries. The transformed Stokes equations in vorticity-stream function form are then solved on a uniform grid using an iterative algorithm. In Part I coefficients of the resistance matrix representing the forces and torque on the cylinder due to its translational motion parallel or perpendicular to the boundaries or due to rotation about its axis are calculated. The solutions are obtained for a wide range of particle radii and positions across the channel. It is found that the force on a particle translating parallel to the boundaries without rotation exhibits a minimum at a position between the channel centerline and the wall and a local maximum on the centerline. The resistance matrix is utilized to calculate translational and angular velocities of a free particle settling under gravity in a vertical channel. It is found that the translational velocity has a maximum at some lateral position and a minimum on the centerline; the particle angular velocity is opposite in sign to that of a particle rolling along the nearer channel wall except when the gap between the particle and the wall is very small. These results are compared with existing analytical solutions for a small cylindrical particle situated on the channel centerline, and with solutions of related 3-D problems for a spherical particle in a circular tube and in a plane channel. It is shown that the behavior of cylindrical and spherical particles in a channel in many cases is qualitatively different. This is attributed to different flow patterns in these two cases. The motion of a spherical particle in a circular tube has qualitative and quantitative features similar to those for a cylindrical particle in a plane channel.

References

  1. J Comput Phys. 1986 Nov;67(1):73-90 - PubMed
  2. Comput Fluids. 1987;15(4):405-419 - PubMed

Publication Types

Grant support