A note on collocation methods for Volterra integro-differential equations with weakly singular kernels

Spline collocation methods can be used to solve Volterra integro-differentialequations with weakly singular kernels. In order to obtain optimalconvergence behavior, collocation on suitably graded mesheswas considered by H. Brunner [1] This work extends his resultsto more practical values of the grading exponent.     

Block boundary value methods for linear weakly singular Volterra integro-differential equations

A class of block boundary value methods (BBVMs) is constructed for linear weakly singular Volterra integro-differential equations (VIDEs). The convergence and stability of these methods is analysed. It is shown that optimal convergence rates can be obtained by using special graded meshes. Numerical examples are given to illustrate the sharpness of our theoretical results and the computational effectiveness of the methods. Moreover, a numerical comparison with piecewise polynomial collocation methods for VIDEs is given, which shows that the BBVMs are comparable in numerical precision.

     

Numerical methods for 2D weakly singular Volterra integral equations of the second kind

The paper is dedicated to the numerical solution of 2D weakly singular Volterra integral equations with two different kernels. In order to weaken a singularity influence on the numerical computations we transform these equations in both cases into equivalent equations. The piecewise polynomial approximation of the exact solution is then applied. For numerical computing of weakly singular integrals we construct the special Gauss-type cubature formula based on a nonuniform grid. Also we suggest the practical mesh which is less nonuniform than standard graded mesh tk = (k/N)^{r /(1–α)} T, and at the same time it gives an equivalent approximation error for the functions from C^r,α (0, T ]. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

High performance parallel numerical methods for Volterra equations with weakly singular kernels

Non-stationary discrete time waveform relaxation methods for Abel systems of Volterra integral equations using fractional linear multistep formulae are introduced. Fully parallel discrete waveform relaxation methods having an optimal convergence rate are constructed. A significant expression of the error is proved, which allows us to estimate the number of iterations needed to satisfy a prescribed tolerance and allows us to identify the problems where the optimal methods offer the best performance. The numerical experiments confirm the theoretical expectations.     

Stability of collocation for weakly singular Volterra equations

For Abel integral equations of the second kind, we investigatethe stability of the collocation method with polynomial splines.We characterize the stability domain in terms of the eigenvaluesof the power series, which arises from the method employed.Furthermore, we prove strong stability and provide a rigorousdiscussion of the boundary of the stability region. A few examplesare considered in more depth and numerical illustrations aregiven.     

Iterated fast multiscale Galerkin methods for Fredholm integral equations of second kind with weakly singular kernels

We propose iterated fast multiscale Galerkin methods for the second kind Fredholm integral equations with mildly weakly singular kernel by combining the advantages of fast methods and iteration post-processing methods. To study the super-convergence of these methods, we develop a theoretical framework for iterated fast multiscale schemes, and apply the scheme to integral equations with weakly singular kernels. We show theoretically that even the computational complexity is almost optimal, our schemes improve the accuracy of numerical solutions greatly, and exhibit the global super-convergence. Numerical examples are presented to illustrate the theoretical results and the efficiency of the methods.     

Product integration for volterra integral equations of the second kind with weakly singular kernels

We introduce a new numerical approach for solving Volterra integral equations of the second kind when the kernel contains a mild singularity. We give a convergence result. We also present numerical examples which show the performance and efficiency of our method.

     

A spectral collocation method for a weakly singular Volterra integral equation of the second kind

The solution y of a weakly singular Volterra equation of the second kind posed on the interval ?1 ≤ t ≤ 1 has in general a certain singular behaviour near t = ?1: typically, \(|y^{\prime }(t)| \sim (1+t)^{-\mu }\) for a parameter μ ∈ (0, 1). Various methods have been proposed for the numerical solution of these problems, but up to now there has been no analysis that takes into account this singularity when a spectral collocation method is applied directly to the problem. This gap in the literature is filled by the present paper.     

A Jacobi-collocation method for solving second kind Fredholm integral equations with weakly singular kernels

In this work, we propose a Jacobi-collocation method to solve the second kind linear Fredholm integral equations with weakly singular kernels. Particularly, we consider the case when the underlying solutions are sufficiently smooth. In this case, the proposed method leads to a fully discrete linear system. We show that the fully discrete integral operator is stable in both infinite and weighted square norms. Furthermore, we establish that the approximate solution arrives at an optimal convergence order under the two norms. Finally, we give some numerical examples, which confirm the theoretical prediction of the exponential rate of convergence.     

Superconvergence of collocation methods for a class of weakly singular Volterra integral equations

We discuss the application of spline collocation methods to a certain class of weakly singular Volterra integral equations. It will be shown that, by a special choice of the collocation parameters, superconvergence properties can be obtained if the exact solution satisfies certain conditions. This is in contrast with the theory of collocation methods for Abel type equations. Several numerical examples are given which illustrate the theoretical results.     

