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



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



  • 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



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



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



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


  • 管理员

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



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


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1