The current analysis focusses on the three-dimensional viscous flow caused by an expanding or contracting porous slider under the impact of heat source\(/\)sink, magnetic field and waste discharge concentration. The amount of fluid injected to levitate the slider is not constant; rather, it varies over time based on the slider's location at any given moment. The fluid flow between the ground and the slider is represented by the unsteady nonlinear equations of motion, which are converted into ordinary differential equations by employing the apt similarity transformations. These simplified equations are solved by applying the hybrid residual power series method (HRPSM). Using graphical representations, the physical behaviour of significant parameters is investigated. Tabular data are used to compare the hybrid residual power series method and the numerical scheme to check the convergence and accuracy of the hybrid residual power series method. A rise in the magnetic parameter results in a very modest rise in the fluid's longitudinal velocity in the regions closest to the upper or lower plates and a clear drop in the liquid's longitudinal velocity in the centre region. An increase in pollutant external-source parameter increases the amount of pollutants introduced into the fluid system.