Spectral Approximation of the Free-Space Heat Kernel

Abstract Many problems in applied mathematics, physics, and engineering require the solution of the heat equation in unbounded domains. Integral equation methods are particularly appropriate in this setting for several reasons: they are unconditionally stable, they are insensitive to the complexity of the geometry, and they do not require the artificial truncation of the computational domain as do finite difference and finite element techniques. Methods of this type, however, have not become widespread due to the high cost of evaluating heat potentials. When mpoints are used in the discretization of the initial data, Mpoints are used in the discretization of the boundary, and Ntime steps are computed, an amount of work of the order O(N2M2+NMm) has traditionally been required. In this paper, we present an algorithm which requires an amount of work of the order O(NMlog M+mlog m) and which is based on the evolution of the continuousspectrum of the solution. The method generalizes an earlier technique developed by Greengard and Strain (1990, Comm. Pure Appl. Math.43, 949) for evaluating layer potentials in bounded domains.

[1]  G. Beylkin.  On the Fast Fourier Transform of Functions with Singularities , 1995 .

[2]  A. Dutt,et al.  On the Rapid Evaluation of Trigonometric Series , 1992 .

[3]  E. Hille.  A Class of Reciprocal Functions , 1926 .

[4]  Kevin Barraclough,et al.  I and i , 2001, BMJ : British Medical Journal.

[5]  J. Strain.  Fast potential theory. II: Layer potentials and discrete sums , 1992 .

[6]  John Strain,et al.  Fast Adaptive Methods for the Free-Space Heat Equation , 1994, SIAM J. Sci. Comput..

[7]  Patrick Po-Yen Lin,et al.  On the numerical solution of the heat equation in unbounded domains , 1993 .

[8]  R. Stephenson.  A and V , 1962, The British journal of ophthalmology.

[9]  D. I. Meiron,et al.  Numerical simulations of unsteady crystal growth , 1992 .

[10]  V. Hutson.  Integral Equations , 1967, Nature.

[11]  A. Friedman.  Partial Differential Equations of Parabolic Type , 1983 .

[12]  J. Mcintyre.  Boundary integral solutions of the heat equation , 1986 .

[13]  Carlos Alberto Brebbia,et al.  Topics in Boundary Element Research , 1985 .

[14]  J. Sethian,et al.  Crystal growth and dendritic solidification , 1992 .

[15]  Petros Koumoutsakos,et al.  Direct numerical simulations of unsteady separated flows using vortex methods , 1993 .

[16]  Witold Pogorzelski,et al.  Integral equations and their applications , 1966 .

[17]  Philip Rabinowitz,et al.  Methods of Numerical Integration , 1985 .

[18]  P. Noon.  The Single Layer Heat Potential and Galerkin Boundary Element Methods for the Heat Equation , 1988 .

[19]  L. Greengard,et al.  A fast algorithm for the evaluation of heat potentials , 1990 .

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

[21]  Harry Dym,et al.  Fourier series and integrals , 1972 .

[22]  Leslie Greengard,et al.  The Fast Gauss Transform , 1991, SIAM J. Sci. Comput..