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

[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函数