initial commit / add all books

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