A meshless local Galerkin method for the numerical solution of Hammerstein integral equations based on the moving least squares technique

In this paper, a computational scheme is proposed to estimate the solution of one- and two-dimensional Fredholm-Hammerstein integral equations of the second kind. The method approximates the solution using the discrete Galerkin method based on the moving least squares (MLS) approach as a locally weighted least squares polynomial fitting. The discrete Galerkin technique for integral equations results from the numerical integration of all integrals in the system corresponding to the Galerkin method. Since the proposed method is constructed on a set of scattered points, it does not require any background meshes and so we can call it as the meshless local discrete Galerkin method. The implication of the scheme for solving two-dimensional integral equations is independent of the geometry of the domain. The new method is simple, efficient and more flexible for most classes of nonlinear integral equations. The error analysis of the method is provided. The convergence accuracy of the new technique is tested over several Hammerstein integral equations and obtained results confirm the theoretical error estimates.     

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.     

Approximate solution of a class of singular integral equations of second kind

A simple method based on polynomial approximation of a function is employed to obtain approximate solution of a class of singular integral equations of the second kind. For a hypersingular integral equation of the second kind, this method avoids the complex function-theoretic method and produces the known exact solution to Prandtl's integral equation as a special case. For a particular singular integro-differential equation of the second kind, this also produces an approximate solution which compares favourably with numerical results obtained by various Galerkin methods. The convergence of the method for both the equations is also established.     

求解一类具有Hibert核的奇异积分方程的小波方法

1　引　　言近年来,用小波方法数值求解积分方程越来越引起人们的注意.文献[1]提出的算法可将一类积分算子所对应的矩阵稀疏化,为小波方法快速求解积分方程开辟了一条新的道路这方面的研究不仅可以深入发展小波理论和应用算法,深入发展小波方法的功效,而且对边界元方法有重要的指导意义.然而研究稳健快速的数值方法,一直是这方面研究的难点问题.本文考虑带Hilbert核的奇异积分方程q(y)=12π∫2π0f(x)ctg12(x-y)dx,y∈[0,2π],(1.1)的小波数值解法;其中f(x)∈H2π,q(y)∈H2π是以2π为周期的Holder类函数;q(y)已知,f(x)待求解;(1.1)式右…     

Explicit expressions for three-dimensional boundary integrals in linear elasticity

On employing isoparametric, piecewise linear shape functions over a flat triangle, exact formulae are derived for all surface potentials involved in the numerical treatment of three-dimensional singular and hyper-singular boundary integral equations in linear elasticity. These formulae are valid for an arbitrary source point in space and are represented as analytical expressions along the edges of the integration triangle. They can be employed to solve integral equations defined on triangulated surfaces via a collocation method or may be utilized as analytical expressions for the inner integrals in a Galerkin technique. A numerical example involving a unit triangle and a source point located at various distances above it, as well as sample problems solved by a collocation boundary element method for the Lamé equation are included to validate the proposed formulae.     

Discretization Methods for Three-Dimensional Singular Integral Equations of Electromagnetism

Theorems providing the convergence of approximate solutions of linear operator equations to the solution of the original equation are proved. The obtained theorems are used to rigorously mathematically justify the possibility of numerical solution of the 3D singular integral equations of electromagnetism by the Galerkin method and the collocation method.     

An augmented galerkin procedure for the boundary integral method applied to two-dimensional screen and crack problems

Here we present a new solution procedure for Helm-holtz and Laplacian Dirichlet screen and crack problems in IR² via boundary integral equations of the first kind having as an unknown the jump of the normal derivative across the screen or a crack curve T. Under the assumption of local finite energy we show the equivalence of the integral equations and the original boundary value problem. Via the method of local Mellin transform in [5]-[lo] and the calculus of pseudodifferential operators we derive existence, uniqueness and regularity results for the solution of our boundary integral equations together with its explicit behaviour near the screen or crack tips.With our integral equations we set up a Galerkin scheme on T and obtain high quasi-optimal convergence rates by using special singular elements besides regular splines as test and trial functions.     

Asymptotic expansions of the error of spline Galerkin boundary element methods

