A convenient technique for solving linear and nonlinear Abel integral equations by the Adomian decomposition method

In this paper, linear and nonlinear Abel integral equations are transformed in such a manner that the Adomian decomposition method can be applied. Some examples with closed-form solutions are studied in detail to further illustrate the proposed technique, and the results obtained indicate this approach is indeed practical and efficient.     

General linear methods for Volterra integral equations

We investigate the class of general linear methods of order p and stage order q=p for the numerical solution of Volterra integral equations of the second kind. Construction of highly stable methods based on the Schur criterion is described and examples of methods of order one and two which have good stability properties with respect to the basic test equation and the convolution one are given.     

A homotopy perturbation algorithm to solve a system of Fredholm–Volterra type integral equations

In this paper, an application of He’s homotopy perturbation (HPM) method is applied to solve the system of Fredholm and Volterra type integral equations, the results revealing that the HPM is very effective and simple.     

On an initial-value method for quickly solving Volterra integral equations: A review

A method of converting nonlinear Volterra equations to systems of ordinary differential equations is compared with a standard technique, themethod of moments, for linear Fredholm equations. The method amounts to constructing a Galerkin approximation when the kernel is either finitely decomposable or approximated by a certain Fourier sum. Numerical experiments from recent work by Bownds and Wood serve to compare several standard approximation methods as they apply to smooth kernels. It is shown that, if the original kernel decomposes exactly, then the method produces a numerical solution which is as accurate as the method used to solve the corresponding differential system. If the kernel requires an approximation, the error is greater, but in examples seems to be around 0.5% for a reasonably small number of approximating terms. In any case, the problem of excessive kernel evaluations is circumvented by the conversion to the system of ordinary differential equations.     

Applications of integral equations in particle-size statistics

We discuss the application of integral equations techniques to two broad areas of particle statistics, namely, stereology and packing. Problems in stereology lead to the inversion of Abel-type integral equations; and we present a brief survey of existing methods, analytical and numerical, for doing this. Packing problems lead to Volterra equations which, in simple cases, can be solved exactly and, in other cases, need to be solved numerically. Methods for doing this are presented along with some numerical results.     

Numerical solution of Volterra integral and integro-differential equations of convolution type by using operational matrices of piecewise constant orthogonal functions

In this paper, we use operational matrices of piecewise constant orthogonal functions on the interval [0,1)$[0,1)$ to solve Volterra integral and integro-differential equations of convolution type without solving any system. We first obtain Laplace transform of the problem and then we find numerical inversion of Laplace transform by operational matrices. Numerical examples show that the approximate solutions have a good degree of accuracy.     

Fast Runge–Kutta methods for nonlinear convolution systems of Volterra integral equations

In this paper fast implicit and explicit Runge–Kutta methods for systems of Volterra integral equations of Hammerstein type are constructed. The coefficients of the methods are expressed in terms of the values of the Laplace transform of the kernel. These methods have been suitably constructed in order to be implemented in an efficient way, thus leading to a very low computational cost both in time and in space. The order of convergence of the constructed methods is studied. The numerical experiments confirm the expected accuracy and computational cost. AMS subject classification (2000)  65R20, 45D05, 44A35, 44A10     

Block by block method for the systems of nonlinear Volterra integral equations

The approach given in this paper leads to numerical methods for solving system of Volterra integral equations which avoid the need for special starting procedures. The method has also the advantages of simplicity of application and at least four order of convergence which is easy to achieve. Also, at each step we get four unknowns simultaneously. A convergence theorem is proved for the described method. Finally numerical examples presented to certify convergence and accuracy of the method.     

An efficient algorithm for the parallel solution of high-dimensional differential equations

The study of high-dimensional differential equations is challenging and difficult due to the analytical and computational intractability. Here, we improve the speed of waveform relaxation (WR), a method to simulate high-dimensional differential-algebraic equations. This new method termed adaptive waveform relaxation (AWR) is tested on a communication network example. Further, we propose different heuristics for computing graph partitions tailored to adaptive waveform relaxation. We find that AWR coupled with appropriate graph partitioning methods provides a speedup by a factor between 3 and 16.     

A stable numerical method for Volterra integral equations with discontinuous kernel

Numerical methods for Volterra integral equations with discontinuous kernel need to be tuned to their peculiar form. Here we propose a version of the trapezoidal direct quadrature method adapted to such a type of equations. In order to delineate its stability properties, we first investigate about the behavior of the solution of a suitable (basic) test equation and then we find out under which hypotheses the trapezoidal direct quadrature method provides numerical solutions which inherit the properties of the continuous problem.     

Boundary value methods for Volterra integral and integro-differential equations

