The steady flow of non-Newtonian power-law fluids across a periodic square array of infinitely long circular cylinders is studied numerically using an unstructured finite volume method. The local and global flow characteristics have extensively been explored by the systematic variations of the pertinent dimensionless parameters as follows: fluid volume fraction (ϕ f = .70-.99), Reynolds number (Re = 1-40), and power-law index (n = .4-1.8). Qualitatively, the dense and curved streamlines are seen with the increasing inertial forces and shear-thinning behavior across all the fluid volume fractions. Further, the pressure coefficient over the surface of periodic cylinders is significantly influenced by the governing parameters and found to be maximum and minimum for the upstream and downstream cylinders, respectively. The individual and total drag coefficients have shown complex dependence on n, ϕ f , and Re. For shear-thinning fluids (n < 1), the pressure drag coefficient dominates over the friction drag coef-