AbstractNatural convection in a bundle of vertical cylinders, arranged in equilateral triangular spacing, has been investigated numerically using a boundary‐fitted coordinate system. Numerical calculations for center‐to‐center distance between cylinders S/D = 1.1 to 1.9, 3.0, 4.0, and 7.0 were made of natural convection of air at modified Grashof numbers Gr* from 10 to 108. Local Nusselt number Nu for uniform wall heat flux indicates the same value at the axial locations except for the thermal entrance region. The region for respective cylinder spacing is noted to diminish with decreasing Grashof number. Numerical values of local Nusselt number Nui are in relatively good agreement with those obtained from the experiment for air. © 2003 Wiley Periodicals, Inc. Heat Trans Asian Res, 32(4): 330–341, 2003; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/htj.10095