使用滑移网格求解旋转问题

zhouqishan

目前在做一个简单的算例，一个尖拱头部——圆柱体结合的导弹绕弹轴旋转，来流马赫数0.4，攻角9.5，转速设为200rad/s，使用moving mesh，设定弹周围小的计算域旋转，弹relative to adjacent  cell zone，速度为0，外层计算域不动，外层远场条件。使用隐式格式，湍流模型为RNG k-e，时间步长取为0.0001s。但是现在计算出的马格努斯力跟理论和实验结果差很多啊，而且马格努斯力是振荡的，貌似时间步取的越大振荡越厉害，是否取足够小的时间步能得到固定的马格努斯力？小计算域的选取大小有没什么要求吗？
希望有经验的各位大侠能不吝赐教，万分感激

[ 本帖最后由 zhouqishan 于 2011-7-14 10:08 编辑 ]

xrs333

在弹体坐标系上来看，弹体附近流动仍然是非定常的，流动现象丰富。
0.0001s的时间步长弹体旋转0.02rad，大约一周的1/300，似乎够小了，但是对于湍流的模拟来说还不够。

zhouqishan

原帖由 xrs333 于 2011-7-14 10:08 发表

                               
登录/注册后可看大图

在弹体坐标系上来看，弹体附近流动仍然是非定常的，流动现象丰富。
0.0001s的时间步长弹体旋转0.02rad，大约一周的1/300，似乎够小了，但是对于湍流的模拟来说还不够。

我感觉最终得到的马格努斯力应该是个固定值，但现在结果一直振荡，而且经常还达到正值。如果时间步继续取小的话那计算时间就有点太长了。可否开始计算时先取一个大的时间步，计算一段时间后再变小，最后得到马格努斯力？

xrs333

非定常，怎么会是定值。先理解流动，然后再作计算。你所谓的固定值，实际上是稳定飞行时的时均值。如果你只关心稳定飞行时的气动力，可以采用任何措施，包括取大一点的时间步长，来加快从计算初始条件到稳定飞行状态的计算过程。

zhouqishan

对，目前只关心飞行时的气动力。那最终的气动力应该怎么得出呢？取一个周期的平均值难道？

xrs333

把难道去掉，就对了。

zhouqishan

谢谢斑竹大大