extend版本应该是类似“A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries”这篇文献的实现方法 ，力源项应该是“An immersed boundary method with direct forcing for the simulation of particulate flows”文章的思想，虽然本质上差不多，但是实现边界条件的方式还是不一样的

@OItoCFD extend版本的IBM是采用插值重构的方式植入的，我们本来想植入直接在动量方程中施加体积力源项的形式的IBM方法，不过最后还是使用了extend版本带的IBM方法

@东岳 谢谢

感谢，最近正好在用它画网格！还想问一下李老师平时用什么软件画网格！

还有这种情况会造成一个owner（或者neighbour） 对应两个neighbour（或者owner）那他们共享的通量有特殊处理么？

谢谢东岳老师及时回复，那关于第二点呢？

最近要使用refineMesh，使用后用paraview打开发现粗细网格交接处有四面体网格过渡，如下图

在网上查找发现和paraview中properties中选项有关，若选择默认的interpolated Volume Fields则会出现上图情况；
若修改该选项为use vtk polyhedorn过渡的四面体消失，如下图

我想知道这是paraview显示的问题，还是网格本身就存在这个过渡四面体？
还有这种情况会造成一个owner（或者neighbour） 对应两个neighbour（或者owner）那他们共享的通量有特殊处理么？

非常感谢李老师，cellLimited Gauss linear 1 格式比Gauss linear好很多，运行程序之前运行renumberMesh -overwrite改善网格质量也有作用

最近发现在openfoam计算压力梯度时，用blockMesh生成的均匀网格，给定压力场为均匀初始场，当压力初始场为零时，全场压力梯度为零，但是当压力初始场为101325时，梯度在x,y（2D）方向上绝对值大小为1e-6(由于网格质量引起的数值误差？)，虽然这个值量级比较小，但是目前计算的系统网格长度量级为1e-5，会对最终的结果产生误差，我分别在openfoam2.2.x以及openfoam6.0上测试都存在该问题，压力梯度计算方法为Gauss linear，不知道有没有人碰到类似的情况，或者有没有推荐的绘制网格质量较高的软件，我试试在不同网格下计算看看
下图为压力初始值为101325的压力梯度分布
（另附测试算例链接https://pan.baidu.com/s/1uFbx_xqbO-5QAbcmTSr-fA）

@东岳 老师，DEM这块不是我编写的，工作主要是在CFD部分，据我了解DEM是计算单个颗粒的曳力，不用指定气泡直径大小吧；那我再看一下CFD部分曳力计算吧，谢谢您

@东岳 确实有聚并破碎