A numerical technique based on B-spline for a class of time-fractional diffusion equation

This paper presents an efficient numerical technique for solving a class of time-fractional diffusion equation. The time-fractional derivative is described in the Caputo form. The L1 scheme is used for discretization of Caputo fractional derivative and a collocation approach based on sextic B-spline basis function is employed for discretization of space variable. The unconditional stability of the fully-discrete scheme is analyzed. Two numerical examples are considered to demonstrate the accuracy and applicability of our scheme. The proposed scheme is shown to be sixth order accuracy with respect to space variable and (2 − α)-th order accuracy with respect to time variable, where α is the order of temporal fractional derivative. The numerical results obtained are compared with other existing numerical methods to justify the advantage of present method. The CPU time for the proposed scheme is provided.     

A new numerical approach for single rational soliton solution of Chen‐Lee‐Liu equation with Riesz fractional derivative in optical fibers

In this paper, time‐splitting spectral approximation technique has been proposed for Chen‐Lee‐Liu (CLL) equation involving Riesz fractional derivative. The proposed numerical technique is efficient, unconditionally stable, and of second‐order accuracy in time and of spectral accuracy in space. Moreover, it conserves the total density in the discretized level. In order to examine the results, with the aid of weighted shifted Grünwald‐Letnikov formula for approximating Riesz fractional derivative, Crank‐Nicolson weighted and shifted Grünwald difference (CN‐WSGD) method has been applied for Riesz fractional CLL equation. The comparison of results reveals that the proposed time‐splitting spectral method is very effective and simple for obtaining single soliton numerical solution of Riesz fractional CLL equation.     

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

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

A Fourier spectral method for fractional-in-space Cahn–Hilliard equation

In this paper, a fractional extension of the Cahn–Hilliard (CH) phase field model is proposed, i.e. the fractional-in-space CH equation. The fractional order controls the thickness and the lifetime of the interface, which is typically diffusive in integer order case. An unconditionally energy stable Fourier spectral scheme is developed to solve the fractional equation with periodic or Neumann boundary conditions. This method is of spectral accuracy in space and of second-order accuracy in time. The main advantages of this method are that it yields high precision and high efficiency. Moreover, an extra stabilizing term is added to obey the energy decay property while maintaining accuracy and simplicity. Numerical experiments are presented to confirm the accuracy and effectiveness of the proposed method.     

A finite difference scheme for the nonlinear time-fractional partial integro-differential equation

In this paper, a finite difference scheme is proposed for solving the nonlinear time-fractional integro-differential equation. This model involves two nonlocal terms in time, ie, a Caputo time-fractional derivative and an integral term with memory. The existence of numerical solutions is shown by the Leray-Schauder theorem. And we obtain the discrete L₂ stability and convergence with second order in time and space by the discrete energy method. Then the uniqueness of numerical solutions is derived. Moreover, an iterative algorithm is designed for solving the derived nonlinear system. Numerical examples are presented to validate the theoretical findings and the efficiency of the proposed algorithm.     

带非线性源项的双侧空间分数阶扩散方程的隐式中点方法

本文用隐式中点方法离散一阶时间偏导数,并用拟紧差分算子逼近Riemann-Liouville空间分数阶偏导数,构造了求解带非线性源项的空间分数阶扩散方程的数值格式.给出了数值方法的稳定性和收敛性分析.数值试验表明数值方法是有效的.     

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.     

分数阶Cahn-Hilliard方程的高效数值算法

给出了时空分数阶Cahn-Hilliard方程的一个高效数值算法.首先,利用Laplace变换将时空分数阶Cahn-Hilliard方程转化为空间分数阶Cahn-Hilliard方程;然后,结合Fourier谱方法和有限差分法得到一个时间二阶、空间谱精度的高效数值格式;最后,通过数值实验验证本文数值算法的有效性,并验证其满足能量耗散性质和质量守恒定律.     

An unstructured mesh finite element method for solving the multi-term time fractional and Riesz space distributed-order wave equation on an irregular convex domain

In this paper, the numerical analysis for a multi-term time fracstional and Riesz space distributed-order wave equation is discussed on an irregular convex domain. Firstly, the equation is transformed into a multi-term time-space fractional wave equation using the mid-point quadrature rule to approximate the distributed-order Riesz space derivative. Next, the equation is solved by discretising in time using a Crank–Nicolson scheme and in space using the finite element method (FEM) with an unstructured mesh, respectively. Furthermore, stability and convergence are investigated by introducing some important lemmas on irregular convex domain. Finally, some examples are provided to show the effectiveness and correctness of the proposed numerical method.     

Finite difference spectral approximation for the time‐space fractional telegraph equation and its parameter estimation

The object of this paper is to present the numerical solution of the time‐space fractional telegraph equation. The proposed method is based on the finite difference scheme in temporal direction and Fourier spectral method in spatial direction. The fast Fourier transform (FFT) technique is applied to practical computation. The stability and convergence analysis are strictly proven, which shows that this method is stable and convergent with (2?α) order accuracy in time and spectral accuracy in space. Moreover, the Levenberg‐Marquardt (L‐M) iterative method is employed for the parameter estimation. Finally, some numerical examples are given to confirm the theoretical analysis.     

Numerical Solution of the Two-Sided Space–Time Fractional Telegraph Equation Via Chebyshev Tau Approximation

