Current collection by positively polarized cylindrical Langmuir probes immersed in flowing plasmas is analyzed using a non-stationary direct Vlasov-Poisson code. A detailed description of plasma density spatial structure as a function of the probe-to-plasma relative velocity U is presented. Within the considered parametric domain, the well-known electron density maximum close to the probe is weakly affected by U. However, in the probe wake side, the electron density minimum becomes deeper as U increases and a rarified plasma region appears. Sheath radius is larger at the wake than at the front side. Electron and ion distribution functions show specific features that are the signature of probe motion. In particular, the ion distribution function at the probe front side exhibits a filament with positive radial velocity. It corresponds to a population of rammed ions that were reflected by the electric field close to the positively biased probe. Numerical simulations reveal that two populations of trapped electrons exist: one orbiting around the probe and the other with trajectories confined at the probe front side. The latter helps to neutralize the reflected ions, thus explaining a paradox in past probe theory. V C 2014 AIP Publishing LLC.