openFoam4.0 学习记录(六):Tutorials之Lid-driven Cavity Flow (后续2)

今晚想将第一个教程后面的部分做完。之前做仿真的时候雷诺数用的是10,非常的低,现在想要增加到100,看看结果如何。

先到run的文件夹里,如何克隆一下cavity的文件,成为新的文件:

run
foamCloneCase -latestTime cavity cavityHighRe
cd cavityHighRe

第一步:pre-processing

在constant文件夹内,打开transportProperties,因为雷诺数从原来的10增加到100倍,所以transportProperties里面的nu(运动粘度kinematic viscosity)改成原来的0.1倍,原理看下面的公式,那个像斜着的v一样的就是nu,运动粘度。

Re=\frac{d\left | U \right |}{\nu }

从0.01改成0.001,大功告成!

另外,我们这个案例是从0.5秒开始的,要修改一下controlDict文件:

openFoam4.0 学习记录(六):Tutorials之Lid-driven Cavity Flow (后续2)_第1张图片其他不变,endTime改成2秒,意思是从0.5s开始,用一开始从cavity的case里面的0.5s的各变量的值作为本次case的初值,直到2秒。

第二步:运行run!

输入icoFoam运行,报错——》忘记将startTime从0改到0.5啦!

再来:

openFoam4.0 学习记录(六):Tutorials之Lid-driven Cavity Flow (后续2)_第2张图片收敛成功!

paraFoam打开看看结果:

openFoam4.0 学习记录(六):Tutorials之Lid-driven Cavity Flow (后续2)_第3张图片这是速度U在0.5s时的状态

openFoam4.0 学习记录(六):Tutorials之Lid-driven Cavity Flow (后续2)_第4张图片两秒后,因为雷诺数增大,速度场不再对称。

openFoam4.0 学习记录(六):Tutorials之Lid-driven Cavity Flow (后续2)_第5张图片用glyph插件将上图改成矢量图。

openFoam4.0 学习记录(六):Tutorials之Lid-driven Cavity Flow (后续2)_第6张图片另外这是压力场

openFoam4.0 学习记录(六):Tutorials之Lid-driven Cavity Flow (后续2)_第7张图片用等压线显示是这样的,可以看出其不对称!

小结:

提高雷诺数是通过修改nu来实现的,对于这个例子而言。

虽然这次的雷诺数不高(100),但也可以看出粘度对流体的影响了!

下一次,要试试更高的雷诺数,并使用RAS方法对其进行湍流模拟!

你可能感兴趣的:(openfoam,userguide)