Adaptive meshless Galerkin boundary node methods for hypersingular integral equations

Adaptive refinement techniques are developed in this paper for the meshless Galerkin boundary node method for hypersingular boundary integral equations. Two types of error estimators are derived. One is a perturbation error estimator that is formulated based on the difference between numerical solutions obtained using two consecutive nodal arrangements. The other is a projection error estimator that is formulated based on the difference between the numerical solution itself and its projection. These error estimators are proven to have an upper and a lower bound by the constant multiples of the exact error in the energy norm. A localization scheme is presented to accomodate the non-local property of hypersingular integral operators for the needed computable local error indicators. The convergence of the adaptive meshless techniques is verified theoretically. To confirm the theoretical results and to show the efficiency of the adaptive techniques, numerical examples in 2D and 3D with high singularities are provided.     

A random integral quadrature method for numerical analysis of the second kind of Volterra integral equations

In this paper, a novel meshless technique termed the random integral quadrature (RIQ) method is developed for the numerical solution of the second kind of the Volterra integral equations. The RIQ method is based on the generalized integral quadrature (GIQ) technique, and associated with the Kriging interpolation function, such that it is regarded as an extension of the GIQ technique. In the GIQ method, the regular computational domain is required, in which the field nodes are scattered along straight lines. In the RIQ method however, the field nodes can be distributed either uniformly or randomly. This is achieved by discretizing the governing integral equation with the GIQ method over a set of virtual nodes that lies along straight lines, and then interpolating the function values at the virtual nodes over all the field nodes which are scattered either randomly or uniformly. In such a way, the governing integral equation is converted approximately into a system of linear algebraic equations, which can be easily solved.     

Fast Fourier-Galerkin methods for first-kind logarithmic-kernel integral equations on open arcs

We propose a fully discrete fast Fourier-Galerkin method for solving an integral equation of the first kind with a logarithmic kernel on a smooth open arc,which is a reformulation of the Dirichlet problem of the Laplace equation in the plane.The optimal convergence order and quasi-linear complexity order of the proposed method are established.A precondition is introduced.Combining this method with an efficient numerical integration algorithm for computing the single-layer potential defined on an open arc,we...     

On Taylor-series expansion methods for the second kind integral equations

In this paper, we comment on the recent papers by Yuhe Ren et al. (1999) [1] and Maleknejad et al. (2006) [7] concerning the use of the Taylor series to approximate a solution of the Fredholm integral equation of the second kind as well as a solution of a system of Fredholm equations. The technique presented in Yuhe Ren et al. (1999) [1] takes advantage of a rapidly decaying convolution kernel k(|s−t|) as |s−t| increases. However, it does not apply to equations having other types of kernels. We present in this paper a more general Taylor expansion method which can be applied to approximate a solution of the Fredholm equation having a smooth kernel. Also, it is shown that when the new method is applied to the Fredholm equation with a rapidly decaying kernel, it provides more accurate results than the method in Yuhe Ren et al. (1999) [1]. We also discuss an application of the new Taylor-series method to a system of Fredholm integral equations of the second kind.     

Extrapolation method for solving weakly singular nonlinear Volterra integral equations of the second kind

Based on a new generalization of discrete Gronwall inequality in [L. Tao, H. Yong, A generalization of discrete Gronwall inequality and its application to weakly singular Volterra integral equality of the second kind, J. Math. Anal. Appl. 282 (2003) 56-62], Navot's quadrature rule for computing integrals with the end point singularity in [I. Navot, A further extension of Euler-Maclaurin summation formula, J. Math. Phys. 41 (1962) 155-184] and a transformation in [P. Baratella, A. Palamara Orsi, A new approach to the numerical solution of weakly singular Volterra integral equations, J. Comput. Appl. Math. 163 (2004) 401-418], a new quadrature method for solving nonlinear weakly singular Volterra integral equations of the second kind is presented. The convergence of the approximation solution and the asymptotic expansion of the error are proved, so by means of the extrapolation technique we not only obtain a higher accuracy order of the approximation but also get a posteriori estimate of the error.     

