个人觉得没必要考虑吧

@Rachel0096 就加在你的所用的求解器的c文件更新粒子附近的地方吧，我只是写了个大概的，具体的你得看你用的求解器里面的具体内容

@Rachel0096 ```
label particleNumber;
particleNumber = 0;
forAllIter( basicKinematicCloud , kinematicCloud, iter)
{
parcelType& p = iter();
// 判断颗粒位置是否在你的slice附近范围内，如果是:
particleNumber+= p.nParticle();
//否则不计数
}
reduce(particleNumber, sumOp<label>());

@李东岳 我在使用过程中，一个网格里多个ib点基本上没问题，问题都出在如果只是固体边界布置一圈浸入边界点，内部是虚拟流体，这样容易透水，但是这方法我算出来的力很光滑。对于很多复杂界面，我为了防止透水，固体表面包络的所有内部流体网格里我也全都每个网格中心布置一个浸入边界点，这样可以避免大部分透水问题，但是这样有时候算出来的力会有小震荡，应该是还是对内部的插值做些特殊判断我没做

@eason 理论上来说是的，用外在的导出得到的散点不一定满足ib点条件，在这种情况下，我会把散点包络的所有内部虚拟流体网格也全部作为浸入边界点了，这样能有效避免透水问题。我原本的方法是只是边界布置点内部虚拟流体，但是确实稍有不慎就透水

@eason 就建模软件建好，导出stl或其它格式，把其中散点数据单独复制出来，然后再openfoam读入进行后续操作

@Rachel0096 r如果只是边界布置进入边界点的IBM，内部是有流场的，他这个应该是有的我好像看过这篇，他只是后处理额外画个圆柱盖住内部虚拟流体

@jasper-0 我觉得不需要

@jasper-0 我的话我是没考虑你说的这个 主要还是先看你问题吧 你的颗粒属于哪类型的 一般需要什么力 你甚至可以不用他自带的力 全关掉 自己循环颗粒自己算你所有应该有的力 加在一起给U方程

@李东岳 他这里应该是自己用其他不带LPT的求解器额外连接openfoam的LPT包，在U方程里直接加比较简单