请教贴：heThermo类中属性he_是内能还是焓的问题

[yingqing]

在rhoCentralFoam求解器中，存在以下代码

solve
        (
            fvm::ddt(rhoE)
          + fvc::div(phiEp)
          - fvc::div(sigmaDotU)
        );
        
		dimensionedScalar Emin("Emin", dimEnergy/dimMass, 71800);
        e =max(rhoE/rho - 0.5*magSqr(U),Emin);
        e.correctBoundaryConditions();

控制方程求解出的是总能，获得内能e。
热物性表constant/thermophysicalProperties中设置的
thermo 是hConst，这意味着he_应该是比焓
energy 设置应该是把he_传递到sensibleEnthalpy类中，使焓来转化为温度感觉这个合理；
但thermo 是hConst
energy设置为sensibleInternalEnergy，这意味着he_应该是焓，传递给sensibleInternalEnergy来转化温度，感觉不是很合理呀。这个地方有理解的老师，能给讲讲嘛，已经懵了

[yingqing]

和一位师兄交流了一下。是我的判断搞错了，在heThermo类中he_属性取决于energy的选择。如果energy 是sensibleEnthalpy，那么he_=h;如果energy是sensibleInternalEnergy，那么he_=e。对于rhoCentralFoam求解的控制方程是总能，计算的是e。所以energy只能使用sensibleInternalEnergy。
我去openFOAM带的rhoCentralFoam案例库查了一下,确实都是sensibleInternalEnergy！！！

[李东岳]

是的，sensibleEnthalpy是h，sensibleInternalEnergy是e
不过这俩应该不影响rhoCentralFoam，他能量方程自己硬植入了