The numerical study reported here deals with the passive flow control around a two-dimensional D-shaped bluff body at a Reynolds number of . A small circular control cylinder located in the near wake behind the main bluff body is employed as a local disturbance of the shear layer and the wake. 3D simulations are carried out using a newly developed very large eddy simulation (VLES) method, based on the standard k − ε turbulence model. The aim of this study is to validate the performance of this method for the complex flow control problem. Numerical results are compared with available experimental data, including global flow parameters and velocity profiles. Good agreements are observed. Numerical results suggest that the bubble recirculation length is increased by about 36% by the local disturbance of the small cylinder, which compares well to the experimental observations in which the length is increased by about 38%. A drag reduction of about 18% is observed in the VLES simulation, which is quite close to the experimental value of 17.5%. It is found that the VLES method is able to predict the flow control problem quite well.