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

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

@李东岳 谢谢老师

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

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

这有点像我用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 好的，感谢大佬建议，我去试一下这个

明白了，谢谢东岳老师，我接下来试试3D模型。老师那您觉得是用层流还是湍流比较好呢

{5531D08E-7EEB-4a66-8334-0712B6413B99}.png
对slice 和clip 看网格用crinkle 看 会显示真实网格 不用的话 切面在网格内布 就会有杂线

@李东岳 李老师，计算过程中没发散啊，第一幅图我是在 paraview 中用了 ISO volume 过滤器出来的液相图，也搞不懂是啥原因提取出来的东西并不连续

我这个8个图，最笨的方法纯自己手搞的。你用宏更好，我没时间研究。有时候研究了半天就过去了。直接开8个算例完事了。

您好，发散问题您解决了吗，我也遇到了这个问题，这个求解器一直发散

可能因为您用的激活函数是ReLU，正好可以很好地拟合甚至泛化y=x，但如果学习一个高次多项式，损失就不一定很小了。