A Chebyshev spectral method based on operational matrix for initial and boundary value problems of fractional order

We are concerned with linear and nonlinear multi-term fractional differential equations (FDEs). The shifted Chebyshev operational matrix (COM) of fractional derivatives is derived and used together with spectral methods for solving FDEs. Our approach was based on the shifted Chebyshev tau and collocation methods. The proposed algorithms are applied to solve two types of FDEs, linear and nonlinear, subject to initial or boundary conditions, and the exact solutions are obtained for some tested problems. Numerical results with comparisons are given to confirm the reliability of the proposed method for some FDEs.     

The BEM for numerical solution of partial fractional differential equations

A numerical method is presented for the solution of partial fractional differential equations (FDEs) arising in engineering applications and in general in mathematical physics. The solution procedure applies to both linear and nonlinear problems described by evolution type equations involving fractional time derivatives in bounded domains of arbitrary shape. The method is based on the concept of the analog equation, which in conjunction with the boundary element method (BEM) enables the spatial discretization and converts a partial FDE into a system of coupled ordinary multi-term FDEs. Then this system is solved using the numerical method for the solution of such equations developed recently by Katsikadelis. The method is illustrated by solving second order partial FDEs and its efficiency and accuracy is validated.     

An efficient method based on the second kind Chebyshev wavelets for solving variable-order fractional convection diffusion equations

In this paper, a class of variable-order fractional convection diffusion equations have been solved with assistance of the second kind Chebyshev wavelets operational matrix. The operational matrix of variable-order fractional derivative is derived for the second kind Chebyshev wavelets. By implementing the second kind Chebyshev wavelets functions and also the associated operational matrix, the considered equations will be reduced to the corresponding Sylvester equation, which can be solved by some appropriate iterative solvers. Also, the convergence analysis of the proposed numerical method to the exact solutions and error estimation are given. A variety of numerical examples are considered to show the efficiency and accuracy of the presented technique.     

Numerical Solutions of Fractional Differential Equations by Using Fractional Taylor Basis

In this paper, a new numerical method for solving fractional differential equations (FDEs) is presented. The method is based upon the fractional Taylor basis approximations. The operational matrix of the fractional integration for the fractional Taylor basis is introduced. This matrix is then utilized to reduce the solution of the fractional differential equations to a system of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of this technique.      

Analysis of inhomogeneous anisotropic viscoelastic bodies described by multi-parameter fractional differential constitutive models

The response to static loads of plane inhomogeneous anisotropic bodies made of linear viscoelastic materials is investigated. Multi-parameter differential viscoelastic constitutive equations are employed, which are generalized using fractional order time derivatives. The governing equations, which are derived by considering the equilibrium of the plane body element, are two coupled linear fractional evolution partial differential equations in terms of the displacement components. Using the Analog Equation Method (AEM) in conjunction with the Boundary Element Method (BEM) these equations are transformed into a system of multi-term ordinary fractional differential equations (FDEs), which are solved using a numerical method for FDEs developed recently by Katsikadelis. Numerical examples are presented, which not only demonstrate the efficiency of the solution procedure and validate its accuracy, but also permit a better understanding of the response of plane bodies described by different viscoelastic models.     

A Chebyshev Finite Difference Method For Solving A Class Of Optimal Control Problems

A new Chebyshev finite difference method for solving class of optimal control problem is proposed. The algorithm is based on Chebyshev approximations of the derivatives arising in system dynamics. In the performance index, we use Chebyshev approximations for integration. The numerical examples illustrate the robustness, accuracy and efficiency of the proposed technique.     

Approximation of Laplace transform of fractional derivatives via Clenshaw–Curtis integration

In this paper, we approximate the Laplace transform of fractional derivatives via Clenshaw–Curtis integration. The idea of applying Chebyshev polynomial to the numerical computation of integrals is extended to Laplace transform of fractional derivatives. The numerical stability of forward recurrence relations is considered, which depends on the asymptotic behaviour of the coefficients. Error estimation for the Laplace approximation of the fractional derivatives is also considered. Finally, from the numerical examples, the method seems to be promising for approximation of the Laplace transform of fractional derivative.     

