Numerical integration schemes for the BEM solution of hypersingular integral equations

In this paper we consider singular and hypersingular integral equations associated with 2D boundary value problems defined on domains whose boundaries have piecewise smooth parametric representations. In particular, given any (polynomial) local basis, we show how to compute efficiently all integrals required by the Galerkin method. The proposed numerical schemes require the user to specify only the local polynomial degrees; therefore they are quite suitable for the construction of p‐ and h‐p versions of Galerkin BEM. Copyright © 1999 John Wiley & Sons, Ltd.     

A pure boundary element method approach for solving hypersingular boundary integral equations

This paper is concerned with discretization and numerical solution of a regularized version of the hypersingular boundary integral equation (HBIE) for the two-dimensional Laplace equation. This HBIE contains the primary unknown, as well as its gradient, on the boundary of a body. Traditionally, this equation has been solved by combining the boundary element method (BEM) together with tangential differentiation of the interpolated primary variable on the boundary. The present paper avoids this tangential differentiation. Instead, a “pure” BEM method is proposed for solving this class of problems. Dirichlet, Neumann and mixed problems are addressed in this paper, and some numerical examples are included in it.     

Thermoelastic fracture mechanics with regularized hypersingular boundary integral equations

This paper is concerned with thermoelastic fracture mechanics, in three-dimensional linear elasticity, using hypersingular boundary integral equations (HBIEs). The HBIEs are regularized by employing modes of deformation (or “simple solutions”). In addition to rigid-body and linear displacement modes, which have been used before for isothermal linear elastic fracture mechanics (LEFM), a new mode is employed here. This new mode is a thermal one in which the body is elevated to a uniform temperature but fully constrained (i.e., zero prescribed displacements) on its bounding surface. An existing isothermal LEFM computer code called BES is extended in this work to include thermoelastic terms. Some numerical results are presented from this new code.     

Traffic flow continuum modeling by hypersingular boundary integral equations

The quantity of data necessary in order to study traffic in dense urban areas through a traffic network, and the large volume of information that is provided as a result, causes managerial difficulties for the said model. A study of this kind is expensive and complex, with many sources of error connected to each step carried out. A simplification like the continuous medium is a reasonable approximation and, for certain dimensions of the actual problem, may be an alternative to be kept in mind. The hypotheses of the continuous model introduce errors comparable to those associated with geometric inaccuracies in the transport network, with the grouping of hundreds of streets in one same type of link and therefore having the same functional characteristics, with the centralization of all journey departure points and destinations in discrete centroids and with the uncertainty produced by a huge origin/destination matrix that is quickly phased out, etc. In the course of this work, a new model for characterizing traffic in dense network cities as a continuous medium, the diffusion–advection model, is put forward. The model is approached by means of the boundary element method, which has the fundamental characteristic of only requiring the contour of the problem to be discretized, thereby reducing the complexity and need for information into one order versus other more widespread methods, such as finite differences and the finite element method. On the other hand, the boundary elements method tends to give a more complex mathematical formulation. In order to validate the proposed technique, three examples in their fullest form are resolved with a known analytic solution. Copyright © 2009 John Wiley & Sons, Ltd.     

Evaluation of gradients on the boundary using fully regularized hypersingular boundary integral equations

Summary The evaluation of the gradient of the primary variable on the boundary for the Laplace problem, and the stress for the elasticity problem, involves hypersingular boundary integrals (HBIEs). To obtain any meaningful results from these integrals, an appropriate regularization scheme needs to be developed. We present an elegant way of calculation of gradients on the boundary of a body, starting from HBIEs regularized by using simple solutions or modes. Our method is currently limited to the calculation of gradients at regular points on the boundary at which the gradients of the primary variables are continous. Theiterative scheme developed in this paper is shown to be extremely robust for the calculations of the gradients. The method is tested on two Laplace problems and two problems in linear elasticity. This method does not involve any limiting process and can be easily extended to 3-dimensions. The approach developed in this paper can also be extended to other problems like acoustics and elastodynamics.     

Regularization of hypersingular boundary integral equations: a new approach for axisymmetric elasticity

