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 这一章。

@vbcwl
我用的 openfoam.com 的版本，这个问题可能是 of9 的bug。

@xiaoyang-luan
RNG k-epsilon 预测不了转捩。
转捩相关的湍流模型，参考：https://turbmodels.larc.nasa.gov/
“Turbulence+Transition Models ” 部分

@vbcwl
哦，不好意思刚才看错了。-0.06788236085 这个数字很奇怪，没看出来是怎么来的。我试了一下也没复现出来。

@xiaoyang-luan
所以如果需要模拟 层流到湍流的转捩过程，需要特殊的湍流模型，一般的 k-epsilon， k-omega 模型确实是假定整个场都是湍流状态。

@vbcwl

那是因为你打开来查看的那个核，恰好没有包含任何outer这个边界的面单元。 nonuniform 0() 表示一个空的list。

@chen_hao 在 关于修改湍流模型中非线性雷诺应力项的问题 中说：

您所指的加上雷诺应力非线性部分就是UEqn.H中turbulence->divDevRhoReff(rho, U)项的修改对吗？

不是，不需要修改这个函数，你仔细看看 ShihQuadraticKE 这个模型的继承关系，它继承的是 nonLinearEddyViscousity 类。这个类计算雷诺应力的方式不一样，除了线性项还有一个 nonLinearStress 项。你只要按照这个框架去写，把你的非线性项赋值给 nonLinearStress 就行了。

@chen_hao
你可以看看 src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/ShihQuadraticKE 这个湍流模型
这是一个非线性的湍流模型而且也只考虑了二阶非线性项，跟你要实现的那个类似。
这类模型除了修改 $\nu_t$ ，还需要把雷诺应力的非线性部分加上。

@vbcwl
哦，如果是这样的话，那只能通过你描述的
$$
y^+ = \frac{y u_\tau}{\nu}
$$
的方式来算 $y^+$ 了。postProcess 这种计算方式调用的是湍流模型的东西，适合用湍流模型的场景。

@vbcwl
请问你计算的时候，nut 的壁面边界条件设置的是什么？
postProcess 算出来的 yPlus 是调用 nut 的边界条件里面定义的 yPlus 函数算出来的。

@hitsc30
看你的离散格式，对流项 div(phi,U) 使用的是 linearUpwind，这个格式不适合用于 LES 仿真，数值耗散太大。

建议尝试一下 LUST，或者 filteredLinear2 这类格式。

还有，如果用 PISO 求解器，p 不要用松弛。你现在用的是 0.3，这个是错的。