Electroosmosis is the predominant mechanism for flow generation in lab-on-chip devices. Since most biofluids encountered in these devices are considered to be non-Newtonian, it is vital to study the flow characteristics of common non-Newtonian models under electroosmotic body force. In this paper, the hydrodynamically fully developed electroosmotic flow of power-law fluids in rectangular microchannels is analyzed. The electrical potential and momentum equations are numerically solved through a finite difference procedure for a non-uniform grid. A thoroughgoing parametric study reveals that the Poiseuille number is an increasing function of the channel aspect ratio, the zeta potential, the flow behavior index, and the dimensionless Debye-Hückel parameter. It is also found that the validity range of the Debye-Hückel approximation for shear-thickening fluids is much wider than that of shear-thinnings. Furthermore, while the dimensionless mean velocity is an increasing function of the channel aspect ratio and the dimensionless Debye-Hückel parameter, it is a decreasing function of the flow behavior index. Moreover, to increase the zeta potential is to increase the dimensionless mean velocity for shear-thinnings, nevertheless, its effect is not significant for shear-thickenings.
