The boundary element-free method for 2D interior and exterior Helmholtz problems

The boundary element-free method (BEFM) is developed in this paper for numerical solutions of 2D interior and exterior Helmholtz problems with mixed boundary conditions of Dirichlet and Neumann types. A unified boundary integral equation is established for both interior and exterior problems. By using the improved interpolating moving least squares method to form meshless shape functions, mixed boundary conditions in the BEFM can be satisfied directly and easily. Detailed computational formulas are derived to compute weakly and strongly singular integrals over linear and higher order integration cells. Three numerical integration procedures are developed for the computation of strongly singular integrals. Numerical examples involving acoustic scattering and radiation problems are presented to show the accuracy and efficiency of the meshless method.     

A semi-analytical algorithm for the evaluation of the nearly singular integrals in three-dimensional boundary element methods

The nearly singular integrals occur in the boundary integral equations when the source point is close to an integration element (as compared to its size) but not on the element. In this paper, the concept of a relative distance from a source point to the boundary element is introduced to describe possible influence of the singularity of the integrals. Then a semi-analytical algorithm is proposed for evaluating the nearly strongly singular and hypersingular integrals in the three-dimensional BEM. By using integration by parts, the nearly singular surface integrals on the elements are transformed to a series of line integrals along the contour of the element. The singular behavior, which appears as factor, is separated from remaining regular integrals. Consequently standard numerical quadrature can provide very accurate evaluation of the resulting line integrals. The semi-analytical algorithm is applied to analyzing the three-dimensional elasticity problems, such as very thin-walled structures. Meanwhile, the displacements and stresses at the interior points very close to its bounding surface are also determined efficiently. The results of the numerical investigation demonstrate the accuracy and effectiveness of the algorithm.     

A node-centric indirect boundary element method: three-dimensional displacement discontinuities

An indirect boundary element formulation based on unknown physical values, defined only at the nodes (vertices) of a boundary discretization of a linear elastic continuum, is introduced. As an adaptation of this general framework, a linear displacement discontinuity density distribution using a flat triangular boundary discretization is considered. A unified element integration methodology based on the continuation principle is introduced to handle regular as well as near-singular and singular integrals. The boundary functions that form the basis of the integration methodology are derived and tabulated in the appendix for linear displacement discontinuity densities. The integration of the boundary functions is performed numerically using an adaptive algorithm which ensures a specified numerical accuracy. The applications include verification examples which have closed-form analytical solutions as well as practical problems arising in rock engineering. The node-centric displacement discontinuity method is shown to be numerically efficient and robust for such problems.     

Numerical evaluation of singular and finite-part integrals on curved surfaces using symbolic manipulation

The numerical integration of all singular surface integrals arising in 3-d boundary element methods is analyzed theoretically and computationally. For all weakly singular integrals arising in BEM, Duffy's triangular or local polar coordinates in conjunction with tensor product Gaussian quadrature are efficient and reliable for bothh-andp-boundary elements. Cauchy- and hypersingular surface integrals are reduced to weakly singular ones by analytic regularization which is done automatically by symbolic manipulation.     

Bending of anisotropic plates by charge simulation method

A boundary element method, called the charge simulation method, is presented for analysis of anisotropic thin-plate bending problems. In this method the singular integrals involved in the other boundary element methods are eliminated and there is no numerical integration involved. Further, the domain integral is replaced by a polynomial particular integral; hence the domain discretization is avoided. This method is conceptually very simple. The results obtained by this method are compared with the available analytical solutions for various anisotropic and symmetric laminates and the results are in good agreement.     

Analytical integral algorithm applied to boundary layer effect and thin body effect in BEM for anisotropic potential problems

A new completely analytical integral algorithm is proposed and applied to the evaluation of nearly singular integrals in boundary element method (BEM) for two-dimensional anisotropic potential problems. The boundary layer effect and thin body effect are dealt with. The completely analytical integral formulas are suitable for the linear and non-isoparametric quadratic elements. The present algorithm applies the analytical formulas to treat nearly singular integrals. The potentials and fluxes at the interior points very close to boundary are evaluated. The unknown potentials and fluxes at boundary nodes for thin body problems with the thickness-to-length ratios from 1E−1 to 1E−8 are accurately calculated by the present algorithm. Numerical examples on heat conduction demonstrate that the present algorithm can effectively handle nearly singular integrals occurring in boundary layer effect and thin body effect in BEM. Furthermore, the present linear BEM is especially accurate and efficient for the numerical analysis of thin body problems.     

Evaluation of singular integrals in the symmetric Galerkin boundary element method

