包含传热传质的两流体算法

[ZY-GONG]

各位老师，前辈，我想学习下multiphaseEuler,或者是reactingTwoPhaseEulerFoam的算法，请问各位前辈有什么相关的文章可以推荐下吗？非常感谢
阅读了xiaopingqiu老师的博客，也看了东岳老师的reactingTwoPhaseEulerFoam的解析，有点疑惑,东岳老师写的连续性方程

ddt(alpha,rho) + div(alphaRhoPhi) = 0

如果考虑传质的话，右侧为什么没有源项呢？

[李东岳]

我那个没有考虑传质。如果需要传质。我可以写一下。但是这个工作排期就是不知道得啥时候了。

[ZY-GONG]

@李东岳 谢谢李老师，我先去找点相关的文章看一看吧

[ZY-GONG]

@李东岳
reactingTwoPhaseEulerFoam求解器的相变求解过程
李老师，我想请问下在org这边的multiphaseEuler，pEqn里既然考虑了传质，但是为啥找不到对能量进行的修正呢，从代码里似乎只能看到在EEqn求解后进行了thermo.correct()，是我理解有问题吗
还有就是能否请您大概提示下这个求解器的求解过程或者pEqn的迭代过程，万分感激！

[李东岳]

pEqn里既然考虑了传质，但是为啥找不到对能量进行的修正呢

这个怎么理解？

[ZY-GONG]

@李东岳 我的理解里，传质的过程应该伴随着能量的变化吧，求出正确的p后才能确定对应的温度这些？
能量方程右侧也有传质的源项，如果传质没求解的话怎么得到正确的温度呢
我的认识非常有限，如果错了请李老师批评指正！

[李东岳]

传质的过程应该伴随着能量的变化吧，求出正确的p后才能确定对应的温度这些？能量方程右侧也有传质的源项，如果传质没求解的话怎么得到正确的温度呢

我看过老版本的openfoam代码。传质是脱离流动方程之外的。不过像你说的这个描述，目前我还不100%确认。我这面算法的获得，都是靠拆openfoam代码来获得的。传质这方面我还没系统的拆解。目前只能透露这些信息了。春节期间没啥事，或许我可以抽时间把openfoam中的算法拔出来。

[ZY-GONG]

@李东岳 您的意思是在EEqn这里求解出的温度，他就默认是最终的正确温度是吗？然后在pEqn这里只是对质量速度压力密度这些的迭代求解，这里就相当于buoyantPimpleFoam里的比较EOS得到的rho和连续性方程的rho这个过程。

春节期间没啥事，或许我可以抽时间把openfoam中的算法拔出来。

期待！感谢李老师！

[李东岳]

是的。

openfoam这面PU的求解框架大体一样。