[pedrobrito]

Hi Marc,

I used the Smagorinsky turbulence model, by calling the SmagorinskyMRTdynamic
and by defining the MTR D2Q9 lattice arrangement.

Is there a limitation in the LBM for the number of refinement levels? (multi-grid approach) Is there an example in the folder where the SuperLatice/BlockLattice is used?

Best,
Pedro.

[pedrobrito]

Hi,

I tried to simulate the cylinder case for Re = 100000 using the BGK Smagorinsky (C=0.15). The size of the inlet was increased to 1m and I applied periodic BC at the walls of the domain (y-direction). However, the results are very different when a put the cylinder at the center line or dislocated some lattice units. When centered we do not see vortices and the Cd are considerably lower than the literature (when not centered the Cd is too high and does not show physical significance). Is there any mistake in my set up?

The resolution of the simulation was N=5 and M=2. Physical time was 20sec.
I tried to insert the MRT, it compiled with no warnings, but it always diverges.

Best,
Pedro.

[pedrobrito]

Hi Marc.

Thank you. It is working well now.

Best,
Pedro.

[pedrobrito]

Hi, Marc.

Thank you for answering me so quickly. I made the modification that you suggested. However, I keep receiving one error message.

cylinder2d.cpp: In function ‘int main(int, char**)’:
cylinder2d.cpp:381:10: error: type ‘class olb::SmagorinskyBGKdynamics<double, olb::descriptors::D2Q9Descriptor>’ argument given to ‘delete’, expected pointer
delete bulkDynamics;
^
Makefile:63: recipe for target ‘/mnt/c/Users/pedro/Documents/olb/examples/cylinder2d/cylinder2d.o’ failed
make: *** [/mnt/c/Users/pedro/Documents/olb/examples/cylinder2d/cylinder2d.o] Error 1
root@DESKTOP-TR9LPLJ:/mnt/c/Users/pedro/Documents/olb/examples/cylinder2d#

What else should I modify to make it run?

Best,
Pedro.

[pedrobrito]

Thanks for answering me, Marc.

Best,
Pedro.

[Author]

Posts