A second-order accurate numerical method for the two-dimensional fractional diffusion equation

Spatially fractional order diffusion equations are generalizations of classical diffusion equations which are used in modeling practical superdiffusive problems in fluid flow, finance and others. In this paper, we present an accurate and efficient numerical method to solve a fractional superdiffusive differential equation. This numerical method combines the alternating directions implicit (ADI) approach with a Crank–Nicolson discretization and a Richardson extrapolation to obtain an unconditionally stable second-order accurate finite difference method. The stability and the consistency of the method are established. Numerical solutions for an example super-diffusion equation with a known analytic solution are obtained and the behavior of the errors are analyzed to demonstrate the order of convergence of the method.     

Numerical Solutions of Fractional Boussinesq Equation

Based upon the Adomian decomposition method, a scheme is developed to obtain numerical solutions of a fractional Boussinesq equation with initial condition, which is introduced by replacing some order time and space derivatives by fractional derivatives. The fractional derivatives are described in the Caputo sense. So the traditional Adomian decomposition method for differential equations of integer order is directly extended to derive explicit and numerical solutions of the fractional differential equations. The solutions of our model equation are calculated in the form of convergent series with easily computable components.     

Numerical solution of the conformable space-time fractional wave equation

In this paper, an efficient numerical method is considered for solving space-time fractional wave equation. The fractional derivatives are described in the conformable sense. The method is based on shifted Chebyshev polynomials of the second kind. Unknown function is written as Chebyshev series with the N term. The space-time fractional wave equation is reduced to a system of ordinary differential equations by using the properties of Chebyshev polynomials. The finite difference method is applied to solve this system of equations. Numerical results are provided to verify the accuracy and efficiency of the proposed approach.     

Finite difference methods for two-dimensional fractional dispersion equation

Fractional order partial differential equations, as generalizations of classical integer order partial differential equations, are increasingly used to model problems in fluid flow, finance and other areas of application. In this paper we discuss a practical alternating directions implicit method to solve a class of two-dimensional initial-boundary value fractional partial differential equations with variable coefficients on a finite domain. First-order consistency, unconditional stability, and (therefore) first-order convergence of the method are proven using a novel shifted version of the classical Grünwald finite difference approximation for the fractional derivatives. A numerical example with known exact solution is also presented, and the behavior of the error is examined to verify the order of convergence.     

Analytical approximate solution for nonlinear space-time fractional Klein Gordon equation

The fractional derivatives in the sense of Caputo and the homotopy analysis method are used to construct an approximate solution for the nonlinear space-time fractional derivatives Klein-Gordon equation. The numerical results show that the approaches are easy to implement and accurate when applied to the nonlinear space-time fractional derivatives KleinGordon equation. This method introduces a promising tool for solving many space-time fractional partial differential equations. This method is efficient and powerful in solving wide classes of nonlinear evolution fractional order equations.     

Numerical approaches to fractional calculus and fractional ordinary differential equation

Nowadays, fractional calculus are used to model various different phenomena in nature, but due to the non-local property of the fractional derivative, it still remains a lot of improvements in the present numerical approaches. In this paper, some new numerical approaches based on piecewise interpolation for fractional calculus, and some new improved approaches based on the Simpson method for the fractional differential equations are proposed. We use higher order piecewise interpolation polynomial to approximate the fractional integral and fractional derivatives, and use the Simpson method to design a higher order algorithm for the fractional differential equations. Error analyses and stability analyses are also given, and the numerical results show that these constructed numerical approaches are efficient.     

Numerical simulation for two-dimensional Riesz space fractional diffusion equations with a nonlinear reaction term

Fractional differential equations have attracted considerable interest because of their ability to model anomalous transport phenomena. Space fractional diffusion equations with a nonlinear reaction term have been presented and used to model many problems of practical interest. In this paper, a two-dimensional Riesz space fractional diffusion equation with a nonlinear reaction term (2D-RSFDE-NRT) is considered. A novel alternating direction implicit method for the 2D-RSFDE-NRT with homogeneous Dirichlet boundary conditions is proposed. The stability and convergence of the alternating direction implicit method are discussed. These numerical techniques are used for simulating a two-dimensional Riesz space fractional Fitzhugh-Nagumo model. Finally, a numerical example of a two-dimensional Riesz space fractional diffusion equation with an exact solution is given. The numerical results demonstrate the effectiveness of the methods. These methods and techniques can be extended in a straightforward method to three spatial dimensions, which will be the topic of our future research.     

