An efficient solution of time domain boundary integral equations for acoustic scattering and its acceleration by Graphics Processing Units

The present paper is aimed at developing a fast numerical solution of the time domain boundary integral equation (TDBIE) reformulated from the convective wave equation for large scale wave scattering and propagation problems. Historically, numerical solutions of boundary integral equation in the time domain have encountered two major difficulties. The first is the intrinsic numerical instability in the early time domain boundary integral equation formulations. And the second is the formidably high computational cost associated with the direct solution of the time-domain boundary integral equation. In this paper, both issues are addressed. A stable Burton-Miller type formulation is proposed for the time domain boundary integral equation in the presence of a mean flow. A justification for stability through the energy equation associated with the convective wave equation is given. A comparison of the current formulation with a previous one in literature is also offered. The boundary integral equation is solved by a time domain boundary element method (TDBEM), using high-order basis functions and unstructured surface elements. To significantly reduce the computational cost, a Time Domain Propagation and Distribution (TDPD) algorithm is proposed, making use of the delayand amplitude-compensated field with a mean flow. Implemented in multi-level interactions, the current algorithm shows a computational cost of O(N) per time step where N is the total number of unknowns on surface elements. Furthermore, GPU computing has been utilized to speedup the computation. Numerical aspects of the GPU computing for boundary element solutions are discussed. Comparison with CPU executions is also given. Numerical examples that demonstrate the capabilities of the proposed method are presented.

[1]  T. W. Wu,et al.  A HYPERSINGULAR INTEGRAL FORMULATION FOR ACOUSTIC RADIATION IN MOVING FLOWS , 1997 .

[2]  Shaojing Li,et al.  Fast evaluation of Helmholtz potential on graphics processing units (GPUs) , 2010, J. Comput. Phys..

[3]  D. L. Hawkings,et al.  Sound generation by turbulence and surfaces in arbitrary motion , 1969, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[4]  Wolfgang Hackbusch,et al.  Sparse convolution quadrature for time domain boundary integral formulations of the wave equation , 2008 .

[5]  Christopher R. Anderson,et al.  An Implementation of the Fast Multipole Method without Multipoles , 1992, SIAM J. Sci. Comput..

[6]  T. Takahashi,et al.  GPU‐accelerated boundary element method for Helmholtz' equation in three dimensions , 2009 .

[7]  David P. Lockard,et al.  AN EFFICIENT, TWO-DIMENSIONAL IMPLEMENTATION OF THE FFOWCS WILLIAMS AND HAWKINGS EQUATION , 2000 .

[8]  Zydrunas Gimbutas,et al.  A wideband fast multipole method for the Helmholtz equation in three dimensions , 2006, J. Comput. Phys..

[9]  Ramani Duraiswami,et al.  A broadband fast multipole accelerated boundary element method for the three dimensional Helmholtz equation. , 2009, The Journal of the Acoustical Society of America.

[10]  M. Myers,et al.  On the application of the Kirchhoff formula for moving surfaces , 1990 .

[11]  Jun Meng,et al.  A multilevel Cartesian non-uniform grid time domain algorithm , 2010, J. Comput. Phys..

[12]  C. Vuik,et al.  Design of Temporal Basis Functions for Time Domain Integral Equation Methods With Predefined Accuracy and Smoothness , 2013, IEEE Transactions on Antennas and Propagation.

[13]  Martin Schanz,et al.  Time Domain BEM: Numerical Aspects of Collocation and Galerkin Formulations , 2009 .

[14]  Shaojing Li,et al.  Graphics Processing Unit Accelerated $O(N)$ Micromagnetic Solver , 2010, IEEE Transactions on Magnetics.

[15]  W. Hackbusch,et al.  Introduction to Hierarchical Matrices with Applications , 2003 .

[16]  Michael Garland,et al.  Understanding throughput-oriented architectures , 2010, Commun. ACM.

[17]  J.A.M. Carrer,et al.  Stability of Galerkin and collocation time domain boundary element methods as applied to the scalar wave equation , 2000 .

[18]  Mike Houston,et al.  A closer look at GPUs , 2008, Commun. ACM.

[19]  J. Ih,et al.  Stabilization of time domain acoustic boundary element method for the exterior problem avoiding the nonuniqueness. , 2013, The Journal of the Acoustical Society of America.

[20]  C. Lubich,et al.  On the multistep time discretization of linear\newline initial-boundary value problems and their boundary integral equations , 1994 .

[21]  T. Ha-Duong,et al.  A Galerkin BEM for transient acoustic scattering by an absorbing obstacle , 2003 .

[22]  W. Hwang.  Boundary spectral method for acoustic scattering and radiation problems , 1997 .

