Numerical study of the nonlinear anomalous reaction–subdiffusion process arising in the electroanalytical chemistry

This paper presents a meshless method based on the finite difference scheme derived from the local radial basis function (RBF-FD). The algorithm is used for finding the approximate solution of nonlinear anomalous reaction–diffusion models. The time discretization procedure is carried out by means of a weighted discrete scheme covering second-order approximation, while the spatial discretization is accomplished using the RBF-FD. The theoretical discussion validates the stability and convergence of the time-discretized formulation which are analyzed in the perspective to the H¹-norm. This approach benefits from a local collocation technique to estimate the differential operators using the weighted differences over local collection nodes through the RBF expansion. Two test problems illustrate the computational efficiency of the approach. Numerical simulations highlight the performance of the method that provides accurate solutions on complex domains with any distribution node type.