The convective heat transfer is an important component in the total energy balance for arrays of buildings immersed in a turbulent boundary layer. This study examines convective heat transfer coefficients (CHTC) at the external windward, leeward, lateral and top surfaces of buildings. This study uses Large Eddy Simulations (LES) with Smagorinsky-Lily model to predict CHTC and then compares the simulation results with experimental data. In addition, steady RANS including the realizable k-epsilon and the shear stress transport k-omega turbulence models are also validated with the experimental data. However, both of these RANS models overestimated CHTC values for arrays of buildings in contrast to LES predictions. Regular arrays of cubic buildings are modeled to investigate CHTC in an urban environment, which are arranged according to different plan area densities (lambda(p) = 0.44, 0.25, 0.16, 0.11, 0.063 and 0.04). This morphological parameter (lambda(p)) represents different urban neighborhoods and it is used to characterize different flow regimes in an urban environment. The CHTC distributions are independent of the Reynolds number based on different incoming wind velocities at a height of 10 m above the ground (U-10) for windward, leeward, lateral and top building surfaces. Furthermore, CHTC distributions for different lambda(p) and U-10 = 5 m/s are compared with flow characteristics in building arrays. Finally, the CHTC correlations as a function of lambda(p) and U-10 were obtained, where lambda(p) varies from 0.04 to 0.25 and the Reynolds number ranges from 7 x 10(5) to 5 x 10(6). With the increase of the plan area densities from 0.04 to 0.25, CHTCs increases 15% for the leeward surface and decreases 16% for the lateral surfaces. Consequently, a total energy balance and the energy consumption predictions for a building need to take into account the urban density of the building surroundings. (C) 2013 Elsevier Ltd. All rights reserved.