A high-order unconditionally stable numerical method for a class of multi-term time-fractional diffusion equation arising in the solute transport models

In this paper, we study a high-order unconditionally stable numerical method to approximate the class of multi-term time-fractional diffusion equations. This type of problem appears in the modelling of transport of certain quantities such as heat, mass, energy, solutes in ground water and soils. The multi-term time-fractional derivative is approximated by using the Crank-Nicolson method for the Caputo's time derivative. The space derivative is approximated by using the collocation method based on quintic B-spline basis functions. We have established the stability and convergence analysis of the proposed numerical scheme thoroughly, and it is shown that the order of convergence in space variable is almost four and in the time variable is O (Delta t(2-max{gamma,gamma i})). To prove the accuracy and efficiency of the developed method, we consider four numerical examples and perform the numerical simulation. The developed algorithm works well andvalidate the theoretical results. The developed method is fourth-order convergent in the space variable, which is almost two orders of magnitude higher than the other spline collocation methods.

Automated summary

In this paper, a high-order numerical method is studied to approximate multi-term time-fractional diffusion equations. The method is shown to be stable and convergent, and its accuracy and efficiency are validated through numerical simulations.