1. 网格质量如何？如果网格太差（checkMesh有严重报错，最大非正交大于85，skewness 超过20，等），建议回炉重造。
2. k 和 epsilon 的对流项离散格式用 linearUpwind，grad(k) 和 grad(epsilon) 用 cellimited，如果还不行，直接上 upwind试试。也可以先用upwind算一段时间再换成 linearUpwind。
3. 边界条件的设置取决于你使用的湍流模型，对于高雷诺数模型，入口 k 和 epsilon 给固定值，根据湍流度来估算；壁面上 k 用 kqrWallFunction，epsilon 用 epsilonWallFunction，nut 用 nutUSpaldingWallFunction。高雷诺数k-epsilon湍流模型还要注意壁面 y+，一般认为需要 y+>30，如果y+太小，结果反而不好，也可能会引起计算不稳定。

@金石为开 本质上是 hashTable 的使用。在基类里调用宏函数 declareRunTimeSelectionTable 创建了一个 hashTable，其 key 的类型为 “word”，value 类型为函数指针，这个函数指针指向一个返回派生类临时对象的智能指针（autoPtr）。

在派生类里，每次调用 addToRunTimeSelectionTable 的时候，相当于往 hashTable 里添加一组元素。

在求解器里，创建的是基类的智能指针，这个指针指向基类中定义的 New 函数的返回值。New 的作用是选择具体的派生类，这里涉及到C++的一个特性，即基类指针可以指向派生类。New 函数根据你在算例里提供的模型的名字（typeName），去 hashTable 里匹配，然后将匹配到的派生类的对象返回。这样，就实现了模型的选择。

具体参考：http://xiaopingqiu.github.io/2016/03/12/RTS1/

@子仲无未 用CAD画一个多面体，坐标要跟你的网格适配，保证这个多面体刚刚好包围了你想要的网格。然后将该多边形导出成 STL 格式，然后用 topoSet 将改多边形所包围的网格提取到一个 cellSet，topoSetDict大致是这样写：

{
 name c0;
type cellSet;
action delete;
source surfaceToCell;
sourceInfo
{
    file "constant/triSurface/block.stl";
   outsidePoints (( 1.012 0.01 0.902));
   //includeCut false;
   includeCut true;
   includeInside true;
   includeOutside false;
   nearDistance 0.0;
   curvature 0;
 }
 }

最后，在 setFieldsDict里，这样写：

defaultFieldValues ( volScalarFieldValue alpha 0 );

regions
( 
    cellToCell 
    { 
        set c0 ;
         fieldValues ( volScalarFieldValue alpha 0.60 ) ; 
    }
 );

用这种方法，可以做到对任意形状的网格区块设置初始场。

@myler

functions
{
    cellSource1
    {
        type         cellSource;
        functionObjectLibs ("libfieldFunctionObjects.so");
        enabled         yes;
        outputControl   outputTime;
        log             yes;
        valueOutput    yes;
        source          all;
        operation       max;

        fields
        (
            T
        );
    }
}

这个加到 controlDict 里就能输出最大值了。

其实还有一种办法，既然你说在终端的输出信息中包括了最大和最小温度，那么如果你保存了终端输出信息，你还可以用 linux 的工具来从那些众多的信息中将你想要的最大和最小信息提取出来，
比如，假设你的文件名为 log

cat log | grep "max(" | sed 's/^\ *//g' | cut -d ' ' -f 5

感兴趣的话可以试试这个命令，可能不能完全适合你的情况，你需要做些修改。

PS：一直不清楚为什么在.bashrc里的那4句注释掉就可以，之前从来没有对这4句太在意，而是直接添加OF的环境变量就好了。。。

这里解释了为什么：http://www.evernote.com/l/AYvYUlXKzfxBIZ5sJbyk-cREUR1SGcAPIto/

@wwzhao 在OpenFOAM-2.1的时候，程序里写的就是“CrankNicholson” ，所以controlDict里也应该写“CrankNicholson”，虽然这个拼写确实是不正确的。但是从OpenFOAM-2.3开始，这个名字改成了正确的拼写：“CrankNicolson”，相差一个字母。所以这就导致了2.1和2.3的算例会在这里不兼容。

常用出入口边界条件：
https://cpp.openfoam.org/v4/a11022.html

壁函数：
https://cpp.openfoam.org/v4/a11020.html

每一个边界条件都附有简短的说明以及使用示例。

@random_ran
是 GAMG 在拖你的后腿，这么多核的情况下，你可以试试用 PCG + DIC 来求解 p 方程，这个方法并行效率会高一些。

@qjh888 在 OpenFOAM中kOmegaSST和yPlusRAS 中说：

yPlus要在30～200之间才可以使用wallFunction,

