CFD中文网

    CFD中文网
    • 登录
    • 搜索
    • 最新

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

    OpenFOAM
    7
    35
    2004
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 五好青年
      五好青年 最后由 编辑

      请教各位，关于浸入边界法的问题。
      我最近在用浸入边界法，做圆柱绕流。发现2个问题请教下

      1. 圆柱附近的网格，是否要比不用浸入边界法（实际建模）模拟时，网格要多很多呢？
      2. 用浸入边界做圆柱绕流，结果后处理速度云图，为啥圆柱位置也有蓝色云图呢？

      QQ截图20210422210133.png

      I am a CFD machine with no emotions. Welcome to browse my Bilibili, search: seeeeeeeeeeer

      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

        @五好青年 在 请教关于浸入边界法的问题，IBM方法的问题 中说：

        用浸入边界做圆柱绕流，结果后处理速度云图，为啥圆柱位置也有蓝色云图呢？

        我不做浸没边界法。不过这个我知道，只是显示问题，实际没速度哈。

        第一个问题我想应该也是的，否则处理出来的边界倾向于是直线

        CFD高性能服务器 http://dyfluid.com/servers.html

        1 条回复 最后回复 回复 引用
        • Z
          zlwei 最后由 zlwei 编辑

          1. 圆柱附近的网格，是否要比不用浸入边界法（实际建模）模拟时，网格要多很多呢？
            不一定。但是IB方法当中，边界稳定性可能会比非IB的差一点，更密的网格可以帮助解决这个问题。

          2. 用浸入边界做圆柱绕流，结果后处理速度云图，为啥圆柱位置也有蓝色云图呢？
            这很正常，IB方法就是用圆柱内的速度去补偿圆柱外的，从而在圆柱表面得到需要的速度（比如no-slip）。所以IB方法里，圆柱内是会有速度的，都是artificial velocity。

          五好青年 R 2 条回复 最后回复 回复 引用
          • 五好青年
            五好青年 @zlwei 最后由 编辑

            @zlwei
            您好

            1. 最近在用foam-extend中的浸入边界法，算圆柱绕流（刚性问题）
            2. 如果我想考虑圆柱的流固耦合fsi，请问设置上的区别在哪呢？
              即：浸入边界法，做流固耦合问题，和刚性问题，设置上区别在哪呢？
              万分感谢！

            I am a CFD machine with no emotions. Welcome to browse my Bilibili, search: seeeeeeeeeeer

            O Z 2 条回复 最后回复 回复 引用
            • O
              OItoCFD @五好青年 最后由 编辑

              @五好青年 一般圆柱内就是有速度 速度反向 刚好内外插值到边界速度为0 满足无滑移

              五好青年 1 条回复 最后回复 回复 引用
              • 五好青年
                五好青年 @OItoCFD 最后由 编辑

                @oitocfd
                您好，老师。我是借这个楼，问你浸入边界下面这个问题呀

                最近在用foam-extend中的浸入边界法，算圆柱绕流（刚性问题）
                如果我想考虑圆柱的流固耦合fsi，请问设置上的区别在哪呢？
                即：浸入边界法，做流固耦合问题，和刚性问题，设置上区别在哪呢？

                万分感谢！

                I am a CFD machine with no emotions. Welcome to browse my Bilibili, search: seeeeeeeeeeer

                1 条回复 最后回复 回复 引用
                • Z
                  zlwei @五好青年 最后由 编辑

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

                  1 条回复 最后回复 回复 引用
                  • S
                    Samuel-Tu 最后由 编辑

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

                    五好青年 1 条回复 最后回复 回复 引用
                    • 五好青年
                      五好青年 @Samuel-Tu 最后由 编辑

                      @samuel-tu

                      您好，

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

                      谢谢

                      I am a CFD machine with no emotions. Welcome to browse my Bilibili, search: seeeeeeeeeeer

                      S 1 条回复 最后回复 回复 引用
                      • S
                        Samuel-Tu @五好青年 最后由 编辑

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

                        五好青年 1 条回复 最后回复 回复 引用
                        • 五好青年
                          五好青年 @Samuel-Tu 最后由 编辑

                          @samuel-tu

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

                          I am a CFD machine with no emotions. Welcome to browse my Bilibili, search: seeeeeeeeeeer

                          S 1 条回复 最后回复 回复 引用
                          • S
                            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就不用加了。。现在不做这块了。。

                            O 1 条回复 最后回复 回复 引用
                            • O
                              OItoCFD @Samuel-Tu 最后由 编辑

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

                              S 1 条回复 最后回复 回复 引用
                              • S
                                Samuel-Tu @OItoCFD 最后由 Samuel-Tu 编辑

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

                                O 1 条回复 最后回复 回复 引用
                                • O
                                  OItoCFD @Samuel-Tu 最后由 编辑

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

                                  S 1 条回复 最后回复 回复 引用
                                  • S
                                    Samuel-Tu @OItoCFD 最后由 编辑

                                    @oitocfd :143:

                                    1 条回复 最后回复 回复 引用
                                    • R
                                      Rachel0096 @zlwei 最后由 编辑

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

                                      O 1 条回复 最后回复 回复 引用
                                      • O
                                        OItoCFD @Rachel0096 最后由 编辑

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

                                        R 1 条回复 最后回复 回复 引用
                                        • R
                                          Rachel0096 @OItoCFD 最后由 编辑

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

                                          1 条回复 最后回复 回复 引用
                                          • E
                                            eason 最后由 编辑

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

                                            O 1 条回复 最后回复 回复 引用
                                            • O
                                              OItoCFD @eason 最后由 编辑

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

                                              E R 2 条回复 最后回复 回复 引用
                                              • E
                                                eason @OItoCFD 最后由 编辑

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

                                                O 1 条回复 最后回复 回复 引用
                                                • 李东岳
                                                  李东岳 管理员 最后由 编辑

                                                  理论上来说1grid1point比较好吧

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

                                                  CFD高性能服务器 http://dyfluid.com/servers.html

                                                  R O 2 条回复 最后回复 回复 引用
                                                  • R
                                                    Rachel0096 @李东岳 最后由 编辑

                                                    @李东岳

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

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

                                                    E 1 条回复 最后回复 回复 引用
                                                    • R
                                                      Rachel0096 @OItoCFD 最后由 编辑

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

                                                      E 1 条回复 最后回复 回复 引用
                                                      • E
                                                        eason @Rachel0096 最后由 编辑

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

                                                        1 条回复 最后回复 回复 引用
                                                        • E
                                                          eason @Rachel0096 最后由 编辑

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

                                                          R 1 条回复 最后回复 回复 引用
                                                          • R
                                                            Rachel0096 @eason 最后由 编辑

                                                            @eason 我看还有做啥复杂机翼的，压根不知道咋生成的拉格朗日点:zoule:
                                                            961591c9-dbff-4af6-81f7-520f17b26cf7-image.png

                                                            E 2 条回复 最后回复 回复 引用
                                                            • O
                                                              OItoCFD @eason 最后由 编辑

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

                                                              1 条回复 最后回复 回复 引用
                                                              • O
                                                                OItoCFD @李东岳 最后由 编辑

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

                                                                1 条回复 最后回复 回复 引用
                                                                • E
                                                                  eason @Rachel0096 最后由 编辑

                                                                  @Rachel0096 是的 脑子疼

                                                                  1 条回复 最后回复 回复 引用
                                                                  • E
                                                                    eason @Rachel0096 最后由 李东岳 编辑

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

                                                                    code_text

                                                                    R 1 条回复 最后回复 回复 引用
                                                                    • R
                                                                      Rachel0096 @eason 最后由 编辑

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

                                                                      1 条回复 最后回复 回复 引用
                                                                      • 李东岳
                                                                        李东岳 管理员 最后由 编辑

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

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

                                                                        替代文字

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

                                                                        CFD高性能服务器 http://dyfluid.com/servers.html

                                                                        E 1 条回复 最后回复 回复 引用
                                                                        • E
                                                                          eason @李东岳 最后由 编辑

                                                                          @李东岳 delta函数

                                                                          1 条回复 最后回复 回复 引用
                                                                          • First post
                                                                            Last post