A hybrid collocation method for a nonlinear Volterra integral equation with weakly singular kernel

This work is concerned with the numerical solution of a nonlinear weakly singular Volterra integral equation. Owing to the singular behavior of the solution near the origin, the global convergence order of product integration and collocation methods is not optimal. In order to recover the optimal orders a hybrid collocation method is used which combines a non-polynomial approximation on the first subinterval followed by piecewise polynomial collocation on a graded mesh. Some numerical examples are presented which illustrate the theoretical results and the performance of the method. A comparison is made with the standard graded collocation method.     

Numerical solution of a Fredholm integro-differential equation modelling neural networks

We compare piecewise linear and polynomial collocation approaches for the numerical solution of a Fredholm integro-differential equations modelling neural networks. Both approaches combine the use of Gaussian quadrature rules on an infinite interval of integration with interpolation to a uniformly distributed grid on a bounded interval. These methods are illustrated by numerical experiments on neural networks equations.     

A multi-interval Chebyshev collocation approach for the stability of periodic delay systems with discontinuities

This paper investigates the stability of periodic delay systems with non-smooth coefficients using a multi-interval Chebyshev collocation approach (MIC). In this approach, each piecewise continuous interval is expanded in a Chebyshev basis of the first order. The boundaries of these intervals are placed at the points of discontinuity to recover the fast convergence properties of spectral methods. Stability is examined for a set of case studies that contain the complexities of periodic coefficients, delays and discontinuities. The new approach is also compared to the conventional Chebyshev collocation method.     

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.     

Spline collocation methods for linear multi-term fractional differential equations

Some regularity properties of the solution of linear multi-term fractional differential equations are derived. Based on these properties, the numerical solution of such equations by piecewise polynomial collocation methods is discussed. The results obtained in this paper extend the results of Pedas and Tamme (2011) [15] where we have assumed that in the fractional differential equation the order of the highest derivative of the unknown function is an integer. In the present paper, we study the attainable order of convergence of spline collocation methods for solving general linear fractional differential equations using Caputo form of the fractional derivatives and show how the convergence rate depends on the choice of the grid and collocation points. Theoretical results are verified by some numerical examples.     

Spline collocation methods for linear multi-term fractional differential equations

Some regularity properties of the solution of linear multi-term fractional differential equations are derived. Based on these properties, the numerical solution of such equations by piecewise polynomial collocation methods is discussed. The results obtained in this paper extend the results of Pedas and Tamme (2011) [15] where we have assumed that in the fractional differential equation the order of the highest derivative of the unknown function is an integer. In the present paper, we study the attainable order of convergence of spline collocation methods for solving general linear fractional differential equations using Caputo form of the fractional derivatives and show how the convergence rate depends on the choice of the grid and collocation points. Theoretical results are verified by some numerical examples.     

Piecewise polynomial collocation for linear boundary value problems of fractional differential equations

We consider a class of boundary value problems for linear multi-term fractional differential equations which involve Caputo-type fractional derivatives. Using an integral equation reformulation of the boundary value problem, some regularity properties of the exact solution are derived. Based on these properties, the numerical solution of boundary value problems by piecewise polynomial collocation methods is discussed. In particular, we study the attainable order of convergence of proposed algorithms and show how the convergence rate depends on the choice of the grid and collocation points. Theoretical results are verified by two numerical examples.     

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.

     

Parallel collocation solution of index-1 BVP-DAEs arising from constrained optimal control problems

The indirect solution of constrained optimal control problems gives rise to two-point boundary value problems (BVPs) that involve index-1 differential-algebraic equations (DAEs) and inequality constraints. This paper presents a parallel collocation algorithm for the solution of these inequality constrained index-1 BVP-DAEs. The numerical algorithm is based on approximating the DAEs using piecewise polynomials on a nonuniform mesh. The collocation method is realized by requiring that the BVP-DAE be satisfied at Lobatto points within each interval of the mesh. A Newton interior-point method is used to solve the collocation equations, and maintain feasibility of the inequality constraints. The implementation of the algorithm involves: (i) parallel evaluation of the collocation equations; (ii) parallel evaluation of the system Jacobian; and (iii) parallel solution of a boarded almost block diagonal (BABD) system to obtain the Newton search direction. Numerical examples show that the parallel implementation provides significant speedup when compared to a sequential version of the algorithm.     

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

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

Spectral collocation and radial basis function methods for one-dimensional interface problems

Differential equations with singular sources or discontinuous coefficients yield non-smooth or even discontinuous solutions. This problem is known as the interface problem. High-order numerical solutions suffer from the Gibbs phenomenon in that the accuracy deteriorates if the discontinuity is not properly treated. In this work, we use the spectral and radial basis function methods and present a least squares collocation method to solve the interface problem for one-dimensional elliptic equations. The domain is decomposed into multiple sub-domains; in each sub-domain, the collocation solution is sought. The solution should satisfy more conditions than the given conditions associated with the differential equations, which makes the problem over-determined. To solve the over-determined system, the least squares method is adopted. For the spectral method, the weighted norm method with different scaling factors and the mixed formulation are used. For the radial basis function method, the weighted shape parameter method is presented. Numerical results show that the least squares collocation method provides an accurate solution with high efficacy and that better accuracy is obtained with the spectral method.     

