Modification of the Discrete Singularity Method for Nonuniform Grids As Applied to One-Dimensional Integral Equations with a Strong Singularity in the Kernel

Differential Equations - Formulas are constructed for the approximate calculation of one-dimensional singular and hypersingular integrals over an interval with an arbitrary partition. On the basis...     

Interpolation-type quadrature formulas for hypersingular integrals in the interval of integration

A hypersingular integral on the interval of integration with a weight function is considered. Spectral ratios for hypersingular integrals on the interval [?1, 1] are proved. Interpolationtype quadrature formulas for certain integrals with weight functions are constructed. The error is estimated.     

Semi-analytic treatment of the three-dimensional Poisson equation via a Galerkin BIE method

A systematic treatment of the three-dimensional Poisson equation via singular and hypersingular boundary integral equation techniques is investigated in the context of a Galerkin approximation. Developed to conveniently deal with domain integrals without a volume-fitted mesh, the proposed method initially converts domain integrals featuring the Newton potential and its gradient into equivalent surface integrals. Then, the resulting boundary integrals are evaluated by means of well-established cubature methods. In this transformation, weakly-singular domain integrals, defined over simply- or multiply-connected domains with Lipschitz boundaries, are rigorously converted into weakly-singular surface integrals. Combined with the semi-analytic integration approach developed for potential problems to accurately calculate singular and hypersingular Galerkin surface integrals, this technique can be employed to effectively deal with mixed boundary-value problems without the need to partition the underlying domain into volume cells. Sample problems are included to validate the proposed approach.     

Semi-analytic treatment of the three-dimensional Poisson equation via a Galerkin BIE method

A systematic treatment of the three-dimensional Poisson equation via singular and hypersingular boundary integral equation techniques is investigated in the context of a Galerkin approximation. Developed to conveniently deal with domain integrals without a volume-fitted mesh, the proposed method initially converts domain integrals featuring the Newton potential and its gradient into equivalent surface integrals. Then, the resulting boundary integrals are evaluated by means of well-established cubature methods. In this transformation, weakly-singular domain integrals, defined over simply- or multiply-connected domains with Lipschitz boundaries, are rigorously converted into weakly-singular surface integrals. Combined with the semi-analytic integration approach developed for potential problems to accurately calculate singular and hypersingular Galerkin surface integrals, this technique can be employed to effectively deal with mixed boundary-value problems without the need to partition the underlying domain into volume cells. Sample problems are included to validate the proposed approach.     

Transformation of hypersingular integrals and black-box cubature

In this paper, we will consider hypersingular integrals as they arise by transforming elliptic boundary value problems into boundary integral equations. First, local representations of these integrals will be derived. These representations contain so-called finite-part integrals. In the second step, these integrals are reformulated as improper integrals. We will show that these integrals can be treated by cubature methods for weakly singular integrals as they exist in the literature.

     

An integral equation of the first kind for a free boundary value problem of the stationary Stokes' equations

We investigate a free boundary value problem of the stationary Stokes' equations. In a previous paper adapted hydrodynamical potentials have been constructed and their jump relations have been discussed. Here we study a direct method to obtain an equivalent boundary integral equations' system of the first kind. Its solution properties are investigated in the framework of strongly elliptic pseudodifferential operators. For numerical purposes a suitable representation formula for the variational equation is given in terms of integro-differential operators which avoids the evaluation of hypersingular integrals.     

Integration of singular Galerkin-type boundary element integrals for 3D elasticity problems

Summary. A Galerkin approximation of both strongly and hypersingular boundary integral equation (BIE) is considered for the solution of a mixed boundary value problem in 3D elasticity leading to a symmetric system of linear equations. The evaluation of Cauchy principal values (v. p.) and finite parts (p. f.) of double integrals is one of the most difficult parts within the implementation of such boundary element methods (BEMs). A new integration method, which is strictly derived for the cases of coincident elements as well as edge-adjacent and vertex-adjacent elements, leads to explicitly given regular integrand functions which can be integrated by the standard Gauss-Legendre and Gauss-Jacobi quadrature rules. Problems of a wide range of integral kernels on curved surfaces can be treated by this integration method. We give estimates of the quadrature errors of the singular four-dimensional integrals. Received June 25, 1995 / Revised version received January 29, 1996     

边界元方法的一些研究进展

本文旨在综述我们小组近二十年来在边界元方法这一领域的一些研究成果,在简要介绍边界元方法的基本思想后,主要介绍了一类非线性界面问题的有限元-边界元耦合方法、求解电磁散射问题的有限元-边界元耦合方法和超奇异积分的一类计算方法.     

THE APPROXIMATE COMPUTATION OF HYPERSINGULAR INTEGRALS ON INTERVAL

In this paper some numerical methods for computing hypersingular integrals on interval are given. Using geometric meshes, these methods lead 10 an exponential convergence in the range of engineering compulation. A numerical example shows their effeclivity and accuracy.     

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.     

Quasiwidths and optimization of methods of mixed approximation of multidimensional singular integrals with kernels of hilbert type

We consider the problem of application of mixed methods to the construction of algorithms, optimal in accuracy, for the calculation of multidimensional singular integrals with Hilbert-type kernels. We propose a method for the optimization of cubature formulas for singular integrals with Hilbert-type kernels based on the theory of quasiwidths.     

位势问题边界元法中几乎奇异积分的正则化

将一种通用算法应用于平面位势问题边界元法中近边界点几乎奇异积分的正则化。对线性单元,位势问题近边界点的几乎强和超奇异积分可归纳为两种形式。通过分部积分,将引起奇异的积分元素变换到积分号之外,从而对这两种积分分别给出了无奇异的正则化计算公式。除了线性元,二次元也应用于该算法。与近边界点临近的二次单元划分为两段线性单元,该算法仍然适用。算例证明了方法的有效性和精确性。对曲线边界问题,联合二次元和线性元可提高计算结果精确度。     

