In this paper two-dimensional Navier-Stokes equations and energy equation governing incompressible laminar flow past a bundle of cylinders were numerically solved by using the finite element method.The velocity correction method was used for time advancement, and spatial discretization was carried out with the Galerkin weighted residual method.Viscous flows past the cylinder banks arranged in in-line cylinder bundles and staggered cylinder bundles, coupled with heat transfer, were investigated for pitch-diameter ratios of 1.5 and 2.0 and the Reynolds numbers from 50 to 500.Flow structures and heat transfer behavior were discussed.The results obtained agree well with available numerical data.