A New Approach for Solving a Class of Delay Fractional Partial Differential Equations

In this article, a new numerical approach has been proposed for solving a class of delay time-fractional partial differential equations. The approximate solutions of these equations are considered as linear combinations of Müntz–Legendre polynomials with unknown coefficients. Operational matrix of fractional differentiation is provided to accelerate computations of the proposed method. Using Padé approximation and two-sided Laplace transformations, the mentioned delay fractional partial differential equations will be transformed to a sequence of fractional partial differential equations without delay. The localization process is based on the space-time collocation in some appropriate points to reduce the fractional partial differential equations into the associated system of algebraic equations which can be solved by some robust iterative solvers. Some numerical examples are also given to confirm the accuracy of the presented numerical scheme. Our results approved decisive preference of the Müntz–Legendre polynomials with respect to the Legendre polynomials.     

A new operational vector approach for time-fractional subdiffusion equations of distributed order based on hybrid functions

This research study deals with the numerical solutions of linear and nonlinear time-fractional subdiffusion equations of distributed order. The main aim of our approach is based on the hybrid of block-pulse functions and shifted Legendre polynomials. We produce a novel and exact operational vector for the fractional Riemann–Liouville integral and use it via the Gauss–Legendre quadrature formula and collocation method. Consequently, we reduce the proposed equations to systems of equations. The convergence and error bounds for the new method are investigated. Six problems are tested to confirm the accuracy of the proposed approach. Comparisons between the obtained numerical results and other existing methods are provided. Numerical experiments illustrate the reliability, applicability, and efficiency of the proposed method.     

Backward substitution method based on Müntz polynomials for solving the nonlinear space fractional partial differential equations

A computational technique based on the Müntz polynomials and meshless method has been presented for the solution of nonlinear and linear space fractional partial differential equations (PDEs). The meshless method that is used in this study is the new version of backward substitution method (BSM). First, the time-derivative term is discretized by the Crank-Nicolson method. Then, the approximate solution is given as the separation of the approximation of the boundary data and the correcting functions using the Müntz polynomials. In general form, this approximate solution that does not necessarily satisfy the original equation is shown as the sum of the system basics in which it consists free parameters. Finally, these parameters are determined by the BSM method inside the domain. The main advantage of the method is efficiency and reliability that is examined by six numerical experiments.     

Numerical solution of multiterm variable‐order fractional differential equations via shifted Legendre polynomials

In this paper, shifted Legendre polynomials will be used for constructing the numerical solution for a class of multiterm variable‐order fractional differential equations. In the proposed method, the shifted Legendre operational matrix of the fractional variable‐order derivatives will be investigated. The fundamental problem is reduced to an algebraic system of equations using the constructed matrix and the collocation technique, which can be solved numerically. The error estimate of the proposed method is investigated. Some numerical examples are presented to prove the applicability, generality, and accuracy of the suggested method.     

A numerical technique for solving fractional variational problems by Müntz–Legendre polynomials

This paper presents an efficient approximate method for solving a class of fractional variational problems (FVPs). The fractional derivatives are described in the Caputo sense. In the our method, we use the Müntz–Legendre orthonormal basis and the properties of Rayleigh–Ritz method for fractional calculus to reduce FVPs to solve a system of algebraic equations which solved using Newton’s iterative method. Also, an estimation of the error is given in the sense of Sobolev norms. The illustrative examples are provided to demonstrate the applicability and simplicity of the new technique.     

The modulus-based matrix splitting iteration methods for second-order cone linear complementarity problems

This paper presents a new computational technique for solving fractional pantograph differential equations. The fractional derivative is described in the Caputo sense. The main idea is to use Müntz-Legendre wavelet and its operational matrix of fractional-order integration. First, the Müntz-Legendre wavelet is presented. Then a family of piecewise functions is proposed, based on which the fractional order integration of the Müntz-Legendre wavelets are easy to calculate. The proposed approach is used this operational matrix with the collocation points to reduce the under study problem to a system of algebraic equations. An estimation of the error is given in the sense of Sobolev norms. The efficiency and accuracy of the proposed method are illustrated by several numerical examples.     

Solution of nonlinear weakly singular Volterra integral equations using the fractional‐order Legendre functions and pseudospectral method

