I’m new in LBM and trying to use SmargorinskyBGKDynamics to simulation a High Reynoalds flow but the result is always diverging.
I consider the material as air crossing a 1m*1m*5m wind tunnel , parameters like belows
I doubt there is any fault in the parameters. If I change the CharNu higher such as 1e-3, the result seemed good.
By the way, I imitate examples outflow as below:
sLattice.defineDynamics(superGeometry, 4, &bulkDynamics);
bc.addPressureBoundary(superGeometry, 4, omega);
lattice.defineRho(superGeometry, 4, rhoF);
I’m not quite understand it ,does it mean that the outflow is a constant pressure boundary?
In my past CFD simulation, I usually define the outflow a free pressure boundary that the pressure gradient is zero. How can I realize it?