Cantera  2.1.2
RedlichKwong.cpp
Go to the documentation of this file.
1 //! @file RedlichKwong.cpp
2 #include "RedlichKwong.h"
3 
4 namespace tpx
5 {
6 double RedlichKwong::up()
7 {
8  return -Pp()/Rho + hresid() + m_energy_offset;
9 }
10 
11 double RedlichKwong::hresid()
12 {
13  double hh = m_b * (Rho/m_mw);
14  double hresid_mol_RT = z() - 1.0
15  - (1.5*m_a/(m_b*8314.3*T*sqrt(T)))*log(1.0 + hh);
16  return 8314.3*T*hresid_mol_RT/m_mw;
17 }
18 
19 double RedlichKwong::sresid()
20 {
21  double hh = m_b * (Rho/m_mw);
22  double sresid_mol_R = log(z()*(1.0 - hh))
23  - (0.5*m_a/(m_b*8314.3*T*sqrt(T)))*log(1.0 + hh);
24  return 8314.3*sresid_mol_R/m_mw;
25 }
26 
27 double RedlichKwong::sp()
28 {
29  const double Pref = 101325.0;
30  double rgas = 8314.3/m_mw;
31  //double ss = rgas*(log(Pref/(Rho*rgas*T)));
32  double sr = sresid();
33  double p = Pp();
34  return rgas*(log(Pref/p)) + sr + m_entropy_offset;
35 }
36 
37 double RedlichKwong::z()
38 {
39  return Pp()*m_mw/(Rho*8314.3*T);
40 }
41 
42 double RedlichKwong::Pp()
43 {
44  double R = 8314.3;
45  double V = m_mw/Rho;
46  return R*T/(V - m_b) - m_a/(sqrt(T)*V*(V+m_b));
47 }
48 
49 double RedlichKwong::Psat()
50 {
51  double tt = m_tcrit/T;
52  double lpr = -0.8734*tt*tt - 3.4522*tt + 4.2918;
53  return m_pcrit*exp(lpr);
54 }
55 
56 double RedlichKwong::ldens()
57 {
58  double c;
59  int i;
60  double sqt = sqrt(T);
61  double v = m_b, vnew;
62  double pp = Psat();
63  double Rhsave = Rho;
64  for (i = 0; i < 50; i++) {
65  c = m_b*m_b + m_b*GasConstant*T/pp - m_a/(pp*sqt);
66  vnew = (1.0/c)*(v*v*v - GasConstant*T*v*v/pp - m_a*m_b/(pp*sqt));
67  v = vnew;
68  }
69  Rho = Rhsave;
70  return m_mw/vnew;
71 }
72 
73 }