气液固三项流，耦合interFoam和DPMFoam可以实现吗？
-
大家好，我研究的方向是海底水流对泥沙的冲刷。其中包括海面上的空气流，海水和海底的沙石颗粒。我个人的理解应该是一种不相容的三项流的模拟。所以我想用interFoam加DPMFoam来完成。但是有个问题：应该在interFoam里加入颗粒还是应该在现有的DPMFoam里把它的单项流换成我想要的双项流？希望各位能多多帮忙，互相学习。谢谢。
-
基金会做了一个mppicinterfoam，基本是是这样的，但是openfoam的拉格朗日追踪还是太落后了，它没法处理粒子比网格大的情况,因此欧拉项界面的精细解析没法耦合在一起，有不少人在开发这样的三相流求解器，我有看到解决了这个问题的论文，但是我还没有看到在openfoam上开源的代码
-
如果你的沙子很小，且不会越过气液界面，我想基金会的mppicinterfoam能够满足你的要求
-
@yfclark 谢谢您的回复。我的项目的确满足您说的这种情况。沙子的尺寸很小，并且不会跃出水面。但是我看了MPPICInterFoam，这个仍然是解决不相容的两相流的问题。怎么添加固体颗粒啊？希望能得到您的再次指导，谢谢。
-
@zhe 你看看官方提供的算例，可以人工入射拉格朗日粒子的，也可以参考其他拉格朗日求解器的设置
-
理论上可以，
同事在做三相bubble column 里面有颗粒。
Bubble部分基于interFoam编写，particle部分使用的就是kinematic parcel。
准确来说是自己修改的interFoam 和 Euler - Lagrangian的 Hard sphere.
DEM和DPM的区别就是Hard sphere 和soft sphere的区别嘛~ps:我到现在一直不理解DEM和 Euler- Lagrangian的区别。
我个人觉得，再interFoam中加入颗粒应该会比较简单
主要修改两个地方，evolve parcel 和 修改动量方程源项（假设不考虑其他影响）
建议看一下ReactingParcelFoam（及更上层的sprayFoam），这个应该是整个Lagrangian model中包含各种模型最全的。
以及下方
https://www.foamacademy.com/wp-content/uploads/2018/03/particles_slides.pdf
-
欧拉拉格朗日倾向于是一种框架。连续相在欧拉框架下算，离散型在拉格朗日框架下。欧拉框架有很多算法，N-S方程是一种算法，矩方法是一种算法。同样拉格朗日框架下有很多算法。DEM属于一种算法。还比如MPPIC，PIC。DPM要比DEM简单些。但都属于拉格朗日跟踪
-
@星星星星晴 谢谢您的回复。但我有一个担心的问题就是，如果在interFoam里加入颗粒的话，需要更多的去考虑液固耦合，和颗粒与颗粒之间的碰撞机制。这样就不只是加入颗粒那么简单了。而如果在DPMFoam中换作两相流的话，是不是自带的就解决了流固耦合，和颗粒之间的碰撞的问题？
-
@yfclark 谢谢，谢谢。我看了它的机制，觉得可以一试。过段时间如果再有问题，会继续向您请教。如果可以成功使用，我也会留言给以后用到的同学。感谢🙏
-
@zhe
情况是这样的，不知道我理解的对不对。
毕竟我不做VOF这一块。。interFoam 直接对于alpha的处理已经足够了。所以你的mixsolver中比较南的部分就已经解决了。
如果你反着来，implement 就可能会更复杂一点。
如果参考下面给你的那个PPT，需要添加源相即可。至于你说的碰撞机制，再sprayFoam中就有 stochastic collision model（O'Rouke）。
以及我同事的论文
Influence of inter‐particle collisions and agglomeration on cyclone performance and collection efficiency（2019）
作者
Oscar Lino Sgrott, Martin Sommerfeld颗粒碰撞的问题就操作比较复杂了，你要深入到src/Lagrangian/intermediate的深层次一点了。
对Parcel进行调整。
-
@星星星星晴 的确像您说的，如果在DPMFoam中加入interFoam的话，alpha的这部分很难。我尝试过，但是加入两相流的话还要在考虑颗粒占的体积分数。需要额外定义一个alpha3。惭愧的是，我一直也没生成这个求解器。可能是从理论上都没过关。我再试试您给的这个ppt上的方法。同时，上面yfclark提到的mppicinterfoam也是值得尝试的一种方法。我会接着在这个论租里更新我的情况。有问题了，还希望您能抽时间帮我看一下。真心感谢。也希望能帮到以后用到这个方法的同学。谢谢。
-
@zhe 哈哈 我不做VOF 但是我做Lagrangian 这部分。
我觉得alpha3应该就是一个scalar 然后写出场即可吧？
另外PPT中只是 固液 coupling，如果涉及到气泡什么的，particle particle就需要自己根据模型写啦。
-
@星星星星晴 是，只是一个scalar。我是新能源和环境工程出身，对code的确很陌生。但是我真心觉得OpenFoam可以用在一些具体工程问题上，所以很想把他学好。谢谢您一直耐心回复，我会多个方法都试试，寻求出较好的方法。多谢！
-
@zhe it's easy 慢慢来，我也是半路出家到多相流，做多相流的模拟不到一年。。
-
@星星星星晴 以后可以多交流交流。我现在还在利物浦大学读博，研究还要很久，并且之后也会接着再做学术研究，希望咱们能多交流多互相借鉴！
-
@zhe :
-
@yfclark 您好，我遇到一个小问题。我用的是openfoam7的版本。MPPICInterFoam是由openfoam.com发布的，只有openfoam-1612v以上的版本才有。我需要换成另外的版本吗？还是说可以下载openfoam.com上的求解器到openfoam 7上？等您回复，谢谢。
-
@东岳 您好，李老师。我下载了MPPICInterFoam后，想要安装。但是一run就出来了错误。希望您能帮我看一下，看看是哪里出了问题，谢谢您。!
-
换基金会版本，两个源代码不通用
-
已解决，路径问题。不好意思。
-
@zhe 您好，我也遇到了和你一样的问题，劳烦问一下，该怎么解决呢
-
@xiaolin 我是因爲安裝了不同版本的OpenFOAM而導致的問題.比如說OpenFOAM6和OpenFOAM18.06,這種不同編號的版本是不同組織創建的.所以裏面有一些東西不兼容.你需要下載相應的版本.
