openFoam4.0 学习记录(七):Tutorials之Lid-driven Cavity Flow (后续3-高雷诺数RAS计算)

5月21日刚刚回家,明天爷爷生日,在此祝他福如东海!

今天按计划完成空腔流教程的最后一项——高雷诺数。在此之前,我所做的都只是低雷诺数下的层流。但在现实中,理想的层流是不存在的,更多的是湍流。这块硬骨头,最是难啃。话不多说,我按着userguide做一遍。

这一篇的关键词:pisoFoam,RAS,k − ε 模型。

首先将教程里面的案例复制到run目录内:

run
cp -r $FOAM_TUTORIALS/incompressible/pisoFoam/ras/cavity cavityRAS
cd cavityRAS

先用blockMesh生成网格。在这里,mesh grading 不是必要的,因为我们用的是k − ε 模型。

paraFoam看一看,就是很普通的网格。

openFoam4.0 学习记录(七):Tutorials之Lid-driven Cavity Flow (后续3-高雷诺数RAS计算)_第1张图片看一看,是很普通的网格。

打开0/nut文件:

openFoam4.0 学习记录(七):Tutorials之Lid-driven Cavity Flow (后续3-高雷诺数RAS计算)_第2张图片

注意边界的type“nutkWallFunction”,还有其他的keyword,这些都要通过深入地研究源代码才能,并且有相关的理论基础才能很好地理解。这些会在日后慢慢写,慢慢学,我还有很长的路要走。

在constant/turbulenceProperties的文件里:

openFoam4.0 学习记录(七):Tutorials之Lid-driven Cavity Flow (后续3-高雷诺数RAS计算)_第3张图片

simulationType有三个选项可选:lamilar,RAS和LES。在RAS的子程序里,默认是选择k − ε 模型,即kEpsilon;

输入pisoFoam,运行计算:

openFoam4.0 学习记录(七):Tutorials之Lid-driven Cavity Flow (后续3-高雷诺数RAS计算)_第4张图片

花了半分钟左右,最后用cat log查看:

openFoam4.0 学习记录(七):Tutorials之Lid-driven Cavity Flow (后续3-高雷诺数RAS计算)_第5张图片

停在了10秒,即我们设定的endtime,收敛成功,看到end字样。

最后看看结果:

openFoam4.0 学习记录(七):Tutorials之Lid-driven Cavity Flow (后续3-高雷诺数RAS计算)_第6张图片速度场

openFoam4.0 学习记录(七):Tutorials之Lid-driven Cavity Flow (后续3-高雷诺数RAS计算)_第7张图片压力场

可以看出,高雷诺数下的结果跟低雷诺数有很大的不同。

小结:

我想通过看源代码来理解程序,但是一方面c++代码基础差,另一方面流体的知识也差,今后会两手抓,通过对icoFoam等代码的学习,体会CFD的美丽。

 

 

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