已经没问题了

各位大佬，我想请教下关于EOS的事，我老师想让我弄一个算液态水密度的函数，我之前根据IAPWS97写了一个，但是发现在临界区算出来的密度跟实际好像差了很多，然后又看了下P-R eos,请问各位有这种根据p T计算水密度的相关资料推荐下吗

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

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

期待！感谢李老师！

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

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

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

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

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

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

@zy-gong 好像确实可以，太惭愧了，总是问这种低级问题，学习之路任重而道远
谢谢李老师了

还是说这种也可以直接用lookupobject呢？

@李东岳
抱歉李老师，我重新组织下语言
以我有限的认识，因为求解器用的是heRhoThermo,无论是new()还是correct()似乎最后都会调用heRhoThermo::calculate()来计算热物性
而我想植入的一个气相的eos，需要的参数是液相的温度和气相的温度，大概就是rho(Tgas,Tliquid)，植入eos的方法我大概已经掌握了，目前的想法是建立一个新的thermoType，newThermo的calculate()函数是用两个温度场作为参数来计算rho的，但是我目前不知道怎么在thermo里获取另一相的温度
假设液相的温度，我可以在多项欧拉求解器里通过fluid.phases()[0].T()来访问，那在newThermo里应该怎么做才能读到液相的温度场呢？
感谢李老师!