Boundary element method and wave equation

In this paper the boundary integral expression for a one-dimensional wave equation with homogeneous boundary conditions is developed. This is done using the time dependent fundamental solution of the corresponding hyperbolic partial differential equation. The integral expression developed is a generalized function with the same form as the well-known D'Alembert formula. The derivatives of the solution and some useful invariants on the characteristics of the partial differential equation are also calculated.The boundary element method is applied to find the numerical solution. The results show excellent agreement with analytical solutions.A multi-step procedure for large time steps which can be used in the boundary element method is also described.In addition, the way in which boundary conditions are introduced during the time dependent process is explained in detail. In the Appendix the main properties of Dirac's delta function and the Heaviside unit step function are described.     

An efficient boundary element method for two-dimensional transient wave propagation problems

The boundary element method (BEM) has been recognized by its unique feature of requiring neither internal cells nor their associated domain integrals in the computation. The method preserves its elegance for transient problems by means of a certain time-stepping scheme that initiates the time integration always from the initial time. Unfortunately, this time-marching scheme becomes rather difficult to apply, because the computation time and storage requirement grow dramatically with the increasing number of time steps. This paper shows that a reduction of one half of the computation time as well as the storage requirement can be achieved by an efficient truncation scheme for two-dimensional transient wave propagation problems. In particular, a guiding parameter for the determination of the truncation limit is proposed, and the overall measure of the error with respect to the truncation guide parameter is established.     

A convolution quadrature Galerkin boundary element method for the exterior Neumann problem of the wave equation

The numerical solution of the Neumann problem of the wave equation on unbounded three‐dimensional domains is calculated using the convolution quadrature method for the time discretization and a Galerkin boundary element method for the spatial discretization. The mathematical analysis that has been built up for the Dirichlet problem is extended and developed for the Neumann problem, which is important for many modelling applications. Numerical examples are then presented for one of these applications, modelling transient acoustic radiation. Copyright © 2009 John Wiley & Sons, Ltd.     

The boundary element method with Lagrangian multipliers

On open surfaces, the energy space of hypersingular operators is a fractional order Sobolev space of order 1/2 with homogeneous Dirichlet boundary condition (along the boundary curve of the surface) in a weak sense. We introduce a boundary element Galerkin method where this boundary condition is incorporated via the use of a Lagrangian multiplier. We prove the quasi‐optimal convergence of this method (it is slightly inferior to the standard conforming method) and underline the theory by a numerical experiment. The approach presented in this article is not meant to be a competitive alternative to the conforming method but rather the basis for nonconforming techniques like the mortar method, to be developed. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

The solution of the non-homogeneous Helmholtz equation by means of the boundary element method

One of the most important advantages of the Boundary Element Method (BEM) is that no internal discretization of the domain is required. This advantage, however, is generally lost when source terms are present in the governing differential equation. It is shown here that for the non-homogeneous Helmholtz equation with a harmonic source term, it is possible to transform the volume integral into a surface integral thus retaining this feature. The transformation is achieved using the Green formula. The technique is applied to solve numerically a test problem with known simple analytical solution.     

A Galerkin boundary element method based on interpolatory Hermite trigonometric wavelets

A Galerkin boundary element method based on interpolatory Hermite trigonometric wavelets is presented for solving 2-D potential problems defined inside or outside of a circular boundary in this paper. In this approach, an equivalent variational form of the corresponding boundary integral equation for the potential problem is used; the trigonometric wavelets are employed as trial and test functions of the variational formulation. The analytical formulae of the matrix entries indicate that most of the matrix entries are naturally zero without any truncation technique and the system matrix is a block diagonal matrix. Each block consists of four circular submatrices. Hence the memory spaces and computational complexity of the system matrix are linear scale. This approach could be easily coupled into domain decomposition method based on variational formulation. Finally, the error estimates of the approximation solutions are given and some test examples are presented.     

Convolution quadrature time-domain boundary element method for 2-D fluid-saturated porous media

For large-scale wave analyses of fluid-saturated porous media, a conventional time-domain boundary element method (BEM) cannot be applied because of the following reasons: (1) no time-domain fundamental solutions are known for some problems, (2) the method sometimes suffers from instability, and (3) the analyses require large amounts of computational time and memory. In this study, an innovative time-domain BEM is developed for a fluid-saturated porous medium. The formulation presented herein overcomes the above disadvantages using a convolution quadrature method (CQM), first proposed by Lubich, and hybrid-parallelization with both MPI and OpenMP. Problems involving the scattering of an incident plane wave by cavities in a 2-D poroelastic medium are solved as a means of validating the proposed method.     

A system of boundary integral equations for the transmission problem in acoustics

In this paper, we reduce the classical two-dimensional transmission problem in acoustic scattering to a system of coupled boundary integral equations (BIEs), and consider the weak formulation of the resulting equations. Uniqueness and existence results for the weak solution of corresponding variational equations are established. In contrast to the coupled system in Costabel and Stephan (1985) [4], we need to take into account exceptional frequencies to obtain the unique solvability. Boundary element methods (BEM) based on both the standard and a two-level fast multipole Galerkin schemes are employed to compute the solution of the variational equation. Numerical results are presented to verify the efficiency and accuracy of the numerical methods.     

