An exponential integrator multicontinuum homogenization method for fractional diffusion problem with multiscale coefficients

In this paper, we present a robust and fully discretized method for solving the time fractional diffusion equation with high-contrast multiscale coefficients. We establish the homogenized equation in a coarse mesh using a multicontinuum approach and employ the exponential integrator method for time discretization. The multicontinuum upscaled model captures the physical characteristics of the solution for the high-contrast multiscale problem, including averages and gradient effects in each continuum at the coarse scale. We use the exponential integration method to address the nonlocality induced by the time fractional derivative and the stiffness from the multiscale coefficients in the semi-discretized problem. Convergence analysis of the numerical scheme is provided, along with illustrative numerical examples. Our results demonstrate the accuracy, efficiency, and improved stability for varying order of fractional derivatives.