In the present work, the wake behavior of wind turbines operating under thermallystratified atmospheric boundary layer (ABL) is numerically investigated. The steady state three dimensional Reynolds-Averaged Navier-Stokes (RANS) equations, combined with actuator disk approach, are used in the simulation. The standard ݇-ߝ turbulence model as well as two modified models namely Crespo model and El kasmi model are adopted. Two different methods are used and compared, to represent the atmospheric stratification conditions. In the first method, the energy equation is considered along with mass, momentum, and turbulence model equations. The stratification in the second method is modeled by means of an additional buoyancy production or dissipation term, which is added to the equation of the turbulent kinetic energy, instead of solving the energy equation. The results obtained from both methods show reasonable agreement with the experimental data available from the literature.