An alternative formulation of the boundary element method

The paper suggests an alternative formulation of the Boundary Element Method, in which singular solutions generated by unit dislocations are required and moreover the stresses at the interior points of the body are directly computed from the boundary quantities, without passing through the displacements. Relationships between the singular solutions for unit dislocation and unit force are derived.     

Developing weak Galerkin finite element methods for the wave equation

In this article, we extend the recently developed weak Galerkin method to solve the second‐order hyperbolic wave equation. Many nice features of the weak Galerkin method have been demonstrated for elliptic, parabolic, and a few other model problems. This is the initial exploration of the weak Galerkin method for solving the wave equation. Here we successfully developed and established the stability and convergence analysis for the weak Galerkin method for solving the wave equation. Numerical experiments further support the theoretical analysis. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 868–884, 2017     

求解二维多区域声波散射问题的边界元方法

对于多散射区域的声波散射问题的外Neumann边值问题,用单层位势来逼近每个散射域上的散射波,再利用位势理论的跳跃关系将问题转换为第二类边界积分方程组的求解问题,然后用Nystrom方法进行了求解.对多个随机散射区域的声波散射问题,数值例子体现了该求解方法的可行性和准确性.     

Analysis of the immersed boundary method for a finite element Stokes problem

Convergence results are presented for the immersed boundary (IB) method applied to a model Stokes problem. As a discretization method, we use the finite element method. First, the immersed force field is approximated using a regularized delta function. Its error in the W^{?1, p} norm is examined for 1 ≤ p < n/(n ? 1), with n representing the space dimension. Subsequently, we consider IB discretization of the Stokes problem and examine the regularization and discretization errors separately. Consequently, error estimate of order h^{1 ? α} in the W^{1, 1} × L¹ norm for the velocity and pressure is derived, where α is an arbitrary small positive number. The validity of those theoretical results is confirmed from numerical examples.     

Approximation of boundary values in the boundary element method

Various techniques may be applied to the approximation of the unknown boundary functions involved in the boundary element method (BEM). Several techniques have been examined numerically to find the most efficient. Techniques considered were: Lagrangian polynomials of the zeroth, first and second orders; spline functions; and the novel weighted minimization technique used successfully in the finite difference method (FDM) for arbitrarily irregular meshes. All these approaches have been used in the BEM for the numerical analysis of plates with various boundary conditions.Both coarse and fine grids on the boundary have been assumed. Maximal errors of the deflections of each plate and the bending moments have been found and the effective computer CPU times determined.Analysis of the results showed that, for the same computer time, the greatest accuracy was obtained by the weighted FDM approach. In the case of the Lagrange approximation, higher order polynomials have proved more efficient. The spline technique yielded more accurate results, but with a higher CPU time.Two discretization approaches have been investigated: the least-squares technique and the collocation method. Despite the fact that the simultaneous algebraic equations obtained were not symmetric, the collocation approach has been confirmed as clearly superior to the least-squares technique, because of the amount of computation time used.     

Simple error estimators for the Galerkin BEM for some hypersingular integral equation in 2D

A posteriori error estimation is an important tool for reliable and efficient Galerkin boundary element computations. For hypersingular integral equations in 2D with a positive-order Sobolev space, we analyse the mathematical relation between the (h???h/2)-error estimator from [S. Ferraz-Leite and D. Praetorius, Simple a posteriori error estimators for the h-version of the boundary element method, Computing 83 (2008), pp. 135–162], the two-level error estimator from [M. Maischak, P. Mund, and E. Stephan, Adaptive multilevel BEM for acoustic scattering, 585 Comput. Methods Appl. Mech. Eng. 150 (1997), pp. 351–367], and the averaging error estimator from [C. Carstensen and D. Praetorius, Averaging techniques for the a posteriori bem error control for a hypersingular integral equation in two dimensions, SIAM J. Sci. Comput. 29 (2007), pp. 782–810]. All of these a posteriori error estimators are simple in the following sense: first, the numerical analysis can be done within the same mathematical framework, namely localization techniques for the energy norm. Second, there is almost no implementational overhead for the realization.     

A boundary element method for homogenization of periodic structures

Homogenized coefficients of periodic structures are calculated via an auxiliary partial differential equation in the periodic cell. Typically, a volume finite element discretization is employed for the numerical solution. In this paper, we reformulate the problem as a boundary integral equation using Steklov–Poincaré operators. The resulting boundary element method only discretizes the boundary of the periodic cell and the interface between the materials within the cell. We prove that the homogenized coefficients converge super-linearly with the mesh size, and we support the theory with examples in two and three dimensions.     

Quasi‐local multitrace boundary integral formulations

In the context of time harmonic wave scattering by piecewise homogeneous penetrable objects, we present a new variant of the multitrace boundary integral formulation, that differs from the local multitrace approach of (Jerez‐Hanckes and Hiptmair, 2012) [13] by the presence of regularization terms involving boundary integral operators and localized at junctions i.e. points where at least three subdomains abut. We prove well‐posedness and quasioptimal convergence of conforming Galerkin discretizations for this new formulation, and present numerical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2043–2062, 2015