The implementation of the symmetric Galerkin boundary element method (SGBEM) involves extensive work on the evaluation of various integrals, ranging from regular integrals to hypersingular integrals. In this paper, the treatments of weak singular integrals in the time domain are reviewed, and analytical evaluations for the spatial double integrals which contain weak singular terms are derived. A special scheme on the allocation of Gaussian integration points for regular double integrals in the SGBEM is developed to improve the efficiency of the Gauss–Legendre rule. The proposed approach is implemented for the two-dimensional elastodynamic problems, and two numerical examples are presented to verify the accuracy of the numerical implementation.     

Quadrature for parabolic Galerkin BEM with moving surfaces

The successful implementation of the Galerkin Boundary Element Method hinges on the accurate and effective quadrature of the influence coefficients. For parabolic boundary integral operators quadrature must be performed in space and time where integrals have singularities when source- and evaluation points coincide. For problems where the surface is fixed, the time integration can be performed analytically, but for moving geometries numerical quadrature in space and time must be used. For this case a set of transformations is derived that render the singular space–time integrals into smooth integrals that can be treated with standard tensor product Gauss quadrature rules. This methodology can be applied to the heat equation and to transient Stokes flow.     

Dynamic inelastic structural analysis by boundary element methods

Summary Boundary element methodologies for the determination of the response of inelastic two-and three-dimensional solids and structures as well as beams and flexural plates to dynamic loads are briefly presented and critically discussed. Elastoplastic and viscoplastic material behaviour in the framework of small deformation theories are considered. These methodologies can be separated into four main categories: those which employ the elastodynamic fundamental solution in their formulation, those which employ the elastostatic fundamental solution in their formulation, those which combine boundary and finite elements for the creation of an efficient hybrid scheme and those representing special boundary element techniques. The first category, in addition to the boundary discretization, requires a discretization of those parts of the interior domain expected to become inelastic, while the second category a discretization of the whole interior domain, unless the inertial domain integrals are transformed by the dual reciprocity technique into boundary ones, in which case only the inelastic parts of the domain have to be discretized. The third category employs finite elements for one part of the structure and boundary elements for its remaining part in an effort to combine the advantages of both methods. Finally, the fourth category includes special boundary element techniques for inelastic beams and plates and symmetric boundary element formulations. The discretized equations of motion in all the above methodologies are solved by efficient step-by-step time integration algorithms. Numerical examples involving two-and three-dimensional solids and structures and flexural plates are presented to illustrate all these methodologies and demonstrate their advantages. Finally, directions for future research in the area are suggested.     

An Implementation of Fast Wavelet Galerkin Methods for Integral Equations of the Second Kind

We present a numerical implementation of the fast Galerkin method for Fredholm integral equations of the second kind using the piecewise polynomial wavelets. We focus on addressing critical issues for the numerical implementation of such a method. They include a choice of practical truncation strategy, numerical integration of weakly singular integrals and the error control of the numerical quadrature. We also implement a multiscale iteration method for solving the resulting compressed linear system. Numerical examples are given to demonstrate the proposed ideas and methods.     

On the numerical integration of certain singular integrals

In this paper we deal with the numerical calculation of singular integrals in the sense of Hadamard by simple integration methods. Error estimates and numerical examples are given. The optimality of the quadrature rules will be examined.     

Numerical quadratures for singular and hypersingular integrals

We present a procedure for the design of high-order quadrature rules for the numerical evaluation of singular and hypersingular integrals; such integrals are frequently encountered in solution of integral equations of potential theory in two dimensions. Unlike integrals of both smooth and weakly singular functions, hypersingular integrals are pseudo-differential operators, being limits of certain integrals; as a result, standard quadrature formulae fail for hypersingular integrals. On the other hand, such expressions are often encountered in mathematical physics (see, for example, [1]), and it is desirable to have simple and efficient “quadrature” formulae for them. The algorithm we present constructs high-order “quadratures” for the evaluation of hypersingular integrals. The additional advantage of the scheme is the fact that each of the quadratures it produces can be used simultaneously for the efficient evaluation of hypersingular integrals, Hilbert transforms, and integrals involving both smooth and logarithmically singular functions; this results in significantly simplified implementations. The performance of the procedure is illustrated with several numerical examples.     

Quadrature methods for periodic singular and weakly singular Fredholm integral equations

High-accuracy numerical quadrature methods for integrals of singular periodic functions are proposed. These methods are based on the appropriate Euler-Maclaurin expansions of trapezoidal rule approximations and their extrapolations. They are subsequently used to obtain accurate quadrature methods for the solution of singular and weakly singular Fredholm integral equations. Throughout the development the periodic nature of the problem plays a crucial role. Such periodic equations are used in the solution of planar elliptic boundary value problems such as those that arise in elasticity, potential theory, conformal mapping, free surface flows, etc. The use of the quadrature methods is demonstrated with numerical examples.     

