The mixed convective momentum and heat transfer phenomena of confined square cylinders in non‐Newtonian nanofluids are numerically investigated. The experimental thermophysical properties of alumina‐water‐based nanofluids are adopted from literature and these nanofluids obey shear‐thinning power‐law type non‐Newtonian behavior. The square cylinder is confined in a vertical channel with a confinement ratio of 0.1333. The flow is assumed to be two‐dimensional and the fluid is allowed to flow in upward direction across the confined square cylinder in the vertical channel. The aiding/opposing buoyancy in the flow is incorporated in terms of Richardson number (Ri) in the range of –2 to 2. The ranges of other dimensionless parameters considered are: Reynolds number, Re: 1 to 40; and volume fraction of nanoparticles, ϕ: 0.005 to 0.045. This range of volume fraction of nanoparticles (i.e., ϕ = 0.005 to 0.045) corresponds to the power‐law index (n) of a non‐Newtonian nanofluid in the range of n = 0.88 to 0.5, respectively. Prior to obtaining new results, the solution methodology is validated with existing literature counterparts. Finally, effects of the Reynolds number, Richardson number, and the rheology of non‐Newtonian nanofluids on streamline patterns, surface pressure, surface vorticity, drag coefficients, isotherm contours, local and average Nusselt numbers are delineated.