High-Order Difference Scheme for Time-Fractional Quasilinear Parabolic Equations

Mathematical modeling of heat and mass transfer processes in porous media using fractional derivative equations is of great practical importance. Within the framework of such models, obtaining analytical solutions to the corresponding initial–boundary value problems is generally difficult. In this work, we numerically investigate quasilinear parabolic problems involving Caputo time-fractional derivatives. First, the well-posedness and existence of weak solutions are discussed. Then, we construct and implement a finite-difference scheme that is fourth-order accurate in space and second-order accurate in time. Convergence in the maximum norm is proven. Numerical experiments confirm the accuracy and efficiency of the proposed approach.