CFD中文网

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

    关于重叠网格的疑问

    OpenFOAM
    3
    14
    3255
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 秋
      秋山结月 最后由 编辑

      本人先是跑了一下tutorials的overInterDyMFoam的算例，一切正常，然后将模型换成了自己的一个算例，然后问题来了。
      openfoam识别不了两套网格的交叉区域，两边好像都是在独自在跑，没有插值，具体问题如下：

      Creating cellMask field to block out hole cells

inverseDistance : detected 2 mesh regions
    zone:0 nCells:1331200  voxels:(35 35 35) bb:(-10.0000079951 -5.00000799512 -6.00000799512) (-4.87499200488 -1.62499200488 -0.874992004884)
    zone:1 nCells:238837  voxels:(35 35 35) bb:(-62.0000295973 -25.0000295973 -26.0000295973) (-35.9999704027 -14.9999704027 -15.9999704027)
Overset analysis : nCells : 1570037
    calculated   : 238837
    interpolated : 0 (interpolated from local:0  mixed local/remote:0  remote:0)
    hole         : 1331200

Creating interpolatedCells field

      从计算开始到结束都没有变化，interpolated的数量始终为0

      主要有几个问题:

      1. 该怎么改.....
      2. 设置重叠网格跟setFields里设置的zoneID有关吗，如果无关，控制重叠网格的字典应该是哪一个？
      1 条回复 最后回复 回复 引用
      • 秋
        秋山结月 最后由 编辑

        @wwzhao 请教一下老师，非常感谢

        1 条回复 最后回复 回复 引用
        • W
          wwzhao 教授 最后由 编辑

          看起来两个 zone 没有 overlap？

          秋 2 条回复 最后回复 回复 引用
          • 秋
            秋山结月 @wwzhao 最后由 编辑

            @wwzhao 对 是这样 两套网格似乎是在独立算没有交互
            我后来又换了个算例，结果又有overlap了，然后今天又做了个，又变成零了。。
            这个规律到底是闹哪样。。。

            1 条回复 最后回复 回复 引用
            • 秋
              秋山结月 @wwzhao 最后由 编辑

              @wwzhao 想请教您一个问题

              现在我在计算中遇到了remote单元不断上升的情况，然后接着就是某一步速度或者压力突然迭代次数剧增，然后计算发散。

              我想请问一下您，一般背景跟overlap的网格差距设定为多大合适?

              另外就是PIMPLE各个设定，看您在别的帖子提到需要提高外循环的次数，一般多少比较合适呢?

              W 1 条回复 最后回复 回复 引用
              • W
                wwzhao 教授 @秋山结月 最后由 编辑

                @秋山结月

                保证插值区域处不同网格块的单元尺寸相当就可以。

                提高PIMPLE外循环可以适当增加时间步长。你的发散可能不是这个问题，网格的可能性大些，不过你可以尝试增加PIMPLE外循环次数，看看能不能解决问题。

                秋 2 条回复 最后回复 回复 引用
                • 秋
                  秋山结月 @wwzhao 最后由 编辑

                  @wwzhao 在 关于重叠网格的疑问 中说：

                  处不同网格块的单元尺寸相当就可以。
                  提高PIMPLE外循环可以适当增加时间步长。你

                  我目前的办法是为了保证插值区相差不大，在背景网格里用snappyHexMesh选了个refineBox来局部加密
                  pimple的话我是设定一个高次数的外循环，再用一个relTol控制次数（差不多4次左右就会满足relTol）
                  但是还是会发散，可能真的是本身网格有些问题
                  总而言之先谢谢老师了 我再试试调整一下网格

                  Y 1 条回复 最后回复 回复 引用
                  • 秋
                    秋山结月 @wwzhao 最后由 编辑

                    @wwzhao 另外想请教您一个问题，一般重叠网格（贴体的那个网格），尺寸应该取多大比较合适呢。
                    我看tutorials里的构件都比较小，可能实体是 0.2×0.3×0.4的大小，重叠网格大概就是0.4×0.6×0.8

                    但是我要计算的物体几何尺寸可能有几m到十几米，这样也重叠网格的尺寸也要翻倍？
                    还是说，只要能够保证重叠区包裹住实体就可以了？

                    W 1 条回复 最后回复 回复 引用
                    • W
                      wwzhao 教授 @秋山结月 最后由 编辑

                      @秋山结月 重叠网格的block size没有特殊要求，包住物体就行。

                      秋 2 条回复 最后回复 回复 引用
                      • 秋
                        秋山结月 @wwzhao 最后由 编辑

                        @wwzhao 好的好的 谢谢老师~~

                        1 条回复 最后回复 回复 引用
                        • 秋
                          秋山结月 @wwzhao 最后由 编辑

                          @wwzhao 老师我想请教一个问题。我看到一些论文把重叠网格分为插值边界单元，贡献单元以及孤点单元，我想问着是否跟OpenFOAM的重叠单元的interpolated from local , mixed local/remote 以及 remoted相对应呢？

                          我计算的时候遇到了remoted单元不断上升的情况（初始0，计算100s后可能上升到几千），但是我的计算确实也没有发散，能够进行下去，如果remoted对应孤点单元，是不是代表我这个计算设定是不准确的呢？

                          W 1 条回复 最后回复 回复 引用
                          • W
                            wwzhao 教授 @秋山结月 最后由 编辑

                            @秋山结月 http://cfd-china.com/post/11197

                            秋 1 条回复 最后回复 回复 引用
                            • 秋
                              秋山结月 @wwzhao 最后由 编辑

                              @wwzhao 看起来显示的只是个并行计算的信息？
                              总而言之谢谢老师的指点了

                              1 条回复 最后回复 回复 引用
                              • Y
                                Yingdong @秋山结月 最后由 编辑

                                @秋山结月 背景网格使用sanppy加密部分区域好像真的会发散，而且p_rgh的求解矩阵也会变为对称，不加密是非对称的，请问您最终怎么解决的这个问题

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