New and effective quadrature rules generated by boundary value methods are introduced. We employ the introduced quadrature rules to construct quadrature methods for the second kind Volterra integral equations and Volterra integro-differential equations. These methods are shown to be effective and possess excellent convergence properties. The nonlinear multigrid method is applied to solve the discrete systems derived from the introduced numerical scheme. Numerical simulations are presented and confirm the efficiency and accuracy of the methods.     

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.     

伴有边界摄动的Volterra型积分微分方程组的奇摄动

讨论了伴有边界摄动的二阶非线性Volterra型积分微分方程组的奇摄动.在适当的条件下,利用对角化技巧证明了解的存在性,构造出解的渐近展开式并给出余项的一致有效估计.     

Constant-sign solutions for systems of singular integral equations of Hammerstein type

We consider the system of Hammerstein integral equations

where T>0 is fixed, ρ_i’s are given functions and the nonlinearities f_i(t,x₁,x₂,…,x_n) can be singular at t=0 and x_j=0 where j{1,2,,n}. Criteria are offered for the existence of constant-sign solutions, i.e., θ_iu_i(t)≥0 for t[0,T] and 1≤i≤n, where θ_i{1,−1} is fixed. The tools used are a nonlinear alternative of Leray–Schauder type, Krasnosel’skii’s fixed point theorem in a cone and Schauder’s fixed point theorem. We also include examples and applications to illustrate the usefulness of the results obtained.     

Spectral methods for weakly singular Volterra integral equations with smooth solutions

We propose and analyze a spectral Jacobi-collocation approximation for the linear Volterra integral equations (VIEs) of the second kind with weakly singular kernels. In this work, we consider the case when the underlying solutions of the VIEs are sufficiently smooth. In this case, we provide a rigorous error analysis for the proposed method, which shows that the numerical errors decay exponentially in the infinity norm and weighted Sobolev space norms. Numerical results are presented to confirm the theoretical prediction of the exponential rate of convergence.     

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.     

On a method of Bownds for solving Volterra integral equations

An initial-value method of Bownds for solving Volterra integral equations is reexamined using a variable-step integrator to solve the differential equations. It is shown that such equations may be easily solved to an accuracy ofO(10^–8), the error depending essentially on that incurred in truncating expansions of the kernel to a degenerate one.This work was sponsored by a University of Nevada at Las Vegas Research Grant.     

The Waveform Relaxation method for systems of differential/algebraic equations

This paper reports efforts towards establishing a parallel numerical algorithm known as Waveform Relaxation (WR) for simulating large systems of differential/algebraic equations. The WR algorithm was established as a relaxation based iterative method for the numerical integration of systems of ODEs over a finite time interval. In the WR approach, the system is broken into subsystems which are solved independently, with each subsystem using the previous iterate waveform as “guesses” about the behavior of the state variables in other subsystems. Waveforms are then exchanged between subsystems, and the subsystems are then resolved repeatedly with this improved information about the other subsystems until convergence is achieved.

In this paper, a WR algorithm is introduced for the simulation of generalized high-index DAE systems. As with ODEs, DAE systems often exhibit a multirate behavior in which the states vary as differing speeds. This can be exploited by partitioning the system into subsystems as in the WR for ODEs. One additional benefit of partitioning the DAE system into subsystems is that some of the resulting subsystems may be of lower index and, therefore, do not suffer from the numerical complications that high-index systems do. These lower index subsystems may therefore be solved by less specialized simulations. This increases the efficiency of the simulation since only a portion of the problem must be solved with specially tailored code. In addition, this paper established solvability requirements and convergence theorems for varying index DAE systems for WR simulation.     

A Volterra integral formulation for determining the periodic solutions of some autonomous, nonlinear, third-order ordinary differential equations

A Volterra integral formulation based on the introduction of a term proportional to the velocity times the square of the (unknown) frequency of oscillation, a new independent variable equal to the original one times the (unknown) frequency of oscillation, the method of variation of parameters and series expansions of both the solution and the frequency of oscillation, is used to determine the periodic solutions to three nonlinear, autonomous, third-order, ordinary differential equations. It is shown that the first term of the series expansion of the frequency of oscillation coincides with that determined from a first-order harmonic balance procedure, whereas the two-term approximation to the frequency of oscillation is shown to be more accurate than that of a parameter perturbation procedure for the second equation considered in this paper. For the third equation, it is shown that the two-term approximation presented in this paper is more accurate than the corresponding one for one of the parameter perturbation methods, and for initial velocities less than one, for the other parameter perturbation approach.     

Existence theorems of global solutions for nonlinear Volterra type integral equations in Banach spaces

In this paper, first a new fixed point theorem is established, and then, by the use of it, the existence theorems of global solutions for nonlinear Volterra type integral equations in Banach spaces are investigated. The results obtained in this paper generalize and improve the results corresponding to those obtained by others.