@zhangxc0223 是的

问题已解决。
在.C文件中少了下面这段代码：

Foam::tmp<Foam::scalarField> Foam::rhotabularThermo::rho(const label patchi) const { return rho_.boundaryField()[patchi]; }

把这个添加到controlDict里面，其中zoneName你可以通过topoSet来设定

volFieldValue1 { type volFieldValue; libs ("libfieldFunctionObjects.so"); log true; writeControl writeTime; writeFields true; regionType cellZone; // Create your cellZone name zoneName; // and set its name operation volAverage; //sum or others fields ( T // List of fields to operate on ); }

@东岳 刚刚解决了，是我的.bashrc里面缺了一句关于mpi命令

请教下各位有没有遇到以下的情况：
我基于Openfoam中floatingBody案例进行修改，网格采用的是fluent meshing划分后倒入OF中的，使用的是重叠网格，多面体网格和六面体网格结合，模拟的是垂荡自由衰减，将结构起始的中心提到水面上1米。在计算的后期包裹结构的网格出现了严重的变形，请问有谁知道是什么原因吗？下图分别是开始是结束时的网格。
0_1546229223114_微信图片_20181231120136.png 0_1546229227768_微信图片_20181231120143.png
经过一番的尝试，基本认为是跟网格的导入方式没有关系，问题可能出现在dynamicMeshDict文件的设定上。目前我的dynamicMeshDict的设置如下图，其中floatingObject 是指中间的结构物，它的边界设定为wall。
0_1546229568766_1546229471(1).jpg
另外，在我修改后的case，overInterDyMFoam运行过程中出现很多如下的warning：no base point for face ***，所以上述问题可不可能就是出在这个warning上。
0_1546229826622_微信图片_20181231121656.png
问题困扰好多天了，亟需帮助，非常感谢！

@benqing 谢谢，等我研究研究，不懂的再向您请教。

@yfclark

Jasak PhD thesis, P112:

In order to improve the solver convergence, it is desirable to increase the diagonal
dominance of the system. Discretisation of the linear part of the source term, Eqn.
(3.36), is closely related to this issue. If Sp < 0, its contribution increases diagonal
dominance and Sp is included into the diagonal. In the case of Sp > 0, diagonal
dominance would be decreased. It is more effective to include this term into the
source and update it when the new solution is available. This measure is, however,
not sufficient to guarantee the diagonal dominance of the matrix.

谢谢，我已经看懂了

ok谢谢老师！

0_1545971746839_捕获2.JPG

0_1545971753428_捕获.JPG

上面两个图分别是我从文献里面找的，和用E-QBMM模拟的气泡云图。因为不太了解气泡减租领域，这个领域比较关心什么数值？文献里面好像重点看Cf。有没有关注相分数的? 好像比较少呢

@wwzhao 好的，谢谢老师

@ibelief 谢谢！

写none的关键词后面应该改成directionalSolidThermo、fluidThermo、lookup、solidThermo中的一个吧？！

OLAFlow?那我假设你是做Ocean/Coastal Engineering的，那这个应该是关于波浪的问题？波浪-结构物相互作用的问题？
个人感觉OLAFlow对interFOAM的改动不算大，一般来说，你的这个可能问题可能是由于MULES算法，自由面出现的wiggles

看看 这个 满足你的需求吗？

这里有一个点火的模型可以使用：Github

把你的inlet-air和inlet-water做成两个patch，然后分别设置不同的alpha值，比如inlet-air的alpha为1，inlet-water的alpha为0

好吧，问题解决了，在其他相同case中找到了k的file，然后发现其他patch都可以用nonuniform List <scalar>,只有文件的第一个patch使用的是uniform 0.005，于是改变了一下两个patch的设定顺序就可以正常运行了。

原本不相信同事说的他在编译的过程修改一下公式的顺序就成功编译了的故事。。。现在我相信了！

这是修改后的文件。
0_1545414489934_d5f799e4-d82d-4a3d-b2fb-782978c9951f-image.png

应该是openfoam读取的要求
0_1545414877346_57f889eb-81d8-444c-8d9e-db3fda972d43-image.png
0_1545414807615_59b3f922-24d0-4d55-83d6-935c0235edc7-image.png

数值上，湍流的出现主要取决于精度。在FVM这种精度不高的情况下须添加随机扰动来制造湍流同时进行维持。我不确定你这种情况湍流会不会维持。你可以试试

各位大神们：
如何提取壁面的数据（尤其是曲面或曲线），比如热通量，压力，速度等值，并且配有轴向坐标值。