Thank you very much for your reply. I will use Smag. model with BGK.
One more clarification. I do not know about Pressure BC in OpenLB package, can it handle the following case: v(u, w), u=dw/dy=0 on the upper boundary?
Or Pressure BC handles it somehow different?