Memory-Driven Anomalous Heat Transport in Heterogeneous Media: A Two-Dimensional Time-Fractional Porous Medium Approach

Heat transport in heterogeneous materials can deviate markedly from classical Fourier behavior when microstructural disorder, trapping effects, nonlinear mobility, and long-range temporal correlations interact across multiple spatial and temporal scales. These mechanisms may produce delayed relaxation, persistent thermal footprints, front deformation, and non-classical spreading patterns that are not adequately represented by conventional integer-order diffusion models. In this study, a modeling and simulation framework is developed for anomalous heat transport in heterogeneous media using a two-dimensional time-fractional porous medium equation. The model combines a Caputo fractional time derivative, which represents thermal memory, with nonlinear degenerate porous-medium diffusion, spatially heterogeneous conductivity, localized volumetric heating, and Robin-type convective boundary exchange. A conservative fully discrete numerical scheme is constructed using flux-based finite differences for the heterogeneous nonlinear diffusion operator and an L1 approximation for the Caputo derivative. The nonlinear algebraic system at each time level is solved using an under-relaxed Picard frozen-coefficient iteration with non-negativity enforcement and sparse direct solution of the resulting linear systems. The numerical implementation is verified through a manufactured-solution convergence study, and additional analyses are performed to examine computational cost, Picard iteration behavior, coefficient-regularization sensitivity, strong-source effects, heterogeneous conductivity structures, and long-time thermal-footprint persistence. The results show that heterogeneous conductivity mainly redirects heat through preferential pathways and enlarges the spatial footprint while producing negligible changes in global heat content. Stronger fractional memory, represented by smaller fractional order, increases the persistence and spatial reach of moderate heating, whereas larger porous-medium exponents confine heat near the source and preserve higher local peaks. Source amplitude increases the thermal burden and footprint monotonically over the tested range, including strong forcing, without producing an abrupt localization-spreading transition. Boundary exchange remains secondary in the short-time interior-heating regime considered. These findings demonstrate that the proposed two-dimensional time-fractional porous medium framework provides a verified and physically interpretable model for non-Fourier heat transport in heterogeneous materials, where local intensity, global heat retention, and spatial thermal exposure must be assessed jointly.