8 files changed, 126 insertions, 0 deletions
diff --git a/779/CH5/EX5.1/5_1.sce b/779/CH5/EX5.1/5_1.sce
new file mode 100755
index 000000000..911ecda20
--- /dev/null
+++ b/779/CH5/EX5.1/5_1.sce
@@ -0,0 +1,16 @@
+// Part(a)
+V1 = 0.95;
+P1 = 100e03;
+v1 = 7;
+V2 = 0.19;
+P2 = 700e03;
+v2 = 5;
+w = 0.5;
+u21 = 90e03; // u21 = u2-u1
+Q = -58e03; // As heat is added Q = dQ/dt
+W =  - w*( u21 + (P2*V2-P1*V1) + ((v2^2-v1^2)/2) ) + Q; // W = dW/dt  
+disp("Watt",W,"The rate of work input is")
+// Part (b)
+A = (v2/v1)*(V1/V2); // A = A1/A2
+d = sqrt(A); // d = d1/d2
+disp(d,"The ratio of the inlet pipe diameter and outer pipe diameter is")
+\ No newline at end of file
diff --git a/779/CH5/EX5.2/5_2.sce b/779/CH5/EX5.2/5_2.sce
new file mode 100755
index 000000000..c58be2540
--- /dev/null
+++ b/779/CH5/EX5.2/5_2.sce
@@ -0,0 +1,13 @@
+V1 = 0.37;
+P1 = 600e03;
+v1 = 16;
+V2 = 0.62; 
+P2 = 100e03;
+v2 = 270;
+Z1 = 32;
+Z2 = 0;
+g = 9.81; 
+Q = -9e03; // heat loss Q = dQ/dt
+W = 135e03; // Work done W = dW/dt
+U12 = (P2*V2-P1*V1) + ((v2^2-v1^2)/2) + (Z2-Z1)*g + W - Q; // U12 = U1-U2
+disp("Joule",U12,"The internal energy decreases by") 
+\ No newline at end of file
diff --git a/779/CH5/EX5.3/5_3.sce b/779/CH5/EX5.3/5_3.sce
new file mode 100755
index 000000000..12c24d67d
--- /dev/null
+++ b/779/CH5/EX5.3/5_3.sce
@@ -0,0 +1,13 @@
+P1 = 4e06;
+t1 = 400;
+h1 = 3213e03;
+V1 = 0.073;
+P2 = 3.5e06;
+t2 = 392;
+h2 = 3202e03;
+V2 = 0.084;
+Q = -8.5e03;
+v1 = sqrt((2*(h1-h2+Q))/(1.15^2-1));
+A1 = (%pi/4)*0.2^2;
+w = (A1*v1)/V1;
+disp("Kg/s",w,"The stean flow rate is")
diff --git a/779/CH5/EX5.4/5_4.sce b/779/CH5/EX5.4/5_4.sce
new file mode 100755
index 000000000..c070928d0
--- /dev/null
+++ b/779/CH5/EX5.4/5_4.sce
@@ -0,0 +1,45 @@
+h1 = 313.93;
+h2 = 2676;
+h3 = 419;
+w1 = 4.2;
+w = poly(0,'w') // w = w2
+P = w1*h1 + w*h2 - h3*(4.2+w)
+function [x] = stress(a,b,f)
+  N = 100;
+  eps = 1e-5;
+  if((f(a)*f(b))>0) then
+    error('no root possible f(a)*f(b)>0');
+    abort;
+  end;
+  if(abs(f(a))<eps) then
+    error('solution at a');
+    abort;
+  end
+  if(abs(f(b))<eps) then
+    error('solution at b');
+    abort;
+  end
+  while(N>0)
+    c = (a+b)/2
+    if(abs(f(c))<eps) then
+      x = c ;
+      x;
+      return;
+    end;
+    if((f(a)*f(c))<0 ) then
+      b = c ;
+    else
+      a = c ;
+    end
+    N = N-1;
+  end
+  error('no convergence');
+  abort;
+endfunction
+
+deff('[y]=p(w)',['y = - 441.294 + 2257*w '])
+w = stress(0.1,0.2,p);
+disp("Kg/h",w*3600,"The amount of heat that should be supplied is")
+
+
+
diff --git a/779/CH5/EX5.5/5_5.sce b/779/CH5/EX5.5/5_5.sce
new file mode 100755
index 000000000..dbc9965c4
--- /dev/null
+++ b/779/CH5/EX5.5/5_5.sce
@@ -0,0 +1,10 @@
+t1 = 15; t2 = 800; t3 = 650; t4 = 500;
+v1 = 30; v2 = 30 ; v3 = 60; 
+w = 2;
+cp = 1005;
+Q1_2 = w*cp*(t2-t1);
+disp("KJ/s",Q1_2/1000,"The rate of heat transfer to the air in the heat exchanger is")
+W_T = w*( ((v2^2-v3^2)/2) + cp*(t2-t3));
+disp("KW",W_T/1000,"The power output from the turbine assuming no heat loss")
+v4 = sqrt( (v3^2) + (2*cp*(t3-t4)) );
+disp("m/s",v4,"The velocity at the exit of the nozzle is")
+\ No newline at end of file
diff --git a/779/CH5/EX5.6/5_6.sce b/779/CH5/EX5.6/5_6.sce
new file mode 100755
index 000000000..a9932db8a
--- /dev/null
+++ b/779/CH5/EX5.6/5_6.sce
@@ -0,0 +1,12 @@
+w = 5;
+h1 = 900e03;
+h2 = 400e03;
+v1 = 50; v2 = 150;
+Q = -25*5; // Q = dQ/dt for w = 5kg
+W = w*( (h1-h2) + ((v1^2-v2^2)/2) ) + Q; // W = dW/dt
+disp("kW",W/1000,"The power output of the turbine is")
+R = 285; T1 = 300; P1 = 100e03;
+V = (w*R*T1)/P1; // V = dV/dt
+A1 = V/v1; 
+D1 = sqrt((4*A1)/%pi);
+disp("m",D1,"The diameter of the inlet pipe is")
+\ No newline at end of file
diff --git a/779/CH5/EX5.7/5_7.sce b/779/CH5/EX5.7/5_7.sce
new file mode 100755
index 000000000..a0d99db2d
--- /dev/null
+++ b/779/CH5/EX5.7/5_7.sce
@@ -0,0 +1,9 @@
+ha = 260; // Enthalpy of air
+hg = 912; // Enthalpy of gas
+Va = 270; // Velocity of air
+f = 0.0190; // Fuel to air ratio wf/wa
+Ef = 44500; // Chemical energy of fuel in kJ/kg
+Q = 21; // Heat loss from the engine
+Eg = 0.05*f*Ef/(1+f); // As 5% of chemical energy is not released in reaction
+Vg = sqrt(2000*(((ha+(Va^2*0.001)/2+(f*Ef)-Q)/(1+f))-hg-Eg));
+disp("m/s",Vg,"Velocity of exhaust gas is")
+\ No newline at end of file
diff --git a/779/CH5/EX5.9/5_9.sce b/779/CH5/EX5.9/5_9.sce
new file mode 100755
index 000000000..e30c3bcac
--- /dev/null
+++ b/779/CH5/EX5.9/5_9.sce
@@ -0,0 +1,8 @@
+u0 = 0.718*273*1e03;
+t = poly(0,'t');
+u = u0+718*t; // in SI unit
+hp = u + 285*(t+273); // ""
+h = horner(hp,150); // h = hp(150)
+W = 100; // W = dW/dt
+m = W/h;
+disp("kg/h",m*3600,"The rate at which air flows out of the tank")
+\ No newline at end of file