A hypersingular boundary integral equation (HBIE) formulation, for axisymmetric linear elasticity, has been recently presented by de Lacerda and Wrobel [Int. J. Numer. Meth. Engng 52 (2001) 1337]. The strongly singular and hypersingular equations in this formulation are regularized by de Lacerda and Wrobel by employing the singularity subtraction technique. The present paper revisits the same problem. The axisymmetric HBIE formulation for linear elasticity is interpreted here in a ‘finite part’ sense and is then regularized by employing a ‘complete exclusion zone’. The resulting regularized equations are shown to be simpler than those by de Lacerda and Wrobel.     

A fast spectral Galerkin method for hypersingular boundary integral equations in potential theory

This research is focused on the development of a fast spectral method to accelerate the solution of three-dimensional hypersingular boundary integral equations of potential theory. Based on a Galerkin approximation, the fast Fourier transform and local interpolation operators, the proposed method is a generalization of the precorrected-FFT technique to deal with double-layer potential kernels, hypersingular kernels and higher-order basis functions. Numerical examples utilizing piecewise linear shape functions are included to illustrate the performance of the method. The US Government retains a nonexclusive royalty-free license to publish or reproduce the published form of this contribution, or allow others to do so, for US Government purposes.     

Volume integration in the hypersingular boundary integral equation

The evaluation of volume integrals that arise in conjunction with a hypersingular boundary integral formulation is considered. In a recent work for the standard (singular) boundary integral equation, the volume term was decomposed into an easily computed boundary integral, plus a remainder volume integral with a modified source function. The key feature of this modified function is that it is everywhere zero on the boundary. In this work it is shown that the same basic approach is successful for the hypersingular equation, despite the stronger singularity in the domain integral. Specifically, the volume term can be directly evaluated without a body-fitted volume mesh, by means of a regular grid of cells that cover the domain. Cells that intersect the boundary are treated by continuously extending the integrand to be zero outside the domain. The method and error results for test problems are presented in terms of the three-dimensional Poisson problem, but the techniques are expected to be generally applicable.     

Explicit evaluation of hypersingular boundary integral equations for acoustic sensitivity analysis based on direct differentiation method

This paper presents a new set of boundary integral equations for three dimensional acoustic shape sensitivity analysis based on the direct differentiation method. A linear combination of the derived equations is used to avoid the fictitious eigenfrequency problem associated with the conventional boundary integral equation method when solving exterior acoustic problems. The strongly singular and hypersingular boundary integrals contained in the equations are evaluated as the Cauchy principal values and Hadamard finite parts for constant element discretization without using any regularization technique in this study. The present boundary integral equations are more efficient to use than the usual ones based on any other singularity subtraction technique and can be applied to the fast multipole boundary element method more readily and efficiently. The effectiveness and accuracy of the present equations are demonstrated through some numerical examples.     

An efficient method to evaluate hypersingular and supersingular integrals in boundary integral equations analysis

An efficient algorithm is employed to evaluated hyper and super singular integral equations encountered in boundary integral equations analysis of engineering problems. The algorithm is based on multiple subtractions and additions to separate singular and regular integral terms in the polar transformation domain, primarily established in Refs. (Guiggiani M, Krishnasamy G, Rudolphi TJ, Rizzo FJ. A general algorithm for the numerical solution of hypersingular boundary integral equations. Trans ASME 1992;59:604–614; Guiggiani M, Casalini P. Direct computation of Cauchy principal value integral in advanced boundary element. Int J Numer Meth Engng 1987;24:1711–1720. Guiggiani M, Gigante A. A general algorithm for multidimensional Cauchy principal value integrals in the boundary element method. J Appl Mech Trans ASME 1990;57:906–915). It can be proved that the regular terms have finite analytical solutions in the range of integration, and the singular terms will be replaced by special periodic kernels in the integral equations. The subtractions involve to multiple derivatives of analytical kernels and the additions require some manipulation to separate the remaining regular terms from singular ones. The regular terms are computed numerically. Three examples on numerical evaluation of singular boundary integrals are presented to show the efficiency and accuracy of the algorithm. In this respect, strongly singular and hypersingular integrals of potential flow problems are considered, followed by a supersingular integral which is extracted from the partial differentiation of a hypersingular integral with respect to the source point.     

