mirror of
https://github.com/Cantera/cantera.git
synced 2025-02-25 18:55:29 -06:00
154 lines
3.8 KiB
C++
154 lines
3.8 KiB
C++
/*
|
|
* $Id$
|
|
*/
|
|
#include "stdio.h"
|
|
#include "ct_defs.h"
|
|
#include "ctexceptions.h"
|
|
#include "global.h"
|
|
#include "xml.h"
|
|
#include "ctml.h"
|
|
#include "WaterPDSS.h"
|
|
#include "ThermoPhase.h"
|
|
#include <new>
|
|
using namespace std;
|
|
using namespace Cantera;
|
|
|
|
int main () {
|
|
|
|
double dens, u, s, h, cv, cp, pres;
|
|
try {
|
|
Cantera::ThermoPhase *nnn = 0;
|
|
Cantera::WaterPDSS *w = new Cantera::WaterPDSS(nnn, 0);
|
|
|
|
double rho;
|
|
|
|
/*
|
|
* Print out the triple point conditions
|
|
*/
|
|
double temp = 273.16;
|
|
pres = w->satPressure(temp);
|
|
printf("psat(%g) = %g\n", temp, pres);
|
|
double presLow = 1.0E-2;
|
|
temp = 298.15;
|
|
double oneBar = 1.0E5;
|
|
|
|
printf("Comparisons to NIST: (see http://webbook.nist.gov):\n\n");
|
|
|
|
w->setDensity(1.0E-8);
|
|
w->setState_TP(temp, presLow);
|
|
double h = w->enthalpy_mole();
|
|
printf("H0(298.15) = %g J/kmol\n", h);
|
|
double h298 = h;
|
|
|
|
double s = w->entropy_mole();
|
|
s -= GasConstant * log(oneBar/presLow);
|
|
printf("S0(298.15) = %g J/kmolK\n", s);
|
|
|
|
|
|
double T[20];
|
|
T[0] = 298.15;
|
|
T[1] = 500.;
|
|
T[2] = 600.;
|
|
T[3] = 1000.;
|
|
|
|
double Cp0, delh0, delg0, g;
|
|
|
|
printf("\nIdeal Gas Standard State:\n");
|
|
printf (" T Cp0 S0 "
|
|
" -(G0-H298)/T H0-H298\n");
|
|
printf (" (K) (J/molK) (J/molK) "
|
|
" (J/molK) (kJ/mol)\n");
|
|
for (int i = 0; i < 4; i++) {
|
|
temp = T[i];
|
|
w->setState_TP(temp, presLow);
|
|
h = w->enthalpy_mole();
|
|
delh0 = h - h298;
|
|
g = w->gibbs_mole();
|
|
delg0 = (g - h298)/temp + GasConstant * log(oneBar/presLow);
|
|
Cp0 = w->cp_mole();
|
|
s = w->entropy_mole();
|
|
s -= GasConstant * log(oneBar/presLow);
|
|
printf("%10g %10g %13g %13g %13g\n", temp, Cp0*1.0E-3, s*1.0E-3,
|
|
-delg0*1.0E-3, delh0*1.0E-6);
|
|
}
|
|
printf("\n\n");
|
|
|
|
temp = 298.15;
|
|
w->setDensity(1000.);
|
|
w->setState_TP(temp, oneBar);
|
|
h = w->enthalpy_mole();
|
|
printf("H_liq(298.15, onebar) = %g J/kmol\n", h);
|
|
double h298l = h;
|
|
s = w->entropy_mole();
|
|
printf("S_liq(298.15, onebar) = %g J/kmolK\n", s);
|
|
|
|
|
|
T[0] = 273.19;
|
|
T[1] = 298.15;
|
|
T[2] = 300.;
|
|
T[3] = 373.15;
|
|
T[4] = 400.;
|
|
T[5] = 500.;
|
|
printf("\nLiquid 1bar or psat Standard State\n");
|
|
printf (" T press psat Cp0 S0 "
|
|
" -(G0-H298)/T H0-H298\n");
|
|
printf (" (K) (bar) (bar) (J/molK) (J/molK)"
|
|
" (J/molK) (kJ/mol)\n");
|
|
|
|
for (int i = 0; i < 6; i++) {
|
|
temp = T[i];
|
|
double psat = w->satPressure(temp);
|
|
double press = oneBar;
|
|
if (psat > press) {
|
|
press = psat*1.002;
|
|
}
|
|
w->setState_TP(temp, press);
|
|
h = w->enthalpy_mole();
|
|
delh0 = h - h298l;
|
|
g = w->gibbs_mole();
|
|
delg0 = (g - h298l)/temp;
|
|
Cp0 = w->cp_mole();
|
|
s = w->entropy_mole();
|
|
printf("%10g %10g %12g %13g %13g %13g %13g\n", temp, press*1.0E-5,
|
|
psat*1.0E-5,
|
|
Cp0*1.0E-3, s*1.0E-3,
|
|
-delg0*1.0E-3, delh0*1.0E-6);
|
|
}
|
|
|
|
printf("\nLiquid Densities:\n");
|
|
printf (" T press psat Density molarVol "
|
|
"\n");
|
|
printf (" (K) (bar) (bar) (kg/m3) (m3/kmol)"
|
|
"\n");
|
|
for (int i = 0; i < 6; i++) {
|
|
temp = T[i];
|
|
double psat = w->satPressure(temp);
|
|
double press = oneBar;
|
|
if (psat > press) {
|
|
press = psat*1.002;
|
|
}
|
|
w->setState_TP(temp, press);
|
|
double d = w->density();
|
|
double mw = w->molecularWeight();
|
|
double vbar = mw/d;
|
|
// not implemented
|
|
//w.getPartialMolarVolumes(&vbar);
|
|
|
|
printf("%10g %10g %12g %13g %13g\n", temp, press*1.0E-5,
|
|
psat*1.0E-5, d, vbar);
|
|
|
|
}
|
|
|
|
|
|
delete w;
|
|
} catch (CanteraError) {
|
|
|
|
showErrors();
|
|
Cantera::appdelete();
|
|
return -1;
|
|
}
|
|
|
|
|
|
return 0;
|
|
}
|