各位朋友：

我想通过计算得到U在切向和轴向的速度分量，所以要用到如下的公式：

Uradial = vector(position) & U / mag( vector(position))

读取U场数据用如下代码：

IOobject Uheader ( "U" runTime.timeName(), mesh, IOobject::MUST_READ );

可是在读取各个网格中心坐标的时候，遇到了问题，我的代码如下：

void Foam::calc(const argList& args, const Time& runTime, const fvMesh& mesh) { //bool writeResults = !args.optionFound("noWrite"); IOobject Uheader ( "U", runTime.timeName(), mesh, IOobject::MUST_READ ); // Check U exists if (Uheader.headerOk() ) { Info<< " Reading U" << endl; volVectorField U(Uheader, mesh); // 读取各个网格中心点的坐标 volVectorField& centres = mesh.C().internalField(); // 由于mesh.C()在网格内部返回的是中心点的坐标，在网格边界处返回的是面心坐标，所以用.internalField()，可是在这儿报错 volScalarField URadial ( IOobject ( "URadial", runTime.timeName(), mesh ), ((centres - centres.component(2) ) & U) //mag ((mesh.C() - (0., 0., mesh.C().component(2))) //这个地方我想得到一个平行于z轴的坐标，方法是用这一点的位置向量减去（0， 0， z）可是系统报错告诉我不能这样减。 ); } else { Info<< " Missing p or T" << endl; } Info<< "\nEnd\n" << endl; }

我想问一下，如果想实现mesh.C()- (0, 0, z)的话，需要怎么做？

非常感谢!
Janry

@qjh888 h好的，谢谢啦

@qjh888
Creating MRF model不代表华谊网格和MRF一起用。你要确保MRF被激活了。

转动网格后就不需要添加速度方程其他源项。MRF和华谊网格类似CFD下面的FVM和FEM。为同时求解问题的两个派系。目前还没有见过MRF叠加华谊网格一起求解。

@zhanghan 1129610784一起交流下、

非常感谢及时回复，对我很有帮助。
问题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!