@李东岳 谢谢

是的 方程13两个都是未知数

那么方程13是不是把所有的网格的体心上的值联立在一起求得的？

对

@capillaryFix 这个流弊，老铁

cefb3e3d-d0f9-4931-b852-e1c290800ee0-image.png 我的openfoam版本是2.4.x，是否可以使用postProcess这个命令？需要自行安装吗，还是只有新版本才能使用。🥺

cfd-online也有人提出了相同的问题，但Jasak表示没有问题，或许真的是安装版本的问题，请问用fe40的兄弟们有遇到相同问题吗？

如果可以正常运行，可以分享一下fe40安装包吗？

@李东岳 好的！十分感谢李老师的帮助和指导！数值异常出现在10s之后fd5fa411-1909-482d-add3-fbfd36d837ca-sloshing.tar.xz

随着某个文章出来的代码，你可以直接联系通讯作者，他们会告诉你版本

刚注册的账号。。为什么中间这段写出来是这个样子。。不好意思给各位带来了观看的麻烦。

@李东岳 谢谢老师

@李东岳 文章中没有具体的流场结果，都是流场和粒子耦合在一起的结果，但是改变流体速度值域是可以对的上的

@瑞水南 您好，请问层流里如何制定的这个问题解决了吗，能不能分享一下

这有点像我用starccm+和abaqus一起计算时候报错的情况，就莫名其妙的报错了...

image.png

@yms00yms 不好意思哈，这我真不好说，那个cyclic边界条件我没用过。

构造UEqu的时候报错，编译可以通过，openfoam7版本，
代码：

fvVectorMatrix UEqn ( fvm::ddt(rhoEff,U) + fvm::div(Rhophi, U) + fvc::div(rhoPhiR_PTRS) // pseudo-turbulent Reynolds-stress - fvm::laplacian(voidfraction*fluidViscosity, U) - fvc::div(voidfraction*fluidViscosity*dev2(Foam::T(fvc::grad(U)))) == - fvm::Sp(Ksl,U) );

运行到这的时候，报错：

[GPU01:121865] *** An error occurred in MPI_Recv [GPU01:121865] *** reported by process [1419640833,1] [01:121865] *** on communicator MPI COMMUNICATOR 3 SPLIT FROM 0 [GPU01:121865] *** MPI_ERR_TRUNCATE: message truncated [GPU01:121865] *** MPI_ERRORS_ARE_FATAL (processes in this communicator will now abort, [GPU01:121865] *** and potentially your MPI job) [GPU01:121859] 2 more processes have sent help message help-mpi-errors.txt / mpi_errors_are_fatal [GPU01:121859] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

大家有没有遇到过这样的问题？

9c713226-c506-438a-a314-69db456707a2-image.png

@mark 请问楼主是如何植入相变的呢，能向您请教一下吗，或者有什么相关文献分享一下

@李东岳 麻了，哭。。。我多么希望是我自己的问题

绝对666

问题找到了，内聚能密度设置太大了，导致受力太大，进而导致粒子消失

@Joker 好的，感谢大佬建议，我去试一下这个