The suitable turbulence model is found to be required in the course of establishing a proper analysis methodology for thermal stripping phenomena which are shown in strong temperature variation area such as reactors and propulsion devices. Three different turbulence models of $kappa$-$varepsilon$ model, modified $kappa$-$varepsilon$ model, and full Reynolds stress(FRS) model, are applied to analyze unsteady turbulent flows with temperature variation. Three test cases are selected for verification. These are vertical jet flows with water and sodium, and parallel jet flow with sodium. Analysis yields the conclusion that 3-D computation with FRS betters others. However, modified modeling is required to improve its heat transfer characteristic analysis. Further analysis is performed to find momentum variation effects on temperature distribution. It is found that the momentum increase results increase of fluid mixing and magnitude of temperature variation.