A quasilinear complexity algorithm for the numerical simulation of scattering from a two-dimensional radially symmetric potential

JOURNAL OF COMPUTATIONAL PHYSICS

卷:410

Article

Bremer, James¹ 

[1] Univ Calif Davis, Dept Math, Davis, CA 95616 USA

关键词: Helmholtz equation;    Scattering theory;    Fast algorithms;    Numerical solution of partial differential equations;   

DOI  :  10.1016/j.jcp.2020.109401

来源: Elsevier

【 摘 要 】

Standard solvers for the variable coefficient Helmholtz equation in two spatial dimensions have running times which grow at least quadratically with the wavenumber k. Here, we describe a solver which applies only when the scattering potential is radially symmetric but whose running time is O(klog(k)) in typical cases. We also present the results of numerical experiments demonstrating the properties of our solver, the code for which is publicly available. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free