MRT圆柱绕流问题

guchaonemo

最近初学LBE，在用MRT方法做圆柱绕流模拟，用C编写的程序，但是不知道怎么回事，每次都是前面运算的时候，相对误差越来越小，可是超过一定的步数的时候相对误差就越来越大（一般是1000步左右）。计算的时候，没用用到矩空间m，而是直接用
Inv (M)*diag（S）*M  代替代替碰撞矩阵。请高手帮忙解答一下（圆柱边界用的是bounce-back，其他边界是用的非平衡外推）边界处理应该没问题，用BGK的时候用的同样的边界处理

guchaonemo

今天把程序弄好了，但是发现MRT方法进行圆柱绕流的时候，松弛参数设置错了，很容易使得得不到正确的结果，有没有大神有关于参数设置的资料，推荐一下

luo@odu.edu

TRT (two-relaxation time) model would be your best bet. You can find some details in the following papers:

L.-S. Luo, W. Liao, X. Chen, Y. Peng, and W. Zhang.
Numerics of the lattice Boltzmann method: Effects of collision models on the lattice Boltzmann simulations.
Physical Review E 83(5):056710 (May 2011).

C. Pan, L.-S. Luo, and C.T. Miller.
An evaluation of lattice Boltzmann schemes for porous medium flow simulation.
Computers & Fluids 35(8/9):898-909 (September, 2006).

Both papers can be downloaded from my website

http://www.lions.odu.edu/~lluo/journal-pubs.html

-- LS Luo

guchaonemo

谢谢您，罗老师，我先看一下您的文章。非常感谢

guchaonemo

罗老师，您好，因为我最近在看您2001年那篇矩形网格的MRT模型，试着用C编写矩形网格的程序，我想问下，当a不为1的时候，声速如何确定，我模拟的是顶盖驱动流

luo@odu.edu

The speed of sound is determined by linear analysis as described in the paper.

I do not recommend the model with rectangular LBE. Instead, you can use interpolations, see the following paper:

A. Mussa, P. Asinari, and L.-S. Luo.
Lattice Boltzmann simulations of 2D laminar flows past two tandem cylinders.
Journal of Computational Physics 228(4):983-999 (March, 2009).

and some refs. therein.