Solving systems of fractional differential equations using differential transform method

This paper presents approximate analytical solutions for systems of fractional differential equations using the differential transform method. The fractional derivatives are described in the Caputo sense. The application of differential transform method, developed for differential equations of integer order, is extended to derive approximate analytical solutions of systems of fractional differential equations. The solutions of our model equations are calculated in the form of convergent series with easily computable components. Some examples are solved as illustrations, using symbolic computation. The numerical results show that the approach is easy to implement and accurate when applied to systems of fractional differential equations. The method introduces a promising tool for solving many linear and nonlinear fractional differential equations.     

Numerical comparison of methods for solving linear differential equations of fractional order

In this article, we implement relatively new analytical techniques, the variational iteration method and the Adomian decomposition method, for solving linear differential equations of fractional order. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of fractional differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. This paper will present a numerical comparison between the two methods and a conventional method such as the fractional difference method for solving linear differential equations of fractional order. The numerical results demonstrates that the new methods are quite accurate and readily implemented.     

Numerical methods for nonlinear partial differential equations of fractional order

In this article, we implement relatively new analytical techniques, the variational iteration method and the Adomian decomposition method, for solving nonlinear partial differential equations of fractional order. The fractional derivatives are described in the Caputo sense. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of fractional differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. Numerical results show that the two approaches are easy to implement and accurate when applied to partial differential equations of fractional order.     

Global error estimation and control in linearly-implicit parallel two-step peer W-methods

The class of linearly-implicit parallel two-step peer W-methods has been designed recently for efficient numerical solutions of stiff ordinary differential equations. Those schemes allow for parallelism across the method, that is an important feature for implementation on modern computational devices. Most importantly, all stage values of those methods possess the same properties in terms of stability and accuracy of numerical integration. This property results in the fact that no order reduction occurs when they are applied to very stiff problems. In this paper, we develop parallel local and global error estimation schemes that allow the numerical solution to be computed for a user-supplied accuracy requirement in automatic mode. An algorithm of such global error control and other technical particulars are also discussed here. Numerical examples confirm efficiency of the presented error estimation and stepsize control algorithm on a number of test problems with known exact solutions, including nonstiff, stiff, very stiff and large-scale differential equations. A comparison with the well-known stiff solver RODAS is also shown.     

On Multistep Methods for Differential Equations of Fractional Order

This paper concerns with numerical methods for the treatment of differential equations of fractional order. Our attention is concentrated on fractional multistep methods of both implicit and explicit type, for which order conditions and stability properties are investigated. Dedicated to the memory of Professor Aldo Cossu     

An improved approximate Newton method for implicit Runge-Kutta formulas

Implicit Runge-Kutta (IRK) methods (such as the s-stage Radau IIA method with s=3,5, or 7) for solving stiff ordinary differential equation systems have excellent stability properties and high solution accuracy orders, but their high computing costs in solving their nonlinear stage equations have seriously limited their applications to large scale problems. To reduce such a cost, several approximate Newton algorithms were developed, including a commonly used one called the simplified Newton method. In this paper, a new approximate Jacobian matrix and two new test rules for controlling the updating of approximate Jacobian matrices are proposed, yielding an improved approximate Newton method. Theoretical and numerical analysis show that the improved approximate Newton method can significantly improve the convergence and performance of the simplified Newton method.     

Numerical solution of delay differential equations by uniform corrections to an implicit Runge-Kutta method

Summary In this paper we develop a class of numerical methods to approximate the solutions of delay differential equations. They are essentially based on a modified version, in a predictor-corrector mode, of the one-step collocation method atn Gaussian points. These methods, applied to ODE's, provide a continuous approximate solution which is accurate of order 2n at the nodes and of ordern+1 uniformly in the whole interval. In order to extend the methods to delay differential equations, the uniform accuracy is raised to the order 2n by some a posteriori corrections. Numerical tests and comparisons with other methods are made on real-life problems.This work was supported by CNR within the Progetto Finalizzato Informatica-Sottopr. P1-SOFMAT     

Solitary wave solutions for a time-fraction generalized Hirota–Satsuma coupled KdV equation by an analytical technique

In this work, we implement a relatively analytical technique, the homotopy perturbation method (HPM), for solving nonlinear partial differential equations of fractional order. The fractional derivatives are described in Caputo derivatives. This method can be used as an alternative to obtain analytic and approximate solutions of different types of fractional differential equations which applied in engineering mathematics. The corresponding solutions of the integer order equations are found to follow as special cases of those of fractional order equations. He’s homotopy perturbation method (HPM) which does not need small parameter is implemented for solving the differential equations. It is predicted that HPM can be found widely applicable in engineering.     

Explicit methods for fractional differential equations and their stability properties

The use of explicit methods in the numerical treatment of differential equations of fractional order is an area not yet widely investigated. In this paper stability properties of some multistep methods of explicit type are investigated and new methods with larger intervals of stability are proposed. Some numerical experiments are presented in order to validate theoretical results.     

Periodic orbits of delay differential equations under discretization

This paper deals with the long-time behaviour of numerical solutions of delay differential equations that have asymptotically stable periodic orbits. It is shown that Runge-Kutta discretizations of such equations have attractive invariant curves which approximate the periodic orbit with the order of the method. The research by this author has been made possible by a fellowship of the Royal Netherlands Academy of Arts and Sciences.     

Homotopy perturbation method for fractional KdV equation

In this paper, the homotopy perturbation method is directly applied to derive approximate solutions of the fractional KdV equation. The results reveal that the proposed method is very effective and simple for solving approximate solutions of fractional differential equations.     