这个说的是标准壁函数，OpenFOAM中的omegaWallfunction是混合的，所以理论上对omega的边界处理是可以适用于较大范围的yPlus的。但是其他的湍流变量并没有这种混合，此外，yPlus太大，边界处网格必然变粗，精度也会相应下降。所以，yPlus取多少还是要根据实际情况来决定

还有一个问题，为什么计算yPlusRAS一定要使用nutWallFunction?

那是因为yPlus的计算是依赖于你所使用的湍流模型以及nut的边界条件处理的，所以，yPlusRAS需要调用nutWallFunction类中的一个函数来计算 yPlus

@lixh 在热物理属性那个文件里设置（thermoPhysicsProperties.air）。是否需要设置密度和年度，取决于你的热物理模型的选择，如果equationOfState设置的是 perfaceGas，那么，流体将遵循理想气体状态方程，这时不需要设置密度，因为密度是根据压力算出来的。如果设为rhoConst，这时需要指定流体的密度。粘度的设置取决于 transport 的设置，如果设置为 const，即粘度为常数，则需要指定粘度的值，如果选择其他粘度模型，则会有不同的设置，具体可以参考UserGuide的 thermoPhysics 这一章。

@Lisiye https://cpp.openfoam.org/v8/ 可以去这个页面搜索。到最下面选择与你最接近的版本，然后右上角的搜索框输出关键字来搜索。

@Samuel-Tu

用 feature edges ，不过要比较新版的 paraview 才有这个，我截图的是 paraview-5.6.1 。

@Samuel-Tu 你去研究一下 foamDictionary 这个命令的用法，绝大部分对openfoam 的dict 文件的修改都可以用这个命令来进行。
另外，sed 也有很多技巧能实现这种替换。你的需求可以这样来分解，首先找到 inlet，然后替换 inlet 下面的 type，如果用sed，则可以

 sed  -i '/inlet/,+2s/wall/patch/' boundary

命令解释：

'/inlet/,+2'

是一个范围，表示从匹配到 inlet 的行开始到它之下的2行，s是命令，表示替换，也就是在前面定义的范围内进行替换。/wall/patch/ 表示实际的替换操作，这里是把 wall 替换成 patch。最后的 boundary是文件名。

WARNING: 注意到 sed 之后的 -i 选项，这个选项用了之后，会直接修改 boundary 文件的内容，所以除非你确信你的替换命令是正确的，才建议直接用 -i 选项，如果仅仅是想看一下替换的效果，去掉 -i，这样命令会把替换之后的效果打印到屏幕上，不会直接修改文件。

抛砖引玉。

@Samuel-Tu 不强求法向指向外面，但是建议做成这样。如果法向是混乱的，有的朝里，有的朝外，有可能会导致奇怪的问题。

paraview 导出的 stl 可以用于 shm。

@Samuel-Tu 可以，用 surfaceMeshTriangulate

@莫晓柒 编译之后得到的库和求解器程序都在 platforms 下面，所以只能重新安装了。

是的。

@DY大世界 某个方向设置为 empty，那么这个方向的流场就不会再计算。以前的版本如果empty方向的网格不是固定层数，是无法计算的，checkMesh 也会报一个提示：

 ***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D.

经测试，现在新的版本可以计算，但是结果我没有验证过。从计算的角度来讲，即使这样算出来结果是合理的，也增加了很多没必要的计算，因为这样处理，网格数量会比单层的2维网格多很多。
另外，假设举个例子，一个方形计算区域，中心部分用snappyHexMesh进行了加密，然后x 方向设置为 empty ，那么你计算完的结果，在x 方向上不同位置的截面上流场可能是不一样的，因为网格数量可能不一样，这样显然跟二维计算的概念是不符合的。

@DY大世界
extrudeMesh 之后，是否运行了 checkMesh 来检查网格质量？
不确定你的流程是否对，可以参考一下 /tutorials/incompressible/pimpleFoam/RAS/wingMotion 这个算例。这个是用 snappyHexMesh 生成 2d 网格的例子。

这两种应该都不对
snappyHexMesh会在三个方向上都加密。对于二维网格，要求是 empty 方向上只有一层网格（或者层数要是一样的），用snappyHexMesh加密，empty方向上的网格层数可能在不同位置会不一样。你可以截图看一下，x，y，z，三个方向分别做一个slice，看看是不是我说的那样。如果empty方向上的网格不是一层，那么二维计算是有问题的。
为了得到empty方向上只有一层网格，有一个办法是先按照你描述的第一种方法用snappyHexMesh生成网格，然后，再用 extrudeMesh 来把某一个 empty patch 拉升一层网格，这样得到的就是常见的二维网格了。