Temporal Second-order Scheme for a Hidden-memory Variable Order Time Fractional Diffusion Equation with an Initial Singularity

In this work,a novel time-stepping L1 formula is developed for a hidden-memory variable-order Caputo’s fractional derivative with an initial singularity.This formula can obtain second-order accuracy and an error estimate is analyzed strictly.As an application,a fully discrete difference scheme is established for the initial-boundary value problem of a hidden-memory variable-order time fractional diffusion model.Numerical experiments are provided to support our theoretical results.

A LOCAL DISCONTINUOUS GALERKIN METHOD FOR TIME-FRACTIONAL DIFFUSION EQUATIONS

In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finite difference method with an accuracy of order3-α,and the space discretization is based on the LDG method.For the finite difference method,we summarize and supplement some previous work by others,and apply it to the analysis of the convergence and stability of the proposed scheme.The optimal error estimate is obtained in the L2norm,indicating that the scheme has temporal(3-α)th-order accuracy and spatial(k+1)th-order accuracy,where k denotes the highest degree of a piecewise polynomial in discontinuous finite element space.The numerical results are also provided to verify the accuracy and efficiency of the considered scheme.

Numerical analysis for viscoelastic fluid flow with distributed/variable order time fractional Maxwell constitutive models

Fractional calculus has been widely used to study the flow of viscoelastic fluids recently,and fractional differential equations have attracted a lot of attention.However,the research has shown that the fractional equation with constant order operators has certain limitations in characterizing some physical phenomena.In this paper,the viscoelastic fluid flow of generalized Maxwell fluids in an infinite straight pipe driven by a periodic pressure gradient is investigated systematically.Consider the complexity of the material structure and multi-scale effects in the viscoelastic fluid flow.The modified time fractional Maxwell models and the corresponding governing equations with distributed/variable order time fractional derivatives are proposed.Based on the L1-approximation formula of Caputo fractional derivatives,the implicit finite difference schemes for the distributed/variable order time fractional governing equations are presented,and the numerical solutions are derived.In order to test the correctness and availability of numerical schemes,two numerical examples are established to give the exact solutions.The comparisons between the numerical solutions and the exact solutions have been made,and their high consistency indicates that the present numerical methods are effective.Then,this paper analyzes the velocity distributions of the distributed/variable order fractional Maxwell governing equations under specific conditions,and discusses the effects of the weight coefficient(α)in distributed order time fractional derivatives,the orderα(r,t)in variable fractional order derivatives,the relaxation timeλ,and the frequencyωof the periodic pressure gradient on the fluid flow velocity.Finally,the flow rates of the distributed/variable order fractional Maxwell governing equations are also studied.

Solving the Nonlinear Variable Order Fractional Differential Equations by Using Euler Wavelets

An Euler wavelets method is proposed to solve a class of nonlinear variable order fractional differential equations in this paper.The properties of Euler wavelets and their operational matrix together with a family of piecewise functions are first presented.Then they are utilized to reduce the problem to the solution of a nonlinear system of algebraic equations.And the convergence of the Euler wavelets basis is given.The method is computationally attractive and some numerical examples are provided to illustrate its high accuracy.

A meshless method based on moving Kriging interpolation for a two-dimensional time-fractional diffusion equation

Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.

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.

A modified Tikhonov regularization method for a Cauchy problem of a time fractional diffusion equation

In this paper,we consider a Cauchy problem of the time fractional diffusion equation(TFDE)in x∈[0,L].This problem is ubiquitous in science and engineering applications.The illposedness of the Cauchy problem is explained by its solution in frequency domain.Furthermore,the problem is formulated into a minimization problem with a modified Tikhonov regularization method.The gradient of the regularization functional based on an adjoint problem is deduced and the standard conjugate gradient method is presented for solving the minimization problem.The error estimates for the regularized solutions are obtained under Hp norm priori bound assumptions.Finally,numerical examples illustrate the effectiveness of the proposed method.

NUMERICAL SIMULATIONS FOR A VARIABLE ORDER FRACTIONAL CABLE EQUATION

