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



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



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


  • OpenFOAM讲师

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



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

    入口液体体积分数为1,

    你这是模拟单相流?



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



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


Log in to reply