research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Fixed errors in implementation of carbon dioxide equation of state

Browse source 
The sign error in computing the entroy term proportional to G[0]
introduced errors in the entropy of about 3 percent.
This commit is contained in:
Ray Speth 2012-05-07 18:20:20 +00:00
parent a4520bfe85
commit eb12c611db
1 changed files with 3 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

10
src/tpx/CarbonDioxide.cpp
View file

@ -24,11 +24,10 @@ static const double Roc=464.00; // [kg/m^3] critical density
static const double To=216.54; // [K] reference Temperature
static const double R=188.918; // [] gas constant for CO2 J/kg/K
static const double Gamma=5.0E-6; // [??]
static const double u0=3.217405E5; // [] internal energy at To
static const double u0=3.2174105E5; // [] internal energy at To
static const double s0=2.1396056E3; // [] entropy at To
static const double Tp=250; // [K] ??
//static const double Pc=7.38350E6; // [Pa] critical pressure
static const double Pc=7.3817589E6; // Adjusted to fit the actual functional form
static const double Pc=7.38350E6; // [Pa] critical pressure
static const double M=44.01; // [kg/kmol] molar density
/*
@ -95,9 +94,6 @@ static const double G[]= {
-1.667825E-3,
7.305950E-7,
-1.255290E-10,
3.2174105E5,
2.1396056E3
};
/*
@ -279,7 +275,7 @@ double CarbonDioxide::sp()
}
sum += G[1]*log(T/To);
sum -= G[0]*(1.0/To - 1.0/T);
sum -= G[0]*(1.0/T - 1.0/To);
for (int i=0; i<=6; i++) {