In the present work, the onset and subsequent development of slug flow in horizontal pipes is investigated by solving the transient one-dimensional version of the two-fluid model in a high resolution mesh using a finite volume technique. The methodology (named slug-capturing) was proposed before in the literature and the present work represents a confirmation of its applicability in predicting this very complex flow regime. Further, different configurations are analyzed here and comparisons are performed against different sets of experimental data. Predictions for mean slug variables were in good agreement with experimental data. Additionally, focus is given to the statistical properties of slug flows such as shapes of probability density functions of slug lengths (which were represented by gamma and log-normal distributions) as well as the evolution of the first statistical moments, which were shown to be well reproduced by the methodology.