Abstract

Time-domain Boundary Element Methods (BEM) have been successfully used in acoustics, optics and elastodynamics to solve transient problems numerically. However, the storage requirements are immense, since the fully populated system matrices have to be computed for a large number of time steps or frequencies. In this article, we propose a new approximation scheme for the Convolution Quadrature Method powered BEM, which we apply to scattering problems governed by the wave equation. We use {mathscr {H}}^2-matrix compression in the spatial domain and employ an adaptive cross approximation algorithm in the frequency domain. In this way, the storage and computational costs are reduced significantly, while the accuracy of the method is preserved.

Highlights

  • The numerical solution of wave propagation problems is a crucial task in computational mathematics and its applications

  • In contrast to Finite Element or Difference Methods, Boundary Element Methods (BEM) are based on boundary integral equations posed in terms of the traces of the solution

  • Since the boundary data is only given on ΓD or ΓN, we derive a system of boundary integral equations for its unknown parts

Read more

Summary

Introduction

The numerical solution of wave propagation problems is a crucial task in computational mathematics and its applications. In this context, BEM play a special role, since they only require the discretisation of the boundary instead of the whole domain. BEM are favourable in situations where the domain is unbounded, as it is often the case for scattering problems. In contrast to Finite Element or Difference Methods, BEM are based on boundary integral equations posed in terms of the traces of the solution. For the classical example of the scalar wave equation, the occurring integral operators take the form of so called “retarded potentials” related to Huygen’s principle. In [2], Bamberger and Ha

Seibel
Formulation of the problem
Boundary integral equations
Galerkin formulation
Convolution quadrature method
Galerkin approximation
Hierarchical matrices The boundary element matrices Vk are of the form
Matrix partitions
H 2-matrices
Adaptive cross approximation
Multivariate adaptive cross approximation
Combined algorithm
Error analysis
Complexity and fast arithmetics
Numerical examples
Tensor approximation
Scattering problem
Findings
Conclusion
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.