4

While deciding what tool (library) to use for a simulation of a hybrid gas turbine plant, I found that different libraries/tools of thermodynamic properties yield different values than the other.

I made a simple test to calculate $c_p$ of air at $300\ K$(same goes for enthalpy by the way) only for simplification and comparison and here are the results.

Source              |    Cp at 300 K (J/kg.K)
--------------------------------------------
Thermopy  (python)  |   1004.89
Cantera   (python)  |   1010.06
CoolProps (python)  |   1006.627
EES Software        |   1005
Engineering Toolbox |   1005
Incropera's book    |   1007
Cengel's book       |   1007

What do you think is the reason for the difference in values?

Apparently the one closest to the properties sheets from Incropera and Cengel is CoolProps but is this level of accuracy sufficient?

The case might not be convincing for $c_p$, but it did occur to me once while calculating adiabatic flame temperature of some fuels and the error were in hundred Celsius range!

Link : used Python code

Algo
  • 2,140
  • 1
  • 20
  • 34

1 Answers1

1

What do you think is the reason for the difference in values?

It's likely that they have calculated the coefficients of whatever model they are using from different data sets (or sub-sets). It's also possible that they are using different models.

is this level of accuracy sufficient?

That depends on what you're using it for. In general you probably have an acceptable level of error in your output. To get an idea of whether your input accuracy is acceptable you can do a sensitivity analysis on your simulation to see what kind of accuracy you need for your inputs.

My guess however, is that a variation of less than 1% will probably not be a significant source of error in your output compared to the simplifications of your model, and thus the sensitivity analysis is unnecessary for the case of errors in $C_p$

However, I do not know the accuracy of the rest of the package, so I'd recommend running a sensitivity analysis on your model, and an error estimation of each of your inputs to give you an idea of how much accuracy you need and if the package you use can deliver that accuracy.

Eph
  • 1,335
  • 8
  • 26