Summary. In this paper we analyse the existence of asymptotic expansions of the error of Galerkin methods with splines of arbitrary degree for the approximate solution of integral equations with logarithmic kernels. These expansions are obtained in terms of an interpolation operator and are useful for the application of Richardson extrapolation and for obtaining sharper error bounds. We also present and analyse a family of fully discrete spline Galerkin methods for the solution of the same equations. Following the analysis of Galerkin methods, we show the existence of asymptotic expansions of the error. Received May 18, 1995 / Revised version received April 11, 1996     

Convolution quadrature and discretized operational calculus. II

Summary Operational quadrature rules are applied to problems in numerical integration and the numerical solution of integral equations: singular integrals (power and logarithmic singularities, finite part integrals), multiple timescale convolution, Volterra integral equations, Wiener-Hopf integral equations. Frequency domain conditions, which determine, the stability of such equations, can be carried over to the discretization.This is Part II to the article with the same title (Part I), which was published in Volume 52, No. 2, pp. 129–145 (1988)     

Finite difference scheme based on proper orthogonal decomposition for the nonstationary Navier-Stokes equations

The proper orthogonal decomposition(POD)and the singular value decomposition(SVD) are used to study the finite difference scheme(FDS)for the nonstationary Navier-Stokes equations. Ensembles of data are compiled from the transient solutions computed from the discrete equation system derived from the FDS for the nonstationary Navier-Stokes equations.The optimal orthogonal bases are reconstructed by the elements of the ensemble with POD and SVD.Combining the above procedures with a Galerkin projection approach yields a new optimizing FDS model with lower dimensions and a high accuracy for the nonstationary Navier-Stokes equations.The errors between POD approximate solutions and FDS solutions are analyzed.It is shown by considering the results obtained for numerical simulations of cavity flows that the error between POD approximate solution and FDS solution is consistent with theoretical results.Moreover,it is also shown that this validates the feasibility and efficiency of POD method.     

带常系数的Cauchy型奇异积分方程的快速方法

Petrov-Galerkin 方法是研究Cauchy型奇异积分方程的最基本的数值方法. 用此方法离散积分方程可得一系数矩阵是稠密的线性方程组. 如果方程组的阶比较大, 则求解此方程组所需要的计算复杂度则会变得很大. 因此, 发展此类方程的快速数值算法就变成了必然. 该文将就对带常系数的Cauchy型奇异积分方程给出一种快速数值方法. 首先用一稀疏矩阵来代替稠密系数矩阵, 其次用数值积分公式离散上述方程组得到其完全离散的形式,然后用多层扩充方法求解此完全离散的线性方程组. 证明经过上述过程得到方程组的逼进解仍然保持了最优阶, 并且整个过程所需要的计算复杂度是拟线性的. 最后通过数值实验证明结论.     

Discrete boundary element methods on general meshes in 3D

Abstract. This paper is concerned with the stability and convergence of fully discrete Galerkin methods for boundary integral equations on bounded piecewise smooth surfaces in . Our theory covers equations with very general operators, provided the associated weak form is bounded and elliptic on , for some . In contrast to other studies on this topic, we do not assume our meshes to be quasiuniform, and therefore the analysis admits locally refined meshes. To achieve such generality, standard inverse estimates for the quasiuniform case are replaced by appropriate generalised estimates which hold even in the locally refined case. Since the approximation of singular integrals on or near the diagonal of the Galerkin matrix has been well-analysed previously, this paper deals only with errors in the integration of the nearly singular and smooth Galerkin integrals which comprise the dominant part of the matrix. Our results show how accurate the quadrature rules must be in order that the resulting discrete Galerkin method enjoys the same stability properties and convergence rates as the true Galerkin method. Although this study considers only continuous piecewise linear basis functions on triangles, our approach is not restricted in principle to this case. As an example, the theory is applied here to conventional “triangle-based” quadrature rules which are commonly used in practice. A subsequent paper [14] introduces a new and much more efficient “node-based” approach and analyses it using the results of the present paper. Received December 10, 1997 / Revised version received May 26, 1999 / Published online April 20, 2000 –? Springer-Verlag 2000     

Hypersingular kernel integration in 3D Galerkin boundary element method

We consider Cauchy singular and Hypersingular boundary integral equations associated with 3D potential problems defined on polygonal domains, whose solutions are approximated with a Galerkin boundary element method, related to a given triangulation of the boundary. In particular, for constant and linear shape functions, the most frequently used basis functions, we give explicit results of the analytical inner integrations and suggest suitable quadrature schemes to evaluate the outer integrals required to form the Galerkin matrix elements. These numerical indications are given after an analysis of the singularities arising in the whole integration process, which is valid also for shape functions of higher degrees.     

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     

