8 files changed, 300 insertions, 0 deletions
diff --git a/2870/CH9/EX9.10/Ex9_10.sce b/2870/CH9/EX9.10/Ex9_10.sce
new file mode 100755
index 000000000..d6029c924
--- /dev/null
+++ b/2870/CH9/EX9.10/Ex9_10.sce
@@ -0,0 +1,36 @@
+clc;clear;
+//Example 9.10
+
+//from 9.2
+r=8;
+T0=290;
+T1=290;
+T2=652.4;
+T3=1575.1;
+P2=1.7997;
+P3=4.345;
+qin=800;
+qout=381.83;
+wnet=418.17;
+Tsource=1700;
+
+//constants used
+R=0.287;//in kPa-m^3/kg-K
+
+//calculations
+//s1=s2 ; s3=s4
+s03=3.5045;
+s02=2.4975;
+s32=(s03-s02)-R*log(P3/P2);//s32 stands for s3-s2
+xdest23=T0*(s32-qin/Tsource);
+Tsink=T1;
+xdest41=T0*(-s32+qout/Tsink);
+xdestcycle=xdest23+xdest41;
+disp(xdestcycle,'exergy destrustion associated with Otto cycle inkJ/kg');
+// X4 = (u4 - u0 )- T0*(s4 - s0) + P0(v4 - v0)
+// s4 - s0 = s4 - s1 = s32
+// u4 - u0 = u4 - u1 = qout
+// v4 - v0 = v4 - v1 = 0
+//hence x4 is
+X4=qout-T0*s32;
+disp(X4,'exergy destruction of purge stream in kJ/kg')
diff --git a/2870/CH9/EX9.2/Ex9_2.sce b/2870/CH9/EX9.2/Ex9_2.sce
new file mode 100755
index 000000000..6c3332d9c
--- /dev/null
+++ b/2870/CH9/EX9.2/Ex9_2.sce
@@ -0,0 +1,47 @@
+clc;clear;
+//Example 9.2
+
+//given data
+T1=17+273;//in K
+P1=100;
+r=8;//compression ratio i.e v1/v2
+qin=800;
+
+//constants used
+R=0.287;//in kPa-m^3/kg-K
+
+//from Table A-17
+//at T1
+u1=206.91;
+vr1=676.1;
+
+//calculations
+//Process 1-2
+vr2=vr1/r;
+//at this vr2
+T2=652.4;
+u2=475.11;
+P2=P1*(T2/T1)*(r);
+//Process 2-3
+u3=qin+u2;
+//at this u3
+T3=1575.1;
+vr3=6.108;
+P3=P2*(T3/T2)*1;//factor of 1 as v3=v2
+disp(T3,'maximum temperature in the cycle in K');
+disp(P3/1000,'maximum pressure in MPa');//factor of 1000 to convert into MPa
+//Process 3-4
+vr4=r*vr3;
+//at this vr4 
+T4=795.6;
+u4=588.74;
+//Process 4-1
+qout=u4-u1;
+Wnet=qin-qout;
+disp(Wnet,'net work output in kJ/kg');
+nth=Wnet/qin;
+disp(nth,'thermal efficiency');
+v1=R*T1/P1;
+MEP=Wnet/(v1*(1-1/r));
+MEP=round(MEP);
+disp(MEP,'mean effective pressure in kPa')
diff --git a/2870/CH9/EX9.3/Ex9_3.sce b/2870/CH9/EX9.3/Ex9_3.sce
new file mode 100755
index 000000000..dc8d4e9cc
--- /dev/null
+++ b/2870/CH9/EX9.3/Ex9_3.sce
@@ -0,0 +1,56 @@
+clc;clear;
+//Example 9.3
+
+//given data
+V1=117;
+T1=80+460;//in R
+P1=14.7;
+r=18;
+rc=2;
+
+//constants used
+R=0.3704;//in psia ft^3/lbm R
+cp=0.240;//in Btu/lbm R
+cv=0.171;//in Btu/lbm R
+
+//from Table A-2Ea
+k=1.4;
+
+//calculations
+V2=V1/r;
+V3=rc*V2;
+V4=V1;
+//Process 1-2
+T2=T1*(V1/V2)^(k-1);
+P2=P1*(V1/V2)^k;
+T2=round(T2);
+P2=round(P2);
+disp('Process 1-2');
+disp(T2,'temperature in R');
+disp(P2,'pressure in psia');
+//Process 2-3
+P3=P2;
+T3=T2*(V3/V2);
+T3=round(T3);
+P3=round(P3);
+disp('Process 2-3');
+disp(T3,'temperature in R');
+disp(P3,'pressure in psia');
+//Process 3-4
+T4=T3*(V3/V4)^(k-1);
+P4=P3*(V3/V4)^k;
+T4=round(T4);
+P4=round(P4);
+disp('Process 3-4');
+disp(T4,'temperature in R');
+disp(P4,'pressure in psia');
+m=P1*V1/(R*T1)/1728;//factor of 1728 to covert to ft^3 from in^3
+Qin=m*cp*(T3-T2);
+Qout=m*cv*(T4-T1);
+Wnet=Qin-Qout ;
+disp(Wnet,'work output in Btu');
+nth=Wnet/Qin;
+disp(nth,'thermal efficiency');
+MEP=Wnet/(V1-V2)*778.17*12;//factor of 778.17 and 12 to convert to lbf ft and in from Btu and ft respectively
+MEP=round(MEP);
+disp(MEP,'mean effective pressure in psia')
diff --git a/2870/CH9/EX9.5/Ex9_5.sce b/2870/CH9/EX9.5/Ex9_5.sce
new file mode 100755
index 000000000..b4e388335
--- /dev/null
+++ b/2870/CH9/EX9.5/Ex9_5.sce
@@ -0,0 +1,35 @@
+clc;clear;
+//Example 9.5
+
+//given data
+T1=300;
+r=8;
+T3=1300;
+
+//calcualtions
+//Process 1-2
+//at T1
+h1=300.19;
+Pr1=1.386;
+Pr2=r*Pr1;
+//at Pr2
+T2=540;
+h2=544.35;
+disp(T2,'temperature at exit of compressor in K');
+//Process 3-4
+//at T3
+h3=1395.97;
+Pr3=330.9;
+Pr4=Pr3/r;
+//at Pr4
+T4=770;
+h4=789.37;
+disp(T4,'temperature at turbine exit in K');
+Win=h2-h1;
+Wout=h3-h4;
+rbw=Win/Wout;
+disp(rbw,'back work ratio');
+qin=h3-h2;
+Wnet=Wout-Win;
+nth=Wnet/qin;
+disp(nth,'thermal efficeincy')
diff --git a/2870/CH9/EX9.6/Ex9_6.sce b/2870/CH9/EX9.6/Ex9_6.sce
new file mode 100755
index 000000000..a520763e9
--- /dev/null
+++ b/2870/CH9/EX9.6/Ex9_6.sce
@@ -0,0 +1,27 @@
+clc;clear;
+//Example 9.6
+
+//from 9.5
+Wsc=244.16;//compressor
+Wst=606.60;//turbine
+h1=300.19;
+h3=1395.17;
+
+//given data
+nC=0.8;
+nT=0.85;
+
+//calculations
+Win=Wsc/nC;
+Wout=nT*Wst;
+rbw=Win/Wout;
+disp(rbw,'back work ratio is');
+h2a=h1+Win;
+qin=h3-h2a;
+Wnet=Wout-Win;
+nth=Wnet/qin;
+disp(nth,'thermal efficency is');
+h4a=h3-Wout;
+//from A-17 at h4a
+T4a=853;
+disp(T4a,'turbine exit temperature in K is')
diff --git a/2870/CH9/EX9.7/Ex9_7.sce b/2870/CH9/EX9.7/Ex9_7.sce
new file mode 100755
index 000000000..e04eac5cd
--- /dev/null
+++ b/2870/CH9/EX9.7/Ex9_7.sce
@@ -0,0 +1,18 @@
+clc;clear;
+//Example 9.7
+
+//from 9.6
+h2a=605.39;
+h4a=880.36;
+h3=1395.97;
+Wnet=210.41;
+
+//given data
+n=0.80;
+
+//calculations
+// n = (h5 - h2a) / (h4a - h2a)
+h5=(h4a - h2a)*n+h2a;
+qin=h3-h5;
+nth=Wnet/qin;
+disp(nth,'thermal efficiency is')
diff --git a/2870/CH9/EX9.8/Ex9_8.sce b/2870/CH9/EX9.8/Ex9_8.sce
new file mode 100755
index 000000000..a6429a115
--- /dev/null
+++ b/2870/CH9/EX9.8/Ex9_8.sce
@@ -0,0 +1,42 @@
+clc;clear;
+//Example 9.8
+
+//given data
+T1=300;
+T6=1300;
+r=8;//overall compression ratio
+
+//calculations
+//as it is case of intercoolimg
+ri=sqrt(r);//ri stands for P2/P1 = P4/P3 = P6/P7 = P8/P9
+//from A-17 at T1
+h1=300.19;
+Pr1=1.386;
+Pr2=ri*Pr1;
+//from A-17 at Pr2
+T2=403.3;
+h2=403.31;
+//from A-17 at T6
+h6=1395.97;
+Pr6=330.9;
+Pr7=Pr6/ri;
+//from A-17 at Pr7
+T7=1006.4;
+h7=1053.33;
+//at inlets
+T3=T1;h3=h1;T8=T6;h8=h6;
+//et exits
+T4=T2;h4=h2;T9=T7;h9=h7;h5=h7;
+Win=2*(h2-h1);
+Wout=2*(h6-h7);
+Wnet=Wout-Win;
+qin=(h6-h4)+(h8-h7);
+rbw=Win/Wout;
+disp(rbw,'back work ratio');
+nth=Wnet/qin;
+disp(nth,'thermal efficiency is')
+//part - b
+disp('part - b');
+qin=(h6-h5)+(h8-h7);
+nth=Wnet/qin;
+disp(nth,'thermal efficiency is')
diff --git a/2870/CH9/EX9.9/Ex9_9.sce b/2870/CH9/EX9.9/Ex9_9.sce
new file mode 100755
index 000000000..19afda678
--- /dev/null
+++ b/2870/CH9/EX9.9/Ex9_9.sce
@@ -0,0 +1,39 @@
+clc;clear;
+//Example 9.9
+
+//given data
+m=100;
+P1=5;
+T1=-40+460;//in R
+T4=2000+460;//in R
+V1=850;
+rp=10;
+
+//constants used
+cp=0.240;//in Btu/lbm F
+k=1.4;
+
+//calculations
+//Process 1-2
+T2=T1+V1^2/(2*cp)/25037;//factor of 25037 to covert to Btu/lbm
+P2=P1*(T2/T1)^(k/(k-1));
+//Process 2-3 
+P3=rp*P2;
+P4=P3;
+T3=T2*(P3/P2)^((k-1)/k);
+//Win=Wout
+T5=T4-T3+T2;
+P5=P4*(T5/T4)^(k/(k-1));
+T5=round(T5);
+disp(T5,'temperature at turbine exit in R');
+disp(P5,'pressure at turbine exit in psia');
+//Process 5-6
+P6=P1;
+T6=T5*(P6/P5)^((k-1)/k);
+T6=floor(T6);//round off
+V6=sqrt(2*cp*(T5-T6)*25037);//factor of 25037 to covert to (ft/s)^2
+disp(round(V6),'the velocity of nozzle exit in ft/s');
+Wp=m*(V6-V1)*V1/25037;//factor of 25037 to covert to Btu/lbm
+Qin=m*cp*(T4-T3);
+nP=Wp/Qin;
+disp(nP*100,'propulsive efficiency % is')