Fast multipole method for the biharmonic equation in three dimensions

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^3 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 L"2 accuracy @e"2=10^-^4, and N=3550 for @e"2=10^-^1^2.

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

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

[3]  L. E. Malvern.  Introduction to the mechanics of a continuous medium , 1969 .

[4]  Naoshi Nishimura,et al.  Application of new fast multipole boundary integral equation method to crack problems in 3D , 2001 .

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

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

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

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

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

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

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

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

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

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

[15]  L. Greengard,et al.  On the numerical solution of the biharmonic equation in the plane , 1992 .

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

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

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

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

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

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

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

[23]  C. Alberini,et al.  Memory , 2006, Cellular and Molecular Life Sciences CMLS.

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

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

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

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

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

[29]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .