有关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，你怎么看，这些可以做到吗？谢谢老师

[金石为开]

@李东岳 老师我正研究两相流代码，看到里面有一个Partial Elimination Algorithm（PEA算法）这是什么东西？第一次听说，老师懂么？可否解释一下？或者提供资料看看呢？谢谢

[CFD中文网]

用来离散曳力的，曳力很重要，搞出来不同的方法，PEM，PIM。See: Spalding D B. Numerical computation of multi-phase fluid flow and heat transfer[J]. Recent advances in numerical methods in fluids, 1980, 1: 139-167.

[CFD中文网]

最近测试了几个数据，semi-conservative以及fully-conservative形式的方程求解的过程中，至少对于实验室规模的鼓泡床数据都还可以哈。

对于PEA，建议看Alberto的这个文章，相对更容易，你知道一般原始文章都比较晦涩:sunglasses: ：
Passalacqua A, Fox R O. Implementation of an iterative solution procedure for multi-fluid gas–particle flow models on unstructured grids[J]. Powder technology, 2011, 213(1): 174-187.

[liuyf]

@cfd-china 请问你的计算结果是多长时间的平均值？鼓泡床的计算网格是几维？ 谢谢

[liuyf]

@cfd-china 感觉到OpenFOAM 4.x，发现还是没有解决bubble column上部无液体区域中仍有液速的问题

[李东岳]

@liuyf 出现相分离的情况，openfoam发布的现存的代码都采用不同的方法来解决这个问题。在采用非守恒形式的时候，相分离气体的速度为气泡最终速度。在采用半守恒或者全守恒的时候。完全的相分离气体的速度等于液体的速度。这也比较符合物理的。当气体相分数非常低，气泡速度约等于液体运动速度。

目前更加困难的是壁面附近的相界面交换力的数值处理防止振荡。