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

@fangyuanaza
$u_{\tau}$ 应该是定义在壁面上的一个量，因为这个量跟壁面剪应力紧密关联，内部网格是没有壁面剪应力的。
首先，壁面剪应力的定义为
$$
\tau_w = \nu\frac{\partial u}{\partial y}
$$
$\nu$ 是分子粘度。
使用有限体积方法进行计算时，所有的计算都是离散的，因此
$$
\frac{\partial u}{\partial y} \approx \frac{u_{internal}-u_w}{\Delta y}
$$
$u_{internal}$ 表示第一层网格的速度，$u_w$ 是壁面速度，一般是0。
如果壁面附近网格画得很密（yPlus < 1），那么 $\frac{\partial u}{\partial y}$ 与 $\frac{u_{internal} - u_w}{\Delta y}$ 足够接近，因此这种情况下，壁面剪应力可以用
$$
\tau_w = \nu \frac{u_{internal}-u_w}{\Delta y}
$$

但是如果壁面网格画得粗，那么 $\frac{\partial u}{\partial y}$ 与 $\frac{u_{internal} - u_w}{\Delta y}$ 之间的差异就会很大。 OpenFOAM 的做法是利用 $\nu_t$ 来修正分子粘度，从而使得
$$
\tau_w \approx (\nu+\nu_t) \frac{u_{internal}-u_w}{\Delta y}
$$

回过头来看，其实不管是近壁网格画得粗还是细，都可以用上一个公式来计算 $\tau_w$ ，只不过壁面网格很细（ yPlus<1)时，可以认为 $\nu_t = 0$（但是 $\nu$ 不会等于0，所以壁面剪应力不可能为0，$u_{\tau}$ 也不可能为0）。而对于比较粗的网格，$\nu_t$ 比 $\nu$ 大得多，所以也可以认为 $(\nu_t + \nu) \approx \nu_t$，这应该也就是李老师写的公式的来源。

@fangyuanaza

    forAll(uTau,patchi)

这里有问题，按照这样遍历，你其实遍历的是每一个体网格，假设你有1000个体网格，那么 patchi 的范围是0-999，但是你接下来的代码是把 patchi 当成边界的编号来用，肯定不行的。至于为什么出现 fluidThermo下面的报错，那个是因为湍流模型有着一个复杂的继承关系，可压缩的计算所使用的湍流模型就会用到 fluidThermo 类相关的东西。

@leilei

1. 试过一个方法可以做到，但是比较麻烦，那就是自己从源码编译 ParaView，然后编译的时候指定一个已经装好你想要的第三方模块的python给ParaView，这样得到的 pvpython 可以 import 你想要的包。https://discourse.paraview.org/t/how-can-i-install-and-import-other-modules-inside-pvpython/3067/6 这个链接提供了一些别的解决办法。
2. 目前知道的是使用 paraview 里面的 python shell ，可以补全。

@fangyuanaza 在 yPlus在openfoam代码里面的实现 中说：

u.boundaryField()[patchi].snGrad()

u.boundaryField()[patchi].snGrad()

加上 mag 试试，也就是改成

mag(u.boundaryField()[patchi].snGrad())

@桎梏
ParaView 里面调用filter-> Cell Size，计算出来网格的 Volume，然后再调用 calculator 来计算速度与 Volume 的乘积就可以得到了。

@chengan-wang
因为 新版本多了一个选项，skip zero time。勾选了这个之后，0时刻就不显示了。

@thomasshi
有两个问题：

1. 多孔介质力是体积力，所以不能用表面积分这种方式来算
2. 即使是用表面积分算曳力，那也不能只算p*normal，还有剪应力，需要根据wallShearStress来算。

@thomasshi
是合力，按照多孔介质阻力的公式，算出来多孔介质所作用的 cellZone 里面每一个网格的多孔介质力，然后按矢量加法全部加起来，最终得到一个合力矢量。

@thomasshi

不需要指定 patch
请看 forces.C 里面的 calcForcesMoment 函数， if (porosity_) 这个分支下面的内容，这部分是计算porous force 的方法。

@thomasshi

可以参考 forces 这个 functionObject，里面有计算多孔介质力的方法，前提是你的多孔介质力是用 porosityModel 来算的，比如 fixedCoeff或者DarcyForchheimer模型。