Blood flow through a narrowed artery, such as those caused by plaque buildup, can lead to various cardiovascular diseases. Therefore, studying this phenomenon is essential for developing better treatments and understanding the underlying mechanisms. In this study, a mathematical model is developed to simulate blood flow through stenotic arteries with porous walls. The familiar Navier-Stokes equations are executed to simulate flow, whereas the Carreau fluid model is employed to quantify blood rheology. This model illustrates both Newtonian and non-Newtonian characteristics depending on the shear rate. Other than that, the blood vessel is assumed to have a moderate level of stenosis with porous saturated walls. The solution of the governed partial differential equations is carried out using a finite difference scheme. Our research addresses the impacts of permeability and rheological parameters, stenosis size, Brinkman number, and Prandtl number on blood velocity, flow rate, resistance to flow, wall shear stress, and temperature distribution. The findings demonstrate that while flow resistance reduces, the velocity, shear stress, and flow rate increase with high permeability values. Conversely, as the stenosis amplitude increases, the magnitude of wall shear stresses, velocity, and flow rates decrease. In addition, blood rheology also has the potential to transfigure the aforementioned factors. We also found that the temperature distribution increases with the Brinkman number, stenosis size, permeability, and power law exponent. It drives down, however, by raising the values of Weissenberg number and Prandtl number. Collectively, our findings contribute to a better understanding of the complicated phenomena of blood flow in stenotic vessels with porous saturated walls. The findings of this study may be valuable in establishing improved therapies for cardiovascular disorders and formulating more precise mathematical models to forecast blood flow in vivo.