A novel method for travelling wave solutions of fractional Whitham–Broer–Kaup,fractional modified Boussinesq and fractional approximate long wave equations in shallow water

In this paper, the analytical approximate traveling wave solutions of Whitham–Broer–Kaup (WBK) equations, which contain blow‐up solutions and periodic solutions, have been obtained by using the coupled fractional reduced differential transform method. By using this method, the solutions were calculated in the form of a generalized Taylor series with easily computable components. The convergence of the method as applied to the WBK equations is illustrated numerically as well as analytically. By using the present method, we can solve many linear and nonlinear coupled fractional differential equations. The results justify that the proposed method is also very efficient, effective and simple for obtaining approximate solutions of fractional coupled modified Boussinesq and fractional approximate long wave equations. Numerical solutions are presented graphically to show the reliability and efficiency of the method. Moreover, the results are compared with those obtained by the Adomian decomposition method (ADM) and variational iteration method (VIM), revealing that the present method is superior to others. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical solution of the fractional-order Vallis systems using multi-step differential transformation method

In this paper, we examine a fractional order Vallis systems. We also fulfil a detailed analysis on the stability of equilibrium. Multi-step differential transform method (MsDTM) extends to give approximate and analytical solutions of a fractional order Vallis systems. Numerical simulations are submitted to verify the validity and reliability results obtained from these methods.     

Approximate analytical solutions of fractional gas dynamic equations

In this article, differential transform method (DTM) has been successfully applied to obtain the approximate analytical solutions of the nonlinear homogeneous and non-homogeneous gas dynamic equations, shock wave equation and shallow water equations with fractional order time derivatives. The true beauty of the article is manifested in its emphatic application of Caputo fractional order time derivative on the classical equations with the achievement of the highly accurate solutions by the known series solutions and even for more complicated nonlinear fractional partial differential equations (PDEs). The method is really capable of reducing the size of the computational work besides being effective and convenient for solving fractional nonlinear equations. Numerical results for different particular cases of the equations are depicted through graphs.     

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.     

Approximate analytical solutions of Schnakenberg systems by homotopy analysis method

In this paper, the homotopy analysis method (HAM) has been employed to obtain analytical solution of a two reaction–diffusion systems of fractional order (fractional Schnakenberg systems) which has been modeling morphogen systems in developmental biology. Different from all other analytic methods, HAM provides us with a simple way to adjust and control the convergence region of solution series by choosing proper values for auxiliary parameter h. The fractional derivative is described in the Caputo sense. The reason of using fractional order differential equations (FOD) is that FOD are naturally related to systems with memory which exists in most biological systems. Also they are closely related to fractals which are abundant in biological systems. The results derived of the fractional system are of a more general nature. Respectively, solutions of FOD spread at a faster rate than the classical differential equations, and may exhibit asymmetry. However, the fundamental solutions of these equations still exhibit useful scaling properties that make them attractive for applications.     

Error analysis and applications of the Fourier-Galerkin Runge-Kutta schemes for high-order stiff PDEs

An integrating factor mixed with Runge-Kutta technique is a time integration method that can be efficiently combined with spatial spectral approximations to provide a very high resolution to the smooth solutions of some linear and nonlinear partial differential equations. In this paper, the novel hybrid Fourier-Galerkin Runge-Kutta scheme, with the aid of an integrating factor, is proposed to solve nonlinear high-order stiff PDEs. Error analysis and properties of the scheme are provided. Application to the approximate solution of the nonlinear stiff Korteweg-de Vries (the 3rd order PDE, dispersive equation), Kuramoto-Sivashinsky (the 4th order PDE, dissipative equation) and Kawahara (the 5th order PDE) equations are presented. Comparisons are made between this proposed scheme and the competing method given by Kassam and Trefethen. It is found that for KdV, KS and Kawahara equations, the proposed method is the best.     

Stability of the Radau IA and Lobatto IIIC methods for delay differential equations

Recently, we have proved that the Radau IA and Lobatto IIIC methods are P-stable, i.e., they have an analogous stability property to A-stability with respect to scalar delay differential equations (DDEs). In this paper, we study stability of those methods applied to multidimensional DDEs. We show that they have a similar property to P-stability with respect to multidimensional equations which satisfy certain conditions for asymptotic stability of the zero solutions. The conditions are closely related to stability criteria for DDEs considered in systems theory. Received October 8, 1996 / Revised version received February 21, 1997     

Numerical solutions and solitary wave solutions of fractional KDV equations using modified fractional reduced differential transform method

In this paper, the modified fractional reduced differential transform method (MFRDTM) has been proposed and it is implemented for solving fractional KdV (Korteweg-de Vries) equations. The fractional derivatives are described in the Caputo sense. In this paper, the reduced differential transform method is modified to be easily employed to solve wide kinds of nonlinear fractional differential equations. In this new approach, the nonlinear term is replaced by its Adomian polynomials. Thus the nonlinear initial-value problem can be easily solved with less computational effort. In order to show the power and effectiveness of the present modified method and to illustrate the pertinent features of the solutions, several fractional KdV equations with different types of nonlinearities are considered. The results reveal that the proposed method is very effective and simple for obtaining approximate solutions of fractional KdV equations.     

Homotopy perturbation method for fractional KdV-Burgers equation

In the paper, we extend the homotopy perturbation method to solve nonlinear fractional partial differential equations. The time- and space-fractional KdV-Burgers equations with initial conditions are chosen to illustrate our method. As a result, we successfully obtain some available approximate solutions of them. The results reveal that the proposed method is very effective and simple for obtaining approximate solutions of nonlinear fractional partial differential equations.