diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /2870/CH15/EX15.10/Ex15_10.sce | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '2870/CH15/EX15.10/Ex15_10.sce')
-rwxr-xr-x | 2870/CH15/EX15.10/Ex15_10.sce | 57 |
1 files changed, 57 insertions, 0 deletions
diff --git a/2870/CH15/EX15.10/Ex15_10.sce b/2870/CH15/EX15.10/Ex15_10.sce new file mode 100755 index 000000000..8b85924fa --- /dev/null +++ b/2870/CH15/EX15.10/Ex15_10.sce @@ -0,0 +1,57 @@ +clc;clear;
+//Example 15.10
+
+//given values
+T0=298;//in K
+
+//contansts used
+Ru=8.314;//in kJ/kmol K
+
+//calculations
+// CH4 + 3(O2 + 3.76N2) = CO2 + 2H2O + O2 + 11.28N2
+//from std. values of heat of formation and ideal gasses in Appendix
+//methane as m
+hfm=-74850;
+//oxygen as o
+hfo=0;
+h298o=8682;
+//nitrogen as n
+hfn=0;
+h298n=8669;
+//water as w
+hfw=-241820;
+h298w=9904;
+//carbondioxide as c
+hfc=-393520;
+h298c=9364;
+//x refers to hCO2 + 2hH2O + 11.28hN2
+xac=1*(hfm)+1*(h298c-hfc)+2*(h298w-hfw)+11.28*(h298n-hfn);
+//from EES the Tprod is determined by trial and error
+Tprod=1789;
+disp(Tprod,'the temperature of the products in K');
+//entropy calculations by using table A-26
+//Si = Ni*(si - Ruln yiPm
+//reactants
+Sm=1*(186.16-Ru*log(1*1));
+So=3*(205.04-Ru*log(0.21*1));
+Sn=11.28*(191.61-Ru*log(.79*1));
+Sreact=Sm+So+Sn;
+//products
+Nt=1+2+1+11.28;//total moles
+yc=1/Nt;
+yw=2/Nt;
+yo=1/Nt;
+yn=11.28/Nt;
+Sc=1*(302.517-Ru*log(yc*1));
+Sw=2*(258.957-Ru*log(yw*1));
+So=1*(264.471-Ru*log(yo*1));
+Sn=11.28*(247.977-Ru*log(yn*1));
+Sprod=Sc+Sw+So+Sn;
+Sgen=Sprod-Sreact;
+disp(Sgen,'exergy destruction in kJ/kmol - K');
+Xdestroyed=T0*Sgen/1000;//factor of 1000 for converting kJ to MJ
+Xdestroyed=ceil(Xdestroyed);
+disp(Xdestroyed,'in MJ/kmol');
+//This process involves no actual work. Therefore, the reversible work and energy destroyed are identical
+Wrev=Xdestroyed;
+disp(Wrev,'the reversible work in MJ/kmol')
|