In this article, Crank-Nicolson method is used to study the variable order fractional cable equation. The variable order fractional derivatives are described in the Riemann- Liouville and the Griinwald-Letnikov sense. The stability analysis of the proposed technique is discussed. Numerical results are provided and compared with exact solutions to show the accuracy of the proposed technique.

Alternating Segment Explicit-Implicit and Implicit-Explicit Parallel Difference Method for Time Fractional Sub-Diffusion Equation

The fractional diffusion equations can accurately describe the migration process of anomalous diffusion, which are widely applied in the field of natural science and engineering calculations. This paper proposed a kind of numerical methods with parallel nature which were the alternating segment explicit-implicit (ASE-I) and implicit-explicit (ASI-E) difference method for the time fractional sub-diffusion equation. It is based on the combination of the explicit scheme, implicit scheme, improved Saul’yev asymmetric scheme and the alternating segment technique. Theoretical analyses have shown that the solution of ASE-I (ASI-E) scheme is uniquely solvable. At the same time the stability and convergence of the two schemes were proved by the mathematical induction. The theoretical analyses are verified by numerical experiments. Meanwhile the ASE-I (ASI-E) scheme has the higher computational efficiency compared with the implicit scheme. Therefore it is feasible to use the parallel difference schemes for solving the time fractional diffusion equation.

A Class of Semi-Implicit Parallel Difference Method for Time Fractional Diffusion Equations

In this paper, we construct a class of semi-implicit difference method for time fractional diffusion equations—the group explicit (GE) difference scheme, which is a difference scheme with good parallelism constructed using Saul’yev asymmetric scheme. The stability and convergence of the GE scheme of time fractional diffusion equation are analyzed by mathematical induction. Then, the theoretical analysis is verified by numerical experiments, which shows that the GE scheme is effective for solving the time fractional diffusion equation.

Space-Fractional Diffusion with Variable Order and Diffusivity:Discretization and Direct Solution Strategies

We consider the multidimensional space-fractional diffusion equations with spatially varying diffusivity and fractional order.Significant computational challenges are encoun-tered when solving these equations due to the kernel singularity in the fractional integral operator and the resulting dense discretized operators,which quickly become prohibitively expensive to handle because of their memory and arithmetic complexities.In this work,we present a singularity-aware discretization scheme that regularizes the singular integrals through a singularity subtraction technique adapted to the spatial variability of diffusiv-ity and fractional order.This regularization strategy is conveniently formulated as a sparse matrix correction that is added to the dense operator,and is applicable to different formula-tions of fractional diffusion equations.We also present a block low rank representation to handle the dense matrix representations,by exploiting the ability to approximate blocks of the resulting formally dense matrix by low rank factorizations.A Cholesky factorization solver operates directly on this representation using the low rank blocks as its atomic com-putational tiles,and achieves high performance on multicore hardware.Numerical results show that the singularity treatment is robust,substantially reduces discretization errors,and attains the first-order convergence rate allowed by the regularity of the solutions.They also show that considerable savings are obtained in storage(O(N^(1.5)))and computational cost(O(N^(2)))compared to dense factorizations.This translates to orders-of-magnitude savings in memory and time on multidimensional problems,and shows that the proposed methods offer practical tools for tackling large nonlocal fractional diffusion simulations.

Green Function of Generalized Time Fractional Diffusion Equation Using Addition Formula of Mittag-Leffler Function

In this paper, we use the Mittag-Leffler addition formula to solve the Green function of generalized time fractional diffusion equation in the whole plane and prove the convergence of the Green function.

Variational iteration method for solving time-fractional diffusion equations in porous the medium

The variational iteration method is successfully extended to the case of solving fractional differential equations, and the Lagrange multiplier of the method is identified in a more accurate way. Some diffusion models with fractional derivatives are investigated analytically, and the results show the efficiency of the new Lagrange multiplier for fractional differential equations of arbitrary order.

Multiple Solutions for a Class of Variable-Order Fractional Laplacian Equations with Concave-Convex Nonlinearity

This paper is concerned with the following variable-order fractional Laplacian equations , where N ≥ 1 and N > 2s(x,y) for (x,y) ∈ Ω × Ω, Ω is a bounded domain in RN, s(⋅) ∈ C (RN × RN, (0,1)), (-Δ)s(⋅) is the variable-order fractional Laplacian operator, λ, μ > 0 are two parameters, V: Ω → [0, ∞) is a continuous function, f ∈ C(Ω × R) and q ∈ C(Ω). Under some suitable conditions on f, we obtain two solutions for this problem by employing the mountain pass theorem and Ekeland’s variational principle. Our result generalizes the related ones in the literature.