Nearly singular approximations of CPV integrals with end- and corner-singularities for the numerical solution of hypersingular boundary integral equations

A local numerical approach to cope with the singular and hypersingular boundary integral equations (BIEs) in non-regularized forms is presented in the paper for 2D elastostatics. The approach is based on the fact that the singular boundary integrals can be represented approximately by the mean values of two nearly singular boundary integrals and on the techniques of distance transformations developed primarily in previous work of the authors. The nearly singular approximations in the present work, including the normal and the tangential distance transformations, are designed for the numerical evaluation of boundary integrals with end-singularities at junctures between two elements, especially at corner points where sufficient continuity requirements are met. The approach is very general, which makes it possible to solve the hypersingular BIE numerically in a non-regularized form by using conforming C⁰ quadratic boundary elements and standard Gaussian quadratures without paying special attention to the corner treatment.With the proposed approach, an infinite tension plate with an elliptical hole and a pressurized thick cylinder were analyzed by using both the formulations of conventional displacement and traction boundary element methods, showing encouragingly the efficiency and the reliability of the proposed approach. The behaviors of boundary integrals with end- and corner-singular kernels were observed and compared by the additional numerical tests. It is considered that the weakly singularities remain but should have been cancelled with each other if used in pairs by the corresponding terms in the neighboring elements, where the corner information is included naturally in the approximations.     

On the use of conjugate gradient-type methods for boundary integral equations

This paper deals with a number of iterative methods for solving matrix equations that result from boundary integral equations. The matrices are non-sparse, and in general neither positive definite nor symmetric. Traditional methods like Gauss-Seidel do not give satisfactory results, therefore the use of conjugate gradient- and Krylov-type methods is investigated based on work of Kleinman and Van den Berg, who presented a general framework for these methods. Eleven of these algorithms are given and their performance (without preconditioning) is compared in a test case involving four different integral operators arising in potential theory. For all four matrix equations the Generalized Minimal Residual method (GMRES) outperforms all other iterative methods in both computation time per iteration and total computation time. For the Fredholm equations of the first kind this method also is the fastest with respect to the number of iterations. The Bi-conjugate gradient method (Bi-CG) and the Quasi-minimal residual method (QMR) are the best alternatives. For the Fredholm equations of the second kind more methods can be used efficiently besides GMRES. The Conjugate Gradient Squared method (CGS), the Bi-conjugate gradient method (Bi-CG), its stabilized version (Bi-CGSTAB) and the Quasi Minimal Residual method (QMR) are efficient alternatives.     

On direct numerical treatment of hypersingular integral equations arising in mechanics and acoustics

Summary.  In this paper, we present a treatment of hypersingular integral equations, which have relevant applications in many problems of wave dynamics, elasticity and fluid mechanics with mixed boundary conditions. The main goal of the present work is the development of an efficient direct numerical collocation method. The paper also includes two examples taken from fracture mechanics and acoustics: a single crack in a linear isotropic elastic medium, and diffraction of a plane acoustic wave by a thin rigid screen. Received July 15, 2002; revised November 13, 2002 Published online: May 8, 2003 The paper has been supported by Italian Ministry for Instructions, University and Research (M.I.U.R.) through its national and local projects.     

Complex hypersingular integrals and integral equations in plane elasticity

Summary Complex hypersingular (finite-part) integrals and integral equations are considered in the functional class of N. Muskhelishvili. The appropriate definition is given. Three regularization (equivalence) formulae follow from this definition. They reduce hypersingular integrals to singular ones and allow to derive hypersingular analogues for Sokhotsky-Plemelj's formulae and for conditions that are necessary and sufficient for the function to be piecewise holomorphic. Two approaches to get and investigate complex hypersingular equations follow from these results: one of them is based on the equivalence formulae; as to the other, it is based on above-mentioned conditions. As an example, authors' equation for plane elasticity is studied. The existence of a unique solution is stated and some advantages over singular equations are outlined. To solve hypersingular equations the quadrature rules are presented. The accuracy of different quadrature formulae is compared, the examples being used. They confirm the need to take into account asymptotics and to carry out a thorough analytical investigation to get safe numerical results.     

