In oder to predict the dispersion of a thermal discharge with strong turbulent and buoyant effects, the development of a numerical model using turbulence model and its application are significantly increased. In this study, a three-dimensional steady-state model for the surface discharge heated water into quiescent water body is developed. For the model closure of turbulent terms the four-equation turbulence model is used. For economic mumerical simulation, the elliptic governing equations are transformed to the partially parabolic equations. In general, the simulated results by the present model agree well with the experimental results by Pande and Rajratnam (1977). The model characteristics are presented in comparison with the predicted results from the two-eqauation turbulence model by McGuirk and Rodi (1979).