QUASI-INTERPOLATING OPERATORS AND THEIR APPLICATIONS IN HYPERSINGULAR INTEGRALS

1.IntroductionSinceP.ZwartobtainedanexpressionofbivariateB--spline[2],R.--HWangandC.K.Chuihavedevelopedaseriesofresults,especially,thequasi--interpolatingoperatorsofSI(AL.)onuniformtype-2triangulationanditsapproximationproperties[']whichhavewidespreadapplicationsinMechanicsandEngineering.Furthermore,R-HWangandC.K.ChuialsoobtainedthefunctionwithminimumsupportinSI(A:)onnon-uniformtype-2triangulationandthebasisofSI(a:)[#].Inthispaperweintroducesomequasi-interpolatingoperatorsofSa(~Z~)o…     

A Novel Time‐Domain BIEM for Wave Propagation Analysis in Anisotropic Solids With Cracks

Analysis of elastic wave propagation in anisotropic solids with cracks is of particular interest to quantitative non‐destructive testing and fracture mechanics. For this purpose, a novel time‐domain boundary integral equation method (BIEM) is presented in this paper. A finite crack in an unbounded elastic solid of general anisotropy subjected to transient elastic wave loading is considered. Two‐dimensional plane strain or plane stress condition is assumed. The initial‐boundary value problem is formulated as a set of hypersingular time‐domain traction boundary integral equations (BIEs) with the crack‐opening‐displacements (CODs) as unknown quantities. A time‐stepping scheme is developed for solving the hypersingular time‐domain BIEs. The scheme uses the convolution quadrature formula of Lubich [1] for temporal convolution and a Galerkin method for spatial discretization of the BIEs. An important feature of the present time‐domain BIEM is that it uses the Laplace‐domain instead of the more complicated time‐domain Green's functions. Fourier integral representations of Laplace‐domain Green's functions are applied. No special technique is needed in the present time‐domain BIEM for evaluating hypersingular integrals.     

Improving the accuracy of multiple integral evaluation by applying Romberg’s method

Romberg’s method, which is used to improve the accuracy of one-dimensional integral evaluation, is extended to multiple integrals if they are evaluated using the product of composite quadrature formulas. Under certain conditions, the coefficients of the Romberg formula are independent of the integral’s multiplicity, which makes it possible to use a simple evaluation algorithm developed for one-dimensional integrals. As examples, integrals of multiplicity two to six are evaluated by Romberg’s method and the results are compared with other methods.     

Numerical solution of a certain hypersingular integral equation of the first kind

In this paper, we first discuss the midpoint rule for evaluating hypersingular integrals with the kernel sin ⁻²[(x−s)/2] defined on a circle, and the key point is placed on its pointwise superconvergence phenomenon. We show that this phenomenon occurs when the singular point s is located at the midpoint of each subinterval and obtain the corresponding supercovergence analysis. Then we apply the rule to construct a collocation scheme for solving the relevant hypersingular integral equation, by choosing the midpoints as the collocation points. It’s interesting that the inverse of coefficient matrix for the resulting linear system has an explicit expression, by which an optimal error estimate is established. At last, some numerical experiments are presented to confirm the theoretical analysis.     

Asymptotic Expansions for Two-Dimensional Hypersingular Integrals

Summary. We define and examine two-dimensional hypersingular integrals on [0,1)² and on [0,)² and relate their Hadamard finite-part (HFP) values to Mellin transforms. These integrands have algebraic singularities of a possibly unintegrable nature on the axes and at the origin. Extending our work on one-dimensional integrals reported in 1998, we obtain variants of the classical Euler-Maclaurin expansion for various two-dimensional integrals. In many cases, the constant term in the expansion (which is not necessarily the leading term) provides the value of the HFP integral. These expansions may be used as the basis for the numerical evaluation of a class of HFP integrals by extrapolation.Mathematics Subject Classification (2000): 65D30, 65B15, 65R10This author was supported in part by the Mathematical, Information, and Computational Sciences Division subprogram of the Office of Advanced Scientific Computing Research, Office of Science, U.S. Department of Energy, under Contract W-31-109-Eng 38. Part of this work was performed while this author was visiting professor at the Politecnico di Torino, under the sponsorship of the Italian C.N.R.-Gruppo Nazionale per lInformatica Matematica.This author was supported by the Ministero dellUniversitá e della Ricerca Scientifica e Tecnologica of Italy, and by the C.N.R.-Comitato n.11.     

Explicit calculation of the partition function of a one-dimensional δ-function bose gas

We present a method to calculate explicitly the partition function of a one-dimensional δ-function bose gas. The cluster integrals are obtained for up to the four particles and are shown to be the same as those by the thermal Bethe ansatz method (TBA).     

Semi-analytic integration of hypersingular Galerkin BIEs for three-dimensional potential problems

An accurate and efficient semi-analytic integration technique is developed for three-dimensional hypersingular boundary integral equations of potential theory. Investigated in the context of a Galerkin approach, surface integrals are defined as limits to the boundary and linear surface elements are employed to approximate the geometry and field variables on the boundary. In the inner integration procedure, all singular and non-singular integrals over a triangular boundary element are expressed exactly as analytic formulae over the edges of the integration triangle. In the outer integration scheme, closed-form expressions are obtained for the coincident case, wherein the divergent terms are identified explicitly and are shown to cancel with corresponding terms from the edge-adjacent case. The remaining surface integrals, containing only weak singularities, are carried out successfully by use of standard numerical cubatures. Sample problems are included to illustrate the performance and validity of the proposed algorithm.