可压缩求解器温度场计算结果疑问

[zwl]

这是温度场云图；
利用reactingTwoPhaseEulerFoam计算两相流动，利用perfectFluid模型计算液体压缩性，thermo选择hConst，能量方程对应sensibleEnthalpy，定入口温度，入口液体体积分数为1，入口压力10MPa，出口压力1MPa，这个计算结果有两个问题：
1.根据所选模型代码，温度不应该随压力变化，计算结果为什么却显示受压力影响？
2.在温度受压力影响的基础上，为什么入口温度会出现那么大变化？给的定温入口边界温度对应的不是入口压力？
请各位做可压缩求解器的大佬赐教，指点一二，这是什么原因导致的呢？先行谢过。

[zwl]

@东岳 好像东岳老师之前用twophaseEulerFoam，不知道有没有遇到过这个问题？能解答一下吗？谢谢。

[cccrrryyy]

没用过两相的可压求解器，只用过单相的。你这里面温度为什么不应该随压力变化呀？刚看了perfectFluid定义，密度是温度和压力的函数，温度是根据焓值求出来的，之后再更新psi mu之类的，然后再解压力方程，更新密度，再重新求解动量方程，如此往复。这一套下来，温度的分布、速度的分布等等都跟压力场是耦合在一起的。

[东岳]

@zwl 在 可压缩求解器温度场计算结果疑问 中说：

入口液体体积分数为1，

你这是模拟单相流?

[zwl]

@东岳 出口是空气场，感觉空气的温度场没问题，所以就没截出来，就是水喷到空气里，空气cpMcv＝R,而液体的为0.所以焓算出来的液体温度随压力变化那么大，我感觉有问题？东岳老师觉得呢？

[东岳]

不清楚。我一般通过phaseProperties里面的type pureIsothermalPhaseModel关掉传热。如果不关掉，肯定温度和压力耦合一起了。
我觉得可能跟你的进出口压力边界条件有关系。