The buoyancy and initial momentum fluxes make near-field dominated by buoyant jet when thermal discharge releases underwater. In order to estimate prediction capabilities of those near-field phenomena, non-hydrostatic RANS applied CFD(Computational Fluid Dynamic) model was used. Condition of model was composed based on past laboratory experiments. Numerical simulations carried out for the horizontal buoyant jet in the stagnant flow and vertical buoyant jet into crossflow. The results of simulation are compared with the terms of trajectory and dilution rate of laboratory experiments and analytic model(CorJET) results. CFD model showed a good agreement with them. CFD model can be appropriate for assessment of submerged thermal discharge effect because CFD model can resolve the limitations of near-field analytic model and far-field quasi 3D hydrodynamic model. The accuracy and capability of the CFD model is reviewed in this study. If the computational efficiency get improved, CFD model can be widely applied for simulation of transport and diffusion of submerged thermal discharge.