A matrix formulation of the Tau Method for Fredholm and Volterra linear integro-differential equations

In this paper we obtain the matrix Tau Method representation of a general boundary value problem for Fredholm and Volterra integro-differential equations of orderv. Some theoretical results are given that simplify the application of the Tau Method. The application of the Tau Method to the numerical solution of such problems is shown. Numerical results and details of the algorithm confirm the high accuracy and userfriendly structure of this numerical approach.     

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.     

The L<Superscript>2</Superscript>-convergence of the Legendre spectral Tau matrix formulation for nonlinear fractional integro differential equations

The operational Tau method, a well known method for solving functional equations is employed to approximate the solution of nonlinear fractional integro-differential equations. The fractional derivatives are described in the Caputo sense. The unique solvability of the linear Tau algebraic system is discussed. In addition, we provide a rigorous convergence analysis for the Legendre Tau method which indicate that the proposed method converges exponentially provided that the data in the given FIDE are smooth. To do so, Sobolev inequality with some properties of Banach algebras are considered. Some numerical results are given to clarify the efficiency of the method.     

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.     

The compound Poisson risk model with multiple thresholds

In this paper we consider a multi-threshold compound Poisson risk model. A piecewise integro-differential equation is derived for the Gerber-Shiu discounted penalty function. We then provide a recursive approach to obtain general solutions to the integro-differential equation and its generalizations. Finally, we use the probability of ruin to illustrate the applicability of the approach.     

Systems of nonlinear Volterra integro-differential equations

An efficient method based on operational Tau matrix is developed, to solve a type of system of nonlinear Volterra integro-differential equations (IDEs). The presented method is also modified for the problems with separable kernel. Error estimation of the new schemes are analyzed and discussed. The advantages of this approach and its modification is that, the solution can be expressed as a truncated Taylor series, and the error function at any stage can be estimated. Methods are applied on the four problems with separable kernel to show the applicability and efficiency of our schemes, specially for those problems at broad intervals.     

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.

     

CONVERGENCE ANALYSIS OF SPECTRAL METHODS FOR INTEGRO-DIFFERENTIAL EQUATIONS WITH VANISHING PROPORTIONAL DELAYS

We describe the application of the spectral method to delay integro-differential equations with proportional delays. It is shown that the resulting numerical solutions exhibit the spectral convergence order. Extensions to equations with more general (nonlinear) vanishing delays are also discussed.     

Linear differential-algebraic equations perturbed by Volterra integral operators

We study linear inhomogeneous vector ordinary differential equations of arbitrary order in which the matrix multiplying the highest derivative of the unknown vector function is singular in the domain where the equations are defined. We also study perturbations (not necessarily small) of such equations, which are linear integro-differential equations with a Volterra operator. We obtain sufficient conditions for the solvability of such equations and give representations of their general solutions; solvability and uniqueness conditions are also given for initial value problems for such equations. The influence of small perturbations of the free term and the initial data on the solution is considered. A numerical method is suggested. The results of numerical experiments are given.     

The adaptive operational Tau method for systems of ODEs

As the Tau method, like many other numerical methods, has the limitation of using a fixed step size with some high degree (order) of approximation for solving initial value problems over long intervals, we introduce here the adaptive operational Tau method. This limitation is very much problem dependent and in such case the fixed step size application of the Tau method loses the true track of the solution. But when we apply this new adaptive method the true solution is recovered with a reasonable number of steps. To illustrate the effectiveness of this method we apply it to some stiff systems of ordinary differential equations (ODEs). The numerical results confirm the efficiency of the method.     

A new numerical method for delay and advanced integro-differential equations

A general formulation is constructed for Jacobi operational matrices of integration, product, and delay on an arbitrary interval. The main purpose of this study is to improve Jacobi operational matrices for solving delay or advanced integro–differential equations. Some theorems are established and utilized to reduce the computational costs. All algorithms can be used for both linear and nonlinear Fredholm and Volterra integro-differential equations with delay. An error estimator is introduced to approximate the absolute error when the exact solution of a given problem is not available. The error of the proposed method is less compared to other common methods such as the Taylor collocation, Chebyshev collocation, hybrid Euler–Taylor matrix, and Boubaker collocation methods. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments.     

New analytical and CWENO numerical solutions for axisymmetric horizontal flows with heat sources

Some new nonlinear analytical solutions are found for axisymmetric horizontal flows dominated by strong heat sources. These flows are common in multiscale atmospheric and oceanic flows such as hurricane embryos and ocean gyres. The analytical solutions are illustrated with several examples. The proposed exact solutions provide analytical support for previous numerical observations and can be also used as benchmark problems for validating numerical models. A central weighted essentially non-oscillatory (CWENO) reconstruction is also employed for numerical simulation of the corresponding integro-differential equations. Due to the use of the same polynomial reconstruction for all derivatives and integral terms, the balance between those terms is well preserved, and the method can precisely reproduce the exact solutions, which are hard to capture by traditional upwind schemes. The developed analytical solutions were employed to evaluate the performance of the numerical method, which showed an excellent performance of the numerical model in terms of numerical diffusion and oscillation.     

Numerical treatment of a well-posed Chebyshev Tau method for Bagley-Torvik equation with high-order of accuracy

