diff options
Diffstat (limited to '632/CH12/EX12.19/example12_19.sce')
-rwxr-xr-x | 632/CH12/EX12.19/example12_19.sce | 21 |
1 files changed, 21 insertions, 0 deletions
diff --git a/632/CH12/EX12.19/example12_19.sce b/632/CH12/EX12.19/example12_19.sce new file mode 100755 index 000000000..7fb2f5f97 --- /dev/null +++ b/632/CH12/EX12.19/example12_19.sce @@ -0,0 +1,21 @@ +//clc()
+Nco = 1;//mol CO
+// CO + 1/2 O2 = CO2
+O2r = 1;//mol
+N2r = 3.76;//mol
+COr = 1;//mol
+O2p = 0.5;//mol
+N2p = 3.76;//mol
+CO2p = 1;//mol
+Hco = 29.23;//J/molK
+Ho2 = 34.83;//J/molK
+Hn2 = 33.03;//J/molK
+Hco2 = 53.59;//J/molK
+Hcomb1 = -282.99;//kJ/mol
+T1 = 298;//K
+T2 = 373;//K
+H1 = (O2r * Ho2 + N2r * Hn2 + COr * Hco) * (T1 - T2);
+//For product at temp T, H2 = (O2p * Ho2 + N2p * Hn2 + CO2p * Hco2) * (T - T1)
+//For adiabatic condition, -(H1 + Hcomb1) = H2
+T = -(H1 + Hcomb1 * 1000) / (O2p * Ho2 + N2p * Hn2 + CO2p * Hco2) + T1;
+disp("K",T,"Theoretical flame temperature = ")
\ No newline at end of file |