On a High-Order Accurate Method for the Numerical Study of Natural Convection in a Vertical Square Cavity