A novel numerical method for the time variable fractional order mobile-immobile advection-dispersion model

  • Authors:
  • H. Zhang;F. Liu;Mantha S. Phanikumar;Mark M. Meerschaert

  • Affiliations:
  • School of Mathematical and Computer Sciences, Fuzhou University, Fuzhou 350108, China;School of Mathematical Sciences, Queensland University of Technology, GPO Box 2434, Brisbane, Qld. 4001, Australia;Department of Civil and Environmental Engineering, Michigan State University, East Lansing, MI, USA;Department of Statistics and Probability, Michigan State University, East Lansing, MI 48824, USA

  • Venue:
  • Computers & Mathematics with Applications
  • Year:
  • 2013

Quantified Score

Hi-index 0.09

Visualization

Abstract

Evolution equations containing fractional derivatives can provide suitable mathematical models for describing anomalous diffusion and transport dynamics in complex systems that cannot be modeled accurately by normal integer order equations. Recently, researchers have found that many physical processes exhibit fractional order behavior that varies with time or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. In this paper, we consider the mobile-immobile advection-dispersion model with the Coimbra variable time fractional derivative which is preferable for modeling dynamical systems and is more efficient from the numerical standpoint. A novel implicit numerical method for the equation is proposed and the stability of the approximation is investigated. As for the convergence of the numerical method, we only consider a special case, i.e., the time fractional derivative is independent of the time variable t. The case where the time fractional derivative depends on both the time variable t and the space variable x will be considered in a future work. Finally, numerical examples are provided to show that the implicit difference approximation is computationally efficient.