Tau approximate solution of fractional partial differential equations

In this study, we improve the algebraic formulation of the fractional partial differential equations (FPDEs) by using the matrix-vector multiplication representation of the problem. This representation allows us to investigate an operational approach of the Tau method for the numerical solution of FPDEs. We introduce a converter matrix for the construction of converted Chebyshev and Legendre polynomials which is applied in the operational approach of the Tau method. We present the advantages of using the method and compare it with several other methods. Some experiments are applied to solve FPDEs including linear and nonlinear terms. By comparing the numerical results obtained from the other methods, we demonstrate the high accuracy and efficiency of the proposed method.     

A wavelet approach for the multi-term time fractional diffusion-wave equation

In this paper, a Galerkin method based on the second kind Chebyshev wavelets (SKCWs) is established for solving the multi-term time fractional diffusion-wave equation. To do this, a new operational matrix of fractional integration for the SKCWs must be derived and in order to improve the computational efficiency, the hat functions are proposed to create a general procedure for constructing this matrix. Implementation of these wavelet basis functions and their operational matrix of fractional integration simplifies the problem under consideration to a system of linear algebraic equations, which greatly decreases the computational cost for finding an approximate solution. The main privilege of the proposed method is adjusting the initial and boundary conditions in the final system automatically. Theoretical error and convergence analysis of the SKCWs expansion approve the reliability of the approach. Also, numerical investigation reveals the applicability and accuracy of the presented method.     

The non-standard finite difference scheme for linear fractional PDEs in fluid mechanics

A non-standard finite difference scheme is developed to solve the linear partial differential equations with time- and space-fractional derivatives. The Grunwald–Letnikov method is used to approximate the fractional derivatives. Numerical illustrations that include the linear inhomogeneous time-fractional equation, linear space-fractional telegraph equation, linear inhomogeneous fractional Burgers equation and the fractional wave equation are investigated to show the pertinent features of the technique. Numerical results are presented graphically and reveal that the non-standard finite difference scheme is very effective and convenient for solving linear partial differential equations of fractional order.     

A posteriori error analysis for a fractional differential equation

Numerical treatment for a fractional differential equation (FDE) is proposed and analysed. The solution of the FDE may be singular near certain domain boundaries, which leads to numerical difficulty. We apply the upwind finite difference method to the FDE. The stability properties and a posteriori error analysis for the discrete scheme are given. Then, a posteriori adapted mesh based on a posteriori error analysis is established by equidistributing arc-length monitor function. Numerical experiments illustrate that the upwind finite difference method on a posteriori adapted mesh is more accurate than the method on uniform mesh.     

An efficient numerical technique based on the extended cubic B-spline functions for solving time fractional Black–Scholes model

Financial theory could introduce a fractional differential equation (FDE) that presents new theoretical research concepts, methods and practical implementations. Due to the memory factor of fractional derivatives, physical pathways with storage and inherited properties can be best represented by FDEs. For that purpose, reliable and effective techniques are required for solving FDEs. Our objective is to generalize the collocation method for solving time fractional Black–Scholes European option pricing model using the extended cubic B-spline. The key feature of the strategy is that it turns these type of problems into a system of algebraic equations which can be appropriate for computer programming. This is not only streamlines the problems but speed up the computations as well. The Fourier stability and convergence analysis of the scheme are examined. A proposed numerical scheme having second-order accuracy via spatial direction is also constructed. The numerical and graphical results indicate that the suggested approach for the European option prices agree well with the analytical solutions.

     

分数阶时变切换系统有限时间异步控制

针对一类时变切换系统,当考虑子系统具有分数阶(Fractional Order)特性时,提出了一种基于模型依赖平均驻留时间方法的有限时间稳定性条件及异步切换控制策略.借助Caputo分数阶导数引理和切换Lyapunov函数,利用矩阵不等式技术提出了分数阶时变切换系统有限时间稳定的充分条件.将有限时间稳定的结果进一步推广到有限时间有界的情形,利用平均驻留时间思想提出了分数阶时变切换系统有限时间有界的充分条件,基于该条件设计了系统的异步切换控制器.所给出的设计方法将系统异步切换控制问题转化为矩阵不等式组的求解问题.通过数值仿真验证了所提控制方法的有效性.     

