The flow around a vertical circular pile exposed to a steady current is studied numerically and experimentally. The numerical model is a three-dimensional model. The model validation was achieved against new experimental data (which include two-component laser-Doppler anemometry (LDA) flow measurements and the hot-film bed shear stress measurements, and reported in the present paper) and the data of others, and a $k$–$\omega $ turbulence model was used for closure. The model does not have a free-surface facility and therefore is applicable only to cases where the Froude number is small ($\hbox{\it Fr}\,{<}\,O(0.2)$). The flow model was used to study the horseshoe vortex and lee-wake vortex flow processes around the pile. The influence on the horseshoe vortex of three parameters, namely the boundary-layer thickness, the Reynolds number and the bed roughness, was investigated. In the latter investigation, the steady solution of the model was chosen. A study of the influence of the unsteady solution on the previously mentioned flow processes was also carried out. The ranges of the parameters covered in the numerical simulations are: The boundary-layer-thickness-to-pile-diameter ratio is varied from 2$\,{\times}\, 10^{-2}$ to 10$^{2}$, the pile Reynolds number from 10$^{2}$ to $2\,{\times}\, 10^{6},$ and the pile diameter-to-roughness ratio from 2 to about 10$^{3}.$ The amplification of the bed shear stress around the pile (including the areas under the horseshoe vortex and the lee-wake region) was obtained for various values of the previously mentioned parameters. The steady-state flow model was coupled with a morphologic model to calculate scour around a vertical circular pile exposed to a steady current in the case of non-cohesive sediment. The morphologic model includes (i) a two-dimensional bed load sediment-transport description, and (ii) a description of surface-layer sand slides for bed slopes exceeding the angle of repose. The results show that the present numerical simulation captures all the main features of the scour process. The equilibrium scour depth obtained from the simulation agrees well with the experiments for the upstream scour hole. Some discrepancy (up to 30%) was observed, however, for the downstream scour hole. The calculations show that the amplification of the bed shear stress around the pile in the equilibrium state of the scour process is reduced considerably with respect to that experienced at the initial stage where the bed is plane.