[23]  Luigi Morino,et al.  Is There a Difference Between Aeroacoustics and Aerodynamics? An Aeroelastician's Viewpoint , 2003 .

[24]  Anthony R. Pilon,et al.  Development of an Improved Kirchhoff Method for Jet Aeroacoustics , 1998 .

[25]  A. Curtis,et al.  Exact wave field simulation for finite-volume scattering problems. , 2007, The Journal of the Acoustical Society of America.

[26]  G. F. Miller,et al.  The application of integral equation methods to the numerical solution of some exterior boundary-value problems , 1971, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[27]  Amir Boag,et al.  Nonuniform grid time domain (NGTD) algorithm for fast evaluation of transient wave fields , 2006 .

[28]  Ana F. Tinetti,et al.  Application of Fast Multipole Methods to the NASA Fast Scattering Code , 2008 .

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

[30]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[31]  R. P. Shaw.  Transient scattering by a circular cylinder , 1975 .

[32]  A. Brandt.  Multilevel computations of integral transforms and particle interactions with oscillatory kernels , 1991 .

[33]  M. K. Myers,et al.  Computation of acoustic scattering from a moving rigid surface , 1992 .

[34]  R. J. Astley,et al.  A three-dimensional boundary element scheme for acoustic radiation in low mach number flows , 1986 .

[35]  Lyle N. Long.  The compressible aerodynamics of rotating blades based on an acoustic formulation , 1983 .

[36]  G. Franceschetti,et al.  On the spatial bandwidth of scattered fields , 1987 .

[37]  E. Michielssen,et al.  A multilevel matrix decomposition algorithm for analyzing scattering from large structures , 1996 .

[38]  Yueping Guo.  Slat noise modeling and prediction , 2012 .

[39]  Dj Chappell,et al.  A stable boundary element method for modeling transient acoustic radiation , 2006 .

[40]  F. Farassat,et al.  Extension of Kirchhoff's formula to radiation from moving surfaces , 1988 .

[42]  Willi Möhring,et al.  Energy conservation, time-reversal invariance and reciprocity in ducts with flow , 2001, Journal of Fluid Mechanics.

[43]  K. Mitzner,et al.  Numerical Solution for Transient Scattering from a Hard Surface of Arbitrary Shape—Retarded Potential Technique , 1967 .

[44]  Leonard V. Lopes,et al.  Design of the Next Generation Aircraft Noise Prediction Program: ANOPP2 , 2011 .

[45]  C. Lubich.  Convolution quadrature and discretized operational calculus. II , 1988 .

[46]  A Spectral Boundary Integral Equation Method for the 2D Helmholtz Equation , 1995 .

[47]  Gordon Erlebacher,et al.  High-order finite-element seismic wave propagation modeling with MPI on a large GPU cluster , 2010, J. Comput. Phys..

[48]  Sean F. Wu.  Nonuniqueness of solutions to extended Kirchhoff integral formulations , 1993 .

[49]  F. Hu,et al.  A Three-Dimensional Time-Domain Boundary Element Method for the Computation of Exact Green's Functions in Acoustic Analogy , 2007 .

[50]  E. Michielssen,et al.  Analysis of transient wave scattering from rigid bodies using a Burton–Miller approach , 1999 .

[51]  W. R. M. M.Sc.  XIV. The kirchhoff formula extended to a moving surface , .

[52]  F. Hu.  A FAST NUMERICAL SOLUTION OF SCATTERING BY A CYLINDER: SPECTRAL METHOD FOR THE BOUNDARY INTEGRAL EQUATIONS , 1994 .

[53]  M. B. Friedman,et al.  Diffraction of Pulses by Cylindrical Obstacles of Arbitrary Cross Section , 1962 .

[54]  A. Boag,et al.  Multilevel nonuniform grid algorithm for acceleration of integral equation-based solvers for acoustic scattering , 2010, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[55]  C. Morfey.  Acoustic energy in non-uniform flows , 1971 .

[56]  Ramani Duraiswami,et al.  Fast multipole methods on graphics processors , 2008, J. Comput. Phys..

[57]  H. A. Schenck.  Improved Integral Formulation for Acoustic Radiation Problems , 1968 .

[58]  Ergin,et al.  Fast analysis of transient acoustic wave scattering from rigid bodies using the multilevel plane wave time domain algorithm , 2000, The Journal of the Acoustical Society of America.

[59]  D. Zorin,et al.  A kernel-independent adaptive fast multipole algorithm in two and three dimensions , 2004 .

[60]  Timothy C. Warburton,et al.  Nodal discontinuous Galerkin methods on graphics processors , 2009, J. Comput. Phys..