A New Finite Element Analysis for Inhomogeneous Boundary-Value Problems of Space Fractional Differential Equations

In this paper the framework and convergence analysis of finite element methods (FEMs) for space fractional differential equations (FDEs) with inhomogeneous boundary conditions are studied. Since the traditional framework of Gakerkin methods for space FDEs with homogeneous boundary conditions is not true any more for the case of inhomogeneous boundary conditions, this paper develops a technique by introducing a new fractional derivative space in which the Galerkin method works and proves the convergence rates of the FEMs.     

The numerical solution of fractional differential equations using the Volterra integral equation method based on thin plate splines

Finding the approximate solution of differential equations, including non-integer order derivatives, is one of the most important problems in numerical fractional calculus. The main idea of the current paper is to obtain a numerical scheme for solving fractional differential equations of the second order. To handle the method, we first convert these types of differential equations to linear fractional Volterra integral equations of the second kind. Afterward, the solutions of the mentioned Volterra integral equations are estimated using the discrete collocation method together with thin plate splines as a type of free-shape parameter radial basis functions. Since the scheme does not need any background meshes, it can be recognized as a meshless method. The proposed approach has a simple and computationally attractive algorithm. Error analysis is also studied for the presented method. Finally, the reliability and efficiency of the new technique are tested over several fractional differential equations and obtained results confirm the theoretical error estimates.

     

A new operational matrix for solving fractional-order differential equations

Fractional calculus has been used to model physical and engineering processes that are found to be best described by fractional differential equations. For that reason we need a reliable and efficient technique for the solution of fractional differential equations. This paper deals with the numerical solution of a class of fractional differential equations. The fractional derivatives are described in the Caputo sense. Our main aim is to generalize the Legendre operational matrix to the fractional calculus. In this approach, a truncated Legendre series together with the Legendre operational matrix of fractional derivatives are used for numerical integration of fractional differential equations. The main characteristic behind the approach using this technique is that it reduces such problems to those of solving a system of algebraic equations thus greatly simplifying the problem. The method is applied to solve two types of fractional differential equations, linear and nonlinear. Illustrative examples are included to demonstrate the validity and applicability of the presented technique.     

Numerical solution of fractional delay differential equation by shifted Jacobi polynomials

In this paper, the fractional delay differential equation (FDDE) is considered for the purpose to develop an approximate scheme for its numerical solutions. The shifted Jacobi polynomial scheme is used to solve the results by deriving operational matrix for the fractional differentiation and integration in the Caputo and Riemann–Liouville sense, respectively. In addition to it, the Jacobi delay coefficient matrix is developed to solve the linear and nonlinear FDDE numerically. The error of the approximate solution of proposed method is discussed by applying the piecewise orthogonal technique. The applicability of this technique is shown by several examples like a mathematical model of houseflies and a model based on the effect of noise on light that reflected from laser to mirror. The obtained numerical results are tabulated and displayed graphically.     

Application of Legendre wavelets for solving fractional differential equations

In this paper, we develop a framework to obtain approximate numerical solutions to ordinary differential equations (ODEs) involving fractional order derivatives using Legendre wavelet approximations. The properties of Legendre wavelets are first presented. These properties are then utilized to reduce the fractional ordinary differential equations (FODEs) to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique. Results show that this technique can solve the linear and nonlinear fractional ordinary differential equations with negligible error compared to the exact solution.     

A finite difference/finite element technique with error estimate for space fractional tempered diffusion-wave equation

An efficient numerical technique is proposed to solve one- and two-dimensional space fractional tempered fractional diffusion-wave equations. The space fractional is based on the Riemann–Liouville fractional derivative. At first, the temporal direction is discretized using a second-order accurate difference scheme. Then a classic Galerkin finite element is employed to obtain a full-discrete scheme. Furthermore, for the time-discrete and the full-discrete schemes error estimate has been presented to show the unconditional stability and convergence of the developed numerical method. Finally, two test problems have been illustrated to verify the efficiency and simplicity of the proposed technique.