A study of natural convective heat transfer in a two-dimensional rectangular enclosure with a cold source is reported. PHOENICS CFD software package is used to solve the conservation of momentum equation and energy equation by using Κ-ε turbulence model associated with wall functions. The cold source side is taken as the low temperature, whereas the overall heat transfer coefficient to be constant is taken as the high temperature boundary condition at the other three sides. The results include the distributions of flow velocity and temperature in the enclosure under different aspect ratios, and a correlation of the relationship of mean Nusselt number and modified Rayleigh number. When 1×10^9≤Ra≤1×l0^11,1.0≤AR≤2.0 and 1W/m^2K≤h≤10W/m^2K, the relation is (average)Nu=0.141Ra^0.212(AR)^0.362(hW/k(subscript f))^0.343.