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.     

Adaptive Hybridized Interior Penalty Discontinuous Galerkin Methods for H(curl)-Elliptic Problems

We develop and analyze an adaptive hybridized Interior Penalty Discontinuous Galerkin (IPDG-H) method for H(curl)-elliptic boundary value problems in 2D or 3D arising from a semi-discretization of the eddy currents equations. The method can be derived from a mixed formulation of the given boundary value problem and involves a Lagrange multiplier that is an approximation of the tangential traces of the primal variable on the interfaces of the underlying triangulation of the computational domain. It is shown that the IPDG-H technique can be equivalently formulated and thus implemented as a mortar method. The mesh adaptation is based on a residual-type a posteriori error estimator consisting of element and face residuals. Within a unified framework for adaptive finite element methods, we prove the reliability of the estimator up to a consistency error. The performance of the adaptive symmetric IPDG-H method is documented by numerical results for representative test examples in 2D.     

A posteriori boundary element error estimation

An a posteriori error estimator is presented for the boundary element method in a general framework. It is obtained by solving local residual problems for which a local concept is introduced to accommodate the fact that integral operators are nonlocal operators. The estimator is shown to have an upper and a lower bound by the constant multiples of the exact error in the energy norm for Symm's and hypersingular integral equations. Numerical results are also given to demonstrate the effectiveness of the estimator for these equations. It can be used for adaptive h,p, and hp methods.     

MLPG method for two-dimensional diffusion equation with Neumann's and non-classical boundary conditions

In this paper, a meshless local Petrov-Galerkin (MLPG) method is presented to treat parabolic partial differential equations with Neumann's and non-classical boundary conditions. A difficulty in implementing the MLPG method is imposing boundary conditions. To overcome this difficulty, two new techniques are presented to use on square domains. These techniques are based on the finite differences and the Moving Least Squares (MLS) approximations. Non-classical integral boundary condition is approximated using Simpson's composite numerical integration rule and the MLS approximation. Two test problems are presented to verify the efficiency and accuracy of the method.     

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

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

Analytical formulation of meshless local integral equation method

In this paper, the exact forms of integrals in the meshless local boundary integral equation method are derived and implemented for elastostatic problems. A weak form for a set of governing equations with a unit test function or polynomial test functions is transformed into local integral equations. Each node has its own support domain and is surrounded by a local integral domain with different shapes of boundaries. The meshless approximation based on the radial basis function (RBF) is employed for the implementation of displacements. A completed set of closed forms of the local boundary integrals are obtained. As there are no numerical integrations to be carried out the computational time is significantly reduced. Three examples are presented to demonstrate the application of this approach in solid mechanics.     

Boundary element methods for potential problems with nonlinear boundary conditions

Galerkin boundary element methods for the solution of novel first kind Steklov-Poincaré and hypersingular operator boundary integral equations with nonlinear perturbations are investigated to solve potential type problems in two- and three-dimensional Lipschitz domains with nonlinear boundary conditions. For the numerical solution of the resulting Newton iterate linear boundary integral equations, we propose practical variants of the Galerkin scheme and give corresponding error estimates. We also discuss the actual implementation process with suitable preconditioners and propose an optimal hybrid solution strategy.

     

Convergence of adaptive 3D BEM for weakly singular integral equations based on isotropic mesh‐refinement

We consider the adaptive lowest‐order boundary element method based on isotropic mesh refinement for the weakly‐singular integral equation for the three‐dimensional Laplacian. The proposed scheme resolves both, possible singularities of the solution as well as of the given data. The implementation thus only deals with discrete integral operators, that is, matrices. We prove that the usual adaptive mesh‐refining algorithm drives the corresponding error estimator to zero. Under an appropriate saturation assumption which is observed empirically, the sequence of discrete solutions thus tends to the exact solution within the energy norm. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Adaptive space-time boundary element methods for the wave equation

For the solution of the wave equation a space-time energetic boundary integral formulation is used. The resulting single layer boundary integral equation is discretised by a conforming ansatz space on the lateral boundary. To derive an adaptive scheme an a posteriori error estimator based on the representation formula is used. Finally, numerical examples for a one-dimensional spatial domain are presented. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Use of radial basis functions for solving the second‐order parabolic equation with nonlocal boundary conditions

