Localized method of fundamental solutions for two- and three-dimensional transient convection-diffusion-reaction equations

This paper investigates the use of the localized method of fundamental solutions (LMFS) for the numerical solution of general transient convection-diffusion-reaction equation in both two-(2D) and three-dimensional (3D) materials. The method is developed as a generalization of the author's earlier work on Laplace's equation to transient convection-diffusion-reaction equation. The popular Crank-Nicolson (CN) time-stepping technology is adopted to perform the temporal simulations. The LMFS approach is then introduced for solving the resulting inhomogeneous boundary value problems, where a pseudo-spectral Chebyshev collocation scheme (CCS) is employed for the approximation of the corresponding particular solutions. As compared with the classical MFS and boundary element method (BEM), the present CN-CCS-LMFS approach produces sparse and banded stiffness matrix which makes the method possible to perform large-scale dynamic simulations. Several benchmark numerical examples are presented to demonstrate the efficiency and feasibility of the present method.