圆柱绕流Cd偏小St偏大,救命啊~~~~

若海

之前发帖问过圆柱绕流网格划分的问题,当时的问题是涡脱落方向不朝向下游,而是随机的向左右偏;涡也不是左右依次脱落的. http://www.cfluid.com/bbs/viewthread.php?tid=113332&extra=page%3D2

重新划分网格后,网格如下图(一个象限的局部):圆柱绕流流场宽度7D长18D,流体水 Re55000. 近壁网格约为1/1400 D(圆柱直径),16层后通过两层1:3过渡层把网格尺寸变成1/140 D,距离壁面1/2 D向外逐步放大,到两侧墙壁约为0.2D.

湍流模型用k-e RNG,近壁处理,选中"differential viscosity model"和"pressure gradient effects",入口湍流参数 k0.00251,e0.00536(照书上的公式算的,入口湍流参数如果按默认的10%湍流强度和湍流粘性比10来算无法形成涡脱落.). report-reference里设了面积为D,密度998,深度1m,流速1m/s,圆柱直径0.055米.计算步长试过涡脱落周期的1/50(0.005s)和1/500(0.0005s).

现在的问题是出来的Cd和Cl都很小,Cd时均大概0.5,Cl最大瞬时值还不到0.5,根据涡脱落频率算出来的St大约0.38.跟别人论文里的数字都差的很远(别人算出来大约1.2,St0.2左右).现在又没方向了,该怎么改进呢?希望大家给些建议,多谢!!!

[ 本帖最后由 若海 于 2012-1-27 17:42 编辑 ]

若海

现在的涡街看起来好像比较正常,但是涡脱落频率与算出来的差别比较大

onesupeng

一般是格式耗散过大，我猜想你是Cd偏大，St偏小

onesupeng

另外，可能你计算的St是用Cd计算的，要用Cl计算，Cl的主频和涡脱落频率一致

若海

我算的是Cd和Cl都偏小，St偏大。我是用Cl的周期算St的，因为我的Cd结果波动不明显。

格式耗散过大是怎么回事？该怎么改善呢？多谢指导

[ 本帖最后由 若海 于 2012-1-21 22:38 编辑 ]

若海

重新划分网格,扩大了近区尺寸,取消了边界层到内区的比例过渡层,用k-w SST来算, Cd较接近理论值, Cl还很小,但是St=0.217,也算比之前接近很多了. 接下来还是得从网格划分着手.

flowinflow

很好的网格收敛的例子