The natural convection heat transfer from two attached cylinders set in a vertical array in air has been studied numerically for the Rayleigh number varying from 10 to 10(5). The comprehensive results consisted of streamlines and isothermal contours, temperature and velocity profiles, as well as local and average Nusselt number are presented with different Rayleigh numbers. They show that the heat transfer for both cylinders is degraded relative to that for a single cylinder owing to the interactions of their plumes. That trend mainly stems from the greater effects of the interaction on the upper part of the lower cylinder and the bottom surface of the upper cylinder. By contrast, it has little influence on the bottom surface of the lower cylinder and top region of the upper cylinder thereby leading to the similar results to that for a single cylinder there, but with a small noticeable difference. Also, their stronger interactions lead to the heat transfer rate of the upper cylinder reduced by 16-53% relative to that at S/D = 2. What's more, the correlating equations for individual cylinder and the whole array have been proposed respectively, which provide an accurate prediction for the engineering calculations applying the present configuration. (C) 2017 Elsevier Ltd. All rights reserved.