Skip to content

Re: Smagorinsky/Turbulence – Modification Code



I tried to simulate the cylinder case for Re = 100000 using the BGK Smagorinsky (C=0.15). The size of the inlet was increased to 1m and I applied periodic BC at the walls of the domain (y-direction). However, the results are very different when a put the cylinder at the center line or dislocated some lattice units. When centered we do not see vortices and the Cd are considerably lower than the literature (when not centered the Cd is too high and does not show physical significance). Is there any mistake in my set up?

The resolution of the simulation was N=5 and M=2. Physical time was 20sec.
I tried to insert the MRT, it compiled with no warnings, but it always diverges.