For the numerical simulation of convection-dominated reacting flow problems governed by convection-reaction equations, grids-based Eulerian methods may cause different degrees of either numerical dissipation or unphysical oscillations. In this paper, a Lagrangian particle algorithm based on the smoothed particle hydrodynamics (SPH) method is proposed for convection-reaction equations and is applied to an autocatalytic reaction model with multicomponent reactants. Four typical Eulerian methods are also presented for comparison, including the high-resolution technique with the Superbee flux limiter, which has been considered to be the most appropriate technique for solving convection-reaction equations. Numerical results demonstrated that when comparing with traditional first- and second-order schemes and the high-resolution technique, the present Lagrangian particle algorithm has better numerical accuracy. It can correctly track the moving steep fronts without suffering from numerical diffusion and spurious oscillations.