Improvement of the Stokesian Dynamics method for systems with finite number of particles

K. Ichiki

Research output: Contribution to journalArticleAcademicpeer-review

30 Citations (Scopus)

Abstract

An improvement of the Stokesian Dynamics method for many-particle systems is presented. A direct calculation of the hydrodynamic interaction is used rather than imposing periodic boundary conditions. The two major diculties concern the accuracy and the speed of calculations. The accuracy discussed in this work is not concerned with the lubrication correction but, rather, focuses on the multipole expansion which until now has only been formulated up to the so-called FTS version or the rst order of force moments. This is improved systematically by a real-space multipole expansion with force moments and velocity moments evaluated at the centre of the particles, where the velocity moments are calculated through the velocity derivatives; the introduction of the velocity derivatives makes the formulation and its extensions straightforward. The reduction of the moments into irreducible form is achieved by the Cartesian irreducible tensor. The reduction is essential to form a well-dened linear set of equations as a generalized mobility problem. The order of truncation is not limited in principle, and explicit calculations of two-body problems are shown with order up to 7. The calculating speed is improved by a conjugate-gradient-type iterative method which consists of a dot-product between the generalized mobility matrix and the force moments as a trial value in each iteration. This provides an O(N2) scheme where N is the number of particles in the system. Further improvement is achieved by the fast multipole method for the calculation of the generalized mobility problem in each iteration, and an O(N) scheme for the non-adaptive version is obtained. Real problems are studied on systems with N = 400 000 particles. For mobility problems the number of iterations is constant and an O(N) performance is achieved; however for resistance problems the number of iterations increases as almost N1=2 with a high accuracy of 10¿6 and the total cost seems to be O(N3=2). 1. Introduction The microstructure of suspensions is governed by the hydrodynamic interactions among particles immersed in a viscous fluid, which is modelled using the Stokes approximation, and have attracted much attention from researchers in physics and chemical engineering. The hydrodynamic interactions have a long-range nature varying as 1=r, where r is distance measured from a particle, and further they have a many-body feature, that is, they must be the solution of a boundary-value problem on the surface of all objects in the system. Therefore, analytical approaches are dicult. In fact, even for rigid spherical particles, the exact solution has been obtained only for two-body problems (Jerey & Onishi 1984); of course, this is partially because the symmetry of the geometry of surfaces for two-body problems is much simpler than that on systems with three or more particles. Therefore, numerical approaches have
Original languageUndefined
Pages (from-to)231-262
Number of pages32
JournalJournal of fluid mechanics
Volume452
DOIs
Publication statusPublished - 2002

Keywords

  • IR-36626
  • METIS-202587

Cite this

@article{b4cdcdb5d35e437c906cb8bf84c11f67,
title = "Improvement of the Stokesian Dynamics method for systems with finite number of particles",
abstract = "An improvement of the Stokesian Dynamics method for many-particle systems is presented. A direct calculation of the hydrodynamic interaction is used rather than imposing periodic boundary conditions. The two major diculties concern the accuracy and the speed of calculations. The accuracy discussed in this work is not concerned with the lubrication correction but, rather, focuses on the multipole expansion which until now has only been formulated up to the so-called FTS version or the rst order of force moments. This is improved systematically by a real-space multipole expansion with force moments and velocity moments evaluated at the centre of the particles, where the velocity moments are calculated through the velocity derivatives; the introduction of the velocity derivatives makes the formulation and its extensions straightforward. The reduction of the moments into irreducible form is achieved by the Cartesian irreducible tensor. The reduction is essential to form a well-dened linear set of equations as a generalized mobility problem. The order of truncation is not limited in principle, and explicit calculations of two-body problems are shown with order up to 7. The calculating speed is improved by a conjugate-gradient-type iterative method which consists of a dot-product between the generalized mobility matrix and the force moments as a trial value in each iteration. This provides an O(N2) scheme where N is the number of particles in the system. Further improvement is achieved by the fast multipole method for the calculation of the generalized mobility problem in each iteration, and an O(N) scheme for the non-adaptive version is obtained. Real problems are studied on systems with N = 400 000 particles. For mobility problems the number of iterations is constant and an O(N) performance is achieved; however for resistance problems the number of iterations increases as almost N1=2 with a high accuracy of 10¿6 and the total cost seems to be O(N3=2). 1. Introduction The microstructure of suspensions is governed by the hydrodynamic interactions among particles immersed in a viscous fluid, which is modelled using the Stokes approximation, and have attracted much attention from researchers in physics and chemical engineering. The hydrodynamic interactions have a long-range nature varying as 1=r, where r is distance measured from a particle, and further they have a many-body feature, that is, they must be the solution of a boundary-value problem on the surface of all objects in the system. Therefore, analytical approaches are dicult. In fact, even for rigid spherical particles, the exact solution has been obtained only for two-body problems (Jerey & Onishi 1984); of course, this is partially because the symmetry of the geometry of surfaces for two-body problems is much simpler than that on systems with three or more particles. Therefore, numerical approaches have",
keywords = "IR-36626, METIS-202587",
author = "K. Ichiki",
year = "2002",
doi = "10.1017/S0022112001006735",
language = "Undefined",
volume = "452",
pages = "231--262",
journal = "Journal of fluid mechanics",
issn = "0022-1120",
publisher = "Cambridge University Press",

}