Callocation methods for volterra integrodifferential equations with singular kernels

A collocation method which uses Hermite cubic elements is proposed for the solution of Volterra integrodifferential equations with singular kernels. Optimum error estimates in the uniform norm are obtained by means of interpolation operators. We also report on results of numerical comparisons with one well established method and another new “modified collocation” scheme.     

Collocation methods for fractional integro-differential equations with weakly singular kernels

In this paper, the piecewise polynomial collocation methods are used for solving the fractional integro-differential equations with weakly singular kernels. We present that a suitable transformation can convert fractional integro-differential equations to one type of second kind Volterra integral equations (VIEs) with weakly singular kernels. Then we solve the VIEs by standard piecewise polynomial collocation methods. It is shown that such kinds of methods are able to yield optimal convergence rate. Finally, some numerical experiments are given to show that the numerical results are consistent with the theoretical results.     

A numerical solution of singular integrodifferential equations with variable coefficients

A method for the numerical solution of singular integrodifferential equations is presented where the integrals are discretized by using a convenient quadrature rule. Then the problem is reduced to a system of linear algebraic equations by applying the discretized functional equation to appropriately selected collocation points. This technique constitutes an extension of an analogous method convenient for solving singular integral equations which was proposed by the authors.     

On Spline Collocation for Singular Integral Equations

This paper is devoted to the approximate solution of one-dimensional singular integral equations on a closed curve by spline collocation methods. As the main result we give conditions which are sufficient and in special cases also necessary for the convergence in SOBOLEV norms. The paper is organized as follows. In chapter 1 we indicate some definitions and some facts about projection methods. In chapter 2, we generalize a technique developed in [1] and study the convergence of collocations using splines of odd degree in periodic SOBOLEV spaces. In chapter 3, we apply our method to collocations by splines of even degree and consider the case of systems of equations. And in the last chapter, 4, the results are applied to singular integral equations and compared with known facts about piecewise linear spline collocation for such equations.     

On the numerical solution of Cauchy type singular integral equations by the collocation method

The collocation method for the numerical solution of Fredholm integral equations of the second kind is applied, properly modified, to the numerical solution of Cauchy type singular integral equations of the first or the second kind but with constant coefficients. This direct method of numerical solution of Cauchy type singular integral equations is compared afterwards with the corresponding method resulting from applying the collocation method to the Fredholm integral equation of the second kind equivalent to the Cauchy type singular integral equation, as well as with another method, based also on the regularization procedure, for the numerical solution of the same class of equations. Finally, the convergence of the method is discussed.     

Two‐dimensional shifted Legendre polynomial collocation method for electromagnetic waves in dielectric media via almost operational matrices

In this paper, a numerical solution of fractional partial differential equations (FPDEs) for electromagnetic waves in dielectric media will be discussed. For the solution of FPDEs, we developed a numerical collocation method using an algorithm based on two‐dimensional shifted Legendre polynomials approximation, which is proposed for electromagnetic waves in dielectric media. By implementing the partial Riemann–Liouville fractional derivative operators, two‐dimensional shifted Legendre polynomials approximation and its operational matrix along with collocation method are used to convert FPDEs first into weakly singular fractional partial integro‐differential equations and then converted weakly singular fractional partial integro‐differential equations into system of algebraic equation. Some results concerning the convergence analysis and error analysis are obtained. Illustrative examples are included to demonstrate the validity and applicability of the technique. Copyright © 2017 John Wiley & Sons, Ltd.     

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.     

Open problems in the discretization of volterra integral equations

The open problems (and conjectures) discussed in this paper arise in the discretization of Volterra integral equations, including equations with weakly singular kernels and delay arguments, by collocation methods in piecewise polynomial spaces. They focus on questions of stability versus accuracy; extrapolation on regular and graded meshes; and equations with certain variable delays     

Collocation Methods and Their Modifications for Cauchy Singular Integral Equations on the Interval

In this paper we present polynomial collocation methods and their modi.cations for the numerical solution of Cauchy singular integral equations over the interval [-1, 1]. More precisely, the operators of the integral equations have the form with piecewise continuous coefficients a and b, and with a Jacobi weight . Using the splitting property of the singular values of the collocation methods, we obtain enough stable approximate methods to .nd the least square solution of our integral equation. Moreover, the modifications of the collocation methods enable us to compute kernel and cokernel dimensions of operators from a C*-algebra, which is generated by operators of the Cauchy singular integral equations.