Analysis of an Implicit Finite Difference Scheme for Time Fractional Diffusion Equation

Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.     

New analysis and application of fractional order Schrödinger equation using with Atangana–Batogna numerical scheme

In this work, an analytical approximation to the solution of Schrodinger equation has been provided. The fractional derivative used in this equation is the Caputo derivative. The existence and uniqueness conditions of solutions for the proposed model are derived based on the power law. While solving the fractional order Schrodinger equation, Atangana–Batogna numerical method is presented for fractional order equation. We obtain an efficient recurrence relation for solving these kinds of equations. To illustrate the usefulness of the numerical scheme, the numerical simulations are presented. The results show that the numerical scheme is very effective and simple.     

A new reliable algorithm based on the sinc function for the time fractional diffusion equation

This paper deals with the numerical solution of time fractional diffusion equation. In this work, we consider the fractional derivative in the sense of Riemann-Liouville. At first, the time fractional derivative is discretized by integrating both sides of the equation with respect to the time variable and we arrive at a semi–discrete scheme. The stability and convergence of time discretized scheme are proven by using the energy method. Also we show that the convergence order of this scheme is O(τ2?α). Then we use the sinc collocation method to approximate the solution of semi–discrete scheme and show that the problem is reduced to a Sylvester matrix equation. Besides by performing some theorems, the exponential convergence rate of sinc method is illustrated. The numerical experiments are presented to show the excellent behavior and high accuracy of the proposed hybrid method in comparison with some other well known methods.     

Finite Difference/Collocation Method for Two-Dimensional Sub-Diffusion Equation with Generalized Time Fractional Derivative

In this paper, we propose a finite difference/collocation method for two-dimensional time fractional diffusion equation with generalized fractional operator. The main purpose of this paper is to design a high order numerical scheme for the new generalized time fractional diffusion equation. First, a finite difference approximation formula is derived for the generalized time fractional derivative, which is verified with order $2-\alpha$ $(0<\alpha<1)$. Then, collocation method is introduced for the two-dimensional space approximation. Unconditional stability of the scheme is proved. To make the method more efficient, the alternating direction implicit method is introduced to reduce the computational cost. At last, numerical experiments are carried out to verify the effectiveness of the scheme.     

Fractional Crank–Nicolson–Galerkin finite element scheme for the time‐fractional nonlinear diffusion equation

This article presents a finite element scheme with Newton's method for solving the time‐fractional nonlinear diffusion equation. For time discretization, we use the fractional Crank–Nicolson scheme based on backward Euler convolution quadrature. We discuss the existence‐uniqueness results for the fully discrete problem. A new discrete fractional Gronwall type inequality for the backward Euler convolution quadrature is established. A priori error estimate for the fully discrete problem in L²(Ω) norm is derived. Numerical results based on finite element scheme are provided to validate theoretical estimates on time‐fractional nonlinear Fisher equation and Huxley equation.     

带非线性源项的Riesz回火分数阶扩散方程的预估校正方法

本文针对带非线性源项的Riesz回火分数阶扩散方程,利用预估校正方法离散时间偏导数,并用修正的二阶Lubich回火差分算子逼近Riesz空间回火的分数阶偏导数,构造出一类新的数值格式.给出了数值格式在一定条件下的稳定性与收敛性分析,且该格式的时间与空间收敛阶均为二阶.数值试验表明数值方法是有效的.     

A novel finite volume method for the Riesz space distributed-order advection–diffusion equation

In this paper, we investigate the finite volume method (FVM) for a distributed-order space-fractional advection–diffusion (AD) equation. The mid-point quadrature rule is used to approximate the distributed-order equation by a multi-term fractional model. Next, the transformed multi-term fractional equation is solved by discretizing in space by the finite volume method and in time using the Crank–Nicolson scheme. We use a novel technique to deal with the convection term, by which the Riesz fractional derivative of order 0 < γ < 1 is transformed into a fractional integral form. An important contribution of our work is the use of nodal basis function to derive the discrete form of our model. The unique solvability of the scheme is also discussed and we prove that the Crank–Nicolson scheme is unconditionally stable and convergent with second-order accuracy. Finally, we give some examples to show the effectiveness of the numerical method.     

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.     

