[Zanco]

Hi guys. I'm writing a code using the NACA0012 airfoil tutorial for rhoSimpleFoam compressible solver as a guide.

Actually I'm not able to figure out how the temperature is calculated. I mean, If I'm right everything should be done inside the hePsiThermo.C and thermoI.H files where thermo.correct() and calculate() functions are defined.

Temperature calculation is coded this way:

Test = Tnew;
        Tnew =
            (this->*limit)
            (Test - ((this->*F)(p, Test) - f)/(this->*dFdT)(p, Test));

I cannot understand what the f stands for. To me this procedure looks like a Newton algorithm.

Inside the thermophysicalProperties these are the selected parameters:

thermoType
{
    type              hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst; 
    equationOfState perfectGas; 
    specie            specie;
    energy           sensibleInternalEnergy; 
}

hConst means Cp = const. sensibleInternalEnergy means I have to use the TEs. Thus f = es and F = Es where

template<class Thermo, template<class> class Type>
inline Foam::scalar
Foam::species::thermo<Thermo, Type>::Es(const scalar p, const scalar T) const
{
    return this->Hs(p, T) - p/this->rho(p, T);
}

and

template<class Thermo, template<class> class Type>
inline Foam::scalar
Foam::species::thermo<Thermo, Type>::es(const scalar p, const scalar T) const
{
    return this->Es(p, T)*this->W();
}

So Es should be the internal energy that is defined as enthalpy minus pressure over density while es (read f) is defined as Es times molecular weight.

Does anyone of you know what it stands for and why it should be there? I cannot understand exactly which is the whole functional I have to use inside the Newton iterations.

Thanks a lot to anyone who can help me.