Original articlesLocal discontinuous Galerkin method for multi-term variable-order time fractional diffusion equation
Introduction
Fractional calculus which involves derivatives and integrals of any order can be applied in many fields to describe different problems. It has received extensive attention from various disciplines, such as physics, engineering science, and mathematics. It was produced due to discussions between L’Hopital and Leibniz, and was promoted by famous scientists such as Euler, Abel, and Laplace [29]. In recent years, many scholars have also paid attention to fractional calculus and used it to describe different types of problems. Fang [6] reduced the number of required parameters by using the fractional calculus to the viscoelastic constitutive model. Chen [4] applied fractional calculus to the supply chain financial system and established a three-dimensional supply chain game model which is closer to reality. Lu et al. [19] analyzed complex biological systems and modeled artificial neural networks by it. There are many other applications, such as diffusion in porous media [20], pollutant transport kinetics [36], and so on.
Samko [23] studied the fractional integration and differentiation of variable-order, and discussed its mathematical analysis. The later research results have shown that many complex physical phenomena can be better described by using variable-order fractional differential equations, such as signature verification, algebraic structure, viscoelastic materials, and noise reduction [7], [9], [24], [41]. Variable order fractional derivative which changes with spatial and temporal variables is a good tool to eliminate the nonphysical singularity of the solutions for constant-order FPDEs and model multi-physics phenomena, and has many advantages in characterizing memory property of systems [25], [26], [27], [30].
Numerical methods play an important role in studying partial differential equations with variable order derivatives. Langlands and Henry proposed a finite-difference scheme for fractional-order diffusion equations and studied the accuracy and stability of the scheme [11]. Cao and his collaborators studied a compact finite-difference scheme for variable-order diffusion equations [1]. Li et al. solved the nonlinear time fractional parabolic problem using the L1-discontinuous finite element method and obtained the optimal error estimate [13]. Du et al. [5] studied the implicit scheme and Crank–Nicolson scheme for the variable-order time fractional diffusion equation over the finite domain. Haq et al. [10] used finite difference method to solve the variable order time fraction (1+1)-dimensional and (1+2)-dimensional convection–diffusion equations. Chen and Liu [3] proposed a fully discrete dual-grid modified characteristic curve method based on two-dimensional nonlinear variable-order time fractional convection–diffusion equations. Sadri and Aminikhah [22] studied an algorithm based on Chebyshev polynomials for varied-order time fractional diffusion wave equations. Chen [2] solved the variable-order anomalous subdiffusion equation with first order temporal accuracy and fourth order spatial accuracy. Lin and Liu [16] studied the explicit finite difference method for variable-order nonlinear fractional differential equations. Wang and Zheng [30] proved the wellposedness of a nonlinear variable-order fractional differential equation and proposed a finite difference scheme for the problem.
Although there are some numerical algorithms for variable-order partial differential equations, it is still challenging to construct high-order numerical schemes for the models. The discontinuous Galerkin method is a hybrid of the finite element and finite volume methods, with many of the benefits of each. The resulting equations are local to the produced elements, which is a significant distinction between the DG method and the traditional finite element method. Without looking at adjoining elements, the solution can be reconstructed in each element. The discontinuous Galerkin method is a very attractive tool for solving partial differential equations [38] and has been proposed to numerically simulate fractional differential equations [8], [14], [15], [17], [18], [31], [34] due to its flexibility and efficiency in terms of meshes and shape functions. Readers may refer to [12], [21], [28], [32], [35], [37], [40] for additional information.
In this paper, we will construct a high-order local discontinuous Galerkin method for the following multi-term variable-order time-fractional diffusion equations where are the orders of the fractional derivatives, are continuous functions of in the interval , , and are smooth functions. The solution in this paper is considered to be either periodic or compactly supported.
The variable-order Caputo fractional derivative in (1.1) is defined by
A brief outline of the paper is as follows. In Section 2, some symbols and auxiliary results are described. Then in Section 3, we present the fully discrete local discontinuous Galerkin method for the multi-term variable-order time fractional diffusion equation. In Section 4, we proved that the scheme is unconditionally stable and convergent. Some numerical examples are given to demonstrate the reliability and efficiency of the method in Section 5, and the last section includes concluding comments.
Section snippets
Notations and auxiliary results
Let be a partition of , denote , for , and , .
Denote by and the values of at , from the right cell , and the left cell , respectively.
The associated discontinuous Galerkin space :
Regarding the error estimation, we will use two projections in the one-dimensional interval , denoted by , which means that for each
The schemes
Next, we introduce the numerical scheme of (1.1) solution.
We divide the interval evenly into time step size , be the mesh points.
For , we have
Stability and convergence
In order to simplify the notation we consider the case of in the numerical analysis.
Theorem 4.1 Under periodic or compactly supported boundary conditions, the fully discrete locally discontinuous Galerkin scheme (3.4) is unconditionally stable, and the numerical solution satisfies
Proof Take the test functions in the scheme (3.4), and use the fluxes (3.5), we get
Numerical experiment
In this section, some numerical examples are presented to illustrate the accuracy and performance of the method. From the results of different numerical examples, we can see that the scheme is convergent and effective.
Example 5.1 Consider the problem (1.1) In the interval , the initial condition . Let the function then the exact solution is
Conclusion
In this paper, we propose a high order numerical method to solve a class of multi-term variable-order fractional diffusion equations. By using finite difference method in time and local discontinuous Galerkin method in space, we get a fully discrete LDG scheme. By taking the appropriate numerical flux, we prove that the method is stable and convergent with . Some numerical experiments are also displayed which is consistent with the theoretical results.
Acknowledgments
This work is supported by the Training Plan of Young Backbone Teachers in Henan University of Technology, PR China (21420049), the Training Plan of Young Backbone Teachers in Colleges and Universities of Henan Province, PR China (2019GGJS094), Scientific and Technological Research Projects in Henan Province (212102210612) and the National Natural Science Foundation of China (12126325, 12126315).
References (41)
- et al.
A compact finite difference scheme for variable order subdiffusion equation
Commun. Nonlinear Sci. Numer. Simul.
(2017) - et al.
A two-grid MMOC finite element method for nonlinear variable-order time-fractional mobile/immobile advection-diffusion equations
Comput. Math. Appl.
(2020) - et al.
Combined application of blockchain technology in fractional calculus model of supply chain financial system
Chaos Solitons Fractals
(2020) - et al.
A meshless method in reproducing kernel space for solving variable-order time fractional advection-diffusion equations on arbitrary domain
Appl. Math. Lett.
(2021) - et al.
Variable-order fractional calculus: A change of perspective
Commun. Nonlinear Sci. Numer. Simul.
(2021) - et al.
On an accurate discretization of a variable-order fractional reaction–diffusion equation
Commun. Nonlinear Sci. Numer. Simul.
(2019) - et al.
Numerical solutions of variable order time fractional (1+1)-dimensional and (1+2)-dimensional advection dispersion and diffusion models
Appl. Math. Comput.
(2019) - et al.
The accuracy and stability of an implicit solution method for the fractional diffusion equation
J. Comput. Phys.
(2005) - et al.
A fast linearized conservative finite element method for the strongly coupled nonlinear fractional Schrödinger equations
J. Comput. Phys.
(2018) - et al.
The local discontinuous Galerkin finite element methods for Caputo-type partial differential equations: Numerical analysis
Appl. Numer. Math.
(2019)
The discontinuous Galerkin finite element method for Caputo-type nonlinear conservation law
Math. Comput. Simulation
Stability and convergence of a new explicit finite-difference approximation for the variable-order nonlinear fractional diffusion equation
Appl. Math. Comput.
Discontinuous Galerkin time stepping method for solving linear space fractional partial differential equations
Appl. Numer. Math.
High-order local discontinuous Galerkin method combined with WSGD-approximation for a fractional subdiffusion equation
Comput. Math. Appl.
A fractional porous medium equation
Adv. Math.
Numerical techniques for the variable order time fractional diffusion equation
Appl. Math. Comput.
An iterative method for the design of variable fractional-order FIR differintegrators
Signal Process.
A fully discrete difference scheme for a diffusion-wave system
Appl. Numer. Math.
A meshless method for solving two-dimensional variable-order time fractional advection–diffusion equation
J. Comput. Phys.
Local discontinuous Galerkin approximations to variable-order time-fractional diffusion model based on the Caputo–Fabrizio fractional derivative
Math. Comput. Simulation
Cited by (4)
A Mixed Finite Element Method for the Multi-Term Time-Fractional Reaction–Diffusion Equations
2024, Fractal and FractionalExplicit Chebyshev–Galerkin scheme for the time-fractional diffusion equation
2024, International Journal of Modern Physics CBarycentric rational interpolation method for solving time-dependent fractional convection-diffusion equation
2023, Electronic Research ArchivePetrov-Galerkin Lucas Polynomials Procedure for the Time-Fractional Diffusion Equation
2023, Contemporary Mathematics (Singapore)