CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    搅拌反应器ICEM结构网格划分方式

    Meshy
    2
    12
    579
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • C
      cyberk 最后由 编辑

      几何模型为四挡板 六直叶Rushton桨搅拌反应器,采用mrf模拟桨叶的转动
      挡板,桨叶,圆盘均处理为无厚度壁面
      我在网上找到的教程是将内外区分开划分,最后再将网格merge,最后导入到openfoam中进行计算
      但是现在遇到一些问题,想来请教一下各位前辈

      1. 应该将内外区分开划分还是一起划分?
        (桨叶部分的网格需要切分成多块,分开划分可以简化工作量)

      2. 根据我的理解,分开划分得到的是一个interface,一起划分得到的是interior,这两种应该都可以应用mrf对吗?

      3. 如果是分开划分的话,在interface界面上的网格节点是否应该对应?
        我划分了两套网格用fluent单相简单计算一下
        第一套:
        内区是25万网格,外区是25万,交界面处的周向划分网格数量是130和128,计算起来还好,但是也能隐约看到内外区的交界面,而且感觉外区网格数量相对偏多
        第二套
        内区是25万网格,外区是10万,交界面处的周向划分网格数量是130和64,外区我是按照第一套网格减半划分的,但是计算结果发现在交界线处出现错误,然后我调整网格合并顺序,错误居然消失了,不清楚是什么原因

      4. 如果一起划分的话,内区该如何加密?
        因为一起划分的话,内外区在周向上的划分就是一致的,可以做的加密也就是对桨叶处做加密,体现出来就是轴向以及径向做一定加密,但是这会使外区的网格也增加不必要的加密,增加网格数量

      5. 轴下方的网格过密,而且质量不高该怎么处理?
        当前搅拌桨下部的圆柱形区域我采用o切分,但是这些地方的网格偏小,可能会影响计算结果

      (我之前采用gambit生成的四面体网格以及fluent meshing生成的多面体-六面体网格,然后导入openfoam进行气液搅拌的计算,发现在桨叶处的k epsilon异常,经过查阅资料,所以决定采用ICEM的六面体网格)

      1 条回复 最后回复 回复 引用
      • C
        cyberk 最后由 编辑

        (我之前采用gambit生成的四面体网格以及fluent meshing生成的多面体-六面体网格,然后导入openfoam进行气液搅拌的计算,发现在桨叶处的k epsilon异常,经过查阅资料,所以决定采用ICEM的六面体网格,之前的网格,导入openfoam中,内区和外区的界面会被自动处理为interior,不需要在boundary以及0文件中给定边界条件,但是目前icem分开划分的方法,内外区交界面导入openfoam中会识别为interface,需要给定边界条件,查阅相关资料似乎是应该给定AMI相关的边界,我暂时还没有深入研究,如果最后一起划分,也就不需要处理这个边界了)

        我晚上回去补上相关的截图

        1 条回复 最后回复 回复 引用
        • C
          cyberk 最后由 编辑

          这个是搅拌槽的几何结构
          几何结构.PNG

          这个是内区的桨叶上的网格
          内区桨叶网格.PNG

          这个是内区的interface界面上的网格
          内区界面网格.PNG

          这是内区的一个剖分形式(深绿色的线圈起来的是桨叶所在的位置,深绿虚线是最开始的O切分,然后对桨叶所在的六个面做切分)
          内区剖分形式.PNG

          这是内区里面质量较差的网格的位置,其实就是轴下面的部分
          内区质量较差网格位置.PNG

          这是内区里面质量较差的网格
          内区质量较差网格.PNG

          这个是外区的interface界面上的网格(十万网格的那套外区)
          外区界面网格.PNG

          这个是外区的剖分形式,就是一个o切,然后在挡板和interface处再做两次切割
          外区剖分形式.PNG

          这是外区的整体网格
          外区网格.PNG

          这是外区整体网格的y轴视角
          外区整个网格y轴视角.PNG

          这是外区整体网格的z轴视角
          外区整体网格z轴视角.PNG

          这是外区质量较差网格的位置,也是轴下面的部分
          外区质量较差网格位置.PNG

          这是外区质量较差的网格
          外区质量较差网格.PNG

          这是李老师文章中的网格示意图,似乎是对桨叶所在位置的外区网格在径向和轴向也做了加密,这个网格是内区外区一起划分的吗?
          东岳老师论文中的图片.PNG

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

            应该将内外区分开划分还是一起划分?
            (桨叶部分的网格需要切分成多块,分开划分可以简化工作量)

            不需要分开画。分开画需要做AMI,还会出现交界面。一起画不会出现这个问题。

            根据我的理解,分开划分得到的是一个interface,一起划分得到的是interior,这两种应该都可以应用mrf对吗?

            是的。不过openfoam这面不存在interior这个概念。ICEM那面的interface可以理解为openfoam的AMI。

            如果是分开划分的话,在interface界面上的网格节点是否应该对应?

            不需要对应。从你问的问题上,感觉你走了分开画的路线。但可以一起画。

            如果一起划分的话,内区该如何加密?

            一起画的话,内部自然而然就回密,因此你肯定不需要加密。有时候反而轴的地方太密了,需要放粗。

            轴下方的网格过密,而且质量不高该怎么处理?

            减少网格数量。icem那面可以做分区加密。手头买有icem。就是在block加密那一页,我记得好像是选择block然后加密。图标我忘了。

            (我之前采用gambit生成的四面体网格以及fluent meshing生成的多面体-六面体网格,然后导入openfoam进行气液搅拌的计算,发现在桨叶处的k epsilon异常,经过查阅资料,所以决定采用ICEM的六面体网格)

            采用六面体网格肯定没问题。

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

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

              旋转计算域和静置计算域完全不需要分开画。你就正体画一个网格就行,然后再openfoam里面通过topoSet创建旋转域的区域。

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

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

                @李东岳 好的好的,谢谢李老师,我这就尝试一下

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

                  @李东岳 在 搅拌反应器ICEM结构网格划分方式 中说:

                  减少网格数量。icem那面可以做分区加密。手头买有icem。就是在block加密那一页,我记得好像是选择block然后加密。图标我忘了。

                  李老师您好,您说的分区加密是这个按钮吗?
                  分区加密.PNG

                  我用一个最简单的空心圆柱的几何模型尝试了一下,对靠外的一个block进行加密
                  这样实现的效果就是在靠近中心的block的周向网格数少一倍,可以降低中间的网格数量,这个是您说的效果吗?
                  分区加密结果.PNG

                  但是这样似乎是有悬挂节点的问题?我导入of中,发现这个存在悬挂节点的面在of中显示的不太正常,我用的icem版本是2022r2,openfoam7,是不兼容还是正常现象?
                  分区加密悬挂节点1.PNG
                  分区加密悬挂节点2.PNG
                  分区加密悬挂节点3.PNG

                  checkmesh的时候发现分区加密以后非正交增大了,最大有16,而没有加密的时候,非正交非常小,这个是正常现象吗?
                  分区加密checkmesh.PNG
                  无分区加密checkmesh.PNG

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

                    你目前这一步都是正常的。不太正常是paraview的显示问题

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

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

                      @李东岳 好的,谢谢李老师,为了使用python的一些功能,我的paraview之前重新编译过,我用其他电脑试一下

                      1 条回复 最后回复 回复 引用
                      • C
                        cyberk 最后由 编辑

                        更新一下
                        关闭paraview里面的DecomposePolyhedra就不显示了
                        decomposePolyhedra.PNG
                        但是paraFoam里面没有这个选项,得用paraview打开case.foam才有这个选项,这是为什么?

                        打开DecomposePolyhedra.PNG 关闭DecomposePolyhedra.PNG

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

                          paraFoam是通过另外一个按钮实现的。具体我也忘了。你用paraview能显示就好

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

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

                            @李东岳 好的,谢谢李老师耐心解答

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