MRT模拟初始化的问题

nanosmile

我用MRT做一个圆柱绕流的模拟

采用D2Q9格子划分

对分步函数初试化的时候，我是用

        
f(0)=f(eq)  （密度=1，u=0，相当于f(0)=w)

m(0)=M*f(0)

得到的m初始值，之后进行演化过程



但结果流动只能出来层流的状态，流体绕过圆柱后不能产生波动

我试着改变几个松弛因子，还是不行。

想问大家是不是我的初始化有问题？

这种问题用什么方式初始化比较合适？

[ 本帖最后由 nanosmile 于 2011-4-9 20:34 编辑 ]

luo@odu.edu

What is your boundary conditions and the Reynolds number? Normally what you need to do in the initial condition is to include some "randomness" (in either density or velocity field) to instigate the instability.

nanosmile

谢谢Luo教授
我用rand（）函数生成的随机数初始化速度场，（我用c编的计算程序）
果然流动状态和初始流场有很大关系，而且现在Re的作用也比较明显。
之前的程序是用的标准反弹边界，出口是充分发展流动出口，入口速度是按Poiseuille流动给的。按rho=1，u=0，初始化分布函数。即使我把Re给的很大，Re=2000的时候，也只能得到层流流动的结果。见图。
若果Re<1000，流动完全是层流的，流线贴着障碍流过。

[ 本帖最后由 nanosmile 于 2011-4-11 12:59 编辑 ]

luo@odu.edu

It is important you read some papers/books in the literature. You won't see vortex street (if that's what you are expecting) with the cylinder placed in a small channel.

luo@odu.edu

PS There is a well known benchmark case: a cylinder asymmetrically placed in the channel:

Sch&#228;fer, M. and Turek, S. (1996) ‘Benchmark computations of laminar flow around a cylinder’, in Hirschel, E.H. (Ed.): Notes in Numerical Fluid Mechanics: Flow Simulations with High Performance Computers II, Vieweg, Braunschweig, Germany, pp.547–566.

You may want to try that. I include a work we did a few years ago for this case.

nanosmile

谢谢Luo教授
    这几天我根据您的建议对程序进行了修改，已经可以得到一些比较理想的结果。
因为刚开始接触LBM，理解上还停留在初级阶段，所以很多看似简单的问题，对我来说理解起来都很费劲。
觉得自己还是应该多加强些基础学习，多看书和相关文献。
    刚来这个论坛的时候，看到几乎每个问题您都会回复，还给一些良好的建议。现在像您这样，愿意跟
初学者分享经验的老师，真的不多啊！
    非常感谢您的意见，谢谢！！！

大名青年奥巴牛

果断MARK