The laminar flow and heat transfer across a triangular periodic array of heated cylinders are simulated computationally and analyzed. The study has been carried out at Reynolds number 10-100 for fluid volume fraction ranging from 0.7 to 0.99 and Prandtl number ranging from 0.7 to 50. The size of the wake region increases continuously with an increase in the Reynolds number for all values of fluid volume fraction. The recirculation bubble from the rear of a cylinder is reaching the front of the next cylinder in the same column of the periodic array for low values of free volume fraction, but this is not the case with the highest free volume fraction, i.e., 0.99. At high Reynolds number, the flow is separating early on the cylinder surfaces. The wake size at higher Reynolds number 75 and 100 for the lowest free volume fraction 0.7 is more in comparison with the wake size at free volume fraction 0.99, which is explained by plotting the location of flow separation against Reynolds number for both the extreme values of free volume fractions, i.e., 0.7 and 0.99. The isovorticity contours are concentrated in the vicinity of the cylinders on increasing the Reynolds number irrespective of free volume fraction and then convected downstream. On increasing free volume fraction, the friction and pressure drags in the array decrease. The increase in Reynolds number also results in the decrease in the values of the individual (friction and pressure drag coefficients) as well as total drag coefficients for all values of free volume fraction. At high values of Reynolds number, the emergence of carbuncle or thermal spike on isotherm near the cylinder's surface is observed where the value of the local Nusselt number is observed low. The heat transfer improves and the Nusselt number increases as the Reynolds number and/ or Prandtl number increases. On the contrary, heat transfer decreases as free volume fraction increases.