Efficient solution of two-dimensional wave propagation problems by Cq-Wavelet BEM: Algorithm and applications

In this paper we consider wave propagation problems in two-dimensional unbounded domains, including dissipative effects, reformulated in terms of space-time boundary integral equa- tions. For their solution, we employ a convolution quadrature (CQ) for the temporal and a Galerkin boundary element method (BEM) for the spatial discretization. It is known that one of the main advantages of the CQ-BEMs is the use of the FFT algorithm to retrieve the discrete time integral operators with an optimal linear complexity in time, up to a logarithmic term. It is also known that a key ingredient for the success of such methods is the efficient and accurate evaluation of all the integrals that define the matrix entries associated to the full space-time discretization. This topic has been successfully addressed when standard Lagrangian basis functions are considered for the space discretization. However, it results that, for such a choice of the basis, the BEM matrices are in general fully populated, a drawback that prevents the application of CQ-BEMs to large-scale problems. In this paper, as a possible remedy to reduce the global complexity of the method, we consider approximant functions of wavelet type. In particular, we propose a numerical procedure that, by taking advantage of the fast wavelet transform, allows us on the one hand to compute the matrix entries associated to the choice of wavelet basis functions by maintaining the accuracy of those associated to the Lagrangian basis ones and, on the other hand, to generate sparse matrices without the need of storing a priori the fully populated ones. Such an approach allows in principle the use of wavelet basis of any type and order, combined with CQ based on any stable ordinary differential equations solver. Several numerical results, showing the accuracy of the solution and the gain in terms of computer memory saving, are presented and discussed

A new boundary element integration strategy for retarded potential boundary integral equations

We consider the retarded potential boundary integral equation, arising from the 3D Dirichlet exterior wave equation problem. For its numerical solution we use compactly supported temporal basis functions in time and a standard collocation method in space. Since the accurate computation of the integrals involved in the numerical scheme is a key issue for the numerical stability, we propose a new efficient and competitive quadrature strategy. We compare this approach with the one that uses the Lubich time convolution quadrature, and show pros and cons of both methods

Efficient Solution of Two-Dimensional Wave Propagation Problems by CQ-Wavelet BEM: Algorithm and Applications

In this paper we consider wave propagation problems in two-dimensional unbounded domains, including dissipative effects, reformulated in terms of space-time boundary integral equations. For their solution, we employ a convolution quadrature (CQ) for the temporal and a Galerkin boundary element method (BEM) for the spatial discretization. It is known that one of the main advantages of the CQ-BEMs is the use of the FFT algorithm to retrieve the discrete time integraloperators with an optimal linear complexity in time, up to a logarithmic term. It is also known that a key ingredient for the success of such methods is the efficient and accurate evaluation of all the integrals that define the matrix entries associated to the full space-time discretization. This topic has been successfully addressed when standard Lagrangian basis functions are considered for the space discretization. However, it results that, for such a choice of the basis, the BEM matrices are in general fully populated, a drawback that prevents the application of CQ-BEMs to large-scale problems. In this paper, as a possible remedy to reduce the global complexity of the method, we consider approximant functions of wavelet type. In particular, we propose a numerical procedure that, by taking advantage of the fast wavelet transform, allows us on the one hand to compute the matrix entries associated to the choice of wavelet basis functions by maintaining the accuracy of those associated to the Lagrangian basis ones and, on the other hand, to generate sparse matrices without the need of storing a priori the fully populated ones. Such an approach allows in principle the use of wavelet basis of any type and order, combined with CQ based on any stable ordinary differential equations solver. Several numerical results, showing the accuracy of the solution and the gain in terms of computer memory saving, are presented and discusse

Convolution spline approximations for time domain boundary integral equations

We introduce a new "convolution spline'' temporal approximation of time domain boundary integral equations (TDBIEs). It shares some properties of convolution quadrature (CQ) but, instead of being based on an underlying ODE solver, the approximation is explicitly constructed in terms of compactly supported basis functions. This results in sparse system matrices and makes it computationally more efficient than using the linear multistep version of CQ for TDBIE time-stepping. We use a Volterra integral equation (VIE) to illustrate the derivation of this new approach: at time step t_n = n\dt the VIE solution is approximated in a backwards-in-time manner in terms of basis functions $\phi_j$ by u(t_n-t) \approx \sum_{j=0}^n u_{n-j}\,\phi_j(t/\dt) for $t \in [0,t_n]$. We show that using isogeometric B-splines of degree $m\ge 1$ on $[0,\infty)$ in this framework gives a second order accurate scheme, but cubic splines with the parabolic runout conditions at $t=0$ are fourth order accurate. We establish a methodology for the stability analysis of VIEs and demonstrate that the new methods are stable for non-smooth kernels which are related to convergence analysis for TDBIEs, including the case of a Bessel function kernel oscillating at frequency \oo(1/\dt). Numerical results for VIEs and for TDBIE problems on both open and closed surfaces confirm the theoretical predictions