9 files changed, 382 insertions, 0 deletions
diff --git a/2870/CH10/EX10.1/Ex10_1.sce b/2870/CH10/EX10.1/Ex10_1.sce
new file mode 100755
index 000000000..e5af001b8
--- /dev/null
+++ b/2870/CH10/EX10.1/Ex10_1.sce
@@ -0,0 +1,33 @@
+clc;clear;
+//Example 10.1
+
+//given data
+P1=75;
+P2=3000;//in kPa
+P3=P2;
+T3=350;
+P4=P1;
+
+//from steam tables
+//at state 1
+v1=0.001037;
+h1=384.44;
+//at state 3
+h3=3116.1;
+s3=6.7450;
+//at state 4
+s4=s3;
+sf=1.2132;
+sfg=6.2426;
+hf=384.44;
+hfg=2278;
+
+//calculations
+win=v1*(P2-P1);
+h2=h1+win;
+x4=(s4-sf)/sfg;
+h4=hf+x4*hfg;
+qin=h3-h2;
+qout=h4-h1;
+nth=1-(qout/qin);
+disp(nth*100,'thermal efficency % is')
diff --git a/2870/CH10/EX10.2/Ex10_2.sce b/2870/CH10/EX10.2/Ex10_2.sce
new file mode 100755
index 000000000..6b977d08b
--- /dev/null
+++ b/2870/CH10/EX10.2/Ex10_2.sce
@@ -0,0 +1,33 @@
+clc;clear;
+//Example 10.2
+
+//given data
+P1=9;
+T1=38;
+P2=16000;
+P3=15.9;
+T3=35;
+P4=15.2;
+T4=625;
+P5=15;
+T5=600;
+nT=0.87;
+nP=0.85;
+m=15;
+
+//from steam tables
+v1=0.001009;
+h5=3583.1;
+h6s=2115.3;
+h4=3647.6;
+h3=160.1;
+
+//calculations
+Win=v1*(P2-P1)/nP;
+Wout=nT*(h5-h6s);
+qin=h4-h3;
+Wnet=Wout-Win;
+nth=Wnet/qin;
+disp(nth,'thermal efficiency is');
+Wnet=m*Wnet;
+disp(Wnet/1000,'power output in MW')
diff --git a/2870/CH10/EX10.3/Ex10_3.sce b/2870/CH10/EX10.3/Ex10_3.sce
new file mode 100755
index 000000000..9c44e9688
--- /dev/null
+++ b/2870/CH10/EX10.3/Ex10_3.sce
@@ -0,0 +1,56 @@
+clc;clear;
+//Example 10.3
+
+//given data
+P1=10;
+P2=3000;
+P3=3000;
+T3=350;
+P4=10;
+
+//from steam tables
+//at state 1
+h1=191.81;
+v1=0.00101;
+//at state 2
+//s2=s1
+//at state 3
+h3=3116.1;
+s3=6.7450;
+//at state 4
+s4=s3;
+sf=0.6492;
+sfg=7.4996;
+hf=191.81;
+hfg=2392.1;
+
+//calculations
+//part - a
+win=v1*(P2-P1);
+h2=h1+win;
+x4=(s4-sf)/sfg;
+h4=hf+x4*hfg;
+qin=h3-h2;
+qout=h4-h1;
+nth=1-(qout/qin);
+disp(nth,'the thermal efficiency of this power plant');
+//part - b
+//States 1 and 2 remain the same in this case, and the enthalpies at state 3 (3 MPa and 600°C) and state 4 (10 kPa and s4=s3) are determined to be
+h3=3682.8;
+h4=2380.3;
+x4=0.915;
+qin=h3-h2;
+qout=h4-h1;
+nth=1-(qout/qin);
+disp(nth,'the thermal efficiency if steam is superheated to 600° instead of 350°C');
+//part - c
+//State 1 remains the same in this case, but the other states change. The enthalpies at state 2 (15 MPa and s2  s1), state 3 (15 MPa and 600°C),and state 4 (10 kPa and s4  s3) are determined in a similar manner to be
+h2=206.95;
+h3=3583.1;
+h4=2115.3;
+x4=0.804;
+qin=h3-h2;
+qout=h4-h1;
+nth=1-(qout/qin);
+disp(nth,'the thermal efficiency if the boiler pressure is raised to 15 MPa while the turbine inlet temperature is maintained at 600°C');
+
diff --git a/2870/CH10/EX10.4/Ex10_4.sce b/2870/CH10/EX10.4/Ex10_4.sce
new file mode 100755
index 000000000..443c95bfe
--- /dev/null
+++ b/2870/CH10/EX10.4/Ex10_4.sce
@@ -0,0 +1,44 @@
+clc;clear;
+//Example 10.4
+
+//given data
+P1=10;
+P2=15000;
+P3=15000;
+T3=600;
+P4=4000;
+T5=600;
+P6=10;
+x6=0.896;
+
+//from steam table
+//at state 1
+h1=191.81;
+v1=0.00101;
+//at state 3
+h3=3593.1;
+s3=6.6796;
+//at state 4
+h4=3155;
+T4=375.5;
+//at state 6
+sf=0.6492;
+sfg=7.4996;
+hf=191.81;
+hfg=2392.1;
+
+//calculations
+s6=sf+x6*sfg;
+h6=hf+x6*hfg;
+//s5 = s6
+//from tables
+P5=4000;//in kPa
+h5=3674.9;
+disp(P5/1000,'the pressure at which the steam should be reheated in MPa');
+//s2 = s1
+win=v1*(P2-P1);
+h2=h1+win;
+qin=(h3-h2)+(h5-h4);
+qout=h6-h1;
+nth=1-(qout/qin);
+disp(nth,'thermal efficiency is')
diff --git a/2870/CH10/EX10.5/Ex10_5.sce b/2870/CH10/EX10.5/Ex10_5.sce
new file mode 100755
index 000000000..4ac9e0da0
--- /dev/null
+++ b/2870/CH10/EX10.5/Ex10_5.sce
@@ -0,0 +1,53 @@
+clc;clear;
+//Example 10.5
+
+//given data
+P1=10;
+P2=1200;
+P3=1200;
+P4=15000;
+P5=15000;
+T5=600;
+P6=1200;
+P7=10;
+
+//from steam table
+//at state 1
+h1=191.81;
+v1=0.00101;
+//at state 3
+h3=798.33;
+v3=0.001138;
+//at state 4
+h4=3155;
+T4=375.5;
+//at state 5
+h5=3583.1;
+s5=6.6796;
+//at state 6
+h6=2860.2;
+T6=218.4;
+//at state 7
+P7=10;
+sf=0.6492;
+sfg=7.4996;
+hf=191.81;
+hfg=2392.1;
+
+//calculations
+//s2 = s1
+win=v1*(P2-P1);
+h2=h1+win;
+//s4 = s3
+win=v3*(P4-P3);
+h4=h3+win;
+s7=s5;
+x7=(s7-sf)/sfg;
+h7=hf+(x7*hfg);
+//y is the fraction of steam extracted from the turbine
+y=(h3-h2)/(h6-h2);
+qin=h5-h4;
+qout=(1-y)*(h7-h1);
+nth=1-(qout/qin);
+disp(y,'fraction of steam extracted');
+disp(nth,'thermal efficiency is')
diff --git a/2870/CH10/EX10.6/Ex10_6.sce b/2870/CH10/EX10.6/Ex10_6.sce
new file mode 100755
index 000000000..78640d962
--- /dev/null
+++ b/2870/CH10/EX10.6/Ex10_6.sce
@@ -0,0 +1,46 @@
+clc;clear;
+//Example 10.6
+
+//given data
+P1=10;
+P2=500;
+P3=500;
+P4=15000;
+P5=P4;
+P6=4000;
+P7=P5;
+P8=P7;
+P9=P7;
+P10=P6;
+P11=P10;
+P12=P3;
+P13=10;
+
+//enthalpies at the various states and the pump work per unit mass of fluid flowing through them are
+h1=191.81;
+h2=192.30;
+h3=640.09;
+h4=643.92;
+h5=1087.4;
+h6=h5;
+h7=1101.2;
+h8=1089.8;
+h9=3583.1;
+h10=3155;
+h11=3679.9;
+h12=3014.8;
+h13=2335.7;
+wIin=0.49;
+wIIin=3.83;
+wIIIin=13.77;
+
+//calculations
+y=(h5-h4)/((h10-h6)+(h5-h4));
+z=(1-y)*(h3-h2)/(h12-h2);
+h8=(1-y)*h5+(y*h7);
+qin=(h9-h8)+(1-y)*(h11-h10);
+qout=(1-y-z)*(h13-h1);
+nth=1-(qout/qin);
+disp(y,'fraction of steam extracted from closed feedwater');
+disp(z,'fraction of steam extracted from open feedwater');
+disp(nth,'thermal efficiency is')
diff --git a/2870/CH10/EX10.7/Ex10_7.sce b/2870/CH10/EX10.7/Ex10_7.sce
new file mode 100755
index 000000000..1dc8c2012
--- /dev/null
+++ b/2870/CH10/EX10.7/Ex10_7.sce
@@ -0,0 +1,34 @@
+clc;clear;
+//Example 10.7
+
+//given data
+T0=290;
+Tsource=1600;
+Tsink=T0;
+//from Ex 10.1
+qin=2728.6;
+qout=2018.6;
+h4=2403;
+
+//from steam tables
+s1=1.2132;
+s3=6.7450;
+
+//calculations
+s2=s1;s4=s3;//isentropic processes
+xdest12=0;
+xdest34=0;
+xdest23=T0*(s3-s2-(qin/Tsource));
+xdest41=T0*(s1-s4+(qout/Tsink));
+disp(xdest12,'exergy destruction in 1-2 in kJ/kg');
+disp(round(xdest23),'exergy destruction in 2-3 in kJ/kg');
+disp(xdest34,'exergy destruction in 3-4 in kJ/kg');
+disp(round(xdest41),'exergy destruction in 4-1 in kJ/kg');
+xdestcy=xdest12+xdest23+xdest34+xdest41;
+disp(round(xdestcy),'exergy destruction in cycle in kJ/kg');
+//from steam tables
+//at 290 K and 100 kPa
+h0=71.355;
+s0=0.2533;
+X4=(h4-h0)-T0*(s4-s0);
+disp(round(X4),'exergy of the leaving steam in kJ/kg')
diff --git a/2870/CH10/EX10.8/Ex10_8.sce b/2870/CH10/EX10.8/Ex10_8.sce
new file mode 100755
index 000000000..852646366
--- /dev/null
+++ b/2870/CH10/EX10.8/Ex10_8.sce
@@ -0,0 +1,49 @@
+clc;clear;
+//Example 10.8
+
+//given data
+m1=15;
+P1=7000;
+P2=P1;
+P3=P1;
+P4=500;
+P5=P4;
+P6=5;
+P7=500;
+P8=5;
+P9=7000;
+P10=7000;
+
+//from steam tables
+v7=0.001005;
+v8=0.001093;
+h1=3411.4;
+h2=h1;
+h3=h1;
+h4=h1;
+h5=2739.3;
+h6=2073.0;
+h7=640.09;
+h8=137.75;
+h11=144.78;
+
+//calculations
+wIin=v8*(P9-P8);
+wIIin=v7*(P10-P7);
+h9=h8+wIin;
+h10=h7+wIIin;
+Qmax=m1*(h1-h7);
+disp(Qmax,'the maximum rate in kW');
+Wtout=m1*(h3-h6);//turbine
+Wpin=m1*wIin;//pump
+Wnet=Wtout-Wpin;
+disp(round(Wnet/1000),'the power produced in MW');
+Qp=0;
+Qin=m1*(h1-h11);
+Eu=(Wnet+Qp)/Qin;
+disp(Eu,'the utilization factor');
+m4=0.1*m1;
+m5=0.7*m1;
+m7=m4+m5;
+Qout=m4*h4+m5*h5-m7*h7;
+disp(Qout/1000,'the rate of process heat supply in MW')
diff --git a/2870/CH10/EX10.9/Ex10_9.sce b/2870/CH10/EX10.9/Ex10_9.sce
new file mode 100755
index 000000000..c4a106351
--- /dev/null
+++ b/2870/CH10/EX10.9/Ex10_9.sce
@@ -0,0 +1,34 @@
+clc;clear;
+//Example 10.9
+
+//given data
+P1=5;
+P2=7000;
+P3=P2;
+T3=500;
+P4=P1;
+
+//gas cycle from Ex9-6
+//d stands for '
+h4d=880.36;
+T4d=853;
+qin=790.58;
+wnetg=210.41;
+nth=0.266
+h5d=451.80;
+//steam cycle
+h2=144.78;
+T2=33;
+h3=3411.4;
+T3=500;
+wnets=1331.4;
+nth=0.408;
+
+//calculations
+//Ein = Eout
+//y is the ratio of ms/mg
+y=(h4d-h5d)/(h3-h2);
+disp(y,'the ratio of the mass flow rates of the steam and the combustion gases');
+wnet=wnetg+y*wnets
+nth=wnet/qin;
+disp(nth,'the thermal efficiency of the combined cycle')