Adaptive meshless Galerkin boundary node methods for hypersingular integral equations

Adaptive refinement techniques are developed in this paper for the meshless Galerkin boundary node method for hypersingular boundary integral equations. Two types of error estimators are derived. One is a perturbation error estimator that is formulated based on the difference between numerical solutions obtained using two consecutive nodal arrangements. The other is a projection error estimator that is formulated based on the difference between the numerical solution itself and its projection. These error estimators are proven to have an upper and a lower bound by the constant multiples of the exact error in the energy norm. A localization scheme is presented to accomodate the non-local property of hypersingular integral operators for the needed computable local error indicators. The convergence of the adaptive meshless techniques is verified theoretically. To confirm the theoretical results and to show the efficiency of the adaptive techniques, numerical examples in 2D and 3D with high singularities are provided.     

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.     

Extrapolation algorithms for solving nonlinear boundary integral equations by mechanical quadrature methods

We study the numerical solution procedure for two-dimensional Laplace’s equation subjecting to non-linear boundary conditions. Based on the potential theory, the problem can be converted into a nonlinear boundary integral equations. Mechanical quadrature methods are presented for solving the equations, which possess high accuracy order O(h ³) and low computing complexities. Moreover, the algorithms of the mechanical quadrature methods are simple without any integration computation. Harnessing the asymptotical compact theory and Stepleman theorem, an asymptotic expansion of the errors with odd powers is shown. Based on the asymptotic expansion, the h ³ −Richardson extrapolation algorithms are used and the accuracy order is improved to O(h ⁵). The efficiency of the algorithms is illustrated by numerical examples.     

Hierarchical basis methods for hypersingular integral equations

In this paper we present hierarchical basis methods for theGalerkin approximation of hypersingular integral equations onthe interval = (–1,1). The condition number of the stiffnessmatrix with respect to the hierarchical basis is shown to behavelike O(|logh|²). The implementations are based on the preconditionedconjugate gradient method using a hierarchical basis (HB) preconditioner.The numerical results are presented with a comparison betweenthe HB preconditioner and the BPX (Bramble, Pasciak and Xu)preconditioner.     

Adaptive boundary element methods for strongly elliptic integral equations

Summary Integral operators are nonlocal operators. The operators defined in boundary integral equations to elliptic boundary value problems, however, are pseudo-differential operators on the boundary and, therefore, provide additional pseudolocal properties. These allow the successful application of adaptive procedures to some boundary element methods. In this paper we analyze these methods for general strongly elliptic integral equations and obtain a-posteriori error estimates for boundary element solutions. We also apply these methods to nodal collocation with odd degree splines. Some numerical examples show that these adaptive procedures are reliable and effective.This work was carried out while Dr. De-hao Yu was an Alexander-von-Humboldt-Stiftung research fellow at the University of Stuttgart in 1987, 1988     

An iterative substructuring method for the p-version of the boundary element method for hypersingular integral operators in three dimensions

Summary. We study preconditioners for the -version of the boundary element method for hypersingular integral equations in three dimensions. The preconditioners are based on iterative substructuring of the underlying ansatz spaces which are constructed by using discretely harmonic basis functions. We consider a so-called wire basket preconditioner and a non-overlapping additive Schwarz method based on the complete natural splitting, i.e. with respect to the nodal, edge and interior functions, as well as an almost diagonal preconditioner. In any case we add the space of piecewise bilinear functions which eliminate the dependence of the condition numbers on the mesh size. For all these methods we prove that the resulting condition numbers are bounded by . Here, is the polynomial degree of the ansatz functions and is a constant which is independent of and the mesh size of the underlying boundary element mesh. Numerical experiments supporting these results are reported. Received July 8, 1996 / Revised version received January 8, 1997     

Boundary element preconditioners for a hypersingular integral equation on an interval

We propose an almost optimal preconditioner for the iterative solution of the Galerkin equations arising from a hypersingular integral equation on an interval. This preconditioning technique, which is based on the single layer potential, was already studied for closed curves [11,14]. For a boundary element trial space, we show that the condition number is of order (1 + | log h _min|)², where h _min is the length of the smallest element. The proof requires only a mild assumption on the mesh, easily satisfied by adaptive refinement algorithms.     

Newton methods for a class of nonlinear hypersingular integral equations

Different iterative schemes based on collocation methods have been well studied and widely applied to the numerical solution of nonlinear hypersingular integral equations (Capobianco et al. 2005). In this paper we apply Newton’s method and its modified version to solve the equations obtained by applying a collocation method to a nonlinear hypersingular integral equation of Prandtl’s type. The corresponding convergence results are derived in suitable Sobolev spaces. Some numerical tests are also presented to validate the theoretical results.     

