twoPhaseEulerFoam中的能量方程
-
按照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给注释掉了,我这种做法对吗?