Numerical algorithm for the time fractional Fokker–Planck equation

Anomalous diffusion is one of the most ubiquitous phenomena in nature, and it is present in a wide variety of physical situations, for instance, transport of fluid in porous media, diffusion of plasma, diffusion at liquid surfaces, etc. The fractional approach proved to be highly effective in a rich variety of scenarios such as continuous time random walk models, generalized Langevin equations, or the generalized master equation. To investigate the subdiffusion of anomalous diffusion, it would be useful to study a time fractional Fokker–Planck equation. In this paper, firstly the time fractional, the sense of Riemann–Liouville derivative, Fokker–Planck equation is transformed into a time fractional ordinary differential equation (FODE) in the sense of Caputo derivative by discretizing the spatial derivatives and using the properties of Riemann–Liouville derivative and Caputo derivative. Then combining the predictor–corrector approach with the method of lines, the algorithm is designed for numerically solving FODE with the numerical error O(k^min{1+2α,2})+O(h²), and the corresponding stability condition is got. The effectiveness of this numerical algorithm is evaluated by comparing its numerical results for α=1.0 with the ones of directly discretizing classical Fokker–Planck equation, some numerical results for time fractional Fokker–Planck equation with several different fractional orders are demonstrated and compared with each other, moreover for α=0.8 the convergent order in space is confirmed and the numerical results with different time step sizes are shown.     

Numerical approximations for fractional diffusion equations via a Chebyshev spectral-tau method

In this paper, a class of fractional diffusion equations with variable coefficients is considered. An accurate and efficient spectral tau technique for solving the fractional diffusion equations numerically is proposed. This method is based upon Chebyshev tau approximation together with Chebyshev operational matrix of Caputo fractional differentiation. Such approach has the advantage of reducing the problem to the solution of a system of algebraic equations, which may then be solved by any standard numerical technique. We apply this general method to solve four specific examples. In each of the examples considered, the numerical results show that the proposed method is of high accuracy and is efficient for solving the time-dependent fractional diffusion equations.     

Numerical solution of fractional differential equations via a Volterra integral equation approach

The main focus of this paper is to present a numerical method for the solution of fractional differential equations. In this method, the properties of the Caputo derivative are used to reduce the given fractional differential equation into a Volterra integral equation. The entire domain is divided into several small domains, and by collocating the integral equation at two adjacent points a system of two algebraic equations in two unknowns is obtained. The method is applied to solve linear and nonlinear fractional differential equations. Also the error analysis is presented. Some examples are given and the numerical simulations are also provided to illustrate the effectiveness of the new method.     

A novel method for nonlinear fractional differential equations using symbolic computation

In this paper, a new approach, namely an ansatz method is applied to find exact solutions for nonlinear fractional differential equations in the sense of modified Riemann–Liouville derivative. Based on a nonlinear fractional complex transformation, a certain fractional partial differential equation can be turned into another ordinary differential equation of integer order. For illustrating the validity of this method, we apply it to solve the fractional-order biological population model and the space–time fractional modified equal width equation, and as a result, some dark soliton solutions for them are established.     

Numerical solution of fractional differential equations by using fractional B-spline

In this paper, we present fractional B-spline collocation method for the numerical solution of fractional differential equations. We consider this method for solving linear fractional differential equations which involve Caputo-type fractional derivatives. The numerical results demonstrate that the method is efficient and quite accurate and it requires relatively less computational work. For this reason one can conclude that this method has advantage on other methods and hence demonstrates the importance of this work.     

Modified Kudryashov method via new exact solutions for some conformable fractional differential equations arising in mathematical biology

In this study, the modified Kudryashov method is used to construct new exact solutions for some conformable fractional differential equations. By implementing the conformable fractional derivative and compatible fractional complex transforms, the fractional generalized reaction duffing (RD) model equation, the fractional biological population model and the fractional diffusion reaction (DR) equation with quadratic and cubic nonlinearity are discussed. As an outcome, some new exact solutions are formally established. All solutions have been verified back into its corresponding equation with the aid of maple package program. We assure that the employed method is simple and robust for the estimation of the new exact solutions, and practically capable for reducing the size of computational work for solving a various class of fractional differential equations arising in applied mathematics, mathematical physics and biology.     

Bright and dark soliton solutions for some nonlinear fractional differential equations