Numerical treatment of retarded boundary integral equations by sparse panel clustering

Stefan SauterInstitute for Mathematics, University of Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland ^{We consider the wave equation in a boundary integral formulation.The discretization in time is done by using convolution quadraturetechniques and a Galerkin boundary element method for the spatialdiscretization. In a previous paper, we have introduced a sparseapproximation of the system matrix by cut-off, in order to reducethe storage costs. In this paper, we extend this approach byintroducing a panel clustering method to further reduce thesecosts.}     

Block boundary value methods for solving Volterra integral and integro-differential equations

Reducible quadrature rules generated by boundary value methods are considered in block version and applied to solve the second kind Volterra integral equations and Volterra integro-differential equations. These extended block boundary value methods are shown to possess both excellent stability properties and high accuracy for Volterra-type equations. Numerical experiments are presented and the efficiency, accuracy and stability of the schemes are confirmed.     

Numerical solution of nonlinear two-dimensional Fredholm integral equations of the second kind using Gauss product quadrature rules

The Gauss product quadrature rules and collocation method are applied to reduce the second-kind nonlinear two-dimensional Fredholm integral equations (FIE) to a nonlinear system of equations. The convergence of the proposed numerical method is proved under certain conditions on the kernel of the integral equation. An iterative method for approximating the solution of the obtained nonlinear system is provided and its convergence is proved. Also, some numerical examples are presented to show the efficiency and accuracy of the proposed method.     

Simple error estimators for the Galerkin BEM for some hypersingular integral equation in 2D

A posteriori error estimation is an important tool for reliable and efficient Galerkin boundary element computations. For hypersingular integral equations in 2D with a positive-order Sobolev space, we analyse the mathematical relation between the (h???h/2)-error estimator from [S. Ferraz-Leite and D. Praetorius, Simple a posteriori error estimators for the h-version of the boundary element method, Computing 83 (2008), pp. 135–162], the two-level error estimator from [M. Maischak, P. Mund, and E. Stephan, Adaptive multilevel BEM for acoustic scattering, 585 Comput. Methods Appl. Mech. Eng. 150 (1997), pp. 351–367], and the averaging error estimator from [C. Carstensen and D. Praetorius, Averaging techniques for the a posteriori bem error control for a hypersingular integral equation in two dimensions, SIAM J. Sci. Comput. 29 (2007), pp. 782–810]. All of these a posteriori error estimators are simple in the following sense: first, the numerical analysis can be done within the same mathematical framework, namely localization techniques for the energy norm. Second, there is almost no implementational overhead for the realization.     

Reducible quadrature methods for Volterra integral equations of the first kind

Quadrature rules, generated by linear multistep methods for ordinary differential equations, are employed to construct a wide class of direct quadrature methods for the numerical solution of first kind Volterra integral equations. Our class covers several methods previously considered in the literature. The methods are convergent provided that both the first and second characteristic polynomial of the linear multistep method satisfy the root condition. Furthermore, the stability behaviour for fixed positive values of the stepsizeh is analyzed, and it turns out that convergence implies (fixedh) stability. The subclass formed by the backward differentiation methods up to order six is discussed and illustrated with numerical examples.     

A system of boundary integral equations for the transmission problem in acoustics

In this paper, we reduce the classical two-dimensional transmission problem in acoustic scattering to a system of coupled boundary integral equations (BIEs), and consider the weak formulation of the resulting equations. Uniqueness and existence results for the weak solution of corresponding variational equations are established. In contrast to the coupled system in Costabel and Stephan (1985) [4], we need to take into account exceptional frequencies to obtain the unique solvability. Boundary element methods (BEM) based on both the standard and a two-level fast multipole Galerkin schemes are employed to compute the solution of the variational equation. Numerical results are presented to verify the efficiency and accuracy of the numerical methods.     

Boundary integral methods for singularly perturbed boundary value problems

