Cell‐based volume integration for boundary integral analysis

The evaluation of volume integrals that arise in boundary integral formulations for non‐homogeneous problems was considered. Using the “Galerkin vector” to represent the Green's function, the volume integral was decomposed into a boundary integral, together with a volume integral wherein the source function was everywhere zero on the boundary. This new volume integral can be evaluated using a regular grid of cells covering the domain, with all cell integrals, including partial cells at the boundary, evaluated by simple linear interpolation of vertex values. For grid vertices that lie close to the boundary, the near‐singular integrals were handled by partial analytic integration. The method employed a Galerkin approximation and was presented in terms of the three‐dimensional Poisson problem. An axisymmetric formulation was also presented, and in this setting, the solution of a nonlinear problem was considered. Copyright © 2012 John Wiley & Sons, Ltd.     

Galerkin formulation and singularity subtraction for spectral solutions of boundary integral equations

A new spectral Galerkin formulation is presented for the solution of boundary integral equations. The formulation is carried out with an exact singularity subtraction procedure based on analytical integrations, which provides a fast and precise way to evaluate the coefficient matrices. The new Galerkin formulation is based on the exact geometry of the problem boundaries and leads to a non-element method that is completely free of mesh generation. The numerical behaviour of the method is very similar to the collocation method; for Dirichlet problems, however, it leads to a symmetric coefficient matrix and therefore requires half the solution time of the collocation method. © 1998 John Wiley & Sons, Ltd.     

Boundary integral equation for tangential derivative of flux in Laplace and Helmholtz equations

In this paper, the boundary integral equations (BIEs) for the tangential derivative of flux in Laplace and Helmholtz equations are presented. These integral representations can be used in order to solve several problems in the boundary element method (BEM): cubic solutions including degrees of freedom in flux's tangential derivative value (Hermitian interpolation), nodal sensitivity, analytic gradients in optimization problems, or tangential derivative evaluation in problems that require the computation of such variable (elasticity problems in BEM). The analysis has been developed for 2D formulation. Kernels for tangential derivative of flux lead to high‐order singularities (O(1/r³)). The limit to the boundary analysis has been carried out. Based on this analysis, regularization formulae have been obtained in order to use such BIE in numerical codes. A set of numerical benchmarks have been carried out in order to validate theoretical and practical aspects, by considering known analytic solutions for the test problems. The results show that the tangential BIEs have been properly developed and implemented. Copyright © 2005 John Wiley & Sons, Ltd.     

The computation of the J integral with the traction boundary integral equation

The solutions of the displacement boundary integral equation (BIE) are not uniquely determined in certain types of boundary conditions. Traction boundary integral equations that have unique solutions in these traction and mixed boundary cases are established. For two‐dimensional linear elasticity problems, the divergence‐free property of the traction boundary integral equation is established. By applying Stokes' theorem, unknown tractions or displacements can be reduced to computation of traction integral potential functions at the boundary points. The same is true of the J integral: it is divergence‐free and the evaluation of the J integral can be inverted into the computation of the J integral potential functions at the boundary points of the cracked body. The J integral can be expressed as the linear combination of the tractions and displacements from the traction BIE on the boundary of the cracked body. Numerical integrals are not needed at all. Selected examples are presented to demonstrate the validity of the traction boundary integral and J integral. Copyright © 2004 John Wiley & Sons, Ltd.     

A Galerkin boundary integral method for multiple circular elastic inclusions

The problem of an infinite, isotropic elastic plane containing an arbitrary number of circular elastic inclusions is considered. The analysis procedure is based on the use of a complex singular integral equation. The unknown tractions at each circular boundary are approximated by a truncated complex Fourier series. A system of linear algebraic equations is obtained by using the classical Galerkin method and the Gauss–Seidel algorithm is used to solve the system. Several numerical examples are considered to demonstrate the effectiveness of the approach. Copyright © 2001 John Wiley & Sons, Ltd.     

Symmetric Galerkin boundary integral formulation for interface and multi-zone problems

Domains containing an ‘internal boundary’, such as a bi-material interface, arise in many applications, e.g. composite materials and geophysical simulations. This paper presents a symmetric Galerkin boundary integral method for this important class of problems. In this situation, the physical quantities are known to satisfy continuity conditions across the interface, but no boundary conditions are specified. The algorithm described herein achieves a symmetric matrix of reduced size. Moreover, the symmetry can also be invoked to lessen the numerical work involved in constructing the system of equations, and thus the method is computationally very efficient. A prototype numerical example, with several variations in the boundary conditions and material properties, is employed to validate the formulation and corresponding numerical procedure. The boundary element results are compared with analytical solutions and with numerical results obtained with the finite element method. © 1997 John Wiley & Sons, Ltd.     

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.     

