13 files changed, 315 insertions, 0 deletions
diff --git a/2870/CH5/EX5.1/Ex5_1.sce b/2870/CH5/EX5.1/Ex5_1.sce
new file mode 100755
index 000000000..8770b06e9
--- /dev/null
+++ b/2870/CH5/EX5.1/Ex5_1.sce
@@ -0,0 +1,17 @@
+clc;clear;
+//Example 5.1
+
+//given data
+V=10;
+t=50;
+p=1;//in kg/L
+re=0.8/2/100;//in m
+
+//calculations
+Vd=V/t*3.7854;//factor 0f 3.7854 for gal to L
+disp(Vd,'volumne flow rate through hose in L/s');
+m=p*Vd;
+disp(m,'mass flow rate through hose in kg/s');
+Ae=%pi*re^2;
+Ve=Vd/Ae/1000;//factor of 1000 for L to m^3
+disp(Ve,'average velocity at the nozzle in m/s');
diff --git a/2870/CH5/EX5.10/Ex5_10.sce b/2870/CH5/EX5.10/Ex5_10.sce
new file mode 100755
index 000000000..aed99a6a0
--- /dev/null
+++ b/2870/CH5/EX5.10/Ex5_10.sce
@@ -0,0 +1,27 @@
+clc;clear;
+//Example 5.10
+
+//given data
+T1=15;
+P1=300;
+T2=25;
+T3=70;
+P3=1000;//in kPa
+T4=35;
+mr=6;
+
+//from Table A-4, A-13 and A-11
+h1=62.982;
+h2=104.83;
+h3=303.85;
+h4=100.87;
+
+//calculations
+//mass balance m1=m2=mw and m3=m4=mr
+//energy balance m1*h1 + m3*h3 =  m2*h2 + m4*h4
+//combining them mw*(h1-h2) = mr*(h4-h3)
+mw= mr*(h4-h3)/(h1-h2);
+disp(mw,'mass flow rate of cooling water in kg/min');
+Qin=mw*(h2-h1);
+Qin=round(Qin);
+disp(Qin,'heat transfer rate in kJ/min')
diff --git a/2870/CH5/EX5.11/Ex5_11.sce b/2870/CH5/EX5.11/Ex5_11.sce
new file mode 100755
index 000000000..0961b288a
--- /dev/null
+++ b/2870/CH5/EX5.11/Ex5_11.sce
@@ -0,0 +1,20 @@
+clc;clear;
+//Example 5.11
+
+//giaven data
+T1=17+273;//in K
+P1=100;
+V1=150;
+Win=15;
+Qout=200/1000;//in kJ/s
+
+//constants used
+R=0.287;//in kPa-m^3/kg-K
+cp=1.005;//in kJ/kg C
+
+//calculations
+v1=R*T1/P1;
+m=V1/v1/60;//factor of 6 to convert to s
+// Win - Qout = m*cp*(T2-T1)
+T2= T1 + (Win - Qout)/(m*cp);
+disp((T2-273),'exit temperature in C')
diff --git a/2870/CH5/EX5.12/Ex5_12.sce b/2870/CH5/EX5.12/Ex5_12.sce
new file mode 100755
index 000000000..cd01dc255
--- /dev/null
+++ b/2870/CH5/EX5.12/Ex5_12.sce
@@ -0,0 +1,20 @@
+clc;clear;
+//Example 5.12
+
+//given data
+Pi=1;
+Ti=300;
+P2=1;
+
+//from Table A-6
+hi=3051.6;
+
+//calculations
+//mass balance mi=m2
+//energy balance mi*hi= m2*u2
+//combining them we get,
+u2=hi;
+//from Table A-6
+//we know P2 and u2, so
+T2=456.1;
+disp(T2,'final temperature in tank in C')
diff --git a/2870/CH5/EX5.13/Ex5_13.sce b/2870/CH5/EX5.13/Ex5_13.sce
new file mode 100755
index 000000000..fefeb7f2c
--- /dev/null
+++ b/2870/CH5/EX5.13/Ex5_13.sce
@@ -0,0 +1,38 @@
+clc;clear;
+//Example 5.13
+
+//given data
+V=6/1000;//in m^3
+Pgage=75;
+Patm=100;
+m1=1;
+Qind=0.5;//d stands for .
+t=30*60;//in s
+
+//calculation
+Pabs=Pgage+Patm;
+//from Table A-5, ths saturation temp 
+T=116.04;
+disp(T,'the temperature at which cooking takes plac in C');
+//mass balance me=(m1-m2)cv
+//energy balance Qin - mehe = (m2u2 - m1u1)cv
+Qin=Qind*t;
+//from Table A-5
+he=2700.2;
+vf=0.001;
+vg=1.004;
+uf=486.82;
+ufg=2037.7;
+v1=V/m1;
+x1=(v1-vf)/(vg-vf);
+u1=uf+x1*ufg;
+U=m1*u1;
+//Qin = (m1 - V/v2)*he + (V/v2*u2 - m1*u1)
+//v2=vf + x2*(vg-vf)
+//u2=uf +  x2*ufg
+//combining these equations we get
+//solved using EES
+x2=0.009;
+v2=vf + x2*(vg-vf);
+m2=V/v2;
+disp(m2,'amount of water left in kg')
diff --git a/2870/CH5/EX5.2/Ex5_2.sce b/2870/CH5/EX5.2/Ex5_2.sce
new file mode 100755
index 000000000..38a484838
--- /dev/null
+++ b/2870/CH5/EX5.2/Ex5_2.sce
@@ -0,0 +1,20 @@
+clc;clear;
+//Example 5.2
+
+//given data
+Dtank=3*12;//in inches
+Djet=0.5;
+h0=2;
+h1=4;
+
+//constants used 
+g=32.2;//in ft/s^2
+
+//calculations
+//min - mout = dmCV/dt
+//mout = p*(2*g*h*Ajet)^2
+//mCV = p*Atank*h
+//from these we get dt = Dtank^2/Djet^2 * (dh/(2*g*h)^2)
+t=integrate('Dtank^2/Djet^2*(1/sqrt(2*g*h))','h',h0,h1);
+t=(t/60);//in min
+disp(t,'time taken to drop to 2ft in min')
diff --git a/2870/CH5/EX5.3/Ex5_3.sce b/2870/CH5/EX5.3/Ex5_3.sce
new file mode 100755
index 000000000..c8fefe9bf
--- /dev/null
+++ b/2870/CH5/EX5.3/Ex5_3.sce
@@ -0,0 +1,28 @@
+clc;clear;
+//Example 5.3
+
+//given data
+P=150;
+Vliquid=0.6/1000;//im m^3
+t=40*60;//in sec
+Ac=8*10^-6;
+
+//from Table A-5
+//from P = 150 kPa
+h=2693.1;
+ug=2519.2;
+vf=0.001053;
+vg=1.1594;
+
+//calculations
+m=Vliquid/vf;
+md=m/t;
+disp(md,'mass flow rate in kg/s');
+V=md*vg/(Ac);
+disp(V,'exit velocity in m/s');
+Eflow=h-ug;
+Et=h;
+disp(Eflow,'flow energy in kJ/kg');
+disp(Et,'total energy in kJ/kg');
+Emass=md*Et;
+disp(Emass,'rate at which energy leaves the cooker in kW')
diff --git a/2870/CH5/EX5.4/Ex5_4.sce b/2870/CH5/EX5.4/Ex5_4.sce
new file mode 100755
index 000000000..8c43ea505
--- /dev/null
+++ b/2870/CH5/EX5.4/Ex5_4.sce
@@ -0,0 +1,24 @@
+clc;clear;
+//Example 5.4
+
+//given data
+T1=283;//in K
+P1=80;
+V1=200;
+A1=0.4;
+
+//constants used
+R=0.287;//in kPa-m^3/kg-K
+
+//calulations
+v1=R*T1/P1;
+m=V1*A1/v1;
+disp(m,'mass flow rate of air in kg/s');
+// Ein - Eout = dEsystem / dt
+//from Table A-17
+h1=283.14;
+V2=0;
+h2=h1-(V2^2 - V1^2)/2/1000;//factor of 1000 to convert to kJ/kg
+//from Table A-17 at this value of h2
+T2=303;
+disp(T2,'the temperature in K is');
diff --git a/2870/CH5/EX5.5/Ex5_5.sce b/2870/CH5/EX5.5/Ex5_5.sce
new file mode 100755
index 000000000..eba140162
--- /dev/null
+++ b/2870/CH5/EX5.5/Ex5_5.sce
@@ -0,0 +1,24 @@
+clc;clear;
+//Example 5.5
+
+//given data
+P1=250;
+T1=700;
+A1=0.2;
+qout=1.2;
+m=10;
+P2=200;
+V2=900;
+
+//from Table A-6E
+v1=2.6883;
+h1=1371.4;
+
+//calculations
+V1=m*v1/A1;
+disp(V1,'the inlet velocity in ft/s');
+// Ein - Eout = dEsystem / dt
+h2=h1-qout-(V2^2 - V1^2)/2/25037;//factor of 25037 to convert to Btu/lbm
+//at this value h2, from Tablw A-6E
+T2=662;
+disp(T2,'exit temperature in F')
diff --git a/2870/CH5/EX5.6/Ex5_6.sce b/2870/CH5/EX5.6/Ex5_6.sce
new file mode 100755
index 000000000..1c9435de4
--- /dev/null
+++ b/2870/CH5/EX5.6/Ex5_6.sce
@@ -0,0 +1,19 @@
+clc;clear;
+//Example 5.6
+
+//given data
+T1=280;
+P1=100;
+m=0.02;
+qout=16;
+P2=600;
+T2=400;
+
+//from Table A-17
+h1=280.13;
+h2=400.98;
+
+//calculations
+// Ein - Eout = dEsystem / dt
+Win=m*qout+m*(h2-h1);
+disp(Win,'the input power of compressor in kW')
diff --git a/2870/CH5/EX5.7/Ex5_7.sce b/2870/CH5/EX5.7/Ex5_7.sce
new file mode 100755
index 000000000..8d097c822
--- /dev/null
+++ b/2870/CH5/EX5.7/Ex5_7.sce
@@ -0,0 +1,32 @@
+clc;clear;
+//Example 5.7
+
+//given data
+P1=2;
+T1=400;
+V1=50;
+z1=10;
+P2=15;
+x2=0.9;
+V2=180;
+z2=6;
+Wout=5*1000;//in kJ
+
+//from Table A-6
+h1=3248.4;
+//similarly for P2
+hf=225.94;
+hfg=2372.3;
+
+//constants used 
+g=9.8;//in m/s^2
+
+//calcualtions
+h2=hf+x2*hfg;
+disp((h2-h1),'difference in enthalpies in kJ/kg');
+disp((V2^2-V1^2)/2/1000,'difference in kinetic energy in kJ/kg');//factor of 1000 to convert to kJ/kg
+disp(g*(z2-z1)/1000,'difference in potential energy in kJ/kg');//factor of 1000 to convert to kJ/kg
+wout=-((h2-h1)+(V2^2-V1^2)/2/1000+g*(z2-z1)/1000);//factor of 1000 to convert to kJ/kg
+disp(wout,'work done per unit of mass in kJ/kg');
+m=Wout/wout;
+disp(m,'mass flow rate in kg/s')
diff --git a/2870/CH5/EX5.8/Ex5_8.sce b/2870/CH5/EX5.8/Ex5_8.sce
new file mode 100755
index 000000000..f4b2047a3
--- /dev/null
+++ b/2870/CH5/EX5.8/Ex5_8.sce
@@ -0,0 +1,22 @@
+clc;clear;
+//Example 5.8
+
+//given data
+P1=0.8;
+P2=0.12;
+
+//from Table A-12
+//sat. liq at P1
+T1=31.31;
+h1=95.47;
+//since process is insentropic and at  P2
+h2=h1;
+hf=22.49;
+hg=236.97;
+T2=-22.32;
+
+//calculations
+x2=(h2-hf)/(hg-hf);
+disp(x2,'the final state is');
+dT=T2-T1;
+disp(dT,'temperature drop in C')
diff --git a/2870/CH5/EX5.9/Ex5_9.sce b/2870/CH5/EX5.9/Ex5_9.sce
new file mode 100755
index 000000000..43aacd0dd
--- /dev/null
+++ b/2870/CH5/EX5.9/Ex5_9.sce
@@ -0,0 +1,24 @@
+clc;clear;
+//Example 5.9
+
+//given data
+T1=140;
+T2=50;
+T3=110;
+P=20;
+
+//for a compressed liq at given temp
+h1=107.99;
+h2=18.07;
+h3=78.02;
+
+//calculations
+//Mass balance min =  mout So, m1+m2 = m3
+//Energy balance Ein = Eout So, m1*h1 + m2*h2 = m3*h3
+//combining realations
+//m1*h1 + m2*h2 = (m1+m2)*h3
+//dividing by m2 and y=m1/m2
+//we get, yh1 + h2 = (y+1)*h3
+y=(h3-h2)/(h1-h3);
+y=round(y);
+disp(y,'the ratio of mass flow rates')