Skip to content

how to add Smagorinsky model to LBM code

Viewing 2 posts - 1 through 2 (of 2 total)
  • Author
  • #7690

    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?

    • This topic was modified 9 months, 3 weeks ago by garcfd.

    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.


Viewing 2 posts - 1 through 2 (of 2 total)
  • You must be logged in to reply to this topic.