Abstract: The turbulence model used in the current numerical simulation for heat transfer to supercritical water is developed based on theory and experimental data at subcritical pressures and the density fluctuation is neglected. However, the thermo-physical properties of supercritical water vary significantly as temperature approaches the pseudo-critical temperature so that the density fluctuation should not be neglected for its important role in the turbulence model. The simulation method considering density fluctuations was developed, in which the thermal expansion coefficient fluctuation was also taken into account. The density fluctuation was implemented in the AKN turbulence model and the model prediction was assessed using experimental data. The validation results show that the predicting results agree better with experimental data. It is recommended that the density fluctuation should be taken into account in future turbulence simulation of heat transfer for supercritical water.