Proper wastewater disposal is crucial in various manufacturing and ecological systems. This study aims to prevent and regulate pollution in the water supply. It examines how the pollutant discharge concentration affects the flow of non-Newtonian nanofluids (NNNFs) over a porous Riga surface. Two different types of NNNFs, namely, Walter’s B and second-grade fluids, have been examined. The fluid flow is conveyed in the form of a system of partial differential equations (PDEs), which are first reduced to a non-dimensional set of ordinary differential equations (ODEs) and then to first-order differential equations. The numerical approach parametric continuation method is employed to solve these ODEs. It has been noticed that the energy curve declines with increasing numbers of TiO2-nanoparticles (NPs). The effect of the external pollutant source variation factor enriches the concentration of pollutants in both fluid cases. Furthermore, the viscoelastic parameter K1 plays a notable role in determining the behavior of the fluids. Particularly in NNNFs, the variation of K1 enhances the fluid flow, whereas the rise of second-grade fluid factor decreases the velocity of the fluid. Our findings indicate a substantial impact of the parameters under consideration on the concentration of pollutant discharge. Significantly, it was observed that an increase in the amount of NPs and the thermal radiation parameter led to an improvement in the thermal conductivity of the nanofluid, consequently decreasing the concentration of pollutants in the discharge. The nanofluid has greater efficiency in boosting the energy transfer rate of the base fluid. In the case of the second-grade fluid, the energy propagation rate increases up to 6.25%, whereas, in the case of Walter’s fluid B, it increases up to 7.85%.