CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    TwoPhaseEulerFOAM的一些小问题

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

      我正在准备把在这个Foam上面植入拉格朗日模型,有一点小问题
      (1)求解器里面似乎没有求解连续性方程,把连续性方程带到了动量方程里面,不知道我的理解是否正确
      (2)因为要植入拉格朗日模型,Euler方法里面的曳力模型,粒子判断什么都要关掉,phaseProperties不设置任何东西,是不是就认为关掉了@东岳
      @zwl 老哥似乎做过空化的,连续性方程有源项是不是需要处理后用到动量方程的求解里面

      Z 1 条回复 最后回复 回复 引用
      • Z
        zwl @yfclark 最后由 编辑

        @yfclark 源相是需要加到压力泊松方程里面,当然由于源相的产生导致的动量源相和质量源相也要加到对应方程中的。不设置任何东西应该是不考虑曳力模型的。不清楚你要怎么加拉格朗日模型,为什么不考虑相间曳力?

        1 条回复 最后回复 回复 引用
        • Y
          yfclark 讲师 最后由 编辑

          @zwl 谢谢老哥,我的设计是连续的两相用Euler-Euler方法,第三相是离散相,用拉格朗日方法,第三相与前面两相有质量和动量交换,Euler和Euler之间没有任何关系,所以我打算关掉这些模型

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

            @yfclark 在 TwoPhaseEulerFOAM的一些小问题 中说:

            Euler和Euler之间没有任何关系,

            这是什么工况行为/物理行为?

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

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

            1 条回复 最后回复 回复 引用
            • Y
              yfclark 讲师 最后由 编辑

              做雾化啊,连续的气相和液相用Euler-Euler,两相没有掺混,就没有质量,动量交换,液相破碎生成的液滴用拉格朗日追踪,更@zwl 老哥的那个空化近场雾化很相似,不过他用的是Euler追踪生成的液滴和气泡,因此Euler-Euler之间需要引入动量和质量交换

              1 条回复 最后回复 回复 引用
              • Y
                yfclark 讲师 最后由 编辑

                这个连续性方程是怎么引入的,质量方程没有源项的话,就应该是0吧。
                Snipaste_2019-06-19_11-06-19.png
                我对比的是OpenFOAM2.3.x的代码:

                        U1Eqn =
                        (
                            fvm::ddt(alpha1, rho1, U1) + fvm::div(alphaRhoPhi1, U1)
                          - fvm::Sp(contErr1, U1)
                          + mrfZones(alpha1*rho1 + virtualMassCoeff, U1)
                          + phase1.turbulence().divDevRhoReff(U1)
                         ==
                          - liftForce
                          - wallLubricationForce
                          - turbulentDispersionForce
                          - virtualMassCoeff
                           *(
                                fvm::ddt(U1)
                              + fvm::div(phi1, U1)
                              - fvm::Sp(fvc::div(phi1), U1)
                              - DDtU2
                            )
                

                @东岳

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

                  @yfclark 在 TwoPhaseEulerFOAM的一些小问题 中说:

                  这个连续性方程是怎么引入的

                  http://dyfluid.com/icoFoam.html 我建议你看看这个,双欧拉模型跟这个是一样的

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

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

                  1 条回复 最后回复 回复 引用
                  • Y
                    yfclark 讲师 最后由 编辑

                    @东岳 非常感谢东岳老师,是我的问题没有问清楚,我只是好奇求解动量方程时引入的fvm::Sp(contErr1, U1)即0244e636-3e3a-48a7-96f7-fa1a1776385f-image.png
                    只求解一个动量方程的求解器都没有做这样的处理,因此有些好奇
                    正在研究OpenFOAM2.3的twoPhaseEulerFoam,有问题还要请教东岳老师,请不要见怪

                    1 条回复 最后回复 回复 引用
                    • Y
                      yfclark 讲师 最后由 编辑

                      @东岳 @zwl 您好,我想问一下您写求解器的时候表面张力是如何处理的,我看multiphaseEulerFoam和reactingmultiphaseFoam都有surfaceTension的计算,但是在momentumtransfer里面却没有表面张力,我现在的思路主要是两种:
                      (1)在动量方程中左侧直接植入表面张力
                      (2)在动量方程中忽略表面张力,使用momentumPredict添加表面张力

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

                        @yfclark 处于有界性的考虑,https://openfoam.org/release/2-2-0/numerics-boundedness/

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

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

                        1 条回复 最后回复 回复 引用
                        • Y
                          yfclark 讲师 最后由 编辑

                          @东岳 明白了,谢谢东岳老师

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