Premium
Stokes flow around cylinders in a bounded two‐dimensional domain using multipole‐accelerated boundary element methods
Author(s) -
Mammoli A. A.,
Ingber M. S.
Publication year - 1999
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/(sici)1097-0207(19990310)44:7<897::aid-nme530>3.0.co;2-s
Subject(s) - boundary element method , discretization , multipole expansion , fast multipole method , mathematics , stokes flow , boundary (topology) , mathematical analysis , flow (mathematics) , finite element method , physics , geometry , quantum mechanics , thermodynamics
The multipole technique has recently received attention in the field of boundary element analysis as a means of reducing the order of data storage and calculation time requirements from O ( N 2 ) (iterative solvers) or O ( N 3 ) (gaussian elimination) to O ( N log N ) or O ( N ), where N is the number of nodes in the discretized system. Such a reduction in the growth of the calculation time and data storage is crucial in applications where N is large, such as when modelling the macroscopic behaviour of suspensions of particles. In such cases, a minimum of 1000 particles is needed to obtain statistically meaningful results, leading to systems with N of the order of 10 000 for the smallest problems. When only boundary velocities are known, the indirect boundary element formulation for Stokes flow results in Fredholm equations of the second kind, which generally produce a well‐posed set of equations when discretized, a necessary requirement for iterative solution methods. The direct boundary element formulation, on the other hand, results in Fredholm equations of the first kind, which, upon discretization, produce ill‐conditioned systems of equations. The model system here is a two‐dimensional wide‐gap couette viscometer, where particles are suspended in the fluid between the cylinders. This is a typical system that is efficiently modelled using boundary element method simulations. The multipolar technique is applied to both direct and indirect formulations. It is found that the indirect approach is sufficiently well‐conditioned to allow the use of fast multipole methods. The direct approach results in severe ill‐conditioning, to a point where application of the multipole method leads to non‐convergence of the solution iteration. Copyright © 1999 John Wiley & Sons, Ltd.