看起来是除以0了。 可能komega里面你设的初始条件有点问题。

这个植入在了湍流模型里面。这个turbulence->divDevRhoReff(rho, U)

@星星星星晴 非常感谢友友的建议！

你这个就是标准的类似6dof的动网格原理。你这个力要推动网格点的移动。你要么看一下SDA、或者6dof那面的处理方法。openfoam这面网格点的运动通过septernion进行。不是传统的力来移动3D坐标。要么你就自己写一个网格边界点受力被移动的在笛卡尔坐标系下的xyz函数。

@吴建民 我前两天看LIGGGHTS已经被替代了吧。我个人经验是 OF5 - OF8 在拉格朗日部分计算速度有很大的提升。

千万级别的颗粒是什么意思，是你要追踪千万级别？还是说一共千万个颗粒，在domain中只有10w个？ 或者你考虑一下 一个parcel代表更多的particle。

如大佬所说，你如果考虑碰撞，什么碰撞模型？deterministic的话根本不可能。

@李东岳 确实没看，主要在网上搜来着🤦‍，感谢老师！我现在研究一下！

@李东岳 嗯嗯，好的，明白了，谢谢老师

动网格一般发散的原因就是网格变形太大了导致网格尺寸太小了

@李东岳 感谢老师的回复。foamListtimes -help 里好像没有cp的选项，但在Allrun里加上两行

{latestTime="$(foamListTimes -latestTime)" cp -r "${latestTime}" newTime}

可以正常运行了。

您好，对于粒子轨迹现在您有什么解决办法吗

@AppleKiller 您好，请问您做气泡上浮使用什么求解器

@杨英狄 老铁这个研究怎么样了。发表没有。

请问还有和这个类似的求解器吗，是可以用于可压缩的，这好像是不可压的

时间步长缩小。这个是不好整。主要网格太大。调试起来都费劲。

可以参考这个写一个 https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/meshTools/sets/cellSources/cylinderToCell/cylinderToCell.C

:146: :146: :146:

好的，李老师，谢谢

github上搜一下cfdem就能找到

是的。

@chapaofa 大佬辛苦啦，我跑就是这样，不知道为啥，感觉命令没问题，我换台电脑试试吧。谢谢啦!