k-omega湍流模型模拟二维圆柱升力系数问题

[ycx1042341217]

大家好，我最近在模拟Re=10000时，二维圆柱的升力系数和阻力系数。采用的湍流模型为k-omega，网格划分采用非结构化网格，边界层采用四边形网格，圆柱壁面y+值在1附近。换过不同的网格，时间步长等，阻力系数为1.0-1.2左右，与文献可以对应上，但是升力系数模拟出来在1.25左右，文献数据为0.4-0.6附近。大家能帮着看看问题出在哪里嘛？fluent里再选择k-omega模型后还可以选择低雷诺数修正或者自由剪切修正，我在计算的时候并没有勾选，大家知道这个一般怎么用吗？

[李东岳]

算不准这个太不好判定了，不过帮你找一个专家 @random_ran

[ycx1042341217]

谢谢东岳老师

[ycx1042341217]

@random_ran 老师可以帮忙看看吗？我后面也换了k-omega SST模型，结果也相差不多，看了一下也并没有漩涡脱落产生。

[random_ran]

升力系数 0.4-0.6

Re 10000 圆柱绕流的流场特征应该还在亚临界区域内。时均 Cl 应该是0才对。我觉得楼主说得是升力系数的波动值 fluctuating lift force（Cl'）， 你看实验中的数据都有很大的不确定性。用3D网格，正交网格，效果会好很多。

M. M. Zdravkovich, “Flow around circular cylinders vol. 1: Fundamentals,” (Oxford University Press, New York, 1997).

[ycx1042341217]

@random_ran 好的谢谢老师