In this work, we propose a new approach, namely ansatz method, for solving fractional differential equations based on a fractional complex transform and apply it to the nonlinear partial space–time fractional modified Benjamin–Bona–Mahoney(m BBM) equation, the time fractional m Kd V equation and the nonlinear fractional Zoomeron equation which gives rise to some new exact solutions. The physical parameters in the soliton solutions: amplitude, inverse width, free parameters and velocity are obtained as functions of the dependent model coefficients. This method is suitable and more powerful for solving other kinds of nonlinear fractional PDEs arising in mathematical physics. Since the fractional derivatives are described in the modified Riemann–Liouville sense.     

Analytical Approach for Nonlinear Partial Differential Equations of Fractional Order

The purpose of the paper is to present analytical and numerical solutions of a degenerate parabolic equation with time-fractional derivatives arising in the spatial diffusion of biological populations. The homotopy-perturbation method is employed for solving this class of equations, and the time-fractional derivatives are described in the sense of Caputo. Comparisons are made with those derived by Adomian's decomposition method, revealing that the homotopy perturbation method is more accurate and convenient than the Adomian's decomposition method. Furthermore, the results reveal that the approximate solution continuously depends on the time-fractional derivative and the proposed method incorporating the Caputo derivatives is a powerful and efficient technique for solving the fractional differential equations without requiring linearization or restrictive assumptions. The basis ideas presented in the paper can be further applied to solve other similar fractional partial differential equations.     

Exact solutions for nonlinear partial fractional differential equations

In this article,we use the fractional complex transformation to convert nonlinear partial fractional differential equations to nonlinear ordinary differential equations.We use the improved(G’/G)-expansion function method to calculate the exact solutions to the time-and space-fractional derivative foam drainage equation and the time-and space-fractional derivative nonlinear KdV equation.This method is efficient and powerful for solving wide classes of nonlinear evolution fractional order equations.     

Explicit and implicit finite difference schemes for fractional Cattaneo equation

In this paper, the numerical solution of fractional (non-integer)-order Cattaneo equation for describing anomalous diffusion has been investigated. Two finite difference schemes namely an explicit predictor–corrector and totally implicit schemes have been developed. In developing each scheme, a separate formulation approach for the governing equations has been considered. The explicit predictor–corrector scheme is the fractional generalization of well-known MacCormack scheme and has been called Generalized MacCormack scheme. This scheme solves two coupled low-order equations and simultaneously computes the flux term with the main variable. Fully implicit scheme however solves a single high-order undecomposed equation. For Generalized MacCormack scheme, stability analysis has been studied through Fourier method. Through a numerical test, the experimental order of convergency of both schemes has been found. Then, the domain of applicability and some numerical properties of each scheme have been discussed.     

Finite difference approximations for the fractional advection-diffusion equation

Fractional order diffusion equations are viewed as generalizations of classical diffusion equations, treating super-diffusive flow processes. In this Letter, in order to solve the two-sided fractional advection-diffusion equation, the fractional Crank-Nicholson method (FCN) is given, which is based on shifted Grünwald-Letnikov formula. It is shown that this method is unconditionally stable, consistent and convergent. The accuracy with respect to the time step is of order ²(Δt). A numerical example is presented to confirm the conclusions.     

分数阶Oldroyd-B黏弹性Poiseuille流的Laplace数值反演分析

采用Laplace数值反演的Stehfest算法研究了分数阶Oldroyd-B粘弹性流体在两平板间非定常的Poiseuille流动问题.首先,通过数值解与近似解析解的比较验证了Stehfest算法的有效性.其次,运用Stehfest算法对平板Poiseuille流动进行了研究,揭示了分数阶黏弹性平板流的速度过冲和应力过冲现象,指出这些现象对分数导数的阶数存在明显的依赖性.同时,数值结果表明,整数阶本构方程仅仅是分数阶本构方程的特例,分数阶本构方程较整数阶本构方程具有更广泛的适用性。     

Compact and noncompact structures for nonlinear fractional evolution equations

This Letter deals with compact and noncompact solutions for nonlinear evolution equations with time-fractional derivatives. We present a reliable approach of the homotopy perturbation method to handle nonlinear fractional evolution equations. The validity of the approach is verified through illustrative examples. New exact solitary wave and compacton solutions are developed. The proposed technique could lead to a promising approach for a wide class of nonlinear fractional evolution equations.