分享CFDEM+OpenFOAM+LIGGGHTS初次编译安装过程

[lyc]

向你们学习

[JokerZhu]

您好，安装编译所需文件时报错：Package libvtk6-dev is not available, but is referred to by another package.
This may mean that the package is missing, has been obsoleted, or
is only available from another source

E: Package 'libvtk6-dev' has no installation candidate
E: Unable to locate package python-numpy
怎么解决呢

[lixin19981013]

@JokerZhu
将“libvtk6-dev”修改成“libvtk9-dev”

[JokerZhu]

@lixin19981013 感谢

[小刘lyw]

@lixin19981013 您好，我按照您的回复将“libvtk6-dev”修改成“libvtk9-dev”，但仍出现以下报错：
E: 无法定位软件包 python-numpy
请问是什么原因呢

[小刘lyw]

@小刘lyw 已经成功解决了，将“python-numpy”修改为“python3-numpy”

[李东岳]

感谢分享

[capillaryFix]

之前用过公开版的CFDEM，目前在尝试将OpenFOAM esi 版本(v2312)跟LIGGGHTS耦合，有感兴趣的同学嘛，一起研究！

[Mina Lee]

@capillaryFix 我之前想干的事，是把openfoam和liggghts通过precice耦合起来，但是这个工作量太大了，一个人很难完成。
现在的OF和liggghts是通过侵入式的耦合，导致了升级OF版本，就需要重新修改求解器。

这是我之前的设想，这个设想，不需要修改liggghts源代码，但是，因为现在苦于毕业没文章，所以这个想法只能暂时搁置了。

[李东岳]

我不太懂liggghts，但我对耦合比较感兴趣。liggghts那面计算颗粒DEM，有没有输出力的接口？老铁能把github上的这一块代码发个链接么

OpenFOAM跟别的软件耦合都是一个路子。只要外挂软件有接口，在OpenFOAM里面调用就行。数据驱动LES是用libtorch计算湍流粘度，然后再OpenFOAM里面调用。那CFDEM就是liggghts算颗粒，然后在openfoam里面调用。要做到这个耦合，要懂liggghts的代码。

[Mina Lee]

@李东岳 DEM软件是没有网格的概念的，因为DEM软件就是对每个颗粒求解牛顿第二定律，计算得到每个颗粒的运动。

所以OF耦合liggghts，OF是提取liggghts颗粒的位置、速度......，然后OF识别颗粒所在网格，最后计算拖曳力、体积分数......。反之，根据OF传递进来对每个颗粒的力，liggghts计算牛顿第二定律，得到颗粒位置。

换而言之，CFD-DEM耦合，大多数的活都需要CFD软件干。

所以我觉得，要做这个耦合，OF的代码修改量可能会远远大于liggghts的代码修改量。

[capillaryFix]

@Mina-Lee 感谢回复，我的开发基于CFDEMcoupling-PUBLIC(https://github.com/CFDEMproject/CFDEMcoupling-PUBLIC)，没有使用您提到的preCICE，我对该程序进行了适当的修改，让其识别OpenFOAM v2312的环境，此后基于OpenFOAM v2312再编译CFDEMcoupling这个库，根据每次的报错，修改相应的位置；主要原因是从OpenFOAM5.x版本到v2312版本，部分函数的定义有所不同。目前已经实现了两者的耦合，也进行了验证，还有部分程序目前无法兼容，但是不影响自己的工作，就暂时先放在To Do里。

[capillaryFix]

@Mina-Lee Lee博士的回复很有道理，我此前的开发也是一直在OpenFOAM里，只在LIGGGHTS里植入了一些接触模型之类的。我的浅显理解是，当时DCS computing开的CFDEMcoupling，其实可以理解成是一个OpenFOAM的求解器，该求解器可以跟LIGGGHTS进行数据互换，计算void fraction, drag forces等等。

@李东岳 李导，代码在这里：https://github.com/CFDEMproject/CFDEMcoupling-PUBLIC