【新手贴】为什么Fluent模拟类二维圆柱绕流问题得不出卡门涡街呀?
浏览:16748 回答:9
我用
Fluent二维双精度求解器,General里面的Time项设置为瞬态,来模拟二维类圆柱绕流问题,Re=10000,可是为什么模拟不出来卡门涡街的图形?(我研究的是潜艇升降桅杆横截面的水动力性能问题)
模型尺寸:特征长度D=0.586m
计算模型:由于Re=10000<20000,是层流范围,所以选择Laminar层流模型,设置保持默认
边界条件:
入口:速度入口条件,入口速度V=0.0171m/s(是按Re=10000反算得到的),其他设置默认
出口:压力出口条件,其他设置默认
上下边界:滑移壁面,滑移速度与来流速度一样,其他设置默认
物面条件:固定壁面,其他设置默认
Solution Method:用的是PISO,压力设为Second Order,动量设为Second Order Upwind
时间步长:通过0.2=Sr=f×D/V,反算出涡脱落频率f=0.0061Hz,脱落周期T=1/f=163s,所以取时间间隔Δt=T/1000=0.163s=163ms
我模拟过一次Re=15000(层流),来流速度V=0.0257m/s,设置全都与上述一样,可是这一次却可以模拟出很漂亮的卡门涡街。见图
可为什么同样是层流的Re=10000,同样的设置,这次却模拟不出涡街?而且其他很个Re数值我也模拟过,湍流的时候湍流模型用过RNGk-ε,全都没出现过卡门涡街(具体模拟结果见下图)。按道理Re>3.5×10^6的时候,不是也会出现湍流涡街吗?
这是怎么回事?请高手帮帮忙!我现在是百思不得其解……
模型尺寸:特征长度D=0.586m

计算模型:由于Re=10000<20000,是层流范围,所以选择Laminar层流模型,设置保持默认
边界条件:
入口:速度入口条件,入口速度V=0.0171m/s(是按Re=10000反算得到的),其他设置默认
出口:压力出口条件,其他设置默认
上下边界:滑移壁面,滑移速度与来流速度一样,其他设置默认
物面条件:固定壁面,其他设置默认
Solution Method:用的是PISO,压力设为Second Order,动量设为Second Order Upwind
时间步长:通过0.2=Sr=f×D/V,反算出涡脱落频率f=0.0061Hz,脱落周期T=1/f=163s,所以取时间间隔Δt=T/1000=0.163s=163ms
我模拟过一次Re=15000(层流),来流速度V=0.0257m/s,设置全都与上述一样,可是这一次却可以模拟出很漂亮的卡门涡街。见图

可为什么同样是层流的Re=10000,同样的设置,这次却模拟不出涡街?而且其他很个Re数值我也模拟过,湍流的时候湍流模型用过RNGk-ε,全都没出现过卡门涡街(具体模拟结果见下图)。按道理Re>3.5×10^6的时候,不是也会出现湍流涡街吗?

这是怎么回事?请高手帮帮忙!我现在是百思不得其解……
这是Re=100000的速度分布图:
第二章图因为我已经看到了后面可以稳定出现卡门涡街的情况了,所以就提前停止计算了。所以可以看到远离物面的尾部的涡比较不规则,散得比较开。如果计算完了,应该会只剩很整齐的卡门涡街了。
解决方法其实很简单,就换成LES模型就好了。其他参数全部不变。