A finite difference approach for the initial-boundary value problem of the fractional Klein-Kramers equation in phase space

Considering the features of the fractional Klein-Kramers equation (FKKE) in phase space, only the unilateral boundary condition in position direction is needed, which is different from the bilateral boundary conditions in [Cartling B., Kinetics of activated processes from nonstationary solutions of the Fokker-Planck equation for a bistable potential, J. Chem. Phys., 1987, 87(5), 2638–2648] and [Deng W., Li C., Finite difference methods and their physical constrains for the fractional Klein-Kramers equation, Numer. Methods Partial Differential Equations, 2011, 27(6), 1561–1583]. In the paper, a finite difference scheme is constructed, where temporal fractional derivatives are approximated using L1 discretization. The advantages of the scheme are: for every temporal level it can be dealt with from one side to the other one in position direction, and for any fixed position only a tri-diagonal system of linear algebraic equations needs to be solved. The computational amount reduces compared with the ADI scheme in [Cartling B., Kinetics of activated processes from nonstationary solutions of the Fokker-Planck equation for a bistable potential, J. Chem. Phys., 1987, 87(5), 2638–2648] and the five-point scheme in [Deng W., Li C., Finite difference methods and their physical constrains for the fractional Klein-Kramers equation, Numer. Methods Partial Differential Equations, 2011, 27(6), 1561–1583]. The stability and convergence are proved and two examples are included to show the accuracy and effectiveness of the method.     

On the analysis of vibration equation involving a fractional derivative with Mittag-Leffler law

The present article deals with a fractional extension of the vibration equation for very large membranes with distinct special cases. The fractional derivative is considered in Atangana-Baleanu sense. A numerical algorithm based on homotopic technique is employed to examine the fractional vibration equation. The stability analysis is conducted for the suggested scheme. The maple software package is utilized for numerical simulation. In order to illustrate the effects of space, time, and order of Atangana-Baleanu derivative on the displacement, the outcomes of this study are demonstrated graphically. The results revel that the Atangana-Baleanu fractional derivative is very efficient in describing vibrations in large membranes.     

A second‐order finite difference method for fractional diffusion equation with Dirichlet and fractional boundary conditions

In this paper, we develop a practical numerical method to approximate a fractional diffusion equation with Dirichlet and fractional boundary conditions. An approach based on the classical Crank–Nicolson method combined with spatial extrapolation is used to obtain temporally and spatially second‐order accurate numerical estimates. The solvability, stability, and convergence of the proposed numerical scheme are proved via the Gershgorin theorem. Numerical experiments are performed to confirm the accuracy and efficiency of our scheme.     

Fast evaluation and high accuracy finite element approximation for the time fractional subdiffusion equation

In this article, an efficient algorithm for the evaluation of the Caputo fractional derivative and the superconvergence property of fully discrete finite element approximation for the time fractional subdiffusion equation are considered. First, the space semidiscrete finite element approximation scheme for the constant coefficient problem is derived and supercloseness result is proved. The time discretization is based on the L1‐type formula, whereas the space discretization is done using, the fully discrete scheme is developed. Under some regularity assumptions, the superconvergence estimate is proposed and analyzed. Then, extension to the case of variable coefficients is also discussed. To reduce the computational cost, the fast evaluation scheme of the Caputo fractional derivative to solve the fractional diffusion equations is designed. Finally, numerical experiments are presented to support the theoretical results.     

[3, 3] Padé approximation method for solving space fractional Fokker–Planck equations

The fractional Fokker–Planck equation has been used in various areas of engineering and physics. In this paper, we proposed a novel numerical scheme for solving the space fractional Fokker–Planck equation with the help of the [3, 3] Padé approximation. It is proved that the numerical method is unconditionally stable in view of the matrix analysis method. Finally, a numerical example is proposed to prove the effectiveness of the numerical scheme.     

Finite difference scheme for time-space fractional diffusion equation with fractional boundary conditions

