对时间和网格无关性的研究

[crispy]

背景：
二维圆柱绕流

网格拓扑结构：

时间无关性的研究：

对于同一套网格，阻力系数随时间步长的变化如下图所示。结果表明，当最大库朗数的值不过分小时，时均阻力系数随着时间步长的减小有收敛的趋势。当最大库朗数的值的量级到1e-2时，时均阻力系数发散。

• 疑问：是否可认为在openfoam的计算中，并不是时间步长越小，计算结果越精确。
  

网格无关性的研究：

• 主要针对y+
  对于相同时间步长，阻力系数随y+的变化如下图所示。结果表明，当y+的值不过分小时，时均阻力系数几乎收敛到同一个值。当y+远小于1时，时均阻力系数发散。
  

结论：

根据目前的结果来看，y+过小或者dt过小都使得阻力系数发散。

• 疑问：是openfoam针对小的y+和dt有什么额外的设置吗？应该如何理解目前的结果呢？谢谢各位老师同学。

[crispy]

@李东岳 李老师，这是我换v2306的结果。

[李东岳]

你的nut用的什么壁面函数

[crispy]

@李东岳

    front     
    {
	type            empty; 
    }
   
    back 
    {
	type            empty;  	
    }

    wall
    {
	type            nutkWallFunction;
        value           uniform 0;
    }
	
    inlet
    {
	type            calculated;
        value           uniform 0;
    }
	
    outlet
    {
	type            calculated;
        value           uniform 0; 
    }

[李东岳]

这个一旦进入湍流那面debug的话，就很麻烦。你能不能做一个层流的，先单独验证时间步长的影响。

[crispy]

@李东岳 老师，层流的话雷诺数定在多少呢？网格还需要换吗？

[李东岳]

你用laminar模型算多少雷诺数都可以，这叫quasi-DNS。不需要换网格，要固定一套网格跑。但是你用laminar跑需要保证y+足够小。你可以先试试你的y+=0.2的网格。

所有参数都固定，单纯的变时间步长来看结果。

[crispy]

@李东岳 好的老师，我算完发上来。

[crispy]

@李东岳 老师我遇到个问题，使用y+=0.2的网格和laminar模型进行计算，目前算的1e-4、5e-5、2e-5、1e-5的流场都没有收敛。这是为什么呢？是网格不合适吗？

[李东岳]

@crispy

库朗数多大
流场没有收敛是什么意思，发散了么

[crispy]

@李东岳 圆柱尾涡不是周期性的。maxco从1.91——0.19。

[李东岳]

@crispy 0.19那个格式发图看一下Cd

[crispy]

@李东岳

[李东岳]

也不能说算的不对。我觉得这个Cd应该把高频过滤掉。

不过从时间平均的结果来看，Cd的值正确么。

[crispy]

@李东岳 不对的，实验结果约为1.2。

[李东岳]

@crispy 那你看其他库朗数的结果怎么样，1以下的，是时间步长越大，cd越好么

[lxz ii]

@crispy 大佬，请问你验证的雷诺数是多少呢？现在有比较好的结果了吗？