Abstract:We consider the development and analysis of local discontinuous Galerkin methods for fractional diffusion problems, characterized by having fractional derivatives, parameterized by β ∈ [1, 2]. We show through analysis that one can construct a numerical flux which results in a scheme that exhibit optimal order of convergence O(hk+1) in the continuous range between pure advection (β = 1) and pure diffusion (β = 2). In these limits, known schemes are recovered. The analysis is confirmed through a few examples.

Methods for fractional diffusion equations