Glitches around cylinder boundary layer
OpenLB – Open Source Lattice Boltzmann Code › Forums › on Lattice Boltzmann Methods › General Topics › Glitches around cylinder boundary layer
- This topic has 5 replies, 2 voices, and was last updated 4 years, 10 months ago by Marc.
July 19, 2018 at 6:07 pm #1989
I met a problem when simulating 2D cylinder case at Re 4000. The resolution is 100 while the cylinder diameter is 1. BGK collision operator is applied. As shown in the vorticity contour,
the boundary layer around the cylinder is correctly presented, with negative vorticity on the top and positive vorticity on the bottom. However, there exist some vorticity strips just off the boundary layer, resulting in these unreasonable glitches.
I also tested the same case with MRT collision operator but got same problem. I didn’t find any literature addressing this problem and would like to ask you for help.
Looking forward to hearing from you.
EricJuly 20, 2018 at 12:00 am #2892
first guess: 100 is not a really fine resolution for a Reynoldys number of 4000.
I need some additional information to evaluate possible error sources:
– boundary condition at the cylinder
– used turbulence modell
– max Mach number
Best MarcJuly 20, 2018 at 12:06 pm #2893
– The boundary condition around the cylinder is bouzidi boundary condition.
– Both BGK and MRT solver are applied and glitches are found in both of them.
– The charLatticeVelocity is set equal to dx = 0.01, which mean the max Mach number should be 0.01/Cs = 0.0173 if I am not mistaken.
I am now testing it with resolution 200 but I am worried that if I wanna to simulate higher Reynolds number, increasing resolution will increase computational cost to a great degree.
EricJuly 24, 2018 at 8:50 am #2901
try to use instead of the bouzidi boundary condition the bounce back condition. Furthermore you can use a LES turbulence model SmagorinskyBGK or SmagorinskyMRT to avoid very high resolutions.
Best MarcJuly 26, 2018 at 11:36 am #2902
Thanks for your reply.
However, the bounce-back condition doesn’t seems to work. Are these strips around the cylinder reasonable in simulation? I mean, how do these strips and glitches arise in the vorticity contour?
I tried SmagorinskyBGK and SmagorinskyMRT as well. It seems these glitches can be tuned somewhat by adjusting Cs but they cannot be fully diminished. Also, at very high Re (50,000 in this case), there are upstream vortex formed in front of the cylinder, as can be show in the following plot:
EricJuly 28, 2018 at 2:49 pm #2903
bounce back have to work in your simulation. If bounce back doesn’t work, your material number configuration is maybe not correct. There is no physical reason for this kind of stripes. The reason for the uptstream vortices could be a reflection of the pressure waves from the cylinder at the inlet. You can try to increase the space between your inlet and your cylinder and to change the boundary condition to velocity bounce back at the inlet.
- You must be logged in to reply to this topic.