A Crank–Nicolson difference scheme for the time variable fractional mobile–immobile advection–dispersion equation

A Crank–Nicolson finite difference scheme to solve a time variable order fractional mobile–immobile advection–dispersion equation is introduced and analyzed. Some a priori estimates of discrete L2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^2$$\end{document}-norm with order of convergence O(τ+h2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$O(\tau +h^2)$$\end{document} are established on uniform grids where τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau $$\end{document} and h are the steps sizes in time and space. Stability and convergence of the numerical solutions are presented in detail. Numerical examples are provided to verify the theoretical analysis.