This paper investigates the effects of introducing wavy cylinders of appropriate wavelength into tube bundles for the passive control of flow-induced vibration. The Reynolds-averaged Navier–Stokes approach RANS was used to study the turbulent flow characteristics within staggered tube bundles having transverse and longitudinal pitch-to-diameter ratios of 3.8 and 2.1, respectively. Two wavy models are tested with different wavelength-to-mean diameter ratios, namely: Model-A with [Formula: see text] and Model-B with [Formula: see text], and their performances are compared with those of a smooth circular tube bundle Model-C ([Formula: see text]), in a range of subcritical Reynolds numbers [Formula: see text]. Mean velocity and pressure fields as well as aerodynamic coefficients are presented and compared. According to the obtained results, the wavy tube bundle with [Formula: see text] exhibits the best performance, with slightly smaller [Formula: see text] values than the smooth tube bundle. However, results show that the wavy model with [Formula: see text] is clearly the worst with an increase of [Formula: see text] up to 40%. Numerical outcomes also reveal that the interaction of neighboring rows and columns affects the fluctuation of [Formula: see text] values relative to each tube position. Moreover, [Formula: see text] evolutions depend significantly on the variation of Re as well, where drag reduction is more substantial by increasing Re. Thus, wavy-shaped tubes with optimum parameters could be more suitable for flow control by enhancing the secondary flow, turbulence and fluid mixing.