A boundary element Galerkin method for a hypersingular integral equation on open surfaces

A hypersingular boundary integral equation of the first kind on an open surface piece Γ is solved approximately using the Galerkin method. As boundary elements on rectangles we use continuous, piecewise bilinear functions which vanish on the boundary of Γ. We show how to compensate for the effect of the edge and corner singularities of the true solution of the integral equation by using an appropriately graded mesh and obtain the same convergence rate as for the case of a smooth solution. We also derive asymptotic error estimates in lower-order Sobolev norms via the Aubin–Nitsche trick. Numerical experiments for the Galerkin method with piecewise linear functions on triangles demonstrate the effect of graded meshes and show experimental rates of convergence which underline the theoretical results.     

The h-p version of the coupling of finite element and boundary element methods for transmission problems in polyhedral domains

Summary. This paper analyzes the rate of convergence of the h-p version of the coupling of the finite element and boundary element method for transmission problems with a linear differential operator with variable coefficients in a bounded polyhedral domain and with constant coefficients in the exterior domain . This procedure uses the variational formulation of the differential equation in and involves integral operators on the interface between and . The finite elements are used to obtain approximate solutions of the differential equation in and the boundary elements are used to obtain approximate solutions of the integral equations. For given piecewise analytic data we show that the Galerkin solution of this coupling procedure converges exponentially fast in the energy norm if the h-p version is used both for finite elements and boundary elements. Received February 10, 1996 / Revised version received April 4, 1997     

Extrapolation for the Approximations to the Solution of a Boundary Integral Equation on Polygonal Domains

In this paper, we consider a boundary integral equation of second kind rising from potential theory. The equation may be solved numerically by Galerkin's method using piecewise constant functions. Because of the singularities produced by the corners, we have to grade the mesh near the corner. In general, Chandler obtained the order 2 superconvergence of the iterated Galerkin solution in the uniform norm. It is proved in this paper that the Richardson extrapolation increases the accuracy from order 2 to order 4.     

Artificial Multilevel Boundary Element Preconditioners

A hierarchical multilevel preconditioner is constructed for an efficient solution of a first kind boundary integral equation with the single layer potential operator discretized by a boundary element method. This technique is based on a hierarchical clustering of all boundary elements as used in fast boundary element methods. This hierarchy is applied to define a sequence of nested boundary element spaces of piecewise constant basis functions as used in the definition of the preconditioning multilevel operator.     

A new frequency‐uniform coercive boundary integral equation for acoustic scattering

A new boundary integral operator is introduced for the solution of the soundsoft acoustic scattering problem, i.e., for the exterior problem for the Helmholtz equation with Dirichlet boundary conditions. We prove that this integral operator is coercive in L²(Γ) (where Γ is the surface of the scatterer) for all Lipschitz star‐shaped domains. Moreover, the coercivity is uniform in the wavenumber k = ω/c, where ω is the frequency and c is the speed of sound. The new boundary integral operator, which we call the “star‐combined” potential operator, is a slight modification of the standard combined potential operator, and is shown to be as easy to implement as the standard one. Additionally, to the authors' knowledge, it is the only second‐kind integral operator for which convergence of the Galerkin method in L²(Γ) is proved without smoothness assumptions on Γ except that it is Lipschitz. The coercivity of the star‐combined operator implies frequency‐explicit error bounds for the Galerkin method for any approximation space. In particular, these error estimates apply to several hybrid asymptoticnumerical methods developed recently that provide robust approximations in the high‐frequency case. The proof of coercivity of the star‐combined operator critically relies on an identity first introduced by Morawetz and Ludwig in 1968, supplemented further by more recent harmonic analysis techniques for Lipschitz domains. © 2011 Wiley Periodicals, Inc.     

Edge asymptotics for the radiosity equation over polyhedral boundaries

In the present paper we consider the radiosity equation over the boundary of a polyhedral domain. Similarly to corresponding results on the double‐layer potential equation, the solution of the second kind integral equation with non‐compact integral operator is piecewise continuous. The partial derivatives, however, are not bounded. In the present paper we derive the first term in the asymptotic expansion of the solution in the vicinity of an edge. Note that, knowing this term, optimal mesh gradings can be designed for the numerical solution of this equation. Copyright © 1999 John Wiley & Sons, Ltd.     

含开边界二维Stokes问题的Galerkin边界元解法

本文推导了含有开边界的二维有限域上Stokes问题的边界积分方程, 得出基于单层位势的第一类间接边界积分方程.对与之等价的边界变分方程用Galerkin边界元求解以得出单层位势的向量密度. 对于含有开边界端点的边界单元,采用特别的插值函数, 以模拟其固有的奇异性.论文用若干数值算例模拟了含有开边界的有限区域上不可压缩粘性流体的绕流.       

