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看一看，就是很普通的网格。

看一看，是很普通的网格。

打开0/nut文件：

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

在constant/turbulenceProperties的文件里：

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

输入pisoFoam，运行计算：

花了半分钟左右，最后用cat log查看：

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

最后看看结果：

速度场

压力场

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

小结：

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