1 files changed, 47 insertions, 0 deletions
diff --git a/1358/CH5/EX5.7/Example57.sce b/1358/CH5/EX5.7/Example57.sce
new file mode 100755
index 000000000..8e62f3e05
--- /dev/null
+++ b/1358/CH5/EX5.7/Example57.sce
@@ -0,0 +1,47 @@
+// Display mode
+mode(0);
+// Display warning for floating point exception
+ieee(1);
+clear;
+clc;
+disp("Turbomachinery Design and Theory,Rama S. R. Gorla and Aijaz A. Khan, Chapter 5, Example 7")
+disp("The following equation provides the relationship between the temperature rise and the desired angles:")
+disp("T02 - T01 = lambda*U*Ca*(tan(beta1) - tan(beta2))/Cp")
+disp("T02-T01 = Rise")
+Rise = 24;
+lambda = 0.93;
+U = 205;
+Ca = 155.5;
+Cp = 1005;
+disp("Dif = tan(beta1) - tan(beta2)")
+Dif = Rise*Cp/(U*lambda*Ca)
+disp("Using the degree of reaction equation:")
+disp("DOF = Ca*(tan(beta1) + tan(beta2))/(2*U)")
+disp("tan(beta1) + tan(beta2) = Add")
+DOF = 0.5;
+Add = DOF*2*U/Ca
+beta1 = atan((Add+Dif)/2)*180/%pi
+alpha2 = beta1
+beta2 = atan(Add - tan(beta1*%pi/180))*180/%pi
+alpha1 = beta2
+disp("The mean radius, rm, is given by: in m")
+N = 152;//rpm
+rm = U/(2*%pi*N)
+disp("The blade height, h in m, is given by: m = rho*ACa, where A is the annular area of the flow.")
+C1 = Ca/cos(alpha1*%pi/180)
+T01 = 290;
+disp("Static temperature in kelvin")
+T1 = T01- C1^2 /(2*Cp)
+disp("Using the isentropic P–T relation:")
+disp("Static pressure: P1 in bars")
+P01 = 1;
+P1 = P01*(T1/T01)^3.5
+R = 287/1000;
+disp("Density rho1 in kg/m3")
+rho1 = P1/(R*T1) * 100
+disp("From the continuity equation:")
+m = 22;
+A = m/(rho1*Ca)
+disp("and the blade height:")
+rm = 0.215;
+h =A/(2*%pi*rm)