@浪迹天大 谢谢，我研究一下，这个file 我就是读了一下文件名字，然后放到IOobject里，原来在of4就完全没问题

positionsFile_(this->coeffDict().lookup("positionsFile")), positions_ ( IOobject ( positionsFile_, owner.db().time().constant(), owner.mesh(), IOobject::MUST_READ, IOobject::NO_WRITE ) ), PDFFileName_(this->coeffDict().lookup("PDFFile")), PDFsamples_ ( IOobject ( PDFFileName_, owner.db().time().constant(), owner.db(), IOobject::MUST_READ, IOobject::NO_WRITE ) ),

东岳老师，我最近在做一个类似的模型，颗粒粘度与其体积分数有关的，用的of8的multiphaseEulerFoam，现在完全没有思路，不知道要将其放在哪里:136: 望您指点一二

查到了一个类似的问题，是fluent用tecplot处理的时候也发现壁面速度不为0：
出现这样问题的原因为：fluent的计算数据是存储在每个单元的中心位置，而tecplot的云图显示是按照节点来显示的。如果是直接导入fluent的.cas文件和.dat文件的话，那么在tecplot中，节点上的值会通过单元中心点的值插值得到，从而造成壁面速度不为0的情况。
不知道OF这面有无特殊处理方式

是我没有解释清楚么？这里的genotype[i]就是输入的字符，因此minus 1.0 2.0 对应三个字符串，从i=0开始，对应minus, i++变为1， 返回第一个递归函数evalGeno对应的i=1,也就是进入else, i++ 为2，numValue为1.0，继续减号，然后到第二个递归evalGeno，此时i为2，所以按程序流程就是1-2，但是程序测试结果是2-1的结果，为什么会这样呢？请各位老师指点一下

问题已解决~因为变量不能从字典中正确读取，修改读取语句后可以正常运行

@bestucan 好的哈！谢谢版主！我仔细看看您发的网址

@bestucan 好的哈，谢谢教授回复！不过您说“但是经过相邻面的速度也反了”，速度的方向是（绝对）唯一的吧，就是说该是多少就是多少。感觉还是从高斯定理出发，法向量永远取外侧为正

@OItoCFD 在 新写的程序单核计算无问题，但无法并行计算？ 中说：

mpi并行那些需要考虑吗？

一般不需要考虑

除非你的surface场进行一些特殊的插值

@Yongbo 你说的两种方式是对的，我验证均通过编译了。东岳老师应该是笔误了

@cresendo 好的好的，我去试试，谢谢！！

@李东岳 不好意思啊 东岳老师，我表述的不是很清楚，我说的exe想说的是编译成application.就是下图中的EXE。谢谢老师回复！47afd1597f7c6cb232f1079f897a298.png

最近在做一个多组分求解器，如题，OF中如何实现单组分到多组分的混合呢？看到一个帖子：
https://www.cfd-china.com/topic/2183/openfoam多组分返回热物性?_=1619261653942
我想我也遇到了和楼主相同的问题：比如说我想求解多组分的定压比热Cp，而OpenFOAM中给出了单组分定压比热Cp的函数，但是多组分的Cp是通过一定的混合规则来求解的，OpenFOAM中从单组分到多组分的这个过程是怎么实现的呢？或者说OpenFOAM中的混合规则如何进行的？望了解的铁子指点一下，谢谢！

@cccrrryyy 后来查了一下GAMA可以算非对称矩阵。忘了当时怎么设置的了，现在改成我记得的边界条件GAMG可以算几步，但是算着算着就报错了，应该是发散报错。我现在感觉可能当时是一开始计算就发散了。。

@cfd_上弦月 谢谢回复，有道理，但因为是非均匀网格，即越靠近壁面越密集，我其实想要获取上图中绿色区域内每个网格的高度，（假设长度都是一致的），用于高度方向的积分运算。

你可以加载两个case.foam，然后显示不同的颜色？

最后一张图是X方向的受力，也就是船的阻力。
谢谢回复:146: ，如您所言，setFields确实有问题，我设置的值正好等于流体域边界值，当时没有检查其实边界处的相分数没有完全赋值。将setFields中的box范围改成稍大于流体域就可以正常求解了。我本来还以为是数值振荡这么严重呢，都想离散格式精度调低看看了。现在看来，还是物理模型导致的振荡才更严重。:136:

@wangfei9088 非常感谢:xinxin3:

@浪迹天大 谢谢您，您的总结让我受益良多。

从 ESI 把这个功能偷过来了，现在在 OF-7下能用了。但是限制比较大，需要完全均一的网格。
https://github.com/ZhangYanTJU/functionObjects
这是我做的结果：energySpectrum.png

都是关于并行计算的包。openmpi和metis都可以手动下载安装，你试试？