A stable Gaussian radial basis function method for solving nonlinear unsteady convection–diffusion–reaction equations

We investigate a novel method for the numerical solution of two-dimensional time-dependent convection–diffusion–reaction equations with nonhomogeneous boundary conditions. We first approximate the equation in space by a stable Gaussian radial basis function (RBF) method and obtain a matrix system of ODEs. The advantage of our method is that, by avoiding Kronecker products, this system can be solved using one of the standard methods for ODEs. For the linear case, we show that the matrix system of ODEs becomes a Sylvester-type equation, and for the nonlinear case we solve it using predictor–corrector schemes such as Adams–Bashforth and implicit–explicit (IMEX) methods. This work is based on the idea proposed in our previous paper (2016), in which we enhanced the expansion approach based on Hermite polynomials for evaluating Gaussian radial basis function interpolants. In the present paper the eigenfunction expansions are rebuilt based on Chebyshev polynomials which are more suitable in numerical computations. The accuracy, robustness and computational efficiency of the method are presented by numerically solving several problems.