A multigrid method combined with defect correction for free convection problems at high Rayleigh numbers