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

@Zhoujq 把 fixedValue 改成 nutLowReWallFunction 试试？因为 omegaWallFunction 这个边界条件里面需要用到 nutWallFunction里面的东西，所以，如果你的 omega 壁面用的是 omegaWallFunction这个边界条件，那么你的nut就需要用一个继承自 nutWallFunction 的边界条件才行。 nutLowReWallFunction 这个边界条件其实就是把边界上的 nut 设置为一个非常小的数，跟你想要的是一样的。

@bestucan 在 fvc::DDt是什么意思？ 中说：

@xpqiu 是啊，我的签名就声明了我就是来搞笑的 滚来滚去，滚去滚来，薛定谔的滚蛋状态 布朗滚蛋运动

像您这样的能答上来的人没有出现，只有先搞笑一下，不至于别人问的问题没人搭理。还能刺激一下大家的回答欲。忘了在哪看的研究，一般不愿意回答的问题要是有个错误回答，就会引来很多纠错。因为纠错的成就感更强一些。

@strive_DUT
根据张量等式
$$
\nabla \cdot ( \mathbf{U} \phi) = \mathbf{U} \cdot \nabla \phi + \phi (\nabla \cdot \mathbf{U})
$$
所以，
$$
\frac{\partial \phi}{\partial t} + \mathbf{U} \cdot \nabla \phi = \frac{\partial \phi}{\partial t} + \nabla \cdot ( \mathbf{U} \phi) - \phi (\nabla \cdot \mathbf{U})
$$
这里等式右边的减号就是二楼贴出来代码里面减号的来源了。

另外，

@bestucan 在 fvc::DDt是什么意思？ 中说：

6.0版的定义。估计之前的难理解，又加了个U。

这个是来搞笑的么？

把 platforms 文件夹整个删了就行了。如果ThirdParty下面的你也编译了，那么把ThirdParty下面的 platforms也删了。

@波流力 是要每个圆柱的受力分开监控吗？

@波流力 找到 src/functionObjects/forces 目录，改造下面的 forces.H 和 forces.C，看懂他是怎么计算受力的，然后把有用的部分拿出来，写成一个新函数。需要处理的问题主要有两个，一个是把跟 functionObjects框架的部分剥离，另一个是用其他方式提供本来是有functionObjects框架提供的变量，主要是湍流接口，U场和p场。

@bestucan 然而很可惜你猜错了。omega_() 返回的是 omega_ 这个 volScalarField 类型的变量的 Internal 部分（也就是只有内部网格的值，不包含边界值），而且是带量纲。

@Tian 初始化一定高度的颗粒堆积在流化床底部试试？现在这个设置，侧面全是周期边界，入口固定速度，假设是单相流情形，因为没有壁面剪切力，速度场会是完全均匀的。针对多相流情形，因为你入口设定的固相体积分率也是固定值，所以这个流动里面一切都是那么均匀，没有任何会引起不均匀的因素，所以算出来是均匀的也可以理解。

@kimy 在 extra term twophasesystem.C 中说：
but the solver neglect it.

Is there any error message, like "cannot find dynamic library..." or "duplicate entry ..." ?

@Sloan 多相流的不太好估计，只能尝试不同的值看看了。