An error analysis is given for convolution quadratures based on strongly A-stable RungeKutta methods, for the non-sectorial case of a convolution kernel with a Laplace transform that is polynomially bounded in a half-plane. The order of approximation depends on the classical order and stage order of the Runge-Kutta method and on the growth exponent of the Laplace transform. Numerical experiment...

We consider the wave equation in a boundary integral formulation. The discretization in time is done by using convolution quadrature techniques and a Galerkin boundary element method for the spatial discretization. In a previous paper, we have introduced a sparse approximation of the system matrix by cutoff, in order to reduce the storage costs. In this paper, we extend this approach by introdu...

Wave propagation in partially saturated porous continua is an interesting subject in Civil Engineering, Petroleum Engineering, Bioengineering, Earthquake Engineering, and Geophysics, etc. For such problems, there exist different theories, e.g., an extension of Biot’s theory, the Theory of Porous Media and the Mixture theory. Based on the Mixture theory, a dynamic three–phase model for partially...

A boundary element approach based on the Green’s functions in integral representation and the convolution quadrature method is presented. Proposed approach is designed for analyzing 3D initial boundary-value problems of the dynamics of general anisotropic elastic and piezoelectric linear homogeneous solids with mixed boundary conditions. Numerical modelling of transient dynamics of elastic anis...

In this paper, we will present a new method for evaluating high order divided differences for certain classes of analytic, possibly, operator values functions. This is a classical problem in numerical mathematics but also arises in new applications such as, e.g., the use of generalized convolution quadrature to solve retarded potential integral equations. The functions which we will consider ar...

In this paper we propose and analyze a new, fast method for the numerical solution of time domain boundary integral formulations of the wave equation. We employ Lubich’s convolution quadrature method for the time discretization and a Galerkin boundary element method for the spatial discretization. The coefficient matrix of the arising system of linear equations is a triangular block Toeplitz ma...

A fractional time derivative is introduced into Burger’s equation to model losses of nonlinear waves. This term amounts to a time convolution product, which greatly penalizes the numerical modeling. A diffusive representation of the fractional derivative is adopted here, replacing this nonlocal operator by a continuum of memory variables that satisfy local-in-time ordinary differential equation...