Transient heat conduction analysis in functionally graded materials by the meshless local boundary integral equation method

Advanced computational method for transient heat conduction analysis in continuously nonhomogeneous functionally graded materials (FGM) is proposed. The method is based on the local boundary integral equations with moving least square approximation of the temperature and heat flux. The initial-boundary value problem is solved by the Laplace transform technique. Both Papoulis and Stehfest algorithms are applied for the numerical Laplace inversion to obtain the time-dependent solutions. Numerical results are presented for a finite strip and a hollow cylinder with an exponential spatial variation of material parameters.     

An ancillary boundary integral equation for magnetostatic analysis

The boundary element method (BEM) has been established as an effective means for magnetostatic analysis. Direct BEM formulations for the magnetic vector potential have been developed over the past 20 years. There is a less well-known direct boundary integral equation (BIE) for the magnetic flux density which can be derived by taking the curl of the BIE for the magnetic vector potential and applying properties of the scalar triple product. On first inspection, the ancillary boundary integral equation for the magnetic flux density appears to be homogeneous, but it can be shown that the equation is well-posed and non-homogeneous using appropriate boundary conditions. In the current research, the use of the ancillary boundary integral equation for the magnetic flux density is investigated as a stand-alone equation and in tandem with the direct formulation for the magnetic vector potential.     

On solvability of a boundary integral equation of the first kind for Dirichlet boundary value problems in plane elasticity

The solution of a Dirichlet boundary value problem of plane isotropic elasticity by the boundary integral equation (BIE) of the first kind obtained from the Somigliana identity is considered. The logarithmic function appearing in the integral kernel leads to the possibility of this operator being non-invertible, the solution of the BIE either being non-unique or not existing. Such a situation occurs if the size of the boundary coincides with the so-called critical (or degenerate) scale for a certain form of the fundamental solution used. Techniques for the evaluation of these critical scales and for the removal of the non-uniqueness appearing in the problems with critical scales solved by the BIE of the first kind are proposed and analysed, and some recommendations for BEM code programmers based on the analysis presented are given.     

On the Cauchy principal value of the surface integral in the boundary integral equation of 3D elasticity

A simple demonstration of the existence of the Cauchy principal value (CPV) of the strongly singular surface integral in the Somigliana Identity at a non-smooth boundary point is presented. First a regularization of the strongly singular integral by analytical integration of the singular term in the radial direction in pre-image planes of smooth surface patches is carried out. Then it is shown that the sum of the angular integrals of the characteristic of the tractions of the Kelvin fundamental solution is zero, a formula for the transformation of angles between the tangent plane of a suface patch and the pre-image plane at smooth mapping of the surface patch being derived for this purpose.     

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.     

A fast and accurate algorithm for a Galerkin boundary integral method

A fast and accurate algorithm is presented to increase the computational efficiency of a Galerkin boundary integral method for solving two-dimensional elastostatics problems involving numerous straight cracks and circular inhomogeneities. The efficiency is improved by computing the combined influences of groups, or blocks, of elements—with each element being an inclusion, a hole, or a crack—using asymptotic expansions, multiple shifts, and Taylor series expansions. The coefficients in the asymptotic and Taylor series expansions are computed analytically. Implementation of this algorithm involves a single- or multi-level grid, a clustering technique, and a tree data structure. An iterative procedure is adopted to solve the coefficients in the series expansions of boundary unknowns block by block. The elastic fields in each block are calculated by superposition of the direct influences from the nearby elements and the grouped far-field influences from all the other elements. This fast multipole algorithm is considerably more efficient for large-scale practical problems than the conventional approach.     

Local BIEM for transient heat conduction analysis in 3-D axisymmetric functionally graded solids

An advanced computational method for transient heat conduction analysis in 3-D axisymmetric continuously nonhomogeneous functionally graded materials (FGM) is proposed. The analysed domain is covered by small circular subdomains. On each subdomain local boundary integral equations for the transient heat conduction problem are derived in the Laplace transform domain. The meshless approximation based on the moving least-squares method is employed for the numerical implementation. The Stehfest algorithm is applied for the numerical Laplace inversion to obtain the temporal variation. Numerical results are presented for finite full and hollow cylinders with an exponential variation of material parameters with spatial coordinates. The authors acknowledge the support by the Slovak Science and Technology Assistance Agency registered under number APVT-51-003702, and the Project for Bilateral Cooperation in Science and Technology supported jointly by the International Bureau of the German BMBF and the Ministry of Education of Slovak Republic under the project number SVK 01/020.     

