请教关于浸入边界法的问题，IBM方法的问题

[eason]

请问一个问题，圆柱的ibm点是好生成的，假如比较复杂的三维体应该如何生成ibm点呢？

[OItoCFD]

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

[eason]

@OItoCFD 理论上来说1grid1point比较好吧，建模软件自带的散点应该没法按照网格分辨率输出吧

[李东岳]

理论上来说1grid1point比较好吧

我刚入这一行，但我怎么感觉反了。为什么不是一个网格cell，多个粒子点更好？望大佬指点

[Rachel0096]

@李东岳

• 如果设置的浸入边界点个数过少，使得浸入点间的长度大于单位网格长度（d>h），则不能保证每个流体质点都能受到浸入边界点上体积力的作用，从而导致部分质点从圆柱体内部渗透过去，这显然与实际情况是不符的，计算结果也会相差很多。

• 如果设置的浸入边界点个数过多，虽然能够保证浸入点间的长度小于单位网格长度（d<h），但浸入边界点个数的增多也会带来计算量的增大，从而影响到整体的计算效率。
  所以一个grid一个point比较好 （转载论文里面的话hhh）

[Rachel0096]

@OItoCFD 想问一下老师，应该用什么建模软件比较好呀
对于二维的圆柱，我是通过matlab生成均匀分布的点的，三维还没有试过

[eason]

@Rachel0096 附议，1cell2-3个points也没事，但是多了有时候会崩，可能force和velocity迭代的时候收敛不了？
我不是用of的 我现在想把复杂的几何构型（如stl）转成1cell 1point的sample文件，找不到很好的方法，请问大佬一般如何操作

[eason]

@Rachel0096 论坛里面有个大佬用gmsh，但是我不知道咋操作，三维简单的圆柱方柱还能用matlab码，复杂的构建就太蛋疼了

[Rachel0096]

@eason 我看还有做啥复杂机翼的，压根不知道咋生成的拉格朗日点

[OItoCFD]

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

[OItoCFD]

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

[eason]

@Rachel0096 是的 脑子疼

[eason]

@Rachel0096 这个可以用但是我用不了，这个点云不均匀，我用了会崩

code_text

[Rachel0096]

@eason 谢谢大佬！我去看看！

[李东岳]

这个帖子里面大佬很多，并且讨论很详细。我之前在这里发过一个帖子

https://www.cfd-china.com/topic/6265/浸没边界法-超音速自定义求解器出现-浮点数例外-核心已转储/9

如何从欧拉场映射到拉格朗日场，以及从拉格朗日场映射到欧拉场，感觉算法很多的样子。你们都是用的什么算法？

[eason]

@李东岳 delta函数