Vertical up-flow of supercritical fluid in the subchannel of a heated rod bundle was numerically simulated using the Computational Fluid Dynamics (CFD) codes ANSYS CFX and ANSYS FLUENT. A total of seven cases from three different sets of experiments were simulated. Three-dimensional steady-state predictions of fluid velocity, pressure, and temperature were made using five versions of two-equation RANS turbulence models with accompanying wall treatments. In addition, the temperature distribution in a solid region comprising a heater and sheathing was also computed in some cases. The k-epsilon turbulence model, implemented using CFX and scalable wall functions, provided the numerical results that have the smallest overall deviation from experimental results for three of the seven cases, and predicts the experimental data of the remaining four cases reasonably well, unlike other turbulence models that severely over-predict the experimental data for wall surface temperature.