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: boag@eng.tau.ac.il). 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

VL - 57

SP - 262

EP - 273

JO - IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control

JF - IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control

SN - 0885-3010

IS - 1

M1 - 5361548

ER -