A numerical analysis of the flow behavior in industrial cooling systems based on arrays of impinging jets has been performed, using several eddy viscosity models to determine their modeling capabilities. For the cooling system studied, and in terms of mean Nusselt number values, the best agreement between experimental results and numerical predictions was obtained with the realizable k-e model. On the other hand, numerical predictions of the local Nusselt number and its spatial variations along the wall are better adjusted to the experiments when using either the standard k-e or the standard k-¿ models. The results obtained also show that the predicted thermal field depends strongly on the combination of near-wall treatment and selected turbulence model.