TY - JOUR
T1 - Multilevel nonuniform grid algorithm for acceleration of integral equation-based solvers for acoustic scattering
AU - Brick, Yaniv
AU - Boag, Amir
N1 - Funding Information:
Manuscript received april 19, 2009; accepted october 5, 2009. This research was supported by The Israel science Foundation under Grant 745/07. The authors are with the school of Electrical Engineering, Tel aviv University, Tel aviv, Israel (e-mail: [email protected]). digital object Identifier 10.1109/TUFFc.2010.1404
PY - 2010/1/1
Y1 - 2010/1/1
N2 - A fast algorithm for the evaluation of acoustic fields produced by given source distributions is developed with the aim of accelerating iterative boundary element method (BEM) solvers. The algorithm is based on field smoothing by phase and amplitude compensation, which allows for sampling of the fields radiated by finite-size source distributions over coarse nonuniform (spherical) grids (NGs). Subsequently, the fields at the desired target points can be obtained by an interpolation and phase and amplitude restoration. Combining this approach with the divide-and-conquer strategy, the total field is computed via a hierarchical decomposition of the source domain. In this computational scheme, the phase and amplitude compensated fields produced by neighboring subdomains are gradually aggregated through a multilevel process involving interpolation between increasingly dense NGs and the scatterer surface. This multilevel NG algorithm is used to reduce the computational cost of applying the field evaluation operator and its adjoint, as required in each iteration of the conjugate gradient solver based on the BEM-discretized integral representation of scattering problems. Accuracy and computational efficiency of the NG algorithm are demonstrated on representative examples of elongated, quasi-planar, and full 3-D scatterers.
AB - A fast algorithm for the evaluation of acoustic fields produced by given source distributions is developed with the aim of accelerating iterative boundary element method (BEM) solvers. The algorithm is based on field smoothing by phase and amplitude compensation, which allows for sampling of the fields radiated by finite-size source distributions over coarse nonuniform (spherical) grids (NGs). Subsequently, the fields at the desired target points can be obtained by an interpolation and phase and amplitude restoration. Combining this approach with the divide-and-conquer strategy, the total field is computed via a hierarchical decomposition of the source domain. In this computational scheme, the phase and amplitude compensated fields produced by neighboring subdomains are gradually aggregated through a multilevel process involving interpolation between increasingly dense NGs and the scatterer surface. This multilevel NG algorithm is used to reduce the computational cost of applying the field evaluation operator and its adjoint, as required in each iteration of the conjugate gradient solver based on the BEM-discretized integral representation of scattering problems. Accuracy and computational efficiency of the NG algorithm are demonstrated on representative examples of elongated, quasi-planar, and full 3-D scatterers.
UR - http://www.scopus.com/inward/record.url?scp=73649112787&partnerID=8YFLogxK
U2 - 10.1109/TUFFC.2010.1404
DO - 10.1109/TUFFC.2010.1404
M3 - Article
C2 - 20040452
AN - SCOPUS:73649112787
SN - 0885-3010
VL - 57
SP - 262
EP - 273
JO - IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control
JF - IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control
IS - 1
M1 - 5361548
ER -