Revisiting turbulence modeling of Rayleigh–Bénard convection by using elliptic-relaxation approach