Defect correction from a Galerkin viewpoint

Summary We consider the numerical solution of systems of nonlinear two point boundary value problems by Galerkin's method. An initial solution is computed with piecewise linear approximating functions and this is then improved by using higher-order piecewise polynomials to compute defect corrections. This technique, including numerical integration, is justified by typical Galerkin arguments and properties of piecewise polynomials rather than the traditional asymptotic error expansions of finite difference methods.     

Error estimates for a discretized Galerkin method for a boundary integral equation in two dimensions

We present a priori and a posteriori estimates for the error between the Galerkin and a discretized Galerkin method for the boundary integral equation for the single layer potential on the square plate. Using piecewise constant finite elements on a rectangular mesh we study the error coming from numerical integration. The crucial point of our analysis is the estimation of some error constants, and we demonstrate that this is necessary if our methods are to be used. After the determination of these constants we are in the position to prove invertibility and quasioptimal convergence results for our numerical scheme, if the chosen numerical integration formulas are sufficiently precise. © 1992 John Wiley & Sons, Inc.     

On spline collocation methods for boundary integral equations in the plane

Nowadays boundary elemen; methods belong to the most popular numerical methods for solving elliptic boundary value problems. They consist in the reduction of the problem to equivalent integral equations (or certain generalizations) on the boundary Γ of the given domain and the approximate solution of these boundary equations. For the numerical treatment the boundary surface is decomposed into a finite number of segments and the unknown functions are approximated by corresponding finite elements and usually determined by collocation and Galerkin procedures. One finds the least difficulties in the theoretical foundation of the convergence of Galerkin methods for certain classes of equations, whereas the convergence of collocation methods, which are mostly used in numerical computations, has yet been proved only for special equations and methods. In the present paper we analyse spline collocation methods on uniform meshes with variable collocation points for one-dimensional pseudodifferential equations on a closed curve with convolutional principal parts, which encompass many classes of boundary integral equations in the plane. We give necessary and sufficient conditions for convergence and prove asymptotic error estimates. In particular we generalize some results on nodal and midpoint collocation obtained in [2], [7] and [8]. The paper is organized as follows. In Section 1 we formulate the problems and the results, Section 2 deals with spline interpolation in periodic Sobolev spaces, and in Section 3 we prove the convergence theorems for the considered collocation methods.     

Coercivity of Combined Boundary Integral Equations in High‐Frequency Scattering

We prove that the standard second‐kind integral equation formulation of the exterior Dirichlet problem for the Helmholtz equation is coercive (i.e., sign‐definite) for all smooth convex domains when the wavenumber k is sufficiently large. (This integral equation involves the so‐called combined potential, or combined field, operator.) This coercivity result yields k‐explicit error estimates when the integral equation is solved using the Galerkin method, regardless of the particular approximation space used (and thus these error estimates apply to several hybrid numerical‐asymptotic methods developed recently). Coercivity also gives k‐explicit bounds on the number of GMRES iterations needed to achieve a prescribed accuracy when the integral equation is solved using the Galerkin method with standard piecewise‐polynomial subspaces. The coercivity result is obtained by using identities for the Helmholtz equation originally introduced by Morawetz in her work on the local energy decay of solutions to the wave equation. © 2015 Wiley Periodicals, Inc.     

Multiple traces boundary integral formulation for Helmholtz transmission problems

We present a novel boundary integral formulation of the Helmholtz transmission problem for bounded composite scatterers (that is, piecewise constant material parameters in ??subdomains??) that directly lends itself to operator preconditioning via Calderón projectors. The method relies on local traces on subdomains and weak enforcement of transmission conditions. The variational formulation is set in Cartesian products of standard Dirichlet and special Neumann trace spaces for which restriction and extension by zero are well defined. In particular, the Neumann trace spaces over each subdomain boundary are built as piecewise $\widetilde{H}^{-1/2}$ -distributions over each associated interface. Through the use of interior Calderón projectors, the problem is cast in variational Galerkin form with an operator matrix whose diagonal is composed of block boundary integral operators associated with the subdomains. We show existence and uniqueness of solutions based on an extension of Lions?? projection lemma for non-closed subspaces. We also investigate asymptotic quasi-optimality of conforming boundary element Galerkin discretization. Numerical experiments in 2-D confirm the efficacy of the method and a performance matching that of another widely used boundary element discretization. They also demonstrate its amenability to different types of preconditioning.     

Discrete Galerkin method for Fredholm integro-differential equations with weakly singular kernels

Approximations to a solution and its derivatives of a boundary value problem of an nth order linear Fredholm integro-differential equation with weakly singular or other nonsmooth kernels are determined. These approximations are piecewise polynomial functions on special graded grids. For their finding a discrete Galerkin method and an integral equation reformulation of the boundary value problem are used. Optimal global convergence estimates are derived and an improvement of the convergence rate of the method for a special choice of parameters is obtained. To illustrate the theoretical results a collection of numerical results of a test problem is presented.     

