Fast Multipole Method for the Biharmonic Equation

The evaluation of sums (matrix-vector products) of the solutions of the three-dimensional biharmonic equation can be accelerated using the fast multipole method, while memory requirements can also be significantly reduced. We develop a complete translation theory for these equations. It is shown that translations of elementary solutions of the biharmonic equation can be achieved by considering the translation of a pair of elementary solutions of the Laplace equations. The extension of the theory to the case of polyharmonic equations in R is also discussed. An efficient way of performing the FMM for biharmonic equations using the solution of a complex valued FMM for the Laplace equation is presented. Compared to previous methods presented for the biharmonic equation our method appears more efficient. The theory is implemented and numerical tests presented that demonstrate the performance of the method for varying problem sizes and accuracy requirements. In our implementation, the FMM for the biharmonic equation is faster than direct matrix vector product for a matrix size of 550 for a relative L2 accuracy 22 = 10−4, and N = 3550 for 22 = 10−12.

[1]  Lexing Ying,et al.  A New Parallel Kernel-Independent Fast Multipole Method , 2003, ACM/IEEE SC 2003 Conference (SC'03).

[2]  G. Rodin,et al.  Fast solution method for three-dimensional Stokesian many-particle problems , 2000 .

[3]  Michael A. Epton,et al.  Multipole Translation Theory for the Three-Dimensional Laplace and Helmholtz Equations , 1995, SIAM J. Sci. Comput..

[4]  Jean Duchon,et al.  Splines minimizing rotation-invariant semi-norms in Sobolev spaces , 1976, Constructive Theory of Functions of Several Variables.

[5]  L. Greengard,et al.  Integral Equation Methods for Stokes Flow and Isotropic Elasticity in the Plane , 1996 .

[6]  Kenichi Yoshida,et al.  Applications of Fast Multipole Method to Boundary Integral Equation Method , 2001 .

[7]  Sangtae Kim,et al.  Stokes flow past three spheres: An analytic solution , 1987 .

[8]  Fang Chen,et al.  Fast evaluation of vector splines in three dimensions , 1998, Computing.

[9]  Richard K. Beatson,et al.  Fast Evaluation of Radial Basis Functions: Methods for Generalized Multiquadrics in Rn , 2001, SIAM J. Sci. Comput..

[10]  J. Happel,et al.  Low Reynolds number hydrodynamics , 1965 .

[11]  A. Sangani,et al.  An O(N) algorithm for Stokes and Laplace interactions of particles , 1996 .

[12]  N. Nishimura.  Fast multipole accelerated boundary integral equation methods , 2002 .

[13]  L. Greengard.  The Rapid Evaluation of Potential Fields in Particle Systems , 1988 .

[14]  Ramani Duraiswami,et al.  Comparison of the efficiency of translation operators used in the fast multipole method for the 3D Laplace equation , 2005 .

[15]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[16]  L. Greengard,et al.  A new version of the Fast Multipole Method for the Laplace equation in three dimensions , 1997, Acta Numerica.

[17]  L. Greengard,et al.  Regular Article: A Fast Adaptive Multipole Algorithm in Three Dimensions , 1999 .

[18]  Richard K. Beatson,et al.  Reconstruction and representation of 3D objects with radial basis functions , 2001, SIGGRAPH.

[19]  R. Duraiswami,et al.  Fast Multipole Methods for the Helmholtz Equation in Three Dimensions , 2005 .

[20]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[21]  A. V. Filippov,et al.  Phoretic Motion of Arbitrary Clusters of N Spheres , 2001 .

[22]  R. V. D. Geijn,et al.  A fast solution method for three‐dimensional many‐particle problems of linear elasticity , 1998 .

[23]  Martin Head-Gordon,et al.  Rotating around the quartic angular momentum barrier in fast multipole method calculations , 1996 .

[24]  V. Popov,et al.  An O(N) Taylor series multipole boundary element method for three-dimensional elasticity problems , 2001 .