Continuous extensions to Nyström methods for second order initial value problems

In this work we consider interpolants for Nyström methods, i.e., methods for solving second order initial value problems. We give a short introduction to the theory behind the discrete methods, and extend some of the work to continuous, explicit Nyström methods. Interpolants for continuous, explicit Runge-Kutta methods have been intensively studied by several authors, but there has not been much effort devoted to continuous Nyström methods. We therefore extend some of the work by Owren.     

Avoiding the order reduction of Runge-Kutta methods for linear initial boundary value problems

A new strategy to avoid the order reduction of Runge-Kutta methods when integrating linear, autonomous, nonhomogeneous initial boundary value problems is presented. The solution is decomposed into two parts. One of them can be computed directly in terms of the data and the other satisfies an initial value problem without any order reduction. A numerical illustration is given. This idea applies to practical problems, where spatial discretization is also required, leading to the full order both in space and time.

     

Novel numerical techniques based on Fokas transforms,for the solution of initial boundary value problems

The unified transform method of A. S. Fokas has led to important new developments, regarding the analysis and solution of various types of linear and nonlinear PDE problems. In this work we use these developments and obtain the solution of time-dependent problems in a straightforward manner and with such high accuracy that cannot be reached within reasonable time by use of the existing numerical methods. More specifically, an integral representation of the solution is obtained by use of the A. S. Fokas approach, which provides the value of the solution at any point, without requiring the solution of linear systems or any other calculation at intermediate time levels and without raising any stability problems. For instance, the solution of the initial boundary value problem with the non-homogeneous heat equation is obtained with accuracy 10⁻¹⁵, while the well-established Crank–Nicholson scheme requires 2048 time steps in order to reach a 10⁻⁸ accuracy.     

Jacobi-Picard iteration method for the numerical solution of nonlinear initial value problems

In this paper, an effective numerical iterative method for solving nonlinear initial value problems (IVPs) is presented. The proposed iterative scheme, called the Jacobi-Picard iteration (JPI) method, is based on the Picard iteration technique, orthogonal shifted Jacobi polynomials, and shifted Jacobi-Gauss quadrature formula. In comparison with traditional methods, the JPI method uses an iterative formula for updating next step approximations and calculating integrals of the shifted Jacobi polynomials are performed via an exact relation. Also, a vector-matrix form of the JPI method is provided in details which reduce the CPU time. The performance of the presented method has been investigated by solving several nonlinear IVPs. Numerical results show the efficiency and the accuracy of the proposed iterative method.     

Symmetric multistep Obrechkoff methods with zero phase-lag for periodic initial value problems of second order differential equations

In this paper, symmetric multistep Obrechkoff methods of orders 8 and 12, involving a parameter p to solve a special class of second order initial value problems in which the first order derivative does not appear explicitly, are discussed. It is shown that the methods have zero phase-lag when p is chosen as 2π times the frequency of the given initial value problem.     

A unique solution of initial value problems for first order impulsive integro-differential equations of mixed type in Banach spaces

In this paper, we investigate initial value problem for first order impulsive integro-differential equation of mixed type in a Banach space. Without any compactness-type assumption, we obtain a unique solution of the equation, an explicit iterative approximation of the solution and an error estimate of the approximation sequence. Our result improves and extends many recent results.     

On the stability of two new two-step explicit methods for the numerical integration of second order initial value problem on a variable mesh

We present two novel two-step explicit methods for the numerical solution of the second order initial value problem on a variable mesh. In the case of a constant mesh the method is superstable in the sense of Chawla (1985). Numerical experimentation is provided to verify the stability analysis.     

Optimization of the Relaxation Parameter for S.S.O.R. and A.D.I. Preconditioning

We present and compare several approaches for the optimization of the relaxation parameter both for A.D.I. and S.S.O.R. basic iteration and preconditioning conjugate gradient method. For each kind of preconditioning a detailed link between estimates of the spectral radius of the iteration matrix and of the condition number resulting from preconditioning is proposed. It allows to choose the best approach in order to obtain the optimal relaxation parameter and the corresponding optimal estimates either of the spectral radius of the iteration matrix and of the resulting condition mumber of the S.S.O.R. and A.D.I. preconditioning.