[求助]LBM方柱绕流问题

nieinhz

    现在在用lbm做方柱绕流，采用郭照立提出的不可压模型，边界条件也是他的非平衡外推方法。
    计算中粘度系数用以下表达式
         niu=(2*tau-1)/6*(c*dx)
    方柱绕流的Re数定义为：
         Re=U*L/niu，计算中U＝1，L＝1（U为来流速度，L为方柱边长）
    计算区域为：  0<=x<=60L     0=<y<=8L
    在计算中我固定Re（100）、tau（0.65）和dx（1/30），然后确定c，开始我用tau＝1.6，然后发现Re最大只能到30，后来我设tau＝0.65，我发现可以计算出Re＝100，但是不知为何，始终得不到涡脱落的结果，仍然只是定常的对涡现象，不知道问题出在哪里。我想请问一下，无论设置tau和c为多少，只要保证Re数一定，是否计算的流场也应该一定？
    谢谢！

wsl2006

我算的是圆柱绕流，同样存在这个问题
即使Re>100,算不出来Votex Shedding现象  [br][br][以下内容由 wsl2006 在 2006年12月07日 07:11pm 时添加] [br]
不知道哪位大侠有高招，帮帮忙呀