A Hybrid Numerical Framework Based on Radial Basis Functions and Finite Difference Method for Solving Advection–Diffusion–Reaction-Type Interface Models

Advection–diffusion–reaction-type interface models have wide-ranging applications in environmental science, chemical engineering, and biological systems, particularly in modeling pollutant transport in groundwater, reactive flows, and drug diffusion across biological membranes. This paper presents a novel numerical method for the solution of these models. The proposed method integrates the meshless collocation technique with the finite difference method. The temporal derivative is approximated using a finite difference scheme, while spatial derivatives are approximated using radial basis functions. The interface across the fixed boundary is treated with discontinuous diffusion, advection, and reaction coefficients. The proposed numerical scheme is applied to both linear and non-linear models. The Gauss elimination method is used for the linear models, while the quasi-Newton linearization method is employed to address the non-linearity in non-linear cases. The L∞ error is computed for varying numbers of collocation points to assess the method’s accuracy. Furthermore, the performance of the method is compared with the Haar wavelet collocation method and the immersed interface method. Numerical results demonstrate that the proposed approach is more efficient, accurate, and easier to implement than existing methods. The technique is implemented in MATLAB R2024b software.