A general numerical framework is designed for the two-dimensional convection–diffusion–reaction (CDR) system. The compatibility of differential quadrature and finite difference methods (FDM) are utilized for the formulation. The idea is to switch one numerical scheme to another numerical scheme without changing the formulation. The only requirement is to input the weighting coefficients associated with the derivative discretizations to the general algorithm. Three numerical schemes comprising combinations of differential quadrature and FDMs are studied using the general algorithm. Properties of numerical schemes and the algorithm are analyzed by using the simulations of two-dimensional linear CDR system, Burgers’ equation, and Brusselator model.