This work presents a study on free convection in a porous square cavity saturated with a Newtonian fluid. Computations for laminar and turbulent flow are performed. Governing equations were time- and volume averaged according to the double-decomposition concept. Discretization of governing equations was obtained with the control-volume approach and the system of algebraic equation was relaxed via the SIMPLE method. Two energy models were employed, namely the one- and two-temperature models. Results indicated that when the ratio of thermal conductivities equals unity, both models give similar results. However, the overall Nusselt number across the cavity is reduced as porosity or the thermal conductivity ratio increases. A critical value for the Rayleigh number, understood as that when laminar and turbulent solution differ by a substantial amount, was found to be a function of the thermal conductivity ratio.