Analysis of multiple steady states for natural convection of Newtonian fluids in a square enclosure