1 #include "combrx.hpp"
2 using namespace std;
3 
combrx(const string & cas)4 combrx::combrx ( const string & cas ) {
5 
6   // combustion.prop :
7   // -----------------
8   CAS = cas;
9 
10   if (CAS=="64-17-5") {
11     nO2  = 3;
12     nCO2 = 2;
13     nH2O = 3;
14   }
15   else if (CAS=="74-82-8") {
16     nO2  = 2;
17     nCO2 = 1;
18     nH2O = 2;
19   }
20   else if (CAS=="1333-74-0") {
21     nO2  = 0.5;
22     nCO2 = 0;
23     nH2O = 1;
24   }
25   else if (CAS=="100-42-5") {
26     nO2  = 10;
27     nCO2 = 8;
28     nH2O = 4;
29   }
30   else if (CAS=="74-85-1") {
31     nO2  = 3;
32     nCO2 = 2;
33     nH2O = 2;
34   }
35   else if (CAS=="108-88-3") {
36     nO2  = 9;
37     nCO2 = 7;
38     nH2O = 4;
39   }
40   else if (CAS=="100-41-4") {
41     nO2  = 10.5;
42     nCO2 = 8;
43     nH2O = 5;
44   }
45   else if (CAS=="71-43-2") {
46     nO2  = 7.5;
47     nCO2 = 6;
48     nH2O = 3;
49   }
50   else {
51     cout << "ERROR 21" << endl;
52     exit(0);
53   }
54 
55   COMB = new chemical(CAS);
56   O2   = new chemical("7782-44-7");
57   N2   = new chemical("7727-37-9");
58   CO2  = new chemical("124-38-9");
59   H2O  = new chemical("7732-18-5");
60   Hro  = CO2->Ho*nCO2 + nH2O*(H2O->Ho - H2O->dHvap) - COMB->Ho;
61   LFLo = -3420.0/Hro + 0.569e-3*Hro + 0.0538e-6*pow(Hro,2) + 1.8;
62   LFLo = LFLo/100.0;
63   UFLo = 0.0063*Hro + 0.567e-6*pow(Hro, 2) + 23.5;
64   UFLo = UFLo/100.0;
65 
66 }
67