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)
    );
    

    a8753d25-8ea0-4cf2-8fae-10a71aa5cdfd-image.png
    那如果不考虑传热的话,只把导热系数设为0和温度设为一样,这样做对吗?
    我为了不考虑传热,是直接把twoPhaseEulerFoam.C中的EEqns.H给注释掉了,我这种做法对吗?


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2