I have added some Smagorinsky turbulence in to my LBM solver (to offer stability for high Re flows) but the problem is that the turbulence seems to maintain the vortex shedding behaviour even at high Reynolds number Re=3million say.
In the image below the flow is only Re=300 top (no turbulence), and bottom (with turbulence). So I am trying to figure out why adding turbulence in this way, causes this laminar type behaviour. The turbulence is added into the relaxation equation for omega as follows:
omega = 1.0/(3.0*(nuLB+nut(i,j,k)) + 0.5)