Constrained Newtonian fluids containing reacting pollutants flowing through porous media were modeled using a continuum theory of mixtures approach, accounting for viscous linear and quadratic drag terms. The work considered very low pollutant concentrations and similar mass densities for the fluid and all the pollutant constituents. These constituents may undergo chemical reactions among themselves, but they do not react with the fluid constituent. In this case, equal velocities for the pollutant constituents and the fluid constituent may be assumed. A solid constituent (the porous matrix, rigid and at rest), an inert gas (to account for the mixture's compressibility), an incompressible fluid (the fluid constituent), and N pollutant constituents compose the mixture. The model encompasses the mass equations for the fluid constituent and the N-pollutant constituents, as well as the fluid constituent's momentum balance. The Glimm scheme, associated with an operator-splitting technique, was used to advance in time. The complete solution to the associated Riemann problem was presented. The numerical results show the evolution of the fluid constituent's velocity and saturation, and the pollutant constituents' concentrations in the mixture for selected initial data, accounting for Darcy and Forchheimer drag terms, and chemical reactions among the pollutants. The novelty lies in the combination of linear and quadratic drag terms with simple chemical reactions among the pollutants, the use of a new constitutive relation for pressure (accounting for kinematic constraints), and a new protocol for the Glimm scheme.