CFDEM：Liggghts和OpenFOAM耦合实现CFD-DEM（四）

        岩石、混凝土材料的颗粒之间是相互胶结的，与流化床模拟时用到的松散的颗粒不同，在DEM-CFD模拟岩石流固耦合时需要添加粘结力，本博文编译了PFM版CFDEMcouping，并对计算结果后处理可视化。

一、CFDEMcoupling-PFM安装

     此处还是参考CFDEM-ACADEMIC安装教程 - 哔哩哔哩 (bilibili.com)这是目前能找到的最详细的中文版安装教程。PFM版本的Liggghts可以实现颗粒间bond，与OpenFOAM 5.x实现对岩石等颗粒间存在胶结的工况的模拟。

二、算例测试

     使用如下路径的算例测试bond的施加：

/home/chunqi/CFDEM/LIGGGHTS/examples/LIGGGHTS/Tutorials_public/bondspackage/sphere

结果会输出*.bond文件，注意在PFM版Liggghts并不能直接输出力链的vtk文件，Public版的输出力链的代码不通用，可以用以下代码输出力链的.dump文件，之后用Paraview（需要插件，后面介绍）打开，生成力链的.dump文件的代码如下：

compute     fc             all pair/gran/local id pos force
dump         forcechain     all local 1000 post/fc*.dump  c_fc[1] c_fc[2] c_fc[3] c_fc[4] c_fc[5] c_fc[6] c_fc[7] c_fc[8] c_fc[9] c_fc[10] c_fc[11] c_fc[12] #x1 y1 z1 x2 y2 z2 id1 id2 periodic fx fy fz

直接使用Paraview是不能直接打开.bond以及.dump文件的，需要安装相应的插件，但是插件所需的Paraview的版本较低（可能因为PFM版本基于的是老版本的Public版开发的，最近几年的Public版可以直接生成vtk），最新的版本是5.8.1的，插件安装方法参考Paraview 批处理绘制力链图 - 知乎 (zhihu.com)

结果如下：

颗粒间bond展示

力链分布