Nodal Spline Integration Rules for Certain 2-D Cauchy Principal Value Integrals

In this paper we consider integration rules based on products of one-dimensional nodal spline quadratures for the numerical evaluation of 2-D singular integrals defined in the Hadamard finite part sense. We derive convergence results and error bounds. Some numerical tests are also presented.     

On the selection of collocation points for the numerical solution of singular integral equations with generalized kernels appearing in elasticity problems

A modification of the collocation method for the numerical solution of Cauchy type singular integral equations with generalized kernels is proposed. In accordance with this modification, although the abscissas and weights used in the numerical integration rule for the approximation of the integrals of the integral equation remain unaltered, yet the collocation points are selected in such a way that the poles of the integrands due not only to the Cauchy principal value part of the kernel, but also to the singularities of the generalized part of the kernel are taken into account. This modification assures the convergence of the method to the correct results since the error terms, usually neglected for the collocation points nearest to the end-points of the integration interval and generally tending to infinity, are now taken into consideration for the selection of the collocation points. The method was applied to the singular integral equations derived for the antiplane and plane elasticity problems of a crack terminating at a bimaterial interface.     

Computational Aspects Of Numerical Integration Based On Optimal Nodal Splines

The numerical construction of nodal spline integration rules is considered for the evaluation both of weakly singular integrals and of certain integrals defined in the Hadamard finite part sense. A bound for the quadrature sum condition number is also given.     

A Sequential Discontinuous Galerkin Method for the Coupling of Flow and Geomechanics

A decoupled and sequential numerical method is proposed and analyzed for solving the linear poroelasticity equations. Unlike other splitting approaches, this method is not iterative, which results in a speed-up of the computational time. The interior penalty discontinuous Galerkin method is employed for the spatial discretization and is combined with the backward Euler method for the time discretization. We provide a convergence analysis of the scheme along with numerical results that confirm the theoretical results.     

Limit Values of Derivatives of the Cauchy Integrals and Computation of the Logarithmic Potentials

We derive limit values of high-order derivatives of the Cauchy integrals, which are extensions of the Plemelj-Sokhotskyi formula. We then use them to develop the Taylor expansion of the logarithmic potentials at the normal direction. Based on the Taylor expansion and numerical integration methods for weekly singular functions using grid points, we design fast algorithms for computing the logarithmic potentials. We prove that these methods have an optimal order of convergence with a linear computational complexity. Numerical examples are included to confirm the theoretical estimates for the methods.     

Shape sensitivity analysis with design-dependent loadings—equivalence between continuum and discrete derivatives

The purpose of this paper is twofold: (1) showing equivalence between continuum and discrete formulations in sensitivity analysis when a linear velocity field is used and (2) presenting shape sensitivity formulations for design-dependent loadings. The equations for structural analysis are often composed of the stiffness part and the applied loading part. The shape sensitivity formulations for the stiffness part were well-developed in the literature, but not for the loading part, especially for body forces and surface tractions. The applied loads are often assumed to be conservative or design-independent. In shape design problems, however, the applied loads are often functions of design variables. In this paper, shape sensitivity formulations are presented when the body forces and surface tractions depend on shape design variables. Especially, the continuum–discrete (C–D) and discrete–discrete (D–D) approaches are compared in detail. It is shown that the two methods are theoretically and numerically equivalent when the same discretization, numerical integration, and linear design velocity fields are used. The accuracy of sensitivity calculation is demonstrated using a cantilevered beam under uniform pressure and an arch dam crown cantilever under gravity and hydrostatic loading at the upstream face of the structure. It is shown that the sensitivity results are consistent with finite difference results, but different from the analytical sensitivity due to discretization and approximation errors of numerical analysis.     

Tailored Finite Point Method for a Singular Perturbation Problem with Variable Coefficients in Two Dimensions

In this paper, we propose a tailored-finite-point method for a type of linear singular perturbation problem in two dimensions. Our finite point method has been tailored to some particular properties of the problem. Therefore, our new method can achieve very high accuracy with very coarse mesh even for very small ε, i.e. the boundary layers and interior layers do not need to be resolved numerically. In our numerical implementation, we study the classification of all the singular points for the corresponding degenerate first order linear dynamic system. We also study some cases with nonlinear coefficients. Our tailored finite point method is very efficient in both linear and nonlinear coefficients cases.