A new meshless regular local boundary integral equation (MRLBIE) approach

A new meshless method based on a regular local integral equation and the moving least‐squares approximation is developed. The present method is a truly meshless one as it does not need a ‘finite element or boundary element mesh’, either for purposes of interpolation of the solution variables, or for the integration of the ‘energy’. All integrals can be easily evaluated over regularly shaped domains (in general, spheres in three‐dimensional problems) and their boundaries. No derivatives of the shape functions are needed in constructing the system stiffness matrix for the internal nodes, as well as for those boundary nodes with no essential‐boundary‐condition‐prescribed sections on their local boundaries. Numerical examples presented in the paper show that high rates of convergence with mesh refinement are achievable, and the computational results for the unknown variable and its derivatives are very accurate. No special post‐processing procedure is required to compute the derivatives of the unknown variable, as the original result, from the moving least‐squares approximation, is smooth enough. Copyright © 1999 John Wiley & Sons, Ltd.     

Modification of the boundary integral equation for the Navier–Lamé equations in 2D elastoplastic problems and its numerical solving

The paper presents the generalization of the modification of classical boundary integral equation and obtaining parametric integral equation system for 2D elastoplastic problems. The modification was made to obtain such equations for which numerical solving does not require application of finite or boundary elements. This was achieved through the use of curves and surfaces for modeling introduced at the stage of analytical modification of the classic boundary integral equation. For approximation of plastic strains the Lagrange polynomials with various number and arrangement of interpolation nodes were used. Reliability of the modification was verified on examples with analytical solutions. Copyright © 2016 John Wiley & Sons, Ltd.     

Meshless analysis of two-dimensional Stokes flows with the Galerkin boundary node method

In this paper, the Galerkin boundary node method (GBNM) is developed for the solution of stationary Stokes problems in two dimensions. The GBNM is a boundary only meshless method that combines a variational form of boundary integral formulations for governing equations with the moving least-squares (MLS) approximations for construction of the trial and test functions. Boundary conditions in this approach are included into the variational form, thus they can be applied directly and easily despite the MLS shape functions lack the property of a delta function. Besides, the GBNM keeps the symmetry and positive definiteness of the variational problems. Convergence analysis results of both the velocity and the pressure are given. Some selected numerical tests are also presented to demonstrate the efficiency of the method.     

A general regularization of the hypersingular integrals in the symmetric Galerkin boundary element method

The symmetric Galerkin boundary element method is used to solve boundary value problems by keeping the symmetric nature of the matrix obtained after discretization. The matrix elements are obtained from a double integral involving the double derivative of Green's operator, which is highly singular. The paper presents a regularization of the hypersingular integrals which depend only on the properties of Green's tensor. The method is presented in the case of Laplace's operator, with an example of application. The case of elasticity is finally addressed theoretically, showing an easy extension to any case of anisotropy. Copyright © 2009 John Wiley & Sons, Ltd.     

Iterative solution of Hermite boundary integral equations

An efficient iterative method for solution of the linear equations arising from a Hermite boundary integral approximation has been developed. Along with equations for the boundary unknowns, the Hermite system incorporates equations for the first‐order surface derivatives (gradient) of the potential, and is therefore substantially larger than the matrix for a corresponding linear approximation. However, by exploiting the structure of the Hermite matrix, a two‐level iterative algorithm has been shown to provide a very efficient solution algorithm. In this approach, the boundary function unknowns are treated separately from the gradient, taking advantage of the sparsity and near‐positive definiteness of the gradient equations. In test problems, the new algorithm significantly reduced computation time compared with iterative solution applied to the full matrix. This approach should prove to be even more effective for the larger systems encountered in three‐dimensional analysis, and increased efficiency should come from pre‐conditioning of the non‐sparse matrix component. Published in 2007 by John Wiley & Sons, Ltd.     

Non‐singular boundary integral formulations for plane interior potential problems

In this article, a non‐singular formulation of the boundary integral equation is developed to solve smooth and non‐smooth interior potential problems in two dimensions. The subtracting and adding‐back technique is used to regularize the singularity of Green's function and to simplify the calculation of the normal derivative of Green's function. After that, a global numerical integration is directly applied at the boundary, and those integration points are also taken as collocation points to simplify the algorithm of computation. The result indicates that this simple method gives the convergence speed of order N ^?3 in the smooth boundary cases for both Dirichlet and mix‐type problems. For the non‐smooth cases, the convergence speed drops at O(N ^?1/2) for the Dirichlet problems. Copyright © 2001 John Wiley & Sons, Ltd.