Rapid solution of integral equations of scattering theory in two dimensions
Journal of Computational Physics
SIAM Journal on Scientific Computing
Handbook of Mathematical Functions, With Formulas, Graphs, and Mathematical Tables,
Handbook of Mathematical Functions, With Formulas, Graphs, and Mathematical Tables,
Journal of Computational Physics
Fast evaluation of Helmholtz potential on graphics processing units (GPUs)
Journal of Computational Physics
Fast evaluation of Helmholtz potential on graphics processing units (GPUs)
Journal of Computational Physics
Hi-index | 31.45 |
A multilevel Cartesian non-uniform grid time domain algorithm (CNGTDA) is introduced to rapidly compute transient wave fields radiated by time dependent three-dimensional source constellations. CNGTDA leverages the observation that transient wave fields generated by temporally bandlimited and spatially confined source constellations can be recovered via interpolation from appropriately delay- and amplitude-compensated field samples. This property is used in conjunction with a multilevel scheme, in which the computational domain is hierarchically decomposed into subdomains with sparse non-uniform grids used to obtain the fields. For both surface and volumetric source distributions, the computational cost of CNGTDA to compute the transient field at N"s observation locations from N"s collocated sources for N"t discrete time instances scales as O(N"tN"slogN"s) and O(N"tN"slog^2N"s) in the low- and high-frequency regimes, respectively. Coupled with marching-on-in-time (MOT) time domain integral equations, CNGTDA can facilitate efficient analysis of large scale time domain electromagnetic and acoustic problems.