RBF-Based Meshless Collocation Method for Time-Fractional Interface Problems with Highly Discontinuous Coefficients

Time-fractional interface problems arise in systems where interacting materials exhibit memory effects or anomalous diffusion. These models provide a more realistic description of physical processes than classical formulations and appear in heat conduction, fluid flow, porous media diffusion, and electromagnetic wave propagation. However, the presence of complex interfaces and the nonlocal nature of fractional derivatives makes their numerical treatment challenging. This article presents a numerical scheme that combines radial basis functions (RBFs) with the finite difference method (FDM) to solve time-fractional partial differential equations involving interfaces. The proposed approach applies to both linear and nonlinear models with constant or variable coefficients. Spatial derivatives are approximated using RBFs, while the Caputo definition is employed for the time-fractional term. First-order time derivatives are discretized using the FDM. Linear systems are solved via Gaussian elimination, and for nonlinear problems, two linearization strategies, a quasi-Newton method and a splitting technique, are implemented to improve efficiency and accuracy. The method’s performance is assessed using maximum absolute and root mean square errors across various grid resolutions. Numerical experiments demonstrate that the scheme effectively resolves sharp gradients and discontinuities while maintaining stability. Overall, the results confirm the robustness, accuracy, and broad applicability of the proposed technique.