CFD中文网

@李东岳 在 OpenFOAM并行计算不同进程数下获取场的积分值不同 中说：

scalar meshVolume(0);
forAll(mesh.V(),cellI)
{
meshVolume += mesh.V()[cellI];
}

Pout << "Mesh volume on this processor: " << meshVolume << endl;
reduce(meshVolume, sumOp<scalar>());
Info << "Total mesh volume on all processors: " << meshVolume << endl;

老师，您的代码我测试了，结果如下：
四核：

[0] Mesh volume on this processor: 829.44 [1] Mesh volume on this processor: 1320.32 [2] Mesh volume on this processor: 718.72 [3] Mesh volume on this processor: 1499.52 Total mesh volume on all processors: 4368

双核：

[0] Mesh volume on this processor: 2177.28 [1] Mesh volume on this processor: 2190.72 Total mesh volume on all processors: 4368

但是domainIntegrate这个函数的作用不就是为了整合各个计算域之间的结果吗：

domainIntegrate(phi) = gSum(fvc::volumeIntegrate(phi));

同时我测试了另一个网格，用24核和36核算了一段时间后发现两者相同时间下差距较小，感觉应该不是某一个核心的结果，而更像是误差？
d35b8301-5575-4f3f-a069-39d31c3892dd-image.png 784bdb80-c837-439c-bf76-0c56f82b271e-image.png

同时还想请教一个问题，就在刚刚我发现前后两次不改任何设置的情况下，进行decompose操作(scotch方法)得到计算域划分结果竟然是不同的：
464b6812-f1ff-440b-91ae-e54c881f8d61-image.png
进而导致一个计算正常进行，一个计算一段时间就发散了，这是什么原因？

我尝试用refineMesh进行多层加密，第一层加密完成后，在其内部的一个子区域进行第二次加密时，出现useHexTopology specified but cell 5018638 on face 7394 of patch bottom is not a hex的问题，我检查发现5018638单元并不在我的第二次加密区域内，为什么仍然要检查他是否为hex？
以下为refineMeshDict：

set expand7; coordinateSystem patchLocal; patchLocalCoeffs { patch bottom; tan1 (0 0 1); } directions ( tan1 tan2 ); useHexTopology yes; geometricCut no; writeMesh no;

如果改用geometricCut会不会导致网格质量下降？

想尝试与一篇文献中实验结果作对比，但根据文中所提供的条件算出的入口速度很小（0.18m/s），以此作为入口的边界条件；同时入口边界的alpha.water也固定为1，计算一段时间后如图所示，发现入口处几乎没有水相的存在了。这是为什么？

ae0a7d49-651d-4636-a27b-c8a0ff4b9334-image.png

入口处的边界条件如下：

U { type fixedValue; value uniform (0.18 0 0); } p_rgh { type zeroGradient; } alpha.water { type fixedValue; value uniform 1; } alpha.air { type zeroGradient; }

@李东岳 LES和laminar都试过，只要把空化核数降一些就不会发散。

@Joe 我也是小白还在折腾这个三相的求解器😂，OF似乎没有现成的三相空化求解器，所以想着只能在两相模型的基础上推广到三相然后自己改一个求解器出来了。

又或者说可以通过一些边界条件将明渠流动化简为管道流然后再用interPhaseChangeFoam计算？

@youhaoyu 好的我试一试

@李东岳 谢谢老师:140:

@李东岳 好的谢谢老师:ok3:

是不是一开始里面就有水然后没给速度

@李东岳 是可以不选择的，但怕设置错了，这玩意用的人确实不多啊，老师还有啥推荐的网格软件吗

@李东岳 好的谢谢老师，课题组只有我一个做这个没有师兄师姐传授经验😢，所以是得多看看文章吗？

@ShaneHEEE 试过了还是不行。。。