补充：

根据网上说的，加入对应代码之后，编译了ThirdParty-$versionNum/Allwmake，已经成功

然后编译src/Pstream文件下面的 Allmake文件，会报错如下：

'/home/zywang2/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/dummy/libPstream.so' is up to date. + case "$WM_MPLIB" in + set +x Note: ignore spurious warnings about missing mpicxx.h headers wmake libso mpi make: mpicc: Command not found SOURCE=UOPwrite.C ; g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -DOMPI_SKIP_MPICXX -IlnInclude -I. -I/home/zywang2/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude -I/home/zywang2/OpenFOAM/OpenFOAM-2.3.0/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOptSYSTEMOPENMPI/UOPwrite.o UOPwrite.C:29:17: fatal error: mpi.h: No such file or directory #include "mpi.h" ^ compilation terminated. make: *** [Make/linux64GccDPOptSYSTEMOPENMPI/UOPwrite.o] Error 1

错误是 没有找到对应的mpi.h文件，，，在以前的论坛上也没有找到有效的解决办法，，，

这个问题应该怎么解决呢？

恳请哪位热心朋友解答一下，谢谢

谢谢教授的夸奖，我的是弄的比较简单:mihu:

@东岳 是的，李老师，感谢李老师！两个版本的定义不一样！OpenFOAM-4中GeometricField字典中是这样定义的：TIM截图20191122141846.png
而foam-extend-4.0中是这样定义的：TIM截图20191122141934.png

@东岳 感谢李老师解答，我再试试

可不可以给速度场给两个入口边界类型呢

@Dingcy 哈哈哈哈，我就瞎看看，我主要做CFDEM方面

谢谢教授

@bestucan 在 openfoam5 安装问题 中说：

sudo apt install paraview

好的谢谢~

ESI版本1906，有个bug，在compressibleInterFoam下，将turbulence models换成twoPhaseTransport时，除了laminar外都不能算，直接报错。bug提交在8个月之前，但是并没有人回应，不知其他人有没有研究。

blockmesh需要在case里运行，这是非常基础的操作。找找东岳流体翻译的openfoam手册看看。

B336D52677F66CD2291AFAA97B39D47C.png
如图所示，不太了解这个错误是由于什么原因引起的，求解器是overOlaDyMFlow，恳请路过的大佬帮忙解答一下！

@东岳 谢谢您，是边界条件的问题。

@Jhhu 你好，我最后没有成功，由于时间很紧张，还是使用了商用软件。

。。。。。。自问自答 ，汗

phiHbyABf不是边界流量。。。。

所以以下代码应该是符合动量方程，猜测可能是进入压力PISO循环之前修改下边界条件的系数

( phiHbyABf[patchi] - rho.boundaryField()[patchi] *MRF.relative(SfBf[patchi] & UBf[patchi], patchi) ) /(magSfBf[patchi]*rhorAUBf[patchi])

找到了一个方法
用动态对象数组

IB **myIB=new IB *[3]; for(int i=0;i<3;i++) { myIB[i]=new IB(i+1); }

注意最后delete。

@nbyjn 界面重构是为了消除界面的数值耗散，vof都需要重构，只不过代数重构不需要重构准确的界面位置

谢谢两位大佬@史浩 @东岳

各位大佬，of里面除了时均值有面均值吗？就是假如有个1m长的水管，网格dx=0.001m，能不能在某个时间步下输出1000个值代表着一共1000个剖面上的某个物理量的面平均值？如果要自己编程解决，那下面这25284个网格怎么和实际坐标联系起来？谢谢

dimensions [0 0 0 0 0 0 0];//整个计算域所有网格的空气相分数 internalField nonuniform List<scalar> 25284 ( 0.878566 0.854193 0.876398 0.860644 0.879329 0.859877 0.878446 0.860776

@东岳 谢谢，改完您建议的就没有那个警告了

把原来流场的小眼睛点开