An additive Schwarz method for the h-p version of the boundary element method for hypersingular integral equations in R3

We study a preconditioner for the h-p version of the boundaryelement method for hypersingular integral equations in threedimensions. The preconditioner is based on a three-level decompositionof the underlying ansatz space, the levels being piecewise bilinearfunctions on a coarse grid, piecewise bilinear functions ona fine grid, and piecewise polynomials of high degree on thefine grid. We prove that the condition number of the preconditionedlinear system is bounded by max_j (1 + log H_jp_j/h_j)² where H_jis the diameter of an element _j of the coarse grid, h_j is thesize of the elements of the fine grid on _j, and p_j is the maximumof the polynomial degrees used in _j. Numerical results supportingour theory are reported. Received 9 March 1999. Accepted 19 July 1999.     

Composition of hypersingular integral operators

By using the notion of the S-asymptotic of distributions, we give several assertions on the ordinary asymptotic behaviour of the distributional Weierstrass transform. The obtained results are ap-plied on the Cauchy problem for the heat equation     

Coupling of boundary integral equation and finite element methods for transmission problems in acoustics

Numerical Algorithms - In this paper, we propose a coupling of finite element method (FEM) and boundary integral equation (BIE) method for solving acoustic transmission problems in two dimensions....     

Multigrid methods for boundary integral equations

Summary Multigrid methods are applied for solving algebraic systems of equations that occur to the numerical treatment of boundary integral equations of the first and second kind. These methods, originally formulated for partial differential equations of elliptic type, combine relaxation schemes and coarse grid corrections. The choice of the relaxation scheme is found to be essential to attain a fast convergent iterative process. Theoretical investigations show that the presented relaxation scheme provides a multigrid algorithm of which the rate of convergence increases with the dimension of the finest grid. This is illustrated for the calculation of potential flow around an aerofoil.     

The coupling of boundary integral and finite element methods for the nonstationary exterior flow

In this article, we represent a new numerical method for solving the nonstationary Navier–Stokes equations in an unbounded domain. The technique consists of coupling the boundary integral and the finite element method. The variational formulation and the well-posedness of the coupling method are obtained. The convergence and optimal error estimates for the approximate solution are provided. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 549–565, 1998     

On the coupling of boundary integral and finite element methods with nonlinear transmission conditions

We apply the coupling of boundary integral and finite element methods to study the weak solvability of certain exterior boundary value problems with nonlinear transmission conditions. As a model we consider a nonlinear second order elliptic equation in divergence form in a bounded inner region of the plane coupled with the Laplace equation in the corresponding exterior domain. The flux jump across the common nonlinear-linear interface is unknown and assumed to depend nonlinearly on the Dirichlet data. We establish the associated variational formulation in an operator equation setting and provide existence, uniqueness and approximation results.     

Effective semi-analytic integration for hypersingular Galerkin boundary integral equations for the Helmholtz equation in 3D

We deal with the Galerkin discretization of the boundary integral equations corresponding to problems with the Helmholtz equation in 3D. Our main result is the semi-analytic integration for the bilinear form induced by the hypersingular operator. Such computations have already been proposed for the bilinear forms induced by the single-layer and the double-layer potential operators in the monograph The Fast Solution of Boundary Integral Equations by O. Steinbach and S. Rjasanow and we base our computations on these results.     

Calculation of line integrals in boundary integral methods

We propose quadrature rules for the approximation of line integrals possessing logarithmic singularities and show their convergence. In some instances a superconvergence rate is demonstrated. This revised version was published online in June 2006 with corrections to the Cover Date.     

Qualocation methods for elliptic boundary integral equations

Summary. The qualocation methods developed in this paper, with spline trial and test spaces, are suitable for classes of boundary integral equations with convolutional principal part, on smooth closed curves in the plane. Some of the methods are suitable for all strongly elliptic equations; that is, for equations in which the even symbol part of the operator dominates. Other methods are suitable when the odd part dominates. Received December 27, 1996 / Revised version received April 14, 1997     

A solution method for hypersingular integral equations

The Neumann problem for the Helmholtz equation is considered. The double-layer potential is used to reduce the problem to a hypersingular integral equation. The properties of the hypersingular operator in a neighborhood lead to a method for approximate solution of the hypersingular equation with an arbitrary contour. Some numerical results are reported.Translated from Matematicheskoe Modelirovanie i Reshenie Obratnykh Zadach. Matematicheskoi Fiziki, pp. 130–136, 1993.