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

@chszkc 在 关于postProcess得到张量的排序 中说：

想请问一下各位OFer，通过命令

postProcess -func "grad(U)"

得到的速度梯度张量共有九个元素，打开数据文件后发现这九个元素是以一维数列的形式保存的，那这九个元素的排列顺序是什么呢？Google了大半天都没找出来.....
之前自己默认顺序是dudx,dudy,dudz,dvdx,dvdy,dvdz,dwdx,dwdy,dwdz，但算出来的结果有点怪，所以怀疑自己是不是在排列顺序上理解错了

应该不是这个顺序，而是转置过来这个顺序

@Micro
好久没看这个了，具体不记得了。这个库在src/finiteVolume 和 src/OpenFOAM 下面做了很多改动，尤其是涉及矩阵的存储和运算的。

@aiweimo
我买过完整版，需要的留邮箱

@aiweimo

UEqn.relax() ，即使松弛因子是1，其实也会对 UEqn 的 系数矩阵进行一些调整，以使它更满足对角占优。fvSolution 里面如果没有定义松弛因子，则默认值为1。

@aiweimo

k-w方程是低雷诺数方程，本身不需要壁面函数

严格来说，应该是 $k-\omega$ 方程可以一直算到 $y^+<1$，但是，$k-\omega$ 方程也可以结合壁函数来算$y^+>30$ 的网格。相比之下，标准$k-\varepsilon$模型的$\varepsilon$ 方程则只在log-law region之外才严格有效，不可直接用来算很小 $y^+$ 的网格。
OF 里面的那些湍流壁面边界条件，有一点容易引起误解的是他所有的名字都带 wallFunction，而大部分人看到 wall function，想到的都是标准壁函数，也就是 log-law region的 $U^+=\frac{1}{\kappa}\ln(Ey^+)$ 以及层流底层的 $U^+=y^+$ 。 实际情况下原比这个复杂，因为很难保证整个region 都满足 $y^+ > 30$ 或者 $y^+<1$，所以，对于实际问题，湍流模型的壁面处理需要考虑如何处理$1<y^+<30$ 的区域。文献中有很多相关的研究，openfoam 的 omegaWallFunction 其实就是一种，里面构造了一个blending function 把适合 $y^+ > 30$和 $y^+<1$的 omega 公式blend一下，来让omega 在任意 $y^+$ 下都能算出来一个相对合理的值。
所以，OF 里面那些 wallFunction 的边界，可以理解为定义多种不同的 wallTreatment 方式，具体怎么用，需要结合着实际情况来。2 楼提到的就是一种常见的用法。文献中壁面处理的方式太多太多了。。。

至于 $\omega_{wall}$ 指定为零梯度或者固定值，这种做法不太合理。因为 $\omega$ 在壁面附近需要满足某种渐近特性（asymptotic），零梯度或者固定值都不满足这个条件，虽然可以可以算出来结果，但是不太符合物理。

@leilei

1. 直接用 p 就可以，OF 里面的 pressure function object 主要作用是用来算总压，静压系数和总压系数。
2. 之所以会出现 p < 0, 是因为 interFoam 求解器是基于不可压缩假设的，也不求解能量方程。因此在 interFoam 的方程里面，起作用的不是压力的绝对值，而是压力梯度。一般情况下，不可亚求解器都至少会有一个边界上的 p 设置为第一类边界条件，也就是指定压力的值。假设，在出口上设置的压力值是0，那么，对于内场，压力小于出口压力的位置，p就是负数，压力大于出口值得位置，p 就是正数。当你需要计算实际作用到某个墙上的压力（也就是正应力）时，你需要先修正interFoam算出来的压力。那么怎么修正呢？还以上面这个例子，计算中你出口边界设置的压力是0，但是实际情况下，出口面的压力则是101325pa，那么所谓修正其实就是把interFoam算出来的p加上101325就可以了。把这个修正后的值再算对墙的正应力，才是合理的。

@calrel 需要显卡支持，如果是 N卡，linux 下面需要安装闭源驱动才行，至少以前是这样的，不知道最新的开源驱动，像Mesa是否支持 surfaceLIC 了。

@Joann

1. volScalarFieldValue zoneID 的意思是有一个场的名字叫 zoneID，默认值是123
2. 两个box 的并集区域设置为水域，之所以弄成两个box，大概是因为水域不是一个规则的 box，所以用两个box。如果是更复杂的形状，那还可能是 3个 box，或者用其他的复杂形状。
3. 不是，而是对属于 c0 这个set 的网格，把场 zoneID 的值设置为0，同理，c1里面的网格，zoneID 设置为 1。

setFieldDict 的逻辑还是比较清楚的，无非是指定一个区域，把这个区域内的某个场设置为某个值。区域的指定方式可以是规则的 box，sphere之类的，也可以是某个 cellSet，至于这个cellSet 则可以使用 topoSet 来生成。topoSet 生成 cellSet 有很多很多选项可以用。

@hery
可以参考这个例子：pisoFoam/LES/pitzDailyMapped

internalField 部分应该是不支持 codeStream 吧