Numerical algorithm for the space-time fractional Fokker-Planck system with two internal states

Abstract

The fractional Fokker-Planck system with multiple internal states is derived in [Xu and Deng, Math. Model. Nat. Phenom., $13$, 10 (2018)], where the space derivative is Laplace operator. If the jump length distribution of the particles is power law instead of Gaussian, the space derivative should be replaced with fractional Laplacian. This paper focuses on solving the two state Fokker-Planck system with fractional Laplacian. We first provide a priori estimate for this system under different regularity assumptions on the initial data. Then we use $L_1$ scheme to discretize the time fractional derivatives and finite element method to approximate the fractional Laplacian operators. Furthermore, we give the error estimates for the space semidiscrete and fully discrete schemes without any assumption on regularity of solutions. Finally, the effectiveness of the designed scheme is verified by numerical experiments.

29 pages

References in corpus (6)

• Numerical Methods for Fractional Diffusion
• Fractional Fokker-Planck dynamics: Numerical algorithm and simulations
• Fast algorithms for convolution quadrature of Riemann-Liouville fractional derivative
• Fractional compound Poisson processes with multiple internal states
• Levy walk with multiple internal states
• Numerical scheme for the Fokker-Planck equations describing anomalous diffusions with two internal states

Cited by in corpus (2)

• Quantifying model uncertainty for the observed non-Gaussian data by the Hellinger distance
• On the bilateral preconditioning for an L2-type all-at-once system arising from time-space fractional Bloch-Torrey equations