收到了,谢谢李老师!!!@李东岳
gengchen
@gengchen
Best posts made by gengchen
Latest posts made by gengchen
-
RE: fluent用LES算水翼边界层选择中心差分引起的震荡如何消除
@李东岳 十分感谢李老师!对OPENfoam不太熟,抱歉耽误了这么长时间回复您,现在使用了三种差分格式 1. filteredlinear 过滤系数取0.5 2. filteredlinear 过滤系数取1 3. LUST格式,25%迎风+75%中心
1.
2.
3.速度场wiggle
可以看到中心差分随着过滤强度从0.5-1,耗散性会变大,导致靠近边壁的涡结构被抹平,逐渐向迎风格式的结果过渡,但是都不能解决震荡问题,(过滤掉了一些高频的震荡(碎涡),但是速度场的wiggle依然存在)。LUST格式下没有震荡,虽然只有25%的迎风,但仍然会造成边界层计算不准确,速度梯度不对,结果和二阶迎风类似。目前我认为应该是这个高雷诺数问题(5e5)本身无法用千万量级的网格算边界层,网格雷诺数太高,(我看有人说cell Re <2才能稳定),中心差分耗散太小,导致震荡难以消除。现在准备暂时放弃边界层定量研究了。。。
-
RE: fluent用LES算水翼边界层选择中心差分引起的震荡如何消除
@李东岳 谢谢李老师!我发现有个filteredlinear格式,准备用这个试试。先研究下您写的中文手册
因为LDV只测到2mm,模拟时取了4mm,其实用不着这么多,对比实验只用了2mm
-
RE: fluent用LES算水翼边界层选择中心差分引起的震荡如何消除
@李东岳 非常感谢李老师!如果方便的话,您可以把Openfoam的文件打包发给我吗,这个文件中转网站速度还行 https://cowtransfer.com/ 我看能不能改改差分格式,不方便也没关系,我这边摸索得比较慢,目前只能跑跑教程算例啥的。
-
RE: fluent用LES算水翼边界层选择中心差分引起的震荡如何消除
@李东岳 谢谢李老师的说明,我这边只试了中心差分Co<1的情况,确定是有震荡的,其他格式结果我还可以再减小时间步长测试下变化,您可以在确保稳定计算的情况增大下时间步长,简单看下foam的中心差分格式结果,有没有类似的不光滑速度等值面,或者破碎的震荡涡,不需要算那么长时间统计平均速度。有改善的话,我就全力研究争取用Openfoam计算。如果太麻烦的话就算啦,还是十分感谢李老师!
关于库朗数可以大于1的原因,我也不甚清楚,看有篇文献说测试平均库朗数到5还能保证计算结果,可能让我没有重视用小时间步计算。 -
RE: fluent用LES算水翼边界层选择中心差分引起的震荡如何消除
@李东岳 辛苦李老师!
1、图中是fluent统计的时间平均速度 2. 最大值为这条线上的最大值
我的计算是这样,图中是中心差分的结果,但是会有震荡。其他结果,比如有界中心差分,会使边界层剖面有整体的改变,二阶迎风能稍微好一点,但是在边界层的粘性底层和过渡层偏差会比较大
Umax差不多11m/s,时间步长的话,因为隐式格式为了快速看结果,我用的5e-5,试了从1e-5到5e-5,结果基本不变。