OpenFOAM-4.x中reactingTwoPhaseEulerFoam中的Partial-elimination phase-flux

[李东岳]

双流体模型中的曳力一般情况下占主导作用，曳力的必要处理有利于数值的稳定性，在OpenFOAM-3.x之前，曳力通semi-implicit方式进行处理，在OpenFOAM比较新的版本中，曳力的处理采用了一种更为稳健的方式。

在链接文本 压力方程中，基金会将植入的算法注释为Partial-elimination phase-flux，在研究过文献之后，这个注释并不正确。植入的方程和Spalding的PEA算法不符。经过推导后，基金会这个处理方法应该也是自创的，通过相对速度强耦合起来。

不明确为什么标注PEA。如果有研究这个的可以注意一下。

The new algorithm implemented in reactingTwoPhaseEulerFoam for drag calculation is not partial elimination method of Spalding.

为了便于google收录，写个英文。

[李东岳]

17年只是粗略看了一下，感觉并不是PEA算法，2年后终于严格的推出来了，这个不是PEA算法！OpenFOAM的document需要好好整整。OpenFOAM官方确实有这种参考文献标识不正确的先例。

这几天整理一下发个reactingTwoPhaseEulerFoam解析

[李东岳]

http://dyfluid.com/reactingTwoPhaseEulerFoam.html 相应的算法写在这里了 我觉得还是非常有意思的