CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    MPPICInterFoam速度公式與interFoam用的一樣嗎?結果很不同.

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

      大家好,我現在在用MPPICInterFoam模擬三相流.首先在不加入顆粒的情況下,對求解器進行了驗證模擬.我的實驗是流體繞柱.圖中第一個是用MPPICInterFoam進行一段時間後的水平速度和垂直速度的圖像.第二個圖是用interFoam做的.比較真實實驗數據,第二個圖更接近.想問問大家,MPPICInterFoam的速度公式跟interFoam的哪裏不一樣了?謝謝大家.!123.png

      1 条回复 最后回复 回复 引用
      • Q
        qingdong_wang 最后由 编辑

        MPPICInterFoam 在左侧加了

        fvm : : ddt ( alphacRho ,  U)+ MRF. DDt( alphacRho ,  U)
        

        应该是 When a particle is transformed into a Eulerian structure, its velocity and volume fraction is added to the continuum fluid

        右侧是the effect of particle motion on the momentum of surrounding fluid。 也即是我们说的粒子力。

        phicForces / rAUcf   //   is   added   as   particle   force 
        

        理论上 , 当粒子为0的时候, 由于alphac为0, 这两项也应为为0才对。理论上是没有影响的。
        不过MPPICInterFoam的alphaEqn改了很多, 说实话,这个我没有看明白。
        不过你这两个case是除了solver不一样外, 其他全部一样吗?比如说fvSchemes或者fvSolutions等等

        你可以看看下面这个文档, 虽然他是用的openfoam3, 但是感觉会对你有帮助

        pg11 http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2016/EbrahimGhahramani/Report_Ghahramani.pdf

        zhe 1 条回复 最后回复 回复 引用
        • zhe
          zhe @qingdong_wang 最后由 编辑

          @qingdong_wang 謝謝你的回覆.除了求解器不同,其他都是相同的,包括fvSchme裏的插值方法.
          我也看到在MPPICInterFoam中添加了對於粒子的方程式.但是因爲我關掉了粒子的添加和耦合,那麼也應該代表粒子力不存在了吧?還有一個,alphac代表的是連續相的體積分數,MPPICInterFoam是將氣體和液體用同一個體積分數代表了,也就是alphac,那麼在我的案例裏,它應該是恆等於1,不是0.除了方程兩邊對於粒子的公式外,如果將alphac等於1,那麼公式跟interFoam是一樣的.
          我會仔細閱讀發來的文獻,並且再檢查一下是不是對於湍流模型的不同要求而有差別的.比如換做LES試試,隨後再更新在這裏.如果有朋友想多討論MPPICInterFoam,也可以互相討論啊,謝謝.

          Q 1 条回复 最后回复 回复 引用
          • Q
            qingdong_wang 最后由 编辑

            谢谢你的指正。
            不过我没看湍流模型, 我用的DNS, 不过看湍流的构造函数的时候,他们的确是不同的。

            1 条回复 最后回复 回复 引用
            • Q
              qingdong_wang @zhe 最后由 编辑

              @zhe 我能加你一下微信吗?我最近也在做这个问题。因为我直接从interfoam改MPPICInterfoam, 所以我觉得还是有些问题。我的微信是 wqdwqdaaa

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

                你那个结果在中间处区别比较大,可以发图看看中间处是什么东西

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

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

                  @东岳 感謝东岳老師的回覆.首先說個題外話,最近在英國這邊一個公衆號上看到您發的一片勸解博後選擇的文章.因爲是轉發的,起初不知道是您發的,後來越看越像您的口吻,最後在下面看到"轉自李东岳"的時候,心裏很是激動.謝謝东岳老師對後輩們的勸解!
                  我又做了一個用MPPICInterFoam在沒有粒子的情況下,繞柱(光滑壁面)的模擬.這次用了LES湍流模型.效果很好.並且並行運算時,LES比K-omega SST運行更快.大概是因爲這種情況更符合實際情況吧.跟實驗數據比較,甚至比之前interFoam做的還要好.圖中我做了MPPICInterFoam和interFoam的對比.上面的是MPPICInterFoam的.在柱後水流形成的湍流旋渦非常明顯.在不同高度的對比,效果也很好.至少說明,在水動力這塊兒,這個求解器配合LES,可以得到不錯的效果.當然,邊界條件也很重要.
                  有問題的大家一起討論啊.
                  3.png
                  3.png

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