The main purpose of this study is to develop and analyze a new high-order operational Tau method based on the Chebyshev polynomials as basis functions for obtaining the numerical solution of Bagley-Torvik equation which has a important role in the fractional calculus. It is shown that some derivatives of the solutions of these equations have a singularity at origin. To overcome this drawback we first change the original equation into a new equation with a better regularity properties by applying a regularization process and thereby the operational Chebyshev Tau method can be applied conveniently. Our proposed method has two main advantages. First, the algebraic form of the Tau discretization of the problem has an upper triangular structure which can be solved by forward substitution method. Second, Tau approximation of the problem converges to the exact ones with a highly rate of convergence under a more general regularity assumptions on the input data in spite of the singularity behavior of the exact solution. Numerical results are presented which confirm the theoretical results obtained and efficiency of the proposed method.     

Numerical solutions of systems of nonlinear integro-differential equations by Homotopy-perturbation method

In this study, the numerical solutions of a system of two nonlinear integro-differential equations, which describes biological species living together, are derived employing the well-known Homotopy-perturbation method. The approximate solutions are in excellent agreement with those obtained by the Adomian decomposition method. Furthermore, we present an analytical approximate solution for a more general form of the system of nonlinear integro-differential equations. The numerical result indicates that the proposed method is straightforward to implement, efficient and accurate for solving nonlinear integro-differential equations.     

分数积分的一种数值计算方法及其应用

提出了一种只需要存储部分历史数据的分数积分的数值计算方法,并给出了误差估计。这种方法可对包含分数积分和分数导数的积分-微分方程进行较长时间的数值计算,克服了存储全部历史数据的困难,并能对计算误差进行控制。作为应用,给出了具有分数导数型本构关系的粘弹性Timoshenko梁的动力学行为研究的控制方程,利用分离变量法讨论梁在简谐激励作用下的动力响应,然后用新提出的数值方法对控制方程进行数值计算,数值计算结果和理论结果进行了比较,它们比较吻合。     

On numerical approximation using differential equations with piecewise-constant arguments

In this paper we give a brief overview of the application of delay differential equations with piecewise constant arguments (EPCAs) for obtaining numerical approximation of delay differential equations, and we show that this method can be used for numerical approximation in a class of age-dependent population models. We also formulate an open problem for stability and oscillation of a class of linear delay equations with continuous and piecewise constant arguments. This research was partially supported by Hungarian NFSR Grant No. T046929.     

水滴清除气溶胶过程的随机算法和数值模拟

气溶胶尺度分布的时间演变可量化气溶胶的湿沉降过程,它在数学上可由考虑湿沉降的通用动力学方程来描述．该方程为一典型的部分积分微分方程,与气溶胶尺度分布和雨滴尺度分布均相关,且由于需要考虑Brown扩散、拦截和惯性碰撞等湿沉降机制而使得清除系数模型非常复杂,普通的数值方法难以求解．为此发展了一种新的多重Monte Carlo算法,以求解考虑湿沉降的通用动力学方程,并用于模拟实际环境中气溶胶的湿沉降．对于对数正态分布的气溶胶尺度分布和雨滴尺度分布, 多重Monte Carlo算法进行的数值模拟表明, 雨滴几何平均尺度越小, 雨滴几何标准偏差越小,越有利于小尺度和中等尺度气溶胶的湿去除,而稍微不利于大尺度气溶胶的湿去除．     

Numerical solution of the system of nonlinear Volterra integro-differential equations with nonlinear differential part by the operational Tau method and error estimation

The operational Tau method, a well-known method for solving functional equations, is employed to solve a system of nonlinear Volterra integro-differential equations with nonlinear differential part. In addition, an error estimation of the method is presented. Some cases of the mentioned equations are solved as examples to illustrate the ability and reliability of the method. The results reveal that the method is very effective and convenient.     

Discrete Galerkin method for Fredholm integro-differential equations with weakly singular kernels

Approximations to a solution and its derivatives of a boundary value problem of an nth order linear Fredholm integro-differential equation with weakly singular or other nonsmooth kernels are determined. These approximations are piecewise polynomial functions on special graded grids. For their finding a discrete Galerkin method and an integral equation reformulation of the boundary value problem are used. Optimal global convergence estimates are derived and an improvement of the convergence rate of the method for a special choice of parameters is obtained. To illustrate the theoretical results a collection of numerical results of a test problem is presented.     

Convergence and Mean-Square Stability of Exponential Euler Method for Semi-Linear Stochastic Delay Integro-Differential Equations

In this paper, the numerical methods for semi-linear stochastic delay integro-differential equations are studied. The uniqueness, existence and stability of analytic solutions of semi-linear stochastic delay integro-differential equations are studied and some suitable conditions for the mean-square stability of the analytic solutions are also obtained. Then the numerical approximation of exponential Euler method for semi-linear stochastic delay integro-differential equations is constructed and the convergence and the stability of the numerical method are studied. It is proved that the exponential Euler method is convergent with strong order $\frac{1}{2}$ and can keep the mean-square exponential stability of the analytical solutions under some restrictions on the step size. In addition, numerical experiments are presented to confirm the theoretical results.