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

[zlwei]

@五好青年
不是很熟悉OpenFoam的IB方法。

[Samuel-Tu]

我是自己在of里写的直接力--浸入边界法。。浸入边界内部的速度都是非物理速度。。一般保证浸入边界上速度为0就行。。OF-EXTEND里的浸入边界法原理好像有点不一样。。

[五好青年]

@samuel-tu

您好，

您是在哪个版本的of上，写的直接力-ibm呢？
另外foam-extend，跟你写的这个直接力-ibm，区别在哪呢，能大概说下么

谢谢

[Samuel-Tu]

@五好青年 of1806,of1912,of2006都写过。。enxtend好像是隐式的浸入边界法，直接力-浸入边界法是属于显式的。。

[五好青年]

@samuel-tu

大佬，您写的浸入边界法，有相关论文么？
您的代码开源么
可以加下您qq？交流下

[Samuel-Tu]

@五好青年 我参考的文献：Ji C, Munjiza A, Williams J J R. A novel iterative direct-forcing immersed boundary method and its finite volume applications[J]. Journal of Computational Physics, 2012, 231(4): 1797-1821.
代码暂不开源。。
QQ就不用加了。。现在不做这块了。。

[OItoCFD]

@samuel-tu 咦 这是我老师的文章哈哈哈 真巧

[Samuel-Tu]

@oitocfd 最近有参加第一届涡激振动前沿国际学术研讨会吗

[OItoCFD]

@samuel-tu 没有没有 及老师的学生应该参加了吧我觉得。 我不做涡激振动 我的老师是JJR Williams 不是及老师哈哈

[Samuel-Tu]

@oitocfd

[Rachel0096]

@zlwei 老师，有一个问题我想请教一下，我计算的圆柱位置也是有蓝色的云图，但是我看文献里面感觉他是实心的，像是人为扣掉了一样（可能我的表述有点不太准确），有点疑惑

[OItoCFD]

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

[Rachel0096]

@OItoCFD 好的好的，谢谢老师！

[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文件，找不到很好的方法，请问大佬一般如何操作