In this paper, the finite difference scheme is developed for the time-space fractional diffusion equation with Dirichlet and fractional boundary conditions. The time and space fractional derivatives are considered in the senses of Caputo and Riemann-Liouville, respectively. The stability and convergence of the proposed numerical scheme are strictly proved, and the convergence order is O(τ^2−α+h²). Numerical experiments are performed to confirm the accuracy and efficiency of our scheme.     

Galerkin-Legendre spectral schemes for nonlinear space fractional Schrödinger equation

In the paper, we first propose a Crank-Nicolson Galerkin-Legendre (CN-GL) spectral scheme for the one-dimensional nonlinear space fractional Schrödinger equation. Convergence with spectral accuracy is proved for the spectral approximation. Further, a Crank-Nicolson ADI Galerkin-Legendre spectral method for the two-dimensional nonlinear space fractional Schrödinger equation is developed. The proposed schemes are shown to be efficient with second-order accuracy in time and spectral accuracy in space which are higher than some recently studied methods. Moreover, some numerical results are demonstrated to justify the theoretical analysis.     

New predictor‐corrector scheme for solving nonlinear differential equations with Caputo‐Fabrizio operator

In this paper, we develop a new, simple, and accurate scheme to obtain approximate solution for nonlinear differential equation in the sense of Caputo‐Fabrizio operator. To derive this new predictor‐corrector scheme, which suits on Caputo‐Fabrizio operator, firstly, we obtain the corresponding initial value problem for the differential equation in the Caputo‐Fabrizio sense. Hence, by fractional Euler method and fractional trapeziodal rule, we obtain the predictor formula as well as corrector formula. Error analysis for this new method is derived. To test the validity and simplicity of this method, some illustrative examples for nonlinear differential equations are solved.     

非光滑时间分数阶齐次扩散方程的修正

当初值不光滑时,时间分数阶齐次扩散方程数值方法的精度会下降.为了得到高阶时间收敛格式,提出加权移位的Grünwald-Letnikov的修正格式,运用Lubich的修正方法,得到非光滑时间分数阶齐次扩散方程的收敛阶仍为O(k²).最后,通过数值算例验证了数值计算结果与理论计算结果一致.     

Numerical solution of fractional oscillator equation

We focus on a numerical scheme applied for a fractional oscillator equation in a finite time interval. This type of equation includes a complex form of left- and right-sided fractional derivatives. Its analytical solution is represented by a series of left and right fractional integrals and therefore is difficult in practical calculations. Here we elaborated two numerical schemes being dependent on a fractional order of the equation. The results of numerical calculations are compared with analytical solutions. Then we illustrate convergence and stability of our schemes.     

Numerical simulation of the nonlinear generalized time‐fractional Klein–Gordon equation using cubic trigonometric B‐spline functions

In this paper, an efficient numerical procedure for the generalized nonlinear time‐fractional Klein–Gordon equation is presented. We make use of the typical finite difference schemes to approximate the Caputo time‐fractional derivative, while the spatial derivatives are discretized by means of the cubic trigonometric B‐splines. Stability and convergence analysis for the numerical scheme are discussed. We apply our scheme to some typical examples and compare the obtained results with the ones found by other numerical methods. The comparison shows that our scheme is quite accurate and can be applied successfully to a variety of problems of applied nature.     

A backward Euler alternating direction implicit difference scheme for the three‐dimensional fractional evolution equation

A backward Euler alternating direction implicit (ADI) difference scheme is formulated and analyzed for the three‐dimensional fractional evolution equation. In our method, the Riemann‐Liouville fractional integral term is treated by means of first order convolution quadrature suggested by Lubich. Meanwhile, an ADI technique is adopted to reduce the multidimensional problem to a series of one‐dimensional problems. A fully discrete difference scheme is constructed with space discretization by finite difference method. Two new inner products and corresponding norms are defined to analyze the scheme. The verification of stability and convergence is based on the nonnegative character of the real quadratic form associated with the convolution quadrature. Numerical experiments are reported to demonstrate the efficiency of our scheme.