twoPhaseEulerFoam中的能量方程

[Dingcy]

按照twoPhaseEulerFoam中的EEqns.H写出能量方程如下（算例文件中能量形式是sensibleInternalEnergy）大神们先看看这么写对吗？

fvScalarMatrix E1Eqn
(
    fvm::ddt(alpha1, rho1, he1) + fvm::div(alphaRhoPhi1, he1)
  - fvm::Sp(contErr1, he1)

  + fvc::ddt(alpha1, rho1, K1) + fvc::div(alphaRhoPhi1, K1)
  - contErr1*K1
  + (
        he1.name() == thermo1.phasePropertyName("e")
      ? fvc::div(fvc::absolute(alphaPhi1, alpha1, U1), p)
      + p*fvc::ddt(alpha1)
      : -alpha1*dpdt
    )

  - fvm::laplacian
    (
        fvc::interpolate(alpha1)
       *fvc::interpolate(thermo1.alphaEff(phase1.turbulence().mut())),
        he1
    )
);

E1Eqn.relax();

E1Eqn -=
(
    Kh*(thermo2.T() - thermo1.T())
  + Kh*he1/Cpv1
  - fvm::Sp(Kh/Cpv1, he1)
  + alpha1*rho1*(U1&g)
  + fvOptions(alpha1, rho1, he1)
);

那如果不考虑传热的话，只把导热系数设为0和温度设为一样，这样做对吗？
我为了不考虑传热，是直接把twoPhaseEulerFoam.C中的EEqns.H给注释掉了，我这种做法对吗？