关于LBM-LES的问题，求高手帮忙

hs6951538

我用Smagorinsky模型做LBM-LES的顶盖方腔流，可不知道为什么雷诺数总上不去，256*256、c=1.0的网格用反弹边界只能算到3500左右，而且所得结果的压力场也很不好，有很多波纹，界限不清晰。我是用LBGK模型做的二维方腔流，看了一些文献，觉得就是把松弛系数改一下就可以了，公式应该是tau_total=sqrt(tau_f*tau_f+18.0*CS*CS*Qave/RHO )/2.0+tau_f/2.0，Qave=sqrt( 2.0*(Qxx*Qxx+Qyy*Qyy+Qxy*Qxy+Qyx*Qyx) )没错吧？可看别人用反弹边界都可以算到Re=7500甚至上万，还是哪里有问题吗？程序就是用何雅玲老师书上的方腔流改的，只是改了边界条件多了计算松弛系数的部分。郁闷啊，搞了很久都不知道哪里的问题 。

wdlxmzd

你可以先用非平衡态外推试一下，然后何老师书后面的程序改到反弹边界时，需要改动一些地方，如果没注意，可能就有问题

hs6951538

谢谢回复啊。不加大涡模拟时反弹边界我也算过，应该没有问题，不论什么边界，只是加了Smagorinsky模型的大涡模拟后Re提高太小了，所得压力场也不理想，所以我觉得是不是自己哪里弄错了，可始终检查不出来。S.Hou的《A Lattice Boltzmann Subgrid Model for High Reynold Number Flows》至少算到上万雷诺数，而且似乎使用的也是反弹边界。请问你加大涡模拟后同样情况下雷诺数能提高多少呢？

luo@odu.edu

You may want to take a look at the paper by H. Yu, L.-S. Luo, and S.S. Girimaji.
"LES of turbulent square jet flow using an MRT lattice Boltzmann model."
Computers & Fluids 35(8/9):957-965 (September, 2006).

hs6951538

谢谢罗老师的提醒！