有关OpenFOAM中的双流体模型

[李东岳]

开源CFD软件OpenFOAM中的双流体模型求解器一直在持续改进。OpenFOAM-1.7.x中的bubbleFoam可谓是双流体模型在OpenFOAM中的小试牛刀，OpenFOAM-4.x中的twoPhaseEulerFoam已经非常完善，其植入了不同的界面动量传递模型、蒸发沸腾模型、湍流模型并支持不同的自定义源项。然而由于双流体模型本身的复杂性，目前的求解器依然需要进行大量的测试研究。

特在此和各位CFDer交流：

1. 是否有人对最新版本twoPhaseEulerFoam中的基于面的动量插值离散技术进行测试，OpenFOAM基金会给出的解释为此技术对于扭曲网格的适应性较好，但官方并未进行测试。基于面的动量插值离散技术

2. 若有CFDer使用过OpenFOAM中的双流体模型求解器发表过相关文章，希望能够拜读一下并学习交流。

谢谢，欢迎交流讨论。

[金石为开]

@李东岳 老师，我目前用的3.0.0版本的openfoam,这个版本的两相欧拉求解器相比于4.X版本怎么样呢？我有没有必要换版本呢？因为之前的算例学习都是在3.0.0版本下的~谢谢

[李东岳]

3.0.x之后的版本大体都差不多了。推荐用twoPhaseEulerFoam，reactingTwoPhaseEulerFoam还需要做测试，官方在收到足够反馈之后会删除twoPhaseEulerFoam，目前她还存在，表明...reactingTwoPhaseEulerFoam还不能代替twoPhaseEulerFoam。有相关对比数据的话，欢迎交流。尤其是基于面的动量插值离散技术

[金石为开]

@李东岳 那我先研究研究3.0.0当下的模型吧，，，研究好了再说。。。

[CS.Hu]

@李东岳 您好，我之前没有用of做过双流体模拟，准备近期试试，先占个楼

[李东岳]

好事是目前把openfoam中所有的双流体模型求解器综述了一遍，从非守恒形式、半守恒形式到全守恒形式，不可压缩到可压缩，低速度和亚音速，常规解法和最新的基于面的动量插值技术，简化和未简化的能量方程，不同的界面动量交换模型，各种不同的湍流模型公式全部搞完。还缺少基于面的动量插值技术验证蒜粒。

[金石为开]

@李东岳 东岳先生目前我就遇到了一个版本的问题，我用3.00版本自己设置一个算例，六自由度，自带算例是sixdofsolve，我新设置算例用velocitysolvr，计算发现输出文件和paraview出入很大，感觉是bug，所以打算用221版本，里面有用volocitysolver算例，我认为是后期开发新的motionsolver导致有点bug，所以想换版本，可是以前版本模型或许不多，算法不够高级，两头为难，怎么办呢？

[金石为开]

@金石为开 东岳先生我看代码这里面的虚拟质量力模型里面的质量系数我看文章是alpha1alpha2，但是代码里面貌似少了alpha2，您看看是不是错了，或者我理解的有问题呀？谢谢老师。

[金石为开]

@李东岳 东岳先生，您有文章讲最新的twophaseEulerFoam求解器的文章吗？老版本的很好理解，新版本的求解器我看了看很复杂，不知道该怎么理解，我也打算用，可是还没理解透彻，老师谢谢你了，还有我提的程序问题，上一条回复，不好意思回复自己了，您看看，非常感谢

[CFD中文网]

我看文章是alpha1alpha2，但是代码里面貌似少了alpha2

alpha2被除掉了吧。

[金石为开]

@cfd-china 为啥除掉呀、。有什么理由或者原因么、？还问大神一个问题，我用twophaseeulerFoam里面的constant diametermodel，常直径模型，我不明白的是为什么要定义直径，我看速度压力方程里面也没有用到粒子直径的公式，想问下哪里用到了呀？我该怎么给他们比如水，空气的直径？谢谢，一直很困扰

[李东岳]

@金石为开 在 有关OpenFOAM中的双流体模型 中说：

东岳先生，您有文章讲最新的twophaseEulerFoam求解器的文章吗？老版本的很好理解，新版本的求解器我看了看很复杂，不知道该怎么理解，我也打算用，可是还没理解透彻，老师谢谢你了，还有我提的程序问题，上一条回复，不好意思回复自己了，您看看，非常感谢

如果你说的老版本是bubbleFoam，其中界面交换项是存在alpha1，alpha2的，但是新版本的求解器没有采取bubbleFoam的求解思路，因此没有alpha1。

[李东岳]

看这个：M. Ishii, T. Hibiki, Thermo-Fluid Dynamics of Two-Phase Flow, Second Edition

[金石为开]

@李东岳 看了看，有点吃力，我没有学过热物理，没有热背景，我气液问题不涉及到传热，我打算不解能量方程了。不过还是谢谢老师

[金石为开]

@李东岳 是的，没有这项，求解思路变了，那有什么依据吗？不是Rusche那篇论文思路了？有论文依据吗？谢谢~

[李东岳]

@金石为开 Rusche中表3.1，右边的是bubbleFoam的方程，左边的是你看到的缺少系数的方程。

[金石为开]

@李东岳 {捂脸}，看来我要把那篇论文打印出来好好看看！谢谢老师。我也开始研究可压求解器了，看到您博客有一步推导，没看懂，推了几次都不对，

19式。老师您看看~~谢谢

[金石为开]

@李东岳 东岳先生，我看多相流计算需要群体平衡模型，可是目前openfoam版本里面两相欧拉貌似没有这个模型，请问下openfoam里面有类似群体模型的实现吗？借鉴借鉴，谢谢

[李东岳]

@金石为开
目前并没有群体平衡模型，官方正在开发，看Henry造化了....

[金石为开]

@李东岳 东岳先生，以你的了解，我要研究气泡流，考虑气泡之间作用，没有群体平衡模型，用其他啥方法可以做到呢？我知道的几种有IATE，DPMFoam，你怎么看，这些可以做到吗？谢谢老师