Spectral collocation methods for nonlinear weakly singular Volterra integro‐differential equations

In this paper, a spectral collocation approximation is proposed for neutral and nonlinear weakly singular Volterra integro‐differential equations (VIDEs) with non‐smooth solutions. We use some suitable variable transformations to change the original equation into a new equation, so that the solution of the resulting equation possesses better regularity, and the the Jacobi orthogonal polynomial theory can be applied conveniently. Under reasonable assumptions on the nonlinearity, we carry out a rigorous error analysis in L^∞ norm and weighted L² norm. To perform the numerical simulations, some test examples (linear and nonlinear) are considered with nonsmooth solutions, and numerical results are presented. Further more, the comparative study of the proposed methods with some existing numerical methods is provided.     

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.     

A fast Fourier-Galerkin method for solving integral equations of second kind with weakly singular kernels

We propose in this paper a convenient way to compress the dense matrix representation of a compact integral operator with a weakly singular kernel under the Fourier basis. This compression leads to a sparse matrix with only ${\mathcal{O}}(n\log n)$ number of nonzero entries, where 2n+1 denotes the order of the matrix. Based on this compression strategy, we develop a fast Fourier-Galerkin method for solving second kind integral equations with weakly singular kernels. We prove that the approximate order of the truncated equation remains optimal and that the spectral condition number of the coefficient matrix of the truncated linear system is uniformly bounded. Furthermore, we develop a fast algorithm for solving the corresponding truncated linear system, which preserves the optimal order of the approximate solution with only ${\mathcal{O}}(n\log^{2}n)$ number of multiplications required. Numerical examples complete the paper.     

On the numerical solution of integral equations of the second kind with weakly singular kernels

The purpose of this paper is to introduce the (Toeplitz) quadrature method for solving Fredholm integral equations of the second kind with mildly singular kernels. We are presented some numerical examples for the computation of the error estimate using the MathCad package.     

Sinc-collocation methods for weakly singular Fredholm integral equations of the second kind

In this paper we propose new numerical methods for linear Fredholm integral equations of the second kind with weakly singular kernels. The methods are developed by means of the Sinc approximation with smoothing transformations, which is an effective technique against the singularities of the equations. Numerical examples show that the methods achieve exponential convergence, and in this sense the methods improve conventional results where only polynomial convergence have been reported so far.     

A discrete collocation method for Fredholm integro-differential equations with weakly singular kernels

A fully discrete version of a piecewise polynomial collocation method is constructed to solve initial or boundary value problems of linear Fredholm integro-differential equations with weakly singular kernels. Using an integral equation reformulation and special graded grids, optimal global convergence estimates are derived. For special values of parameters an improvement of the convergence rate of elaborated numerical schemes is established. Some of our theoretical results are illustrated by numerical experiments.     

Collocation and iterated collocation methods for a class of weakly singular Volterra integral equations

We discuss the convergence properties of spline collocation and iterated collocation methods for a weakly singular Volterra integral equation associated with certain heat conduction problems. This work completes the previous studies of numerical methods for this type of equations with noncompact kernel. In particular, a global convergence result is obtained and it is shown that discrete superconvergence can be achieved with the iterated collocation if the exact solution belongs to some appropriate spaces. Some numerical examples illustrate the theoretical results.     

Superconvergence of collocation methods for Volterra and Abel integral equations of the second kind

Summary This paper deals with the question of the attainable order of convergence in the numerical solution of Volterra and Abel integral equations by collocation methods in certain piecewise polynomial spaces and which are based on suitable interpolatory quadrature for the resulting moment integrals. The use of a (nonlinear) variation of constants formula for the representation of the error function in terms of the defect allows for a unified treatment of equations with continuous and weakly singular kernels.     

The piecewise polynomial collocation method for nonlinear weakly singular Volterra equations

Second-kind Volterra integral equations with weakly singular kernels typically have solutions which are nonsmooth near the initial point of the interval of integration. Using an adaptation of the analysis originally developed for nonlinear weakly singular Fredholm integral equations, we present a complete discussion of the optimal (global and local) order of convergence of piecewise polynomial collocation methods on graded grids for nonlinear Volterra integral equations with algebraic or logarithmic singularities in their kernels.

     

第二类Volterra型积分方程的Chebyshev-Legendre谱配置方法

提出了一种新的求解第二类线性Volterra型积分方程的Chebyshev谱配置方法.该方法分别对方程中积分部分的核函数和未知函数在Chebyshev-Gauss-Lobatto点上进行插值,通过Chebyshev-Legendre变换,把插值多项式表示成Legendre级数形式,从而将积分转换为内积的形式,再利用Legendre多项式的正交性进行计算.利用Chebyshev插值算子在不带权范数意义下的逼近结果,对该方法在理论上给出了L∞范数意义下的误差估计,并通过数值算例验证了算法的有效性和理论分析的正确性.