In this article, our main goal is to render an idea to convert a nonlinear weakly singular Volterra integral equation to a non‐singular one by new fractional‐order Legendre functions. The fractional‐order Legendre functions are generated by change of variable on well‐known shifted Legendre polynomials. We consider a general form of singular Volterra integral equation of the second kind. Then the fractional Legendre–Gauss–Lobatto quadratures formula eliminates the singularity of the kernel of the integral equation. Finally, the Legendre pseudospectral method reduces the solution of this problem to the solution of a system of algebraic equations. This method also can be utilized on fractional differential equations as well. The comparison of results of the presented method and other numerical solutions shows the efficiency and accuracy of this method. Also, the obtained maximum error between the results and exact solutions shows that using the present method leads to accurate results and fast convergence for solving nonlinear weakly singular Volterra integral equations. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of the modified operational matrices in multiterm variable‐order time‐fractional partial differential equations

In this paper, we present a novel discrete scheme based on Genocchi polynomials and fractional Laguerre functions to solve multiterm variable‐order time‐fractional partial differential equations (M‐V‐TFPDEs) in the large interval. In this purpose, the accurate modified operational matrices are constructed to reduce the problems into a system of algebraic equations. Also, the computational algorithm based on the method and modified operational matrices in the large interval is easily implemented. Furthermore, we discuss the error estimation of the proposed method. Ultimately, to confirm our theoretical analysis and accuracy of numerical approach, several examples are presented.     

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.     

A new approach by two‐dimensional wavelets operational matrix method for solving variable‐order fractional partial integro‐differential equations

In this paper, a new computational scheme based on operational matrices (OMs) of two‐dimensional wavelets is proposed for the solution of variable‐order (VO) fractional partial integro‐differential equations (PIDEs). To accomplish this method, first OMs of integration and VO fractional derivative (FD) have been derived using two‐dimensional Legendre wavelets. By implementing two‐dimensional wavelets approximations and the OMs of integration and variable‐order fractional derivative (VO‐FD) along with collocation points, the VO fractional partial PIDEs are reduced into the system of algebraic equations. In addition to this, some useful theorems are discussed to establish the convergence analysis and error estimate of the proposed numerical technique. Furthermore, computational efficiency and applicability are examined through some illustrative examples.     

Numerical solution of nonlinear partial quadratic integro‐differential equations of fractional order via hybrid of block‐pulse and parabolic functions

In this paper, an effective numerical approach based on a new two‐dimensional hybrid of parabolic and block‐pulse functions (2D‐PBPFs) is presented for solving nonlinear partial quadratic integro‐differential equations of fractional order. Our approach is based on 2D‐PBPFs operational matrix method together with the fractional integral operator, described in the Riemann–Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations, which greatly simplifies the problem. By using Newton's iterative method, this system is solved, and the solution of fractional nonlinear partial quadratic integro‐differential equations is achieved. Convergence analysis and an error estimate associated with the proposed method is obtained, and it is proved that the numerical convergence order of the suggested numerical method is O(h³) . The validity and applicability of the method are demonstrated by solving three numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the exact solutions much easier.     

On the applicability of Genocchi wavelet method for different kinds of fractional‐order differential equations with delay

A novel collocation method based on Genocchi wavelet is presented for the numerical solution of fractional differential equations and time‐fractional partial differential equations with delay. In this work, to achieve the approximate solution with height accuracy, we employed the operational matrix of integer derivative and the pseudo‐operational matrix of fractional derivative in Caputo sense. Also, based on Genocchi function properties, we presented delay and pantograph operational matrices of Genocchi wavelet functions (GWFs). Due to operational and pseudo‐operational matrices, the equations under this study can be turned into nonlinear algebraic equations with the unknown GWF coefficients. For illustrating the upper bound of error for the proposed method, we estimate the error in the sense of Sobolev space. In addition, to demonstrate the efficacy of the pseudo‐operational matrix of fractional derivative, we investigate the upper bound of error for the mentioned matrix. Finally, the algorithm based on the proposed approach is implemented for some numerical experiments to confirm accuracy and applicability.     

Numerical solution of two-dimensional time fractional cable equation with Mittag-Leffler kernel

The main motive of this article is to study the recently developed Atangana-Baleanu Caputo (ABC) fractional operator that is obtained by replacing the classical singular kernel by Mittag-Leffler kernel in the definition of the fractional differential operator. We investigate a novel numerical method for the nonlinear two-dimensional cable equation in which time-fractional derivative is of Mittag-Leffler kernel type. First, we derive an approximation formula of the fractional-order ABC derivative of a function t^k using a numerical integration scheme. Using this approximation formula and some properties of shifted Legendre polynomials, we derived the operational matrix of ABC derivative. In the author of knowledge, this operational matrix of ABC derivative is derived the first time. We have shown the efficiency of this newly derived operational matrix by taking one example. Then we solved a new class of fractional partial differential equations (FPDEs) by the implementation of this ABC operational matrix. The two-dimensional model of the time-fractional model of the cable equation is solved and investigated by this method. We have shown the effectiveness and validity of our proposed method by giving the solution of some numerical examples of the two-dimensional fractional cable equation. We compare our obtained numerical results with the analytical results, and we conclude that our proposed numerical method is feasible and the accuracy can be seen by error tables. We see that the accuracy is so good. This method will be very useful to investigate a different type of model that have Mittag-Leffler fractional derivative.     