A boundary integral equation method for three-dimensional crack problems in elasticity

This paper presents a solution procedure for three-dimensional crack problems via first kind boundary integral equations on the crack surface. The Dirichlet (Neumann) problem is reduced to a system of integral equations for the jump of the traction (of the field) across the crack surface. The calculus of pseudodifferential operators is used to derive existence and regularity of the solutions of the integral equations. With the concept of the principal symbol and the Wiener-Hopf technique we derive the explicit behavior of the densities of the integral equations near the edge of the crack surface. Based on the detailed regularity results we show how to improve the boundary element Galerkin method for our integral equations. Quasi-optimal asymptotic estimates for the Galerkin error are given.     

Wavelet approximations for first kind boundary integral equations on polygons

Summary. An elliptic boundary value problem in the interior or exterior of a polygon is transformed into an equivalent first kind boundary integral equation. Its Galerkin discretization with degrees of freedom on the boundary with spline wavelets as basis functions is analyzed. A truncation strategy is presented which allows to reduce the number of nonzero elements in the stiffness matrix from to entries. The condition numbers are bounded independently of the meshwidth. It is proved that the compressed scheme thus obtained yields in operations approximate solutions with the same asymptotic convergence rates as the full Galerkin scheme in the boundary energy norm as well as in interior points. Numerical examples show the asymptotic error analysis to be valid already for moderate values of . Received March 12, 1994 / Revised version received January 9, 1995     

A finite element method for the solution of a potential theory integral equation

This paper discusses a finite element approximation for an integral equation of the second kind deduced from a potential theory boundary value problem in two variables. The equation is shown to admit a unique solution, to be variational and coercive in the Hilbert space of functions σ ε H^1/2(Γ), f_rd γ = 0. The Galerkin method with finite elements as trial functions is shown to lead to an optimal rate of convergence.     

Boundary element approximation for Maxwell's eigenvalue problem

We introduce a new method for computing eigenvalues of the Maxwell operator with boundary finite elements. On bounded domains with piecewise constant material coefficients, the Maxwell solution for fixed wave number can be represented by boundary integrals, which allows to reduce the eigenvalue problem to a nonlinear problem for determining the wave number along with boundary and interface traces. A Galerkin discretization yields a smooth nonlinear matrix eigenvalue problem that is solved by Newton's method or, alternatively, the contour integral method. Several numerical results including an application to the band structure computation of a photonic crystal illustrate the efficiency of this approach. Copyright © 2013 John Wiley & Sons, Ltd.     

Richardson extrapolation of iterated discrete Galerkin solution for two-dimensional Fredholm integral equations

In this paper, we study the numerical solution of two-dimensional Fredholm integral equation by discrete Galerkin and iterated discrete Galerkin method. We are able to derive an asymptotic error expansion of the iterated discrete Galerkin solution. This expansion covers arbitrarily high powers of the discretization parameters if the solution of the integral equation is smooth. The expansion gives rise to Richardson-type extrapolation schemes which rapidly improve the original rate of the convergence. Numerical experiments confirm our theoretical results.     

关于薄板的无网格局部边界积分方程方法中的友解

无网格局部边界积分方程方法是最近发展起来的一种新的数值方法,这种方法综合了伽辽金有限元、边界元和无单元伽辽金法的优点,是一种具有广阔应用前景的、真正的无网格方法．把无网格局部边界积分方程方法应用于求解薄板问题,给出了薄板无网格局部边界积分方程方法所需要的友解及其全部公式．     

Hybrid Galerkin boundary elements: theory and implementation

Summary. In this paper we present a new quadrature method for computing Galerkin stiffness matrices arising from the discretisation of 3D boundary integral equations using continuous piecewise linear boundary elements. This rule takes as points some subset of the nodes of the mesh and can be used for computing non-singular Galerkin integrals corresponding to pairs of basis functions with non-intersecting supports. When this new rule is combined with standard methods for the singular Galerkin integrals we obtain a “hybrid” Galerkin method which has the same stability and asymptotic convergence properties as the true Galerkin method but a complexity more akin to that of a collocation or Nystr?m method. The method can be applied to a wide range of singular and weakly-singular first- and second-kind equations, including many for which the classical Nystr?m method is not even defined. The results apply to equations on piecewise-smooth Lipschitz boundaries, and to non-quasiuniform (but shape-regular) meshes. A by-product of the analysis is a stability theory for quadrature rules of precision 1 and 2 based on arbitrary points in the plane. Numerical experiments demonstrate that the new method realises the performance expected from the theory. Received January 22, 1998 / Revised version received May 26, 1999 / Published online April 20, 2000 –? Springer-Verlag 2000