非常感谢及时回复，对我很有帮助。
问题1：也就是说公式（4）（5）中的1/f是表示nabla算子的下标f，而不应该显示为1/f，是吧？
问题2：后来参考“浅谈计算流体力学中的通量”一文，可以看到图中标有d。
明白了，谢谢！

http://blog.sina.com.cn/s/blog_97cadb820102w2h0.html
在这里找到一些回答

可以参考用户指南4.4，里面还是有不错的概括的。

用rhoCentralFoam计算非定常问题，是每个时间步内计算各个方程，我怎么知道它的每个时间步的计算是收敛的？计算是正确的呢？像pimple那样可以通过每个时间步内迭代残差的下降可以判断，不知道这个求解器是怎么保证的

@金石为开 Foam::fvPatchField::evaluate(const Pstream::commsType) 函数里面有调用 updateCoeffs()

@zhanghan 客气:sunglasses:

@jzheng
relax是松弛的意思。查看此处
另外各大CFD教材都有关于松弛的说明。低松弛...

@cfd-china 对，速度入口造波

参考 $FOAM_ETC/caseDicts/postProcessing/flowRate/flowRatePatch。

OpenFOAM-4.0 可以直接在controlDict中使用 #includeFunc 语句 [1]，用法参考user-guide [2]。

[1] https://github.com/OpenFOAM/OpenFOAM-dev/commit/dd20d9086889ccbf993e8648d8e8f070f1b16b79
[2] http://cfd.direct/openfoam/user-guide/post-processing-cli/

@cfd-china ![其中F的方向指向X轴，在x轴方向的分量为0即Fx=0
Y轴方向的分量为Fy=5y
Z轴方向的分量为Fz=5z

@wwzhao 不是，我程序中有一句写的不对

const vectorField& Cf = patch.Cf()/delta1VD; vectorField meanVelocity(Cf.size());

我改成：

vectorField Cf = patch.Cf()/delta1VD; vectorField meanVelocity(Cf.size());

或者

tmp<vectorField> tCf = patch.Cf()/delta1VD; const vectorField& Cf = tCf(); vectorField meanVelocity(Cf.size());

就可以了

@cfd-china 谢谢啦，:happy: 刚查了下工具书，是在fvschemes文件里设置

当中断后继续从当前时间开始后：

*.oldTime（）就是创建场时的值（相对于本时间为上一时间步）；

*.oldTime（）.oldTime（）的值OpenFOAM会自动存在输出文件中，命名为 *_0

@liujiaming
把你的截图发出来看看？和网格关系非常大。我感觉你的网格变形过大。

@cfd-china :happy: THANK YOU!

@桥下卧龙 也没帮上忙，实在不好意思，我最近也在做类似的东西，不过没出现你这样的情况。

0_1469972545332_foamToTecplot360编译.JPG 下载OF-2.3.0对应的tecplot360版本进行编译即可，我在OF-2.4.0版本安装没问题，你可以试试，我的理解是，OF-2.3.0后tecplot没有变

@administrators 现在可以自动识别代码了？

多谢楼主了，可以计算了，赞一个！