CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    动量方程的2/3μ这一项

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

      我看到东岳老师在buoyantPimpleFOAM里面写的动量方程为:

      动量方程.png

      按可压缩方程来讲,应该还有一项−2/3µ(∇ • U)。

      为何方程中没有体现出来,是归并到某一项中了吗?

      另外我正在参考的一篇论文里面的动量方程的形式也是没有那一项:

      动量方程2.png

      从方程上看,这应该是可压缩的吧?如果省略掉−2/3µ(∇ • U),那么方程应该是不可压缩的,所以我认为这一项应该是归并到某一项中,只是方程形式上没有体现出来。

      不知道大家怎么看?谢谢大家。

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

        多谢指正。你说的是对的,目前我把文档简单改了一下,你看看,若推不明白我再更新 :duang: 这一项在+ turbulence->divDevRhoReff(U)这里处理的

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

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

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

          在代码中是有体现这一项的。

          只不过您给的方程中没有体现这一项,所以是想说归到某一项里面还是其他的操作?谢谢。

          1 条回复 最后回复 回复 引用
          • 队长别开枪
            队长别开枪 教授 最后由 李东岳 编辑

            不可压流动连续方程是$\nabla \cdot \bfU = 0$,因此好多公式里把这项省略了,代码里加入这一项可能是为了统一,可压不可压都可以调用,而且加入这一项有助于加速收敛,因为只有收敛后$\nabla \cdot \bfU = 0$才成立。

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

              关于RANS方程中有效雷诺应力的计算,可以参考这篇文章 。

              C 1 条回复 最后回复 回复 引用
              • C
                Calf.Z-DNS @队长别开枪 最后由 编辑

                @队长别开枪 在 动量方程的2/3μ这一项 中说:

                不可压流动连续方程是$\nabla \cdot \bfU = 0$,因此好多公式里把这项省略了,代码里加入这一项可能是为了统一,可压不可压都可以调用,而且加入这一项有助于加速收敛,因为只有收敛后$\nabla \cdot \bfU = 0$才成立。

                可压与不可压的代码对这一项的代码不一样吧,一个是divDevRhoReff(U),一个是divDevReff(U),它们具体包含的内容不一样。

                对于我给的这个方程,应该是可压的,但是这一项没有在方程中体现,所以感到奇怪。

                队长别开枪 1 条回复 最后回复 回复 引用
                • C
                  Calf.Z-DNS @wwzhao 最后由 编辑

                  @wwzhao 在 动量方程的2/3μ这一项 中说:

                  关于RANS方程中有效雷诺应力的计算,可以参考这篇文章 。

                  谢谢,我参考一下。

                  1 条回复 最后回复 回复 引用
                  • 队长别开枪
                    队长别开枪 教授 @Calf.Z-DNS 最后由 编辑

                    @Calf-Z-DNS 字面不一样,参考wwzhao贴出来的文章,对于线性涡粘模型,divDevReff(U) = divDevRhoReff(U)

                    template<class TransportModel>
                    Foam::tmp<Foam::fvVectorMatrix>
                    Foam::IncompressibleTurbulenceModel<TransportModel>::divDevReff
                    (
                        volVectorField& U
                    ) const
                    {
                        return divDevRhoReff(U);
                    }
                    
                    C 1 条回复 最后回复 回复 引用
                    • C
                      Calf.Z-DNS @队长别开枪 最后由 编辑

                      @队长别开枪 谢谢指导。
                      这么看来,OpenFOAM对可压缩跟不可压缩的问题实现了代码的复用。
                      如果是不可压的,收敛后∇⋅U=0,2/3这项自动为0,也就省略了,不知道这么理解对不对。

                      那么对于动量方程2.png
                      这种形式的动量方程,没有体现2/3μ这一项,是因为不够严谨,还是因为这一项因为∇⋅U=0而省略掉?

                      如果是后者的话,那么这就变成不可压缩的动量方程?那为何还有密度这一项?

                      还是说这是一个统称的写法,具体还要根据问题判断。谢谢。

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