In this paper, the Large Eddy Simulation (LES) was used to study the free-surface turbulent channel flow with passive heat transfer. The three-dimensional filtered incompressible Navier-Stokes equations and energy equation were numerically solved with dynamic Subgrid Scale (SGS) models for modeling turbulent stresses and heat flux. To compare the turbulent behavior of the free-surface and two-walled channel flows, the LES of two-walled turbulent channel flow was performed. The statistical quantities and flow structures of the free-surface turbulence with heat transfer in the vicinity of the free-surface were investigated. The results are also in good agreement with theoretical analysis and available results by Direct Numerical Simulation (DNS).