Improvement of the Stokesian Dynamics method for systems with finite number of particles. / Ichiki, K.

In: Journal of fluid mechanics, Vol. 452, 2002, p. 231-262.

Research output: Contribution to journalArticleAcademicpeer-review

TY - JOUR

T1 - Improvement of the Stokesian Dynamics method for systems with finite number of particles

AU - Ichiki, K.

PY - 2002

Y1 - 2002

N2 - An improvement of the Stokesian Dynamics method for many-particle systems is presented. A direct calculation of the hydrodynamic interaction is used rather than imposing periodic boundary conditions. The two major diculties concern the accuracy and the speed of calculations. The accuracy discussed in this work is not concerned with the lubrication correction but, rather, focuses on the multipole expansion which until now has only been formulated up to the so-called FTS version or the rst order of force moments. This is improved systematically by a real-space multipole expansion with force moments and velocity moments evaluated at the centre of the particles, where the velocity moments are calculated through the velocity derivatives; the introduction of the velocity derivatives makes the formulation and its extensions straightforward. The reduction of the moments into irreducible form is achieved by the Cartesian irreducible tensor. The reduction is essential to form a well-dened linear set of equations as a generalized mobility problem. The order of truncation is not limited in principle, and explicit calculations of two-body problems are shown with order up to 7. The calculating speed is improved by a conjugate-gradient-type iterative method which consists of a dot-product between the generalized mobility matrix and the force moments as a trial value in each iteration. This provides an O(N2) scheme where N is the number of particles in the system. Further improvement is achieved by the fast multipole method for the calculation of the generalized mobility problem in each iteration, and an O(N) scheme for the non-adaptive version is obtained. Real problems are studied on systems with N = 400 000 particles. For mobility problems the number of iterations is constant and an O(N) performance is achieved; however for resistance problems the number of iterations increases as almost N1=2 with a high accuracy of 10¿6 and the total cost seems to be O(N3=2). 1. Introduction The microstructure of suspensions is governed by the hydrodynamic interactions among particles immersed in a viscous fluid, which is modelled using the Stokes approximation, and have attracted much attention from researchers in physics and chemical engineering. The hydrodynamic interactions have a long-range nature varying as 1=r, where r is distance measured from a particle, and further they have a many-body feature, that is, they must be the solution of a boundary-value problem on the surface of all objects in the system. Therefore, analytical approaches are dicult. In fact, even for rigid spherical particles, the exact solution has been obtained only for two-body problems (Jerey & Onishi 1984); of course, this is partially because the symmetry of the geometry of surfaces for two-body problems is much simpler than that on systems with three or more particles. Therefore, numerical approaches have

AB - An improvement of the Stokesian Dynamics method for many-particle systems is presented. A direct calculation of the hydrodynamic interaction is used rather than imposing periodic boundary conditions. The two major diculties concern the accuracy and the speed of calculations. The accuracy discussed in this work is not concerned with the lubrication correction but, rather, focuses on the multipole expansion which until now has only been formulated up to the so-called FTS version or the rst order of force moments. This is improved systematically by a real-space multipole expansion with force moments and velocity moments evaluated at the centre of the particles, where the velocity moments are calculated through the velocity derivatives; the introduction of the velocity derivatives makes the formulation and its extensions straightforward. The reduction of the moments into irreducible form is achieved by the Cartesian irreducible tensor. The reduction is essential to form a well-dened linear set of equations as a generalized mobility problem. The order of truncation is not limited in principle, and explicit calculations of two-body problems are shown with order up to 7. The calculating speed is improved by a conjugate-gradient-type iterative method which consists of a dot-product between the generalized mobility matrix and the force moments as a trial value in each iteration. This provides an O(N2) scheme where N is the number of particles in the system. Further improvement is achieved by the fast multipole method for the calculation of the generalized mobility problem in each iteration, and an O(N) scheme for the non-adaptive version is obtained. Real problems are studied on systems with N = 400 000 particles. For mobility problems the number of iterations is constant and an O(N) performance is achieved; however for resistance problems the number of iterations increases as almost N1=2 with a high accuracy of 10¿6 and the total cost seems to be O(N3=2). 1. Introduction The microstructure of suspensions is governed by the hydrodynamic interactions among particles immersed in a viscous fluid, which is modelled using the Stokes approximation, and have attracted much attention from researchers in physics and chemical engineering. The hydrodynamic interactions have a long-range nature varying as 1=r, where r is distance measured from a particle, and further they have a many-body feature, that is, they must be the solution of a boundary-value problem on the surface of all objects in the system. Therefore, analytical approaches are dicult. In fact, even for rigid spherical particles, the exact solution has been obtained only for two-body problems (Jerey & Onishi 1984); of course, this is partially because the symmetry of the geometry of surfaces for two-body problems is much simpler than that on systems with three or more particles. Therefore, numerical approaches have

KW - IR-36626

KW - METIS-202587

U2 - 10.1017/S0022112001006735

DO - 10.1017/S0022112001006735

M3 - Article

VL - 452

SP - 231

EP - 262

JO - Journal of fluid mechanics

JF - Journal of fluid mechanics

SN - 0022-1120

ER -