In this paper we consider boundary integral methods appliedto boundary value problems for the positive definite Helmholtz-typeproblem –U + ²U = 0 in a bounded or unbounded domain,with the parameter real and possibly large. Applications arisein the implementation of space–time boundary integralmethods for the heat equation, where is proportional to 1/(t),and t is the time step. The corresponding layer potentials arisingfrom this problem depend nonlinearly on the parameter and havekernels which become highly peaked as , causing standard discretizationschemes to fail. We propose a new collocation method with arobust convergence rate as . Numerical experiments on a modelproblem verify the theoretical results.     

Regularized collocation method for Fredholm integral equations of the first kind

In this paper we consider a collocation method for solving Fredholm integral equations of the first kind, which is known to be an ill-posed problem. An “unregularized” use of this method can give reliable results in the case when the rate at which smallest singular values of the collocation matrices decrease is known a priori. In this case the number of collocation points plays the role of a regularization parameter. If the a priori information mentioned above is not available, then a combination of collocation with Tikhonov regularization can be the method of choice. We analyze such regularized collocation in a rather general setting, when a solution smoothness is given as a source condition with an operator monotone index function. This setting covers all types of smoothness studied so far in the theory of Tikhonov regularization. One more issue discussed in this paper is an a posteriori choice of the regularization parameter, which allows us to reach an optimal order of accuracy for deterministic noise model without any knowledge of solution smoothness.     

Sparse convolution quadrature for time domain boundary integral formulations of the wave equation

S. A. Sauter Institut für Mathematik, Universität Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland ^{Many important physical applications are governed by the waveequation. The formulation as time domain boundary integral equationsinvolves retarded potentials. For the numerical solution ofthis problem, we employ the convolution quadrature method forthe discretization in time and the Galerkin boundary elementmethod for the space discretization. We introduce a simple apriori cut-off strategy where small entries of the system matricesare replaced by zero. The threshold for the cut-off is determinedby an a priori analysis which will be developed in this paper.This analysis will also allow to estimate the effect of additionalperturbations such as panel clustering and numerical integrationon the overall discretization error. This method reduces thestorage complexity for time domain integral equations from O(M2N)to O(M2N logM), where N denotes the number of time steps andM is the dimension of the boundary element space.}     

An algorithm based on the regularization and integral mean value methods for the Fredholm integral equations of the first kind

In this paper, an algorithm based on the regularization and integral mean value methods, to handle the ill-posed multi-dimensional Fredholm equations, is introduced. The application of this algorithm is based on the transforming the first kind equation to a second kind equation by the regularization method. Then, by converting the first kind to a second kind, the integral mean value method is employed to handle the resulting Fredholm integral equations of the second kind. The efficiency of the approach will be shown by applying the procedure on some examples.     

Convergence in the integral metric of the general projective method for solving singular integral equations of the first kind with the Cauchy kernel

We study a projective method for solving singular integral equations of the first kind with the Cauchy kernel. Depending on the index of the equation, we introduce pairs of weight spaces which represent a restriction of the space of summable functions. We prove the correctness of the stated problem. We obtain sufficient conditions for the convergence of the projective method in the integral metric.     

含积分边界条件的分数阶微分方程边值问题的正解的存在性

研究了含积分边界条件的分数阶微分方程的边值问题,首先给出格林函数及性质,其次将问题转化为一个等价的积分方程,最后应用Krasnoselkii及Leggett-Williams不动点定理得到了一个及多个正解的存在性,推广了以往的结果．     

Solution methods for the generalized integral volterra equation of the first kind

This paper is devoted to exact and approximate methods (first of all, direct ones) for the solution of integro-operational equations. Themost attention is paid to the theoretical substantiation of the collocation method for the solution of the mentioned equations within the general theory of approximate methods developed by L. V. Kantorovich.     

Numerical treatment of sth order boundary value problems by solving second kind Fredholm integral equations

A Nyström method is proposed for solving Fredholm integral equations equivalent to boundary value problems of order s with complete differential equations. The stability and the convergence of the proposed procedure are proved. Some numerical examples are provided in order to illustrate the accuracy of the method and to compare the procedure with some other ones given in the literature.