Boundary element methods for the wave equation based on hierarchical matrices and adaptive cross approximation

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$$ 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.

Transient elastodynamic analysis with a combination of convolution quadrature method and pseudo-initial condition method

Purpose The Convolution Quadrature Method (CQM) has been widely applied to solve transient elastodynamic problems because of its stability and generality. However, the CQM suffers from the problems of huge memory requirement in case of direct implementation in time domain or CPU time in case of its reformulation in Laplace domain. The purpose of this paper is to combine the CQM with the pseudo-initial condition method (PICM) to achieve a good balance between memory requirement and CPU time. Design/methodology/approach The combined methods first subdivide the whole analysis into a few sub-analyses, which is dealt with the PICM, namely, the results obtained by previous sub-analysis are used as the initial conditions for the next sub-analysis. In each sub-analysis, the time interval is further discretized into a number of sub-steps and dealt with the CQM. For non-zero initial conditions, the pseudo-force method is used to transform them into equivalent body forces. The boundary face method is employed in the numerical implementation. Three examples are analyzed. Results are compared with analytical solutions or FEM results and the results of reformulated CQM. Findings Results demonstrate that the computation time and the storage requirement can be reduced significantly as compared to the CQM, by using the combined approach. Originality/value The combined methods can be successfully applied to the problems of long-time dynamic response, which requires a large amount of computer memory when CQM is applied, while preserving the CQM stability. If the number of time steps is high, then the accuracy of the proposed approach can be deteriorated because of the pseudo-force method.

Boundary Element Method with Runge-Kutta Convolution Quadrature for Three-Dimensional Dynamic Poroelasticity

The paper contains a brief introduction to the state of the art in poroelasticity models, in BIE & BEM methods application to solve dynamic problems in Laplace domain. Convolution Quadrature Method is formulated, as well as Runge-Kutta convolution quadrature modification and scheme with a key based on the highly oscillatory quadrature principles. Several approaches to Laplace transform inversion, including based on traditional Euler stepping scheme and Runge-Kutta stepping schemes, are numerically compared. A BIE system of direct approach in Laplace domain is used together with the discretization technique based on the collocation method. The boundary is discretized with the quadrilateral 8-node biquadratic elements. Generalized boundary functions are approximated with the help of the Goldshteyn’s displacement-stress matched model. The time-stepping scheme can rely on the application of convolution theorem as well as integration theorem. By means of the developed software the following 3d poroelastodynamic problem were numerically treated: a Heaviside-shaped longitudinal load acting on the face of a column.

Calculation of Hydrodynamic Forces for Unsteady Stokes Flows by Singularity Integral Equations Based on Fundamental Solutions

ABSTRACTThe attractive feature of the singularity method for steady Stokes flows is that the hydrodynamic forces acting on the particle can be calculated by the total strength of distributed singularities. For unsteady Stokes flows, however we have to derive hydrodynamic forces acting on a solid body in terms of the strengths of both unsteady Stokeslets as well as unsteady potential dipoles if mass and force sources are both taken into consideration. Since the hydrodynamic force formulation results in a Volterra integral equation of the first kind, and the strengths are numerically approximated by means of the Lubich convolution quadrature method (CQM) in this study. As far as the numerical solutions of time-domain integral formulations of the unsteady Stokes equations are concerned, this paper requires only the Laplace-domain instead of the time- domain fundamental solutions of the governing equations. The stability and accuracy of the proposed method are verified through some well selected numerical examples. In total we include two examples presenting the accuracy of Lubich CQM, and another two examples for calculating general hydrody-namic forces of a sphere in oscillating or non-oscillating unsteady Stokes flows. It is concluded that this study is able to extend the unsteady Stokes flow theory to more general transient motions instead to limit to the oscillating flow assumption.

On a hybrid Laplace-time domain approach to dynamic interaction problems

Dynamic interaction problems between two subdomains lead in time to convolution integrals on the interface when one (linear) subdomain is modelled by an impedance operator and the other exhibits non-linear behaviour. In the present work, a convolution quadrature method is used to address a Laplace transform-based approach to evaluate these convolution products. Its properties are discussed on some numerical soil–structure interaction applications: one that is fully linear and the other showing material non-linearities.