1 #include "combrx.hpp" 2 using namespace std; 3 combrx(const string & cas)4combrx::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