Dear all, I’m using Palabos to simulate a 3D turbulent flow with noslip BC on bottom, freeslip on top, periodic on lateral sides, a velocity profile on inlet and outflow BC on outlet. The Re number in my simulation will up to 1e4 so I use smagorinsky model.
But when I combine these BC together, the result always seems unreasonable. Here is one of my codes I tried:
Box3D top = Box3D(0, nx-1, ny-1, ny-1, 0, nz-1);
Box3D inlet = Box3D(0, 0, 1, ny-1, 0, nz-1);
Box3D outlet = Box3D(nx-1, nx-1, 1, ny-2, 0, nz-1);
boundaryCondition.addVelocityBoundary0N (inlet, lattice );
boundaryCondition.addVelocityBoundary0P (outlet, lattice, boundary::outflow );
boundaryCondition.addVelocityBoundary1P (top, lattice, boundary::freeslip );
T cSmago = 0.12;
setBoundaryVelocity (lattice, inlet, TurbulentVelocity<T>(parameters) );
initializeAtEquilibrium(lattice, lattice.getBoundingBox(), TurbulentDensityAndZeroVelocity<T>(parameters) );
MultiBlockLattice3D<T, DESCRIPTOR> lattice (
parameters.getNx(), parameters.getNy(), parameters.getNz(),
new SmagorinskyRegularizedDynamics<T,DESCRIPTOR>(parameters.getOmega(), cSmago) );
OnLatticeBoundaryCondition3D<T,DESCRIPTOR>* boundaryCondition
// = createLocalBoundaryCondition3D<T,DESCRIPTOR>();
= createInterpBoundaryCondition3D<T,DESCRIPTOR>();
// lattice.periodicity().toggle(0, true); // periodic along the x-axis
lattice.periodicity().toggle(2, true); // periodic along the z-axis
I think there must be something wrong with my definition of BCs. I read some of the topic in this forum but I still can’t make it go right, so any help is highly appreciated.
Another question is when I increase the Re number to around 1e4 or less in simulation of turbulent flow, the velocity field near the newmann BC will become unstable. I tried InterpBC in palabos and the unstability still happen. I wonder is this due to my wrong settings or something else.
Best regard.