TY - CHAP

T1 - Numerical solution of the three-dimensional Stokes' equations in the presence of suspended particles.

AU - Fogelson, A. L.

AU - Peskin, C. S.

PY - 1986

Y1 - 1986

N2 - A new fast numerical method for solving the three dimensional Stokes' equations in the presence of discrete suspended particles is presented. The fluid dynamics equations are solved on a grid. The particles are represented by configurations of localized forces which are not constrained to move on the grid. These forces contribute to the force density term in the Stokes' equations. As a result, a single set of fluid dynamics equations holds at all points of the domain and there are no internal boundaries. Computational work increases only linearly with the number of particles, so large numbers (500-1000) of particles may be studied efficiently. Particle size, shape, and deformability may be prescribed. The method was implemented to run on all CRAY computers: the implementation exploits the CRAY's vectorized arithmetic, and on machines with insufficient memory, it performs efficient disk I/O while storing most of the data on the disk. Applications of the method to sedimentation of one, two, and many particle systems are described. Trajectories and settling speeds for two particle sedimentation, and settling speed for multiparticle sedimentation from initial distributions on a cubic lattice or at random give good quantitative agreement with existing theories.

AB - A new fast numerical method for solving the three dimensional Stokes' equations in the presence of discrete suspended particles is presented. The fluid dynamics equations are solved on a grid. The particles are represented by configurations of localized forces which are not constrained to move on the grid. These forces contribute to the force density term in the Stokes' equations. As a result, a single set of fluid dynamics equations holds at all points of the domain and there are no internal boundaries. Computational work increases only linearly with the number of particles, so large numbers (500-1000) of particles may be studied efficiently. Particle size, shape, and deformability may be prescribed. The method was implemented to run on all CRAY computers: the implementation exploits the CRAY's vectorized arithmetic, and on machines with insufficient memory, it performs efficient disk I/O while storing most of the data on the disk. Applications of the method to sedimentation of one, two, and many particle systems are described. Trajectories and settling speeds for two particle sedimentation, and settling speed for multiparticle sedimentation from initial distributions on a cubic lattice or at random give good quantitative agreement with existing theories.

UR - http://www.scopus.com/inward/record.url?scp=85040886935&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85040886935&partnerID=8YFLogxK

M3 - Chapter

AN - SCOPUS:85040886935

SN - 0898712122

SN - 9780898712124

BT - Unknown Host Publication Title

PB - Soc. Ind. & Appl. Math

ER -