Density-based solver is because of the dependence of density on both temperature and pressure near supercritical temperature and pressure . I have tried using the pressure-based coupled scheme but the probe values are oscillating. One more thing, when I used constant temperature condition on the cold wall convergence resulted. But after using convective boundary condition the convergence rate is extremely slow. I have used steady-state solver