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

[zhe]

大家好，我研究的方向是海底水流对泥沙的冲刷。其中包括海面上的空气流，海水和海底的沙石颗粒。我个人的理解应该是一种不相容的三项流的模拟。所以我想用interFoam加DPMFoam来完成。但是有个问题：应该在interFoam里加入颗粒还是应该在现有的DPMFoam里把它的单项流换成我想要的双项流？希望各位能多多帮忙，互相学习。谢谢。

[yfclark]

基金会做了一个mppicinterfoam，基本是是这样的，但是openfoam的拉格朗日追踪还是太落后了，它没法处理粒子比网格大的情况,因此欧拉项界面的精细解析没法耦合在一起，有不少人在开发这样的三相流求解器，我有看到解决了这个问题的论文，但是我还没有看到在openfoam上开源的代码

[yfclark]

如果你的沙子很小，且不会越过气液界面，我想基金会的mppicinterfoam能够满足你的要求

[zhe]

@yfclark 谢谢您的回复。我的项目的确满足您说的这种情况。沙子的尺寸很小，并且不会跃出水面。但是我看了MPPICInterFoam，这个仍然是解决不相容的两相流的问题。怎么添加固体颗粒啊？希望能得到您的再次指导，谢谢。

[yfclark]

@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简单些。但都属于拉格朗日跟踪

[zhe]

@星星星星晴 谢谢您的回复。但我有一个担心的问题就是，如果在interFoam里加入颗粒的话，需要更多的去考虑液固耦合，和颗粒与颗粒之间的碰撞机制。这样就不只是加入颗粒那么简单了。而如果在DPMFoam中换作两相流的话，是不是自带的就解决了流固耦合，和颗粒之间的碰撞的问题？

[zhe]

@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进行调整。

[zhe]

@星星星星晴 的确像您说的，如果在DPMFoam中加入interFoam的话，alpha的这部分很难。我尝试过，但是加入两相流的话还要在考虑颗粒占的体积分数。需要额外定义一个alpha3。惭愧的是，我一直也没生成这个求解器。可能是从理论上都没过关。我再试试您给的这个ppt上的方法。同时，上面yfclark提到的mppicinterfoam也是值得尝试的一种方法。我会接着在这个论租里更新我的情况。有问题了，还希望您能抽时间帮我看一下。真心感谢。也希望能帮到以后用到这个方法的同学。谢谢。

[星星星星晴]

@zhe 哈哈 我不做VOF 但是我做Lagrangian 这部分。
我觉得alpha3应该就是一个scalar 然后写出场即可吧？
另外PPT中只是 固液 coupling，如果涉及到气泡什么的，particle particle就需要自己根据模型写啦。

[zhe]

@星星星星晴 是，只是一个scalar。我是新能源和环境工程出身，对code的确很陌生。但是我真心觉得OpenFoam可以用在一些具体工程问题上，所以很想把他学好。谢谢您一直耐心回复，我会多个方法都试试，寻求出较好的方法。多谢！

[星星星星晴]

@zhe it's easy 慢慢来，我也是半路出家到多相流，做多相流的模拟不到一年。。

[zhe]

@星星星星晴 以后可以多交流交流。我现在还在利物浦大学读博，研究还要很久，并且之后也会接着再做学术研究，希望咱们能多交流多互相借鉴！

[星星星星晴]

@zhe :

[zhe]

@yfclark 您好，我遇到一个小问题。我用的是openfoam7的版本。MPPICInterFoam是由openfoam.com发布的，只有openfoam-1612v以上的版本才有。我需要换成另外的版本吗？还是说可以下载openfoam.com上的求解器到openfoam 7上？等您回复，谢谢。

[zhe]

@东岳 您好，李老师。我下载了MPPICInterFoam后，想要安装。但是一run就出来了错误。希望您能帮我看一下，看看是哪里出了问题，谢谢您。!

[yfclark]

换基金会版本，两个源代码不通用

[zhe]

已解决，路径问题。不好意思。

[xiaolin]

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

[zhe]

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