请教：为什么没有涡脱落

glandetian

我用间断有限元方法计算可压缩NS方程，圆柱绕流问题，来流马赫数为0.38，雷诺数为200，计算了很久，为什么只能在圆柱背面得到一对对称的涡，而没有涡脱落呢？

    请高手指点啊，谢谢了！

onesupeng

原帖由 glandetian 于 2013-3-15 06:04 发表

                               
登录/注册后可看大图

我用间断有限元方法计算可压缩NS方程，圆柱绕流问题，来流马赫数为0.38，雷诺数为200，计算了很久，为什么只能在圆柱背面得到一对对称的涡，而没有涡脱落呢？

    请高手指点啊，谢谢了！

如果你程序正确的话，耗散过大，或者网格过粗

通流

也许你的程序和网格做的太仔细了，连误差都是对称的？

glandetian

谢谢上面两位指点！

我按下面的做法，得到了脱落涡：首先，将网格逆时针旋转0.5°，来流条件不变，还是沿着x轴方向；计算一定时间后，停止程序，将计算结果写成txt文件；再次启动程序，网格不旋转了，并读入刚才的文件，继续计算；计算一段时间后，输出计算结果，可以得到脱落涡。

我的理解是，我上面的做法相当于是先有来流吹过，过一段时间后，来流方向偏转了0.5°，给了点扰动，所以就有涡脱落了。

不知道我这样做行不行。

glandetian

如果壁面处网格没有细化，会得到几乎跟求解Euler方程一样的结果。

glandetian

我一开始也是这么想的，呵呵