A Matrix Transform Technique for Distributed-Order Time-Fractional Advection–Dispersion Problems

Invoking the matrix transfer technique, we propose a novel numerical scheme to solve the time-fractional advection–dispersion equation (ADE) with distributed-order Riesz-space fractional derivatives (FDs). The method adopts the midpoint rule to reformulate the distributed-order Riesz-space FDs by means of a second-order linear combination of Riesz-space FDs. Then, a central difference approximation is used side by side with the matrix transform technique for approximating the Riesz-space FDs. Based on this, the distributed-order time-fractional ADE is transformed into a time-fractional ordinary differential equation in the Caputo sense, which has an equivalent Volterra integral form. The Simpson method is used to discretize the weakly singular kernel of the resulting Volterra integral equation. Stability, convergence, and error analysis are presented. Finally, simulations are performed to substantiate the theoretical findings.