This work contributes to the study of flow over a circular cylinder at Reynolds number Re= 3900. Although this classical flow is widely documented in the literature, especially for this precise Reynolds number that leads to a subcritical flow regime, there is no consensus about the turbulence statistics immediately just behind the obstacle. Here, the flow is investigated both numerically with large eddy simulation and experimentally with hot-wire anemometry and particle image velocimetry. The numerical simulation is performed using high-order schemes and a specific immersed boundary method. The present study focuses on turbulence statistics and power spectra in the near wake up to ten diameters. Statistical estimation is shown to need large integration times increasing the computational cost and leading to an uncertainty of about 10% for most flow characteristics considered in this study. The present numerical and experimental results are found to be in good agreement with previous large eddy simulation data. Contrary to this, the present results show differences compared to the experimental data found in the literature, the differences being larger than the estimated uncertainty range. Therefore, previous numerical-experimental controversy for this flow seems to be reduced with the data presented in this article.

Three-dimensional direct numerical simulations of vortex shedding behind cylinders have been performed when the body diameter and the incoming flow involved spanwise linear nonuniformity. Four configurations were considered: the shear flow, the tapered cylinder and their combination which gave rise to namely the adverse and aiding cases. In contrast with the observations of other investigators, these computations highlighted distinct vortical features between the shear case and the tapered case. In addition, it was observed that the shear case and the adverse case (respectively tapered case and aiding case), yielded similarities in flow topology. This phenomenon was explained by the spanwise variations of U/D which seemed to govern these flows. Indeed, it was observed that large spanwise variations of U/D seemed to enhance three dimensionality, through the appearance of vortex adhesions and dislocations. Spanwise cellular pattern of vortex shedding was identified. Their modifications in cell size, junction position and number were correlated with the variation of U/D. In the lee side of the obstacle a wavy secondary motion was identified. Induced secondary flow due to the bending of Karman vortices in the vicinity of vortex adhesion and dislocations was suggested to explain this result.