A symmetric Galerkin boundary element method for 3d linear poroelasticity

In this paper, a symmetric Galerkin boundary element formulation for 3D linear poroelasticity is presented. By means of the convolution quadrature method, the time domain problem is decoupled into a set of Laplace domain problems. Regularizing their kernel functions via integration by parts, it is possible to compute all operators for rather general discretizations, only requiring the evaluation of weakly singular integrals. At the end, some numerical results are presented and compared with a collocation BEM. Throughout these studies, the symmetric Galerkin BEM performs better than the collocation method, especially for not optimal discretizations parameters, i.e. a bad relation of mesh to time-step size. The most obvious advantages can be observed in the fluid flux results. However, these advantages are obtained at a higher numerical cost.