Nonlocal mathematical models appear in various problems of physics and engineering. In these models the integral term may appear in the boundary conditions. In this paper the problem of solving the one‐dimensional parabolic partial differential equation subject to given initial and nonlocal boundary conditions is considered. These kinds of problems have certainly been one of the fastest growing areas in various application fields. The presence of an integral term in a boundary condition can greatly complicate the application of standard numerical techniques. As a well‐known class of meshless methods, the radial basis functions are used for finding an approximation of the solution of the present problem. Numerical examples are given at the end of the paper to compare the efficiency of the radial basis functions with famous finite‐difference methods. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A posteriori error analysis for a boundary element method with nonconforming domain decomposition

We present and analyze an a posteriori error estimator based on mesh refinement for the solution of the hypersingular boundary integral equation governing the Laplacian in three dimensions. The discretization under consideration is a nonconforming domain decomposition method based on the Nitsche technique. Assuming a saturation property, we establish quasireliability and efficiency of the error estimator in comparison with the error in a natural (nonconforming) norm. Numerical experiments with uniform and adaptively refined meshes confirm our theoretical results. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 947–963, 2014     

Shape-explicit constants for some boundary integral operators

Among the well-known constants in the theory of boundary integral equations are the coercivity constants of the single-layer potential and the hypersingular boundary integral operator, and the contraction constant of the double-layer potential. Whereas there have been rigorous studies how these constants depend on the size and aspect ratio of the underlying domain, only little is known on their dependency on the shape of the boundary. In this article, we consider the homogeneous Laplace equation and derive explicit estimates for the above-mentioned constants in three dimensions. Using an alternative trace norm, we make the dependency explicit in two geometric parameters, the so-called Jones parameter and the constant in Poincaré's inequality. The latter one can be tracked back to the constant in an isoperimetric inequality. There are many domains with quite irregular boundaries, where these parameters stay bounded. Our results provide a new tool in the analysis of numerical methods for boundary integral equations and in particular for boundary element based domain decomposition methods.     

Energetic boundary element method analysis of wave propagation in 2D multilayered media

The analysis of scalar wave propagation in 2D zonewise homogeneous media with vanishing initial and mixed boundary conditions is carried out. The problem is formulated in terms of time‐dependent boundary integral equations, and then it is set in a weak form, based on a natural energy identity satisfied by the differential problem solution. Several numerical results have been obtained by means of the related energetic Galerkin boundary element method showing accuracy and stability of the method. Copyright © 2012 John Wiley & Sons, Ltd.     

Boundary elements solution of stokes flow between curved surfaces with nonlinear slip boundary condition

This work presents a boundary integral equation formulation for Stokes nonlinear slip flows based on the normal and tangential projection of the Green's integral representational formulae for the velocity field. By imposing the surface tangential velocity discontinuity (slip velocity) in terms of the nonlinear slip flow boundary condition, a system of nonlinear boundary integral equations for the unknown normal and tangential components of the surface traction is obtained. The Boundary Element Method is used to solve the resulting system of integral equations using a direct Picard iteration scheme to deal with the resulting nonlinear terms. The formulation is used to study flows between curved rotating geometries: i.e., concentric and eccentric Couette flows and single rotor mixers, under nonlinear slip boundary conditions. The numerical results obtained for the concentric Couette flow is validated again a semianalytical solution of the problem, showing excellent agreements. The other two cases, eccentric Couette and single rotor mixers, are considered to study the effect of different nonlinear slip conditions in these flow configurations. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

An hp‐adaptive refinement strategy for hypersingular operators on surfaces

An adaptive refinement strategy for the hp‐version of the boundary element method with hypersingular operators on surfaces is presented. The error indicators are based on local projections provided by two‐level decompositions of ansatz spaces with additional bubble functions. Assuming a saturation property and locally quasi‐uniform meshes, efficiency and reliability of the resulting error estimator is proved. A second error estimator based on mesh refinement and overlapping decompositions that better fulfills the saturation property is presented. The performance of the algorithm and the estimators is demonstrated for a model problem. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 396–419, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10011     

