Development of POD reduced-order model and its closure scheme for 2D Rayleigh–Bénard convection