The turbulent vortex shedding flow behind a square cube is numerically analyzed by the LES approach. The eddy viscosity and Smagorinsky models have been used to simulate the SGS Reynolds stress and eddy viscosity, respectively. The performance of the computations is evaluated by comparing with available measurements, as well as with those results obtained from RANs. It is found that the LES approach predicts the Strouhal number and the mean stream‐wise velocity profiles with a better accuracy than RANs.