CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    rhoSimpleFoam计算翼型发散

    OpenFOAM
    4
    13
    2481
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • D
      dzw05 教授 最后由 编辑

      理论上rhoSimpleFoam可以计算可压缩流动问题,官方给的tutorial中也有计算超音速流动的算例。但是我在计算一个三段翼的亚音速流动(U=125m/s)时,却无论如何都不能收敛,但是我把速度降到60之后就很容易收敛。不知道各位大佬有何高见?附件是计算的网格和设置文件,各位有兴趣可以试一试。rhoSimple.zip

      自主匠心,普惠仿真。

      D 2 条回复 最后回复 回复 引用
      • D
        dzw05 教授 @dzw05 最后由 编辑

        @dzw05 rhoSimple.z01.zip

        自主匠心,普惠仿真。

        1 条回复 最后回复 回复 引用
        • D
          dzw05 教授 @dzw05 最后由 编辑

          @dzw05 抱歉网站上传附件有大小和格式限制。所以我把压缩文件分成三份。rhoSimple.z01和rhoSimple.z01文件,下载后请去掉zip后缀。rhoSimple.z02.zip

          自主匠心,普惠仿真。

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

            问题已经解决?下载这俩个文件可以?

            CFD高性能服务器 http://dyfluid.com/servers.html
            2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

            D 1 条回复 最后回复 回复 引用
            • D
              dzw05 教授 @李东岳 最后由 编辑

              @东岳 并没有解决,还发现了新问题。simpleC和simple算的升力系数大概差20%

              自主匠心,普惠仿真。

              1 条回复 最后回复 回复 引用
              • J
                Jacobian 最后由 编辑

                非气动专业,尝试和大佬讨论,想问一下:
                (1)div(phid,p) 的格式为什么不用 bounded Gauss upwind 呢?经验来看,bounded格式稳定很多;
                (2)grad用最小二乘,以及interpolation用迎风也是为了数值稳定性吧?但貌似当前算例都用中心差分也可以,所以为了精度给改了。
                按照以上两项修改试着跑了几万步,不过不是特别收敛,连续性方程不满足,升阻力系数还在震荡。是不是这个物理问题本身是非定常的,或者是网格太粗糙呢?
                还有就是挺好奇这个网格使用sHM画的吗?像这类混合网格如何处理trailing edge这种特变尖锐的地方?

                D 1 条回复 最后回复 回复 引用
                • D
                  dzw05 教授 @Jacobian 最后由 编辑

                  @Jacobian 我也是初用OpenFOAM算气动问题。
                  (1)对于bounded格式,实际上看代码的话,会减去phid对应的场量的散度再乘以p。如果散度为0,那么没有问题。但是phid的散度实际上不为0,所以我倾向于不用bounded。
                  (2)最小二乘其实是为了保证精度。前面的cellLimited是为了数值稳定性。
                  我用Simple-C算没问题,但Simple一直发散,目前还没找到原因。
                  (3)网格其实实用Pointwise做的,SHM不太容易做。

                  自主匠心,普惠仿真。

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

                    下载rhoSimple.z02.zip rhoSimple.z01.zip 打不开呢

                    CFD高性能服务器 http://dyfluid.com/servers.html
                    2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

                    1 条回复 最后回复 回复 引用
                    • J
                      Jacobian 最后由 编辑

                      @dzw05 学习了,div(phid,p) 确实不用 bounded 为好
                      @东岳 把 z01~02 后面的zip去掉,然后对着 rhoSimple.zip 选择解压

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

                        zip去掉之后不久成了rhoSimple.z02和rhoSimple.z01了么?怎么解压

                        CFD高性能服务器 http://dyfluid.com/servers.html
                        2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

                        1 条回复 最后回复 回复 引用
                        • S
                          shampoo 最后由 编辑

                          气动方向菜鸡选手,也想和大佬讨论。

                          1. 大佬要不要试一下OpenFOAM里的SA模型?这个一方程模型我看他们在算气动时还挺常用的,可能会好收敛一些,不过kwsst精度好像会好一些。
                          2. 网格可能也有影响。边界条件或许可以参考OpenFOAM里算翼型的case。
                          3. 还有可能是这个翼型不一定适合这么高的流速,看到的一些文献里马赫数在0.2左右,所以收敛可能就是挺难的。
                          4. 对了还想问楼主SIMPLEC和SIMPLE算法有什么区别呀,看您的讨论似乎会好收敛一些?
                            airfoil.png
                          1 条回复 最后回复 回复 引用
                          • 李东岳
                            李东岳 管理员 最后由 编辑

                            可算解压缩出来了。。。先看网格,大佬这网格画的很有意思 :140: :140: :140:

                            替代文字

                            CFD高性能服务器 http://dyfluid.com/servers.html
                            2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

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

                              https://www.jianguoyun.com/p/DecMC2kQ9s3ZBhiw1P0EIAA 我粗糙的改了一下,目前可以收敛,不过总感觉网格好像有点稀,网格渐变很粗暴,:zoule:

                              上U下p
                              捕获.JPG

                              CFD高性能服务器 http://dyfluid.com/servers.html
                              2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

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