Energetic boundary element method for accurate solution of damped waves hard scattering problems

The paper deals with the numerical solution of 2D wave propagation exterior problems including viscous and material damping coefficients and equipped by Neumann boundary condition, hence modeling the hard scattering of damped waves. The differential problem, which includes, besides diffusion, advection and reaction terms, is written as a space–time boundary integral equation (BIE) whose kernel is given by the hypersingular fundamental solution of the 2D damped waves operator. The resulting BIE is solved by a modified Energetic Boundary Element Method, where a suitable kernel treatment is introduced for the evaluation of the discretization linear system matrix entries represented by space–time quadruple integrals with hypersingular kernel in space variables. A wide variety of numerical results, obtained varying both damping coefficients and discretization parameters, is presented and shows accuracy and stability of the proposed technique, confirming what was theoretically proved for the simpler undamped case. Post-processing phase is also taken into account, giving the approximate solution of the exterior differential problem involving damped waves propagation around disconnected obstacles and bounded domains.