封闭曲线上的奇异积分方程的样条间接逼近解法

利用复插值样条函数，给出了定义于光滑封闭曲线上一般的正则型奇异积分方程的样条间接逼近解法，证明了一致收敛性．对于其中的一类奇异积分方程，还给出了近似解和误差估计．     

非线性积分方程迭代配置法的渐近展开及其外推

非线性积分方程迭代配置法的渐近展开及其外推韩国强（华南理工大学计算机工程与科学系）ＡＳＹＭＰＴＯＴＩＣＥＲＲＯＲＥＸＭＮＳＩＯＮＳＡＮＤＥＸＴＲＡＰＯＬＡＴＩＯＮＦＯＲＴＨＥＩＴＥＲＡＴＥＤＣＯＬＬＯＣＡＴＩＯＮＭＥＴＨＯＤＳＯＦＮＯＮＬＩＮＥＡＲＩ...     

On the numerical solution of the generalized Prandtl equation using variation-diminishing splines

Convergence results are proved for Cauchy principal value integrals of the Schoenberg variation-diminishing splines and its first derivative. The use of such splines in the numerical solution of the Prandtl and generalized Prandtl integral equations is proposed. A Nyström-type method and a modified Nyström method are used and compared computationally.     

A volume integral equation method for periodic scattering problems for anisotropic Maxwell's equations

This paper presents a volume integral equation method for an electromagnetic scattering problem for three-dimensional Maxwell's equations in the presence of a biperiodic, anisotropic, and possibly discontinuous dielectric scatterer. Such scattering problem can be reformulated as a strongly singular volume integral equation (i.e., integral operators that fail to be weakly singular). In this paper, we firstly prove that the strongly singular volume integral equation satisfies a Gårding-type estimate in standard Sobolev spaces. Secondly, we rigorously analyze a spectral Galerkin method for solving the scattering problem. This method relies on the periodization technique of Gennadi Vainikko that allows us to efficiently evaluate the periodized integral operators on trigonometric polynomials using the fast Fourier transform (FFT). The main advantage of the method is its simple implementation that avoids for instance the need to compute quasiperiodic Green's functions. We prove that the numerical solution of the spectral Galerkin method applied to the periodized integral equation converges quasioptimally to the solution of the scattering problem. Some numerical examples are provided for examining the performance of the method.     

Numerical integration in Galerkin meshless methods, applied to elliptic Neumann problem with non-constant coefficients

In this paper, we explore the effect of numerical integration on the Galerkin meshless method used to approximate the solution of an elliptic partial differential equation with non-constant coefficients with Neumann boundary conditions. We considered Galerkin meshless methods with shape functions that reproduce polynomials of degree k?≥?1. We have obtained an estimate for the energy norm of the error in the approximate solution under the presence of numerical integration. This result has been established under the assumption that the numerical integration rule satisfies a certain discrete Green’s formula, which is not problem dependent, i.e., does not depend on the non-constant coefficients of the problem. We have also derived numerical integration rules satisfying the discrete Green’s formula.     

Fully discrete spectral boundary integral methods for Helmholtz problems on smooth closed surfaces in ${\mathbb R}^3$

Summary. In this paper we describe and analyse a class of spectral methods, based on spherical polynomial approximation, for second-kind weakly singular boundary integral equations arising from the Helmholtz equation on smooth closed 3D surfaces diffeomorphic to the sphere. Our methods are fully discrete Galerkin methods, based on the application of special quadrature rules for computing the outer and inner integrals arising in the Galerkin matrix entries. For the outer integrals we use, for example, product-Gauss rules. For the inner integrals, a variant of the classical product integration procedure is employed to remove the singularity arising in the kernel. The key to the analysis is a recent result of Sloan and Womersley on the norm of discrete orthogonal projection operators on the sphere. We prove that our methods are stable for continuous data and superalgebraically convergent for smooth data. Our theory includes as a special case a method closely related to one of those proposed by Wienert (1990) for the fast computation of direct and inverse acoustic scattering in 3D. Received May 29, 2000 / Revised version received March 26, 2001/ Published online December 18, 2001