A weighted finite difference method for the fractional diffusion equation based on the Riemann-Liouville derivative

Ercília Sousa,Can Li

DOI: https://doi.org/10.48550/arXiv.1109.2345

2011-09-12

Abstract:A one dimensional fractional diffusion model with the Riemann-Liouville fractional derivative is studied. First, a second order discretization for this derivative is presented and then an unconditionally stable weighted average finite difference method is derived. The stability of this scheme is established by von Neumann analysis. Some numerical results are shown, which demonstrate the efficiency and convergence of the method. Additionally, some physical properties of this fractional diffusion system are simulated, which further confirm the effectiveness of our method.

Numerical Analysis,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop an effective and stable numerical method for solving one - dimensional fractional - order diffusion equations based on the Riemann - Liouville derivative. Specifically, the authors aim to:

1. **Propose a second - order discretization method**: For the Riemann - Liouville fractional - order derivative (\(1 < \alpha < 2\)), a discretization method with second - order accuracy is proposed.

2. **Design an unconditionally stable weighted finite - difference method**: A weighted average finite - difference method is introduced, which is unconditionally stable for the time - step weight \(\tau\in[1/2, 1]\) and can achieve second - order accuracy when \(\tau = 1/2\).

3. **Verify the effectiveness and convergence of the method**: The stability and convergence of the method are proved by von Neumann stability analysis, and the efficiency and accuracy of the method are demonstrated by numerical experiments.

4. **Simulate physical properties**: The physical properties of the fractional - order diffusion system are studied by numerical simulation, further verifying the effectiveness of the proposed method.

### Formula representation

- The form of the fractional - order diffusion equation is:
\[
\frac{\partial u}{\partial t}(x,t)=d(x)\frac{\partial^\alpha u}{\partial x^\alpha}(x,t)+p(x,t)
\]
where \(1 < \alpha\leq2\), and satisfies the initial condition \(u(x,0)=f(x)\) and the boundary condition \(u(x,t) = 0\) when \(|x|\to\infty\).

- The Riemann - Liouville fractional - order derivative is defined as:
\[
\frac{\partial^\alpha u}{\partial x^\alpha}(x,t)=\frac{1}{\Gamma(n-\alpha)}\frac{\partial^n}{\partial x^n}\int_{a}^{x}u(\xi,t)(x - \xi)^{n-\alpha-1}\,d\xi
\]
where \(n - 1 < \alpha < n\), and \(\Gamma(\cdot)\) is the gamma function.

Through these works, the authors have addressed the challenges in the numerical solution of fractional - order diffusion equations, especially in the case of superdiffusion.