A Compact Finite Volume Scheme for the Multi-Term Time Fractional Sub-Diffusion Equation

In this paper, we introduce high-order finite volume methods for the multi-term time fractional sub-diffusion equation. The time fractional derivatives are described in Caputo’s sense. By using some operators, we obtain the compact finite volume scheme have high order accuracy. We use a compact operator to deal with spatial direction;then we can get the compact finite volume scheme. It is proved that the finite volume scheme is unconditionally stable and convergent in L∞-norm. The convergence order is O(τ2-α + h4). Finally, two numerical examples are given to confirm the theoretical results. Some tables listed also can explain the stability and convergence of the scheme.

Enriched reproducing kernel particle method for fractional advection–diffusion equation

The reproducing kernel particle method （RKPM） has been efficiently applied to problems with large deformations, high gradients and high modal density. In this paper, it is extended to solve a nonlocal problem modeled by a fractional advectiondiffusion equation （FADE）, which exhibits a boundary layer with low regularity. We formulate this method on a moving least-square approach. Via the enrichment of fractional-order power functions to the traditional integer-order basis for RKPM, leading terms of the solution to the FADE can be exactly reproduced, which guarantees a good approximation to the boundary layer. Numerical tests are performed to verify the proposed approach.

A DIRECT DISCONTINUOUS GALERKIN METHOD FOR TIME FRACTIONAL DIFFUSION EQUATIONS WITH FRACTIONAL DYNAMIC BOUNDARY CONDITIONS

This paper deals with the numerical approximation for the time fractional diffusion problem with fractional dynamic boundary conditions.The well-posedness for the weak solutions is studied.A direct discontinuous Galerkin approach is used in spatial direction under the uniform meshes,together with a second-order Alikhanov scheme is utilized in temporal direction on the graded mesh,and then the fully discrete scheme is constructed.Furthermore,the stability and the error estimate for the full scheme are analyzed in detail.Numerical experiments are also given to illustrate the effectiveness of the proposed method.

Fractal and Fractional Diffusion Equations of Price Changing of Commodity

In this paper, three types of modeling of diffusion equations for price changing of commodity are studied. In which, the partial derivatives of price of commodity respected to time on the left hand side are integer-derivative, fractal derivative, and fractional derivative respectively;while just a second order derivative respected to space is considered on the right hand side. The solutions of these diffusion equations are obtained by method of departing variables and initial boundary conditions, by translation of variables, and by translation of operators. The definitions of order of commodity x and the distance between commodity?xi and xj are defined as [1]. Examples of calculation of price of pork, beef and mutton mainly due to price raising of pork in 2007-07 to 2008-02 inChina are given with same market data as [1]. Conclusion is made.

Local Discontinuous Galerkin Methods with Novel Basis for Fractional Diffusion Equations with Non-smooth Solutions

In this paper,we develop novel local discontinuous Galerkin(LDG)methods for fractional diffusion equations with non-smooth solutions.We consider such problems,for which the solutions are not smooth at boundary,and therefore the traditional LDG methods with piecewise polynomial solutions suffer accuracy degeneracy.The novel LDG methods utilize a solution information enriched basis,simulate the problem on a paired special mesh,and achieve optimal order of accuracy.We analyze the L2 stability and optimal error estimate in L2-norm.Finally,numerical examples are presented for validating the theoretical conclusions.

A specific state variable for a class of 3D continuous fractional-order chaotic systems

A specific state variable in a class of 3D continuous fractional-order chaotic systems is presented. All state variables of fractional-order chaotic systems of this class can be obtained via a specific state variable and its （q-order and 2q-order） time derivatives. This idea is demonstrated by using several well-known fractional-order chaotic systems. Finally, a synchronization scheme is investigated for this fractional-order chaotic system via a specific state variable and its （q-order and 2q-order） time derivatives. Some examples are used to illustrate the effectiveness of the proposed synchronization method.