Many designs for natural-convection heat sinks and semi-empirical correlations have been proposed in the recent years, but they are only valid in a limited range of Elenbaas numbers El and were mostly tested for laminar flows. To alleviate those limits, parametric studies with 2D and quasi-3D models were carried out, in ranges of Grashof numbers up to 1.55×1011 and Elenbaas numbers up to 3.42×107. Ansys Fluent’s laminar, transition-SST, SST k-ω and k-ϵ models were applied. In addition, when used in this valid range, i.e., mean Elenbaas numbers, with the simplified quasi-3D model, the transition-SST model could predict better results, overestimating the heat flux by 10 to 15% compared to semi-empirical correlations. The 2D model was not deemed satisfying, regarding turbulence models. Consequently, a quasi-3D model was developed: it appeared to be an efficient trade-off between computational time and prediction accuracy, in particular for turbulence models. New grouping factors were also found, to ensure proper dimensioning of natural-convection heat sinks. They corresponded to non-dimensional parameters that dictated the physical behaviour of the heat sink with respect to the semi-empirical correlations. Typically, the ratio of the spacing to the optimal spacing predicted by Bar-Cohen’s correlation turned out to be an appropriate grouping factor with a threshold of 1, above which the fins could safely be considered as isolated, thus greatly simplifying all further calculations.