A fast numerical method for solving the three-dimensional stokes' equations in the presence of suspended particles

Aaron L. Fogelson, Charles Peskin

Research output: Contribution to journalArticle

Abstract

A new fast numerical method for solving the three-dimensional Stokes' equations in the presence of suspended particles is presented. The fluid dynamics equations are solved on a lattice. A particle is represented by a set of points each of which moves at the local fluid velocity and is not constrained to lie on the lattice. These points are coupled by forces which resist deformation of the particle. These forces contribute to the force density 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. Particle size, shape, and deformability may be prescribed. Computational work increases only linearly with the number of particles, so large numbers (500-1000) of particles may be studied efficiently. The numerical method involves implicit calculation of the particle forces by minimizing an energy function and solution of a finite-difference approximation to the Stokes' equations using the Fourier-Toeplitz method. The numerical method has been implemented to run on all CRAY computers: the implementation exploits the CRAY's vectorized arithmetic, and on machines with insufficient central memory, it performs efficient disk 1/O while storing most of the data on 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.

Original languageEnglish (US)
Pages (from-to)50-69
Number of pages20
JournalJournal of Computational Physics
Volume79
Issue number1
DOIs
StatePublished - 1988

Fingerprint

Sedimentation
Numerical methods
Fluid dynamics
Formability
Particle size
Trajectories
settling
fluid dynamics
Data storage equipment
Fluids
cubic lattices
trajectories
fluids
approximation

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this

A fast numerical method for solving the three-dimensional stokes' equations in the presence of suspended particles. / Fogelson, Aaron L.; Peskin, Charles.

In: Journal of Computational Physics, Vol. 79, No. 1, 1988, p. 50-69.

Research output: Contribution to journalArticle

@article{f02d4d9ef5bb4cdea5d42e0307d8dfcc,
title = "A fast numerical method for solving the three-dimensional stokes' equations in the presence of suspended particles",
abstract = "A new fast numerical method for solving the three-dimensional Stokes' equations in the presence of suspended particles is presented. The fluid dynamics equations are solved on a lattice. A particle is represented by a set of points each of which moves at the local fluid velocity and is not constrained to lie on the lattice. These points are coupled by forces which resist deformation of the particle. These forces contribute to the force density 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. Particle size, shape, and deformability may be prescribed. Computational work increases only linearly with the number of particles, so large numbers (500-1000) of particles may be studied efficiently. The numerical method involves implicit calculation of the particle forces by minimizing an energy function and solution of a finite-difference approximation to the Stokes' equations using the Fourier-Toeplitz method. The numerical method has been implemented to run on all CRAY computers: the implementation exploits the CRAY's vectorized arithmetic, and on machines with insufficient central memory, it performs efficient disk 1/O while storing most of the data on 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.",
author = "Fogelson, {Aaron L.} and Charles Peskin",
year = "1988",
doi = "10.1016/0021-9991(88)90003-4",
language = "English (US)",
volume = "79",
pages = "50--69",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "1",

}

TY - JOUR

T1 - A fast numerical method for solving the three-dimensional stokes' equations in the presence of suspended particles

AU - Fogelson, Aaron L.

AU - Peskin, Charles

PY - 1988

Y1 - 1988

N2 - A new fast numerical method for solving the three-dimensional Stokes' equations in the presence of suspended particles is presented. The fluid dynamics equations are solved on a lattice. A particle is represented by a set of points each of which moves at the local fluid velocity and is not constrained to lie on the lattice. These points are coupled by forces which resist deformation of the particle. These forces contribute to the force density 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. Particle size, shape, and deformability may be prescribed. Computational work increases only linearly with the number of particles, so large numbers (500-1000) of particles may be studied efficiently. The numerical method involves implicit calculation of the particle forces by minimizing an energy function and solution of a finite-difference approximation to the Stokes' equations using the Fourier-Toeplitz method. The numerical method has been implemented to run on all CRAY computers: the implementation exploits the CRAY's vectorized arithmetic, and on machines with insufficient central memory, it performs efficient disk 1/O while storing most of the data on 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 suspended particles is presented. The fluid dynamics equations are solved on a lattice. A particle is represented by a set of points each of which moves at the local fluid velocity and is not constrained to lie on the lattice. These points are coupled by forces which resist deformation of the particle. These forces contribute to the force density 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. Particle size, shape, and deformability may be prescribed. Computational work increases only linearly with the number of particles, so large numbers (500-1000) of particles may be studied efficiently. The numerical method involves implicit calculation of the particle forces by minimizing an energy function and solution of a finite-difference approximation to the Stokes' equations using the Fourier-Toeplitz method. The numerical method has been implemented to run on all CRAY computers: the implementation exploits the CRAY's vectorized arithmetic, and on machines with insufficient central memory, it performs efficient disk 1/O while storing most of the data on 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=0002869426&partnerID=8YFLogxK

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

U2 - 10.1016/0021-9991(88)90003-4

DO - 10.1016/0021-9991(88)90003-4

M3 - Article

VL - 79

SP - 50

EP - 69

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 1

ER -