The object of this study is to establish CFD simulation code of LES so as to calculate a thermally-stratified turbulent boundary layer which is often observed and affects the climate. First, to evaluate a prediction precision of various SGS models of LES under the modeled condition of flow field in urban environment, various thermally-stratified turbulent boundary layers are calculated by the evaluated LES models. Also, LES is the effective calculation technique for an analysis of a turbulent heat and mass transfer in large Reynolds and Prandtl numbers fluids, because LES can calculate turbulent heat and mass transfer phenomena using fewer grid points as compared with DNS. In this study, these SGS models are evaluated in comparison with DNS database, in which problems of LES in the prediction of various thermally-stratified boundary layers are clarified.