A Chebyshev-based rectangular-polar integral solver for scattering by geometries described by non-overlapping patches

This paper introduces a high-order-accurate strategy for integration of singular kernels and edge-singular integral densities that appear in the context of boundary integral equation formulations of the problem of acoustic scattering. In particular, the proposed method is designed for use in conjunction with geometry descriptions given by a set of arbitrary non-overlapping logically-quadrilateral patches---which makes the algorithm particularly well suited for treatment of CAD-generated geometries. Fejer's first quadrature rule is incorporated in the algorithm, to provide a spectrally accurate method for evaluation of contributions from far integration regions, while highly-accurate precomputations of singular and near-singular integrals over certain "surface patches" together with two-dimensional Chebyshev transforms and suitable surface-varying "rectangular-polar" changes of variables, are used to obtain the contributions for singular and near-singular interactions. The overall integration method is then used in conjunction with the linear-algebra solver GMRES to produce solutions for sound-soft open- and closed-surface scattering obstacles, including an application to an aircraft described by means of a CAD representation. The approach is robust, fast, and highly accurate: use of a few points per wavelength suffices for the algorithm to produce far-field accuracies of a fraction of a percent, and slight increases in the discretization densities give rise to significant accuracy improvements.

[1]  Hank Childs,et al.  VisIt: An End-User Tool for Visualizing and Analyzing Very Large Data , 2011 .

[2]  J. Nédélec.  Acoustic and Electromagnetic Equations : Integral Representations for Harmonic Problems , 2001 .

[3]  R. Kress,et al.  Inverse Acoustic and Electromagnetic Scattering Theory , 1992 .

[4]  O. Bruno,et al.  A fast, high-order algorithm for the solution of surface scattering problems: basic implementation, tests, and applications , 2001 .

[5]  Catalin Turc,et al.  Planewave Density Interpolation Methods for 3D Helmholtz Boundary Integral Equations , 2019, SIAM J. Sci. Comput..

[6]  Erich Martensen,et al.  Über eine Methode zum räumlichen Neumannschen Problem mit einer Anwendung für torusartige Berandungen , 1963 .

[7]  J. Waldvogel.  Fast Construction of the Fejér and Clenshaw–Curtis Quadrature Rules , 2006 .

[8]  Oscar P. Bruno,et al.  A high-order integral solver for scalar problems of diffraction by screens and apertures in three-dimensional space , 2012, J. Comput. Phys..

[9]  Ivan G. Graham,et al.  Product integration-collocation methods for noncompact integral operator equations , 1988 .

[10]  Leslie Greengard,et al.  Quadrature by expansion: A new method for the evaluation of layer potentials , 2012, J. Comput. Phys..

[11]  Pérez Arancibia,et al.  Windowed Integral Equation Methods for Problems of Scattering by Defects and Obstacles in Layered Media , 2017 .

[12]  Oscar P. Bruno,et al.  Windowed Green Function Method for Nonuniform Open-Waveguide Problems , 2016, IEEE Transactions on Antennas and Propagation.

[13]  Martin Costabel,et al.  General edge asymptotics of solutions of second-order elliptic boundary value problems I , 1993 .

[14]  Johannes Elschner,et al.  The double-layer potential operator over polyhedral domains II: Spline Galerkin methods , 1992 .

[15]  David Elliott,et al.  A sinh transformation for evaluating nearly singular boundary element integrals , 2005 .

[16]  Catalin Turc,et al.  Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D , 2019, J. Comput. Phys..

[17]  Akash Anand,et al.  Efficient Solution of Three-Dimensional Problems of Acoustic and Electromagnetic Scattering by Open Surfaces , 2011 .

[18]  R. Duraiswami,et al.  Fast Multipole Methods for the Helmholtz Equation in Three Dimensions , 2005 .

[19]  Barbara M. Johnston,et al.  A sinh transformation for evaluating two‐dimensional nearly singular boundary element integrals , 2007 .

[20]  Mark Lyon,et al.  Windowed Green Function Method for Layered-Media Scattering , 2015, SIAM J. Appl. Math..

[21]  Ari Sihvola,et al.  Surface Integral Equation Method for Scattering by DB Objects with Sharp Wedges , 2011 .

[22]  O. Bruno,et al.  Windowed Green function method for the Helmholtz equation in the presence of multiply layered media , 2017, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[23]  Rainer Kußmaul,et al.  Ein numerisches Verfahren zur Lösung des Neumannschen Außenraumproblems für die Helmholtzsche Schwingungsgleichung , 1969, Computing.

[24]  Oscar P. Bruno,et al.  A high-order integral algorithm for highly singular PDE solutions in Lipschitz domains , 2009, Computing.

[25]  James Bremer,et al.  A Nyström method for weakly singular integral operators on surfaces , 2012, J. Comput. Phys..

[26]  W. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[27]  Francisco-Javier Sayas,et al.  Convergence analysis of a high-order Nyström integral-equation method for surface scattering problems , 2011, Numerische Mathematik.

[28]  Oscar P. Bruno,et al.  Regularized integral equations and fast high‐order solvers for sound‐hard acoustic scattering problems , 2012 .

[29]  A. Rathsfeld.  Nystro¨m's method and iterative solvers for the solution of the double-layer potential equation over polyhedral boundaries , 1995 .

[30]  Ivan G. Graham,et al.  A high-order algorithm for obstacle scattering in three dimensions , 2004 .

[31]  Oscar P. Bruno,et al.  Surface scattering in three dimensions: an accelerated high–order solver , 2001, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.