气液固三项流，耦合interFoam和DPMFoam可以实现吗？

[xiaolin]

@zhe 您好，我也遇到了和你一样的问题，劳烦问一下，该怎么解决呢

[zhe]

@xiaolin 我是因爲安裝了不同版本的OpenFOAM而導致的問題.比如說OpenFOAM6和OpenFOAM18.06,這種不同編號的版本是不同組織創建的.所以裏面有一些東西不兼容.你需要下載相應的版本.

[wangtao121]

@zhe 您好，你问题解决了吗？我想问一下后续你是interFoam加入颗粒，还是把DPM单项流给换掉了？谢谢

[zhe]

@wangtao121 之后发现，在openfoam中multiphase的solver中，有一个MPPICInterFoam的solver，把其中的MPPIC换成DPM就可以解决这个问题。实际上是在interfoam里加入了粒子效果。

[马乔]

像这样子嘛？

[OItoCFD]

@马乔 您好 请问 这种是直接把颗粒的lib加到interfoam里就可以吗？粒子从水里到空气中需要做些特殊处理吗？感觉不做也可以呀？

[马乔]

@oitocfd 直接用MPPICInterFoam就可以了

[OItoCFD]

@马乔 好的谢谢

[星星星星晴]

@zhe 雷猴啊，同事要做三相的模拟，请问你做的咋样啦～

[小狗狗]

@星星星星晴 你好，我想请教一下，欧拉-拉格朗日方法中，DPM方法描述的拉格朗日颗粒是不是假设半径为0的点？

[小狗狗]

@李东岳 你好，我想请教一下，欧拉-拉格朗日方法中，DPM方法描述的拉格朗日颗粒是不是假设半径为0的点？

[星星星星晴]

@小狗狗 怎么说呢，parcel是point mass，有质量，有体积的particle的集合。
请问你说的假设半径为0指的是什么？
DPMFoam的parcel是有半径的，但是不会像DEM那样真正的占据一定体积。就是一个数，所以euler-Lagrange才对parcel占网格体积比有要求

[小狗狗]

@星星星星晴 我的意思是在计算粒子轨迹过程中，比如轨迹和网格的几何相交关系，是不是把粒子当作质点来处理，不考虑粒子的实际的体积

[星星星星晴]

@小狗狗 我没处理过这样的问题，无法给你准确答案。但是我认为是不可以的，就类似你计算drag的时候同样也要考虑半径，当你考虑空间关系的话，应该也需要将其考虑为真正的小球。比如用拉格朗日模拟气泡壁面的话，肯定也要考虑气泡与避免的距离一样。

[小狗狗]

@星星星星晴 你说的应该是物理模型层面的东西，算拽力之类的，会给定体积来算的。我说的仅仅是颗粒追踪过程粒子与网格的几何关系，举个例子，一个粒子从一个单元内穿过单元面进入到相邻单元，如何判定这个粒子是否穿过面？是不是就把粒子当作质点来计算？还是说把粒子当作是有实际半径的球来计算？

[星星星星晴]

@小狗狗 你这么说的话 就是质点， 不考虑半径

[小狗狗]

@星星星星晴 嗯嗯，谢谢

[OItoCFD]

@小狗狗 是的 openfoam里这种颗粒，算力是拖曳力系数用类似莫里森方程来算，算是概化模型吧，运动那些都是当为一个支点，不是grid-resolved。如果是精细考虑颗粒，每个颗粒表面也要画网格，受力是根据颗粒表面附近流体节点压力积分来算，再算总的这个颗粒平动转动。

[小狗狗]

@oitocfd 嗯嗯，明白了，谢谢

[zhe]

@星星星星晴 不好意思，一直在工作和学习两边忙，没有怎么来论坛。都还好，谢谢。至于openfoam的粒子追踪模型，都是不考虑粒子体积的，只是一个点。但在欧拉下计算NS方程，是考虑了粒子占网格的比值（这也不是考虑了体积）这只是考虑了由于粒子存在在网格，占据了原本流体的一部分体积，对计算流体造成了影响。总的来说，还是算的流，粒子只是一个源项。