Cd计算不准的问题-LES-Re3900-三维圆柱绕流-pisoFoam



  • @random_ran foam-extend/3.1,广州天河2,用了864核算的,10s左右才能更新一次时间步,deltaT=0.001,效率不是很高,我画的是矩形区域O型网格,网格数量5千多万。计算300s,需要50多天时间,资源消耗太多,看别人的文献DNS圆柱绕流Re=3900资源消耗并不是很大,不知问题出在哪里;我现在准备换成圆形区域网格计算


  • OpenFOAM讲师

    @李东岳
    我是免费用户,那个网站很赞的地方是很多用户把模型都共享了。

    @bingningmeng45
    Hi,多谢分享细节,这个计算套装很酷。感觉你的网格可能画得太细了。

    稍稍满足满足y+的要求就可以。

    我也更倾向用圆形计算域。

    另外可以试试 OpenFOAM-dev。 大规模的计算,过多的processor文件会让你头痛的。

    另外分解计算域的方式,不同scheme的选择,不同求解器的选择,不同的硬件,不同的节点选择,… 太多东西可以影响计算时间。

    目前,我看到有人专门研究这些的是这篇文献:
    S. Keough,2014,“Optimising the Parallelisation of OpenFOAM Simulations”。

    但是计算量也不是很大,供你参考。



  • @random_ran
    Mean streamwise velocity along the centreline of the cylinder这幅图怎么进行z轴(圆柱轴向方向)平均?还有对比其他的量也需要z轴平均。需要在算例怎么设置


  • OpenFOAM讲师

    @bingningmeng45

    方法有很多。

    我自己的做法是 用 perl 生成 python 脚本来作轴向平均。

    再配合MATLAB做最后的平均以及画图处理。



  • @random_ran 你这个做法我不是很懂,能详细说明一下吗,谢谢


  • OpenFOAM讲师

    @bingningmeng45

    就是把ParaView中的操作脚本化。然后对生成的脚本做一些“胶水”处理。我这里用的Perl,你可以用自己熟悉的编程语言。

    获得*.py 脚本的内容这样操作,在ParaView GUI 里面找到:

    Tools - Start Trace

    进行你的图形操作,直到输出数据文件为止。
    

    然后 Stop Trace,就生成了 *.py 的脚本。

    接下来对这个脚本文件处理就好。这个脚本文件中绝大部分都是一样的,只要你在感兴趣的地方稍作修改就好,比如储存地址,切片位置。

    用Perl的目的就是少了做了几次“复制”“粘贴”。

    最后生成一个“真正”的脚本,执行

    pvpython *.py
    


  • @random_ran 0_1514879810615_QQ图片20180102155602.png
    我现在正在运行一个圆柱绕流Re=3900的算例,在处理一维波数谱图(1D wavenumber spectra)中遇到问题,不太清楚下图横纵坐标如何处理得到。
    纵坐标是u’u’么?
    横坐标是波数k_z,如何得到的?大部分的文献说的很笼统,并不能很清楚知道他们的计算处理过程。
    我现在的做法是沿圆柱轴线z方向,截取一条线上的流向(x轴)速度波动,横纵坐标不太会处理,一直纠结在这。


  • OpenFOAM讲师

    我从维基里面扣下来的一段话:

    “从随着时间而变的函数萃取出的一组数据,经过傅里叶变换,会得到一个频率谱;而从随着位置而变的函数萃取出的一组数据,经过傅里叶变换,会得到一个波数谱。”

    相信对于速度在时间历程上进行快速傅里叶变换(FFT),应该很熟悉。时间域内(时间, 某个方向上的速度大小) --> (频率,某个方向上的速度大小); 对应的 位置域内(位置,某个方向上的速度)–> (波数,某个方向上的速度大小)。 [位置域内是我瞎造的], 不过量纲上是讲的过去。 时间FFT后 横轴变换成 s-1; 位置(横座标)FFT后变成 m-1 也就是波数的单位。

    我的感觉是作者取了4条直线。 每条直线上搜集的是沿流向方向上的速度的波动值 (u’u’) [对应雷诺应力],然后以 (纵座标,u’u’) 做FFT。

    但是为什么这样做呢?


  • 网格教授 OpenFOAM教授 管理员

    @random_ran
    你这个文章是不是已经投出去了?接受没呢



  • @random_ran 你好,请教nut和nuTilda边界怎么设置



  • @random_ran @yhdthu 我对于那条蓝色的cd曲线还存在问题,

    这个Cd是平均的整个Cylinder的表面。

    这个是如何平局得到的呢?
    我的理解是计算输出每一步的结果,从某个时刻开始,开始对瞬时值做时间平均,是这样的吗?