[avk1985]

I was going through requirements to input Cp polynomial for a real gas modelled with peng robinson model I came across following information

Source hPolynomialThermoI.H

template<class EquationOfState, int PolySize>
inline Foam::scalar Foam::hPolynomialThermo<EquationOfState, PolySize>::Cp
(
const scalar p, const scalar T
) const
{
return CpCoeffs_.value(T) + EquationOfState::Cp(p, T);
}

Source hPolynomialThermo.H

The polynomial expression is evaluated as so:

\f[
Cp = 1000 - 0.05 T + 0.003 T^2
\f]

Note
- Heat of formation is inputted in [J/kg], but internally uses [J/kmol]
- Standard entropy is inputted in [J/kg/K], but internally uses [J/kmol/K]
- Specific heat at constant pressure polynomial coefficients evaluate to an
expression in [J/(kg.K)].

in hPolynomialThermoI Cp seems to be calculated as summation of 2 components polynomial expansion and EquationOfState::Cp(p, T) so whether polynomial needs to be evaluated subtracting the value of EquationOfState::Cp(p, T) from actual Cp value (test data, refprop data etc)
eg Polynominal (T) = Cpactual(T) - EOSCP(p,T)

but hPolynomialThermo doesnt says so

Can some one please help clarify this ?

Regards
Amod