An improved form of the hypersingular boundary integral equation for exterior acoustic problems

An improved form of the hypersingular boundary integral equation (BIE) for acoustic problems is developed in this paper. One popular method for overcoming non-unique problems that occur at characteristic frequencies is the well-known Burton and Miller (1971) method [7], which consists of a linear combination of the Helmholtz equation and its normal derivative equation. The crucial part in implementing this formulation is dealing with the hypersingular integrals. This paper proposes an improved reformulation of the Burton–Miller method and is used to regularize the hypersingular integrals using a new singularity subtraction technique and properties from the associated Laplace equations. It contains only weakly singular integrals and is directly valid for acoustic problems with arbitrary boundary conditions. This work is expected to lead to considerable progress in subsequent developments of the fast multipole boundary element method (FMBEM) for acoustic problems. Numerical examples of both radiation and scattering problems clearly demonstrate that the improved BIE can provide efficient, accurate, and reliable results for 3-D acoustics.     

Smoothness–relaxation strategies for singular and hypersingular integral equations

Three stages are involved in the formulation of a typical direct boundary element method: derivation of an integral representation; taking a Limit To the Boundary (LTB) so as to obtain an integral equation; and discretization. We examine the second and third stages, focussing on strategies that are intended to permit the relaxation of standard smoothness assumptions. Two such strategies are indicated. The first is the introduction of various apparent or ‘pseudo-LTBs’. The second is ‘relaxed regularization’, in which a regularized integral equation, derived rigorously under certain smoothness assumptions, is used when less smoothness is available. Both strategies are shown to be based on inconsistent reasoning. Nevertheless, reasons are offered for having some confidence in numerical results obtained with certain strategies. Our work is done in two physical contexts, namely two-dimensional potential theory (using functions of a complex variable) and three-dimensional elastostatics. © 1998 John Wiley & Sons, Ltd.     

Analytical regularization of hypersingular integral for Helmholtz equation in boundary element method

This paper presents a gradient field representation using an analytical regularization of a hypersingular boundary integral equation for a two-dimensional time harmonic wave equation called the Helmholtz equation. The regularization is based on cancelation of the hypersingularity by considering properties of hypersingular elements that are adjacent to a singular node. Advantages to this regularization include applicability to evaluate corner nodes, no limitation for element size, and reduced computational cost compared to other methods. To demonstrate capability and accuracy, regularization is estimated for a problem about plane wave propagation. As a result, it is found that even at a corner node the most significant error in the proposed method is due to truncation error of non-singular elements in discretization, and error from hypersingular elements is negligibly small.     

On the implementation of 3D Galerkin boundary integral equations

In this article, a reverse contribution technique is proposed to accelerate the construction of the dense influence matrices associated with a Galerkin approximation of hypersingular boundary integral equations of mixed-type in potential theory. In addition, a general-purpose sparse preconditioner for boundary element methods has also been developed to successfully deal with ill-conditioned linear systems arising from the discretization of mixed boundary-value problems on non-smooth surfaces. The proposed preconditioner, which originates from the precorrected-FFT method, is sparse, easy to generate and apply in a Krylov subspace iterative solution of discretized boundary integral equations. Moreover, an approximate inverse of the preconditioner is implicitly built by employing an incomplete LU factorization. Numerical experiments involving mixed boundary-value problems for the Laplace equation are included to illustrate the performance and validity of the proposed techniques.     

Differentiability of strongly singular and hypersingular boundary integral formulations with respect to boundary perturbations

In this paper, we establish that the Lagrangian-type material differentiation formulas, that allow to express the first-order derivative of a (regular) surface integral with respect to a geometrical domain perturbation, still hold true for the strongly singular and hypersingular surface integrals usually encountered in boundary integral formulations. As a consequence, this work supports previous investigations where shape sensitivities are computed using the so-called direct differentiation approach in connection with singular boundary integral equation formulations. Communicated by T. Cruse, 6 September 1996