Flow of generalized Newtonian fluids in porous media can be modeled as a bundle of capillary tubes or a pore-scale network. In general, both approaches rely on the solution of Hagen–Poiseuille equation using power law to estimate the variations in the fluid viscosity due to the applied shear rate. Despite the effectiveness and simplicity, power law tends to provide unrealistic values for the effective viscosity especially in the limits of zero and infinite shear rates. Here, instead of using power law, Carreau model (bubbles, drops, and particles in non-Newtonian fluids. Taylor & Francis Group, New York, 2007) is used to determine the effective viscosity as a function of the shear strain rate. Carreau model can predict accurately the variation in the viscosity at all shear rates and provide more accurate solution for the flow physics in a single pore. Using the results for a single pore, normalized Fanning friction coefficient has been calculated and plotted as a function of the newly defined Reynolds number based on pressure gradient. For laminar flow, the variation in the friction coefficient with Reynolds number has been plotted and scaled. It is observed that generalized Newtonian fluid flows show Newtonian nature up to a certain Reynolds number. At high Reynolds number, deviation from the Newtonian behavior is observed. The main contribution of this paper is to present a closed-form solution for the flow in a single pore using Carreau model, which allows for fast evaluation of the relationship between flux and pressure gradient in an arbitrary pore diameter. In this way, we believe that our development will open the perspectives for using Carreau models in pore-network simulations at low computational costs to obtain more accurate prediction for generalized Newtonian fluid flows in porous media.