how to add Smagorinsky model to LBM code

[garcfd]

I have added Smagorinsky turbulence model into my LBM code, via the relaxation term
such that omega = 1/(3(nu_l+nu_t)+0.5), where nu_t is the contribution from
Smagorinsky model. However it doesn’t work because the turbulent contribution,
pushes the omega value back into the laminar regime for low Reynolds numbers,
giving lots of vortex shedding even at high Reynolds number. So how do you add
turbulence into LBM code correctly?

https://drive.google.com/file/d/1oJ7Ltr48jzt-YfZyVFS45jsVKt8qybWl/view?usp=drive_link

• This topic was modified 1 year, 4 months ago by garcfd.

[mathias]

Dear garcfd,

there might be a bug with the right scaling. You need to scale it to lattice units, which is sometime the problem. Our OpenLB implementation is working fine and is validated.

Best
Mathias

[Author]

Posts