CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    亚临界雷诺数下二维圆柱涡激振动出现负体积

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

      各位老师,我在做亚临界雷诺数下二维圆柱的涡激振动(y向),目前模拟的工况是Re=5000,圆柱直径0.0381m,速度0.132m/s,约化速度为8,质量比2.6,阻尼因子0.0036,振动方程通过编写四阶龙格库塔的UDF导入到fluent中,振动方程进行了无量纲化。
      微信截图_20200712212314.png
      模型如下,上游10D,下游20D,宽度上下各15D,圆柱附近采用四边形网格,其余为三角形网格。第一层网格高度0.0023D,时间步长0.0014s,湍流模型k-wSST。
      图片1.jpg
      在低雷诺数时基于编好的UDF进行过计算(此时的模型不一样,直径0.0016m),与实验数据接近,但在亚临界时,计算开始一段时间后圆柱位移一直增大,直到出现了负体积,此时的位移已接近10D。自己更改过时间步长,但是问题也是一样,想问下有老师做过类似的模拟吗?这个问题排查了很久还是没能解决,谢谢大家!
      编写的UDF如下。
      微信截图_20200712215323.png

      1 条回复 最后回复 回复 引用
      • Y
        yxg2018 最后由 编辑

        你的动网格更新怎么设置的

        Y 1 条回复 最后回复 回复 引用
        • Y
          ycx1042341217 @yxg2018 最后由 编辑

          @yxg2018 我现在动网格是采用弹簧更新和网格重构,将圆柱边界设置为刚体运动

          Y 1 条回复 最后回复 回复 引用
          • Y
            yxg2018 @ycx1042341217 最后由 编辑

            @ycx1042341217 网格重构应该是只更新三角形网格,你可以将圆柱和边界层网格一起做刚体运动

            Y 1 条回复 最后回复 回复 引用
            • Y
              ycx1042341217 @yxg2018 最后由 编辑

              @yxg2018 谢谢,我将网格按照单元类型分离并将边界层网格也设定为刚体运动区域,计算依旧会报错,提示Divergence detected in AMG solver: x-momentum,我将迭代次数增加并改小了动量的松弛因子,计算还是发散。
              微信截图_20200713155718.png

              L 1 条回复 最后回复 回复 引用
              • L
                LGOHVIV @ycx1042341217 最后由 编辑

                此回复已被删除!
                1 条回复 最后回复 回复 引用
                • First post
                  Last post