The operational matrices of left Caputo fractional derivative, right Caputo fractional derivative, and Riemann–Liouville fractional integral, for shifted Chebyshev polynomials, are presented and derived. We propose an accurate and efficient spectral algorithm for the numerical solution of the two-sided space–time Caputo fractional-order telegraph equation with three types of non-homogeneous boundary conditions, namely, Dirichlet, Robin, and non-local conditions. The proposed algorithm is based on shifted Chebyshev tau technique combined with the derived shifted Chebyshev operational matrices. We focus primarily on implementing the novel algorithm both in temporal and spatial discretizations. This algorithm reduces the problem to a system of algebraic equations greatly simplifying the problem. This system can be solved by any standard iteration method. For confirming the efficiency and accuracy of the proposed scheme, we introduce some numerical examples with their approximate solutions and compare our results with those achieved using other methods.     

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.     

The Numerical Simulation of Space-Time Variable Fractional Order Diffusion Equation

Many physical processes appear to exhibit fractional order behavior that may vary with time or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. Numerical methods and analysis of stability and convergence of numerical scheme for the variable fractional order partial differential equations are quite limited and difficult to derive. This motivates us to develop efficient numerical methods as well as stability and convergence of the implicit numerical methods for the space-time variable fractional order diffusion equation on a finite domain. It is worth mentioning that here we use the Coimbra-definition variable time fractional derivative which is more efficient from the numerical standpoint and is preferable for modeling dynamical systems. An implicit Euler approximation is proposed and then the stability and convergence of the numerical scheme are investigated. Finally, numerical examples are provided to show that the implicit Euler approximation is computationally efficient.     

A Linearized Spectral-Galerkin Method for Three-Dimensional Riesz-Like Space Fractional Nonlinear Coupled Reaction-Diffusion Equations

In this paper, we establish a novel fractional model arising in the chemical reaction and develop an efficient spectral method for the three-dimensional Riesz-like space fractional nonlinear coupled reaction-diffusion equations. Based on the backward difference method for time stepping and the Legendre-Galerkin spectral method for space discretization, we construct a fully discrete numerical scheme which leads to a linear algebraic system. Then a direct method based on the matrix diagonalization approach is proposed to solve the linear algebraic system, where the cost of the algorithm is of a small multiple of $N^4$ ($N$ is the polynomial degree in each spatial coordinate) flops for each time level. In addition, the stability and convergence analysis are rigorously established. We obtain the optimal error estimate in space, and the results also show that the fully discrete scheme is unconditionally stable and convergent of order one in time. Furthermore, numerical experiments are presented to confirm the theoretical claims. As the applications of the proposed method, the fractional Gray-Scott model is solved to capture the pattern formation with an analysis of the properties of the fractional powers.     

一类非线性Schrdinger方程的守恒差分法与Fourier谱方法

考察了一类带导数项的非线性Schrdinger方程的周期边值问题,提出了一种守恒的差分格式,在空间方向上采用Fourier谱方法,证明了格式的稳定性和收敛性.数值试验得到了与理论分析一致的结果.     

一类非线性Schr(o)dinger方程的守恒差分法与Fourier谱方法

考察了一类带导数项的非线性Schrodinger方程的周期边值问题，提出了一种守恒的差分格式，在空间方向上采用Fourier谱方法，证明了格式的稳定性和收敛性．数值试验得到了与理论分析一致的结果．     

A high‐order B‐spline collocation scheme for solving a nonhomogeneous time‐fractional diffusion equation

A high‐accuracy numerical approach for a nonhomogeneous time‐fractional diffusion equation with Neumann and Dirichlet boundary conditions is described in this paper. The time‐fractional derivative is described in the sense of Riemann‐Liouville and discretized by the backward Euler scheme. A fourth‐order optimal cubic B‐spline collocation (OCBSC) method is used to discretize the space variable. The stability analysis with respect to time discretization is carried out, and it is shown that the method is unconditionally stable. Convergence analysis of the method is performed. Two numerical examples are considered to demonstrate the performance of the method and validate the theoretical results. It is shown that the proposed method is of order O(Δx⁴ + Δt^{2 ? α}) convergence, where α ∈ (0,1) . Moreover, the impact of fractional‐order derivative on the solution profile is investigated. Numerical results obtained by the present method are compared with those obtained by the method based on standard cubic B‐spline collocation method. The CPU time for present numerical method and the method based on cubic B‐spline collocation method are provided.     

A fast numerical solution of the 3D nonlinear tempered fractional integrodifferential equation

In this paper, we investigate the numerical solution of the three-dimensional (3D) nonlinear tempered fractional integrodifferential equation which is subject to the initial and boundary conditions. The backward Euler (BE) method in association with the first-order convolution quadrature rule is employed to discretize this equation for time, and the Galerkin finite element method is applied for space, which is combined with an alternating direction implicit (ADI) algorithm, in order to reduce the computational cost for solving the three-dimensional nonlocal problem. Then a fully discrete BE ADI Galerkin finite element scheme can be obtained by linearizing the non-linear term. Thereafter we prove a positive-type lemma, from which the stability and convergence of the proposed numerical scheme are derived based on the energy method. Numerical experiments are performed to verify the effectiveness of the proposed approach.     

Finite difference/$H^1$-Galerkin MFE procedure for a fractional water wave model

In this article, an $H^1$-Galerkin mixed finite element (MFE) method for solving the time fractional water wave model is presented. First-order backward Euler difference method and $L1$ formula are applied to approximate integer derivative and Caputo fractional derivative with order $1/2$, respectively, and $H^1$-Galerkin mixed finite element method is used to approximate the spatial direction. The analysis of stability for fully discrete mixed finite element scheme is made and the optimal space-time orders of convergence for two unknown variables in both $H^1$-norm and $L^2$-norm are derived. Further, some computing results for a priori analysis and numerical figures based on four changed parameters in the studied problem are given to illustrate the effectiveness of the current method     

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.