CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    东岳老师最近做的DNS验证算例疑惑点

    OpenFOAM
    3
    13
    2946
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • C
      Calf.Z-DNS 最后由 编辑

      参考了东岳老师复现的DNS算例(在东岳流体网站的icoFoam解析中)

      里面说到:

      本算例初始场通过boxTurb给定湍流涡结构之后,运行了200个周期,作为稳定的初始场。

      周期是如何计算的呢?这个是要在粗网格上面使用boxTurb并且跑一段时间,然后把结果映射到目标网格上作为初始场吗?还是直接在目标网格上跑?如果是后者,为何有初始场这一说法?

      我目前打算计算pipe flow 的DNS,用的是perturbU小工具产生扰动场,不知道如何检验它的效果呢?
      如果我画一个差不多大小的方形管道(网格数2的n次方),然后用boxTurb,之后把结果映射到我的pipe里面作为初始场,同时进出口使用cyclic,这样可以吗?

      另外我算的是加热可压缩管道流,用piso或者pimple算法来算DNS都是可以的吧?添加meanVelocityForce是必须的吗?

      DNS接触不久,望解答,谢谢!

      李东岳 1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 @Calf.Z-DNS 最后由 编辑

        @Calf-Z-DNS 在 东岳老师最近做的DNS验证算例疑惑点 中说:

        周期是如何计算的呢?

        速度0.1335米,按照这个计算的,就是跑了几圈。不是用boxTurb跑几圈..

        我没还测试perturbU,好像不是OpenFOAm内置的?

        加热可压缩管道流,用piso或者pimple算法来算DNS都是可以的吧?添加meanVelocityForce是必须的吗?

        可以用rhoPimpleFoam,压力驱动流需要meanVelocityForce

        线上CFD课程开始报名:http://www.dyfluid.com/class.html

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

        C sibo 2 条回复 最后回复 回复 引用
        • C
          Calf.Z-DNS @李东岳 最后由 编辑

          @东岳 谢谢!

          perturbU是我在cfdonline下载下来的:https://www.cfd-online.com/OpenFOAM_Discus/messages/1/2946.html
          不知道怎么验证效果。

          我的例子是空气竖直向上流过加热圆管,有重力,可压缩问题。
          进口采用mapped映射速度,速度出口采用advective boundary;用perturbU施加扰动。不知道还要不要添加meanVelocityForce?

          目前正在用buoyantPimpleFoam来跑,但是pimple循环迭代几次就报错了,我就把outcorrecter设为1,改为piso mode。
          但残差一直在规律地震荡,且进出口的质量流量不守恒。还望指点,谢谢!

          李东岳 1 条回复 最后回复 回复 引用
          • 李东岳
            李东岳 管理员 @Calf.Z-DNS 最后由 编辑

            @Calf-Z-DNS 个人感觉你这个是浮力驱动,不需要加。meanVelocityForce用于压力驱动

            线上CFD课程开始报名:http://www.dyfluid.com/class.html

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

            C 1 条回复 最后回复 回复 引用
            • C
              Calf.Z-DNS @李东岳 最后由 编辑

              @东岳 谢谢!
              如果我进出口都采用cyclic边界,也不需要通过fvOptions添加源项吧?因为是浮力驱动

              文件里ncorrect是设成3,不是一般设成2就可以了吗?

              另外我之前用gauss linear感觉收敛性有点差,后来一直用linearupwind。QUICK格式好像也比较难收敛,不知道有什么意见呢?

              谢谢!

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

                @Calf-Z-DNS 在 东岳老师最近做的DNS验证算例疑惑点 中说:

                如果我进出口都采用cyclic边界,也不需要通过fvOptions添加源项吧?因为是浮力驱动

                因为没看过相关文献,个人感觉是这样的,浮力可以保持稳定的速度

                纯六面体+O的网格?不应该收敛性差啊,贴网格瞅瞅

                线上CFD课程开始报名:http://www.dyfluid.com/class.html

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

                C 1 条回复 最后回复 回复 引用
                • C
                  Calf.Z-DNS @李东岳 最后由 编辑

                  @东岳
                  7f4161cd-5bb9-4558-b220-93193a1d9934-图片.png

                  之前用湍流模型模拟超临界CO2的时候会有收敛性的问题,用linearupwind会比较好一些,它也是二阶的。
                  我猜测是超临界状态的CO2物性剧烈变化引起的收敛性问题。

                  李东岳 1 条回复 最后回复 回复 引用
                  • 李东岳
                    李东岳 管理员 @Calf.Z-DNS 最后由 编辑

                    @Calf-Z-DNS :wolaile:

                    个人感觉你这个是浮力驱动,不需要加。meanVelocityForce用于压力驱动

                    验证没,浮力驱动用加meanVelocityForc么

                    线上CFD课程开始报名:http://www.dyfluid.com/class.html

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

                    1 条回复 最后回复 回复 引用
                    • C
                      Calf.Z-DNS 最后由 编辑

                      算例还在跑,一个月了:zoule:
                      现在我没有加meanVelocityForce,就只有重力场,浮力驱动。

                      结果发现壁面温度朝着想要的方向发展,但计算的流场逐渐趋于均匀,感觉是变成层流。。我开始计算前没有加扰动,直接开始计算,不知道这会不会有影响?雷诺数是6000,在理想状态下,流动会变成层流吗?不知道最终结果应该是啥状态,望指点,谢谢!

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

                        @Calf-Z-DNS 在 东岳老师最近做的DNS验证算例疑惑点 中说:

                        我开始计算前没有加扰动,直接开始计算,不知道这会不会有影响?

                        多谢反馈

                        没有浮力驱动的时候,会有很大影响。二阶精度下的壁面不足以创造涡旋发展成湍流,因此需要添加扰动。你有浮力驱动和密度差,不确定会不会从均一场发展成湍流,我个人感觉最好添加个扰动,对比下

                        线上CFD课程开始报名:http://www.dyfluid.com/class.html

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

                        C 1 条回复 最后回复 回复 引用
                        • C
                          Calf.Z-DNS @李东岳 最后由 编辑

                          @东岳 好的谢谢!
                          没加扰动的时候,我看圆管后半段是湍流状态,但随着时间的进行,紊乱的部分越来越少,感觉像是维持不住湍流状态。但我发现壁温变均匀的区域,流场也变均匀,这样子壁温呈现均匀上升时,流场也全是均匀的了。

                          目前我想到添加扰动的方式是用perturbU小工具扰动初始场,然后开始跑;perturbU是cfdonline上的工具,不是官方的。
                          另外一种是画一个等长,出口入口边长等于圆管直径的方腔,用boxturb加扰动,然后把速度场映射到圆管上。
                          不知道哪种方法更好呢?谢谢!

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

                            目前没对比过,也不清楚哪一种好呢 :tishizi:

                            线上CFD课程开始报名:http://www.dyfluid.com/class.html

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

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

                              @东岳 请问东岳老师,可以用速度基求解器比如rhoCentralFoam求解这个算例么?还得加个温度去解压力?
                              没算过DNS,望指点!谢谢!

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