In this paper, numerical computations for analyzing the properties of three-dimensional indoor thermal environment under two different cases, which include with and without fresh-air supply, are conducted based on a Viollet type k-ε turbulence model and a proposed hybrid scheme. The buoyancy effects on the indoor airflow distribution are investigated. Experimental results show the feasibility of the approach.