Skip to content

how to add Smagorinsky model to LBM code

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

    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 3 months, 3 weeks ago by garcfd.
    #7692
    mathias
    Keymaster

    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

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