An a posteriori error estimate for a first-kind integral equation

In this paper we present a new a posteriori error estimate for the boundary element method applied to an integral equation of the first kind. The estimate is local and sharp for quasi-uniform meshes and so improves earlier work of ours. The mesh-dependence of the constants is analyzed and shown to be weaker than expected from our previous work. Besides the Galerkin boundary element method, the collocation method and the qualocation method are considered. A numerical example is given involving an adaptive feedback algorithm.

     

On the numerical method for solving a hypersingular integral equation with the computation of the solution gradient

We consider a linear integral equation, which arises when solving the Neumann boundary value problem for the Laplace equation with the representation of the solution in the form of a double layer potential, with a hypersingular integral treated in the sense of Hadamard finite value. We consider the case in which the exterior or interior problem is solved in a domain whose boundary is a closed smooth surface and the integral equation is written over that surface. A numerical scheme for solving the integral equation is constructed with the use of quadrature formulas of the type of the method of discrete singularities with a regularization for the use of an irregular grid. We prove the convergence, uniform over the grid points, of the numerical solutions to the exact solution of the hypersingular equation and, in addition, the uniform convergence of the values of the approximate finite-difference derivative operator on the numerical solution to the values on the projection of the exact solution onto the subspace of grid functions with nodes at the collocation points.     

Convergence of the piecewise linear approximation and collocation method for a hypersingular integral equation on a closed surface

We study the numerical solution of a linear hypersingular integral equation arising when solving the Neumann boundary value problem for the Laplace equation by the boundary integral equation method with the solution represented in the form of a double layer potential. The integral in this equation is understood in the sense of Hadamard finite value. We construct quadrature formulas for the integral occurring in this equation based on a triangulation of the surface and an application of the linear approximation to the unknown function on each of the triangles approximating the surface. We prove the uniform convergence of the quadrature formulas at the interpolation nodes as the triangulation size tends to zero. A numerical solution scheme for this integral equation based on the suggested quadrature formulas and the collocation method is constructed. Under additional assumptions about the shape of the surface, we prove a uniform estimate for the error in the numerical solution at the interpolation nodes.     

The fast multipole method for the symmetric boundary integral formulation

^** Email: of{at}mathematik.uni-stuttgart.de^*** Email: o.steinbach{at}tugraz.at^**** Email: wendland{at}mathematik.uni-stuttgart.de A symmetric Galerkin boundary-element method is used for thesolution of boundary-value problems with mixed boundary conditionsof Dirichlet and Neumann type. As a model problem we considerthe Laplace equation. When an iterative scheme is employed forsolving the resulting linear system, the discrete boundary integraloperators are realized by the fast multipole method. While thesingle-layer potential can be implemented straightforwardlyas in the original algorithm for particle simulation, the double-layerpotential and its adjoint operator are approximated by the applicationof normal derivatives to the multipole series for the kernelof the single-layer potential. The Galerkin discretization ofthe hypersingular integral operator is reduced to the single-layerpotential via integration by parts. We finally present a correspondingstability and error analysis for these approximations by thefast multipole method of the boundary integral operators. Itis shown that the use of the fast multipole method does notharm the optimal asymptotic convergence. The resulting linearsystem is solved by a GMRES scheme which is preconditioned bythe use of hierarchical strategies as already employed in thefast multipole method. Our numerical examples are in agreementwith the theoretical results.     

A Burton-Miller boundary element-free method for Helmholtz problems

A Burton-Miller boundary element-free method is developed by using the Burton-Miller formulation for meshless and boundary-only analysis of Helmholtz problems. The method can produce a unique solution at all wavenumbers and is valid for Dirichlet, Neumann and mixed problems simultaneously. An efficient numerical integration procedure is presented to handle both strongly singular and hypersingular boundary integrals directly and uniformly. Numerical results reveal that this direct meshless method only involves boundary nodes and can deal with Helmholtz problems at extremely large wavenumbers.