Skip to content

Re: Smagorinsky/Turbulence – Modification Code

Due to recent bot attacks we have changed the sign-up process. If you want to participate in our forum, first register on this website and then send a message via our contact form.

› Forums › OpenLB › General Topics › Smagorinsky/Turbulence – Modification Code › Re: Smagorinsky/Turbulence – Modification Code

October 9, 2017 at 4:48 pm #2733
pedrobrito
Member

Hi,

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.

Best,
Pedro.