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

[zhe]

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

[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]

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

[qingdong_wang]

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

[qingdong_wang]

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

[李东岳]

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

[zhe]

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