This work presents buoyancy-driven mixed convective flow and heat transfer phenomena of an isothermally heated horizontal elliptic cylinder in vertically upward unbounded flow of power-law type non-Newtonian nanofluids using ANSYS Fluent. The governing continuity, momentum and energy equations for the shear-thinning power-law nanofluids along with suitable boundary conditions are simultaneously solved within the limitations of Boussinesq approximation. The semi implicit method for pressure-linked equations algorithm along with the quadratic upstream interpolation for convective kinematics scheme for discretizing the convective terms in both momentum and energy equations are adopted. The ranges of parameters considered for this study are: volume fraction of nanoparticles, 0.005–0.045; aspect ratio of elliptic cylinder, 0.5–2.5; and Richardson number, 0–40; and a representative Reynolds number of 20. The streamline patterns, surface pressure coefficient distributions, total drag coefficients, isotherm contours, and Nusselt numbers are presented for better understanding of heat transfer and flow phenomena around elliptic cylinders. Briefly results indicate that the total drag coefficient is found to increase with the increasing Richardson number whereas it decreases with the increasing volume fraction of nanoparticles. The average Nusselt numbers are found to increase with increasing Richardson number and increasing volume fraction of nanoparticles.