Rapid evaluation of nonreflecting boundary kernels or time-domain wave propagation

被引:162
作者
Alpert, B
Greengard, L
Hagstrom, T
机构
[1] Natl Inst Stand & Technol, Boulder, CO 80303 USA
[2] NYU, Courant Inst Math Sci, New York, NY 10012 USA
[3] Univ New Mexico, Dept Math & Stat, Albuquerque, NM 87131 USA
关键词
Bessel function; approximation; high-order convergence; wave equation; Maxwell's equations; nonreflecting boundary condition; radiation boundary condition; absorbing boundary condition;
D O I
10.1137/S0036142998336916
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We present a systematic approach to the computation of exact nonreflecting boundary conditions for the wave equation. In both two and three dimensions, the critical step in our analysis involves convolution with the inverse Laplace transform of the logarithmic derivative of a Hankel function. The main technical result in this paper is that the logarithmic derivative of the Hankel function H-nu((1))(z) of real order nu can be approximated in the upper half z-plane with relative error epsilon by a rational function of degree d similar to O(log\nu\ log 1/epsilon + log(2) \nu\ + \nu\(-1) log(2) 1/epsilon) as \nu\ --> infinity, epsilon --> 0, with slightly more complicated bounds for nu = 0. If N is the number of points used in the discretization of a cylindrical (circular) boundary in two dimensions, then, assuming that epsilon < 1/N, O(N log N log 1/epsilon) work is required at each time step. This is comparable to the work required for the Fourier transform on the boundary. In three dimensions, the cost is proportional to N-2 log(2) N + N-2 log N log 1/epsilon for a spherical boundary with N-2 points, the first term coming from the calculation of a spherical harmonic transform at each time step. In short, nonreflecting boundary conditions can be imposed to any desired accuracy, at a cost dominated by the interior grid work, which scales like N-2 in two dimensions and N-3 in three dimensions.
引用
收藏
页码:1138 / 1164
页数:27
相关论文
共 30 条