The comparison of two reliable methods for accurate solution of time‐fractional Kaup–Kupershmidt equation arising in capillary gravity waves

In this paper, we compared two different methods, one numerical technique, viz Legendre multiwavelet method, and the other analytical technique, viz optimal homotopy asymptotic method (OHAM), for solving fractional‐order Kaup–Kupershmidt (KK) equation. Two‐dimensional Legendre multiwavelet expansion together with operational matrices of fractional integration and derivative of wavelet functions is used to compute the numerical solution of nonlinear time‐fractional KK equation. The approximate solutions of time fractional Kaup–Kupershmidt equation thus obtained by Legendre multiwavelet method are compared with the exact solutions as well as with OHAM. The present numerical scheme is quite simple, effective, and expedient for obtaining numerical solution of fractional KK equation in comparison to analytical approach of OHAM. Copyright © 2016 John Wiley & Sons, Ltd.     

基于拟Legendre多项式求解一类分数阶微分方程

本文基于移位的Legendre多项式构造一类新的正交拟Legendre多项式求解一类分数阶微分方程.用阶数随所求未知函数的微分的阶数而变化的拟Legendre多项式逼近未知函数;利用分数阶积分的性质推导拟Legendre多项式的积分算子阵,结合算子矩阵的思想和Tau方法,将问题转化为求解代数方程组的问题.最后,给出数值算例证明该方法的有效性.     

一类变系数分数阶微分方程的数值解法

本文研究了一类变系数分数阶微分方程的数值解法问题. 利用Cheyshev小波推导出的分数阶微分方程的算子矩阵把分数阶微分方程转换为代数方程组. 同时给出了Cheyshev小波基的收敛性和误差估计表达式, 并给出数值算例说明所提方法的精确性和有效性     

The transport dynamics in complex systems governing by anomalous diffusion modelled with Riesz fractional partial differential equations

In this paper, numerical solutions of fractional Fokker–Planck equations with Riesz space fractional derivatives have been developed. Here, the fractional Fokker–Planck equations have been considered in a finite domain. In order to deal with the Riesz fractional derivative operator, shifted Grünwald approximation and fractional centred difference approaches have been used. The explicit finite difference method and Crank–Nicolson implicit method have been applied to obtain the numerical solutions of fractional diffusion equation and fractional Fokker–Planck equations, respectively. Numerical results are presented to demonstrate the accuracy and effectiveness of the proposed numerical solution techniques. Copyright © 2016 John Wiley & Sons, Ltd.     

Fractional-Order Legendre Functions and Their Application to Solve Fractional Optimal Control of Systems Described by Integro-differential Equations

In this paper, we introduce a set of functions called fractional-order Legendre functions (FLFs) to obtain the numerical solution of optimal control problems subject to the linear and nonlinear fractional integro-differential equations. We consider the properties of these functions to construct the operational matrix of the fractional integration. Also, we achieved a general formulation for operational matrix of multiplication of these functions to solve the nonlinear problems for the first time. Then by using these matrices the mentioned fractional optimal control problem is reduced to a system of algebraic equations. In fact the functions of the problem are approximated by fractional-order Legendre functions with unknown coefficients in the constraint equations, performance index and conditions. Thus, a fractional optimal control problem converts to an optimization problem, which can then be solved numerically. The convergence of the method is discussed and finally, some numerical examples are presented to show the efficiency and accuracy of the method.     

A new Jacobi operational matrix: An application for solving fractional differential equations

In this paper, we derived the shifted Jacobi operational matrix (JOM) of fractional derivatives which is applied together with spectral tau method for numerical solution of general linear multi-term fractional differential equations (FDEs). A new approach implementing shifted Jacobi operational matrix in combination with the shifted Jacobi collocation technique is introduced for the numerical solution of nonlinear multi-term FDEs. The main characteristic behind this approach is that it reduces such problems to those of solving a system of algebraic equations which greatly simplifying the problem. The proposed methods are applied for solving linear and nonlinear multi-term FDEs subject to initial or boundary conditions, and the exact solutions are obtained for some tested problems. Special attention is given to the comparison of the numerical results obtained by the new algorithm with those found by other known methods.