Heat transfer often occurs effectively from horizontal elements of relatively complex shapes in natural convective cooling of electronic and electrical devices used in industrial applications. The effect of complex surface shapes on laminar natural convective heat transfer from horizontal isothermal polygons of hexagonal and octagonal flat surfaces facing upward and downward of different aspect ratios has been numerically investigated. The polygons’ surface is embedded in a large surrounding plane adiabatic surface, where the adiabatic surface is in the same plane as the surface of the heated element. For the Boussinesq approach used in this work, the density of the fluid varies with temperature, which causes the buoyancy force, while other fluid properties are assumed constants. The numerical solution of the full three-dimensional form of governing equations is obtained by using the finite volume method-based computational fluid dynamics (CFD) code, FLUENT14.5. The solution parameters include surface shape, dimensionless surface width, different characteristic lengths, the Rayleigh number, and the Prandtl number. These parameters are considered as follows: the Prandtl number is 0.7, the Rayleigh numbers are between 10 3 and 10 8 , and for various surface shapes the width-to-height ratios are between 0 and 1. The effect of different characteristic lengths has been investigated in defining the Nusselt and Rayleigh numbers for such complex shapes. The effect of these parameters on the mean Nusselt number has been studied, and correlation equations for the mean heat transfer rate have been derived.