From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 1358/CH7/EX7.10/Example710.sce | 34 ++++++++++++++++++++++++++++++++++
 1 file changed, 34 insertions(+)
 create mode 100755 1358/CH7/EX7.10/Example710.sce

(limited to '1358/CH7/EX7.10/Example710.sce')

diff --git a/1358/CH7/EX7.10/Example710.sce b/1358/CH7/EX7.10/Example710.sce
new file mode 100755
index 000000000..06008189f
--- /dev/null
+++ b/1358/CH7/EX7.10/Example710.sce
@@ -0,0 +1,34 @@
+// 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 7, Example 10")
+disp("For no loss up to throat Ps in bars")
+P01 = 4;//bar
+gam = 1.33;
+Ps = P01*(2/(gam+1))^(gam/(gam-1))
+T01 = 1100;//K
+Ts = 944;//K
+Cpg = 1147;
+U = 300;//m/s
+C = (2*Cpg*(T01-Ts))^0.5//m/s
+R = 0.287;
+rhos = Ps*100/(R*Ts)//kg/m3
+disp("Throat area in m2")
+m=20;//kg/s
+A = m/(rhos*C)
+disp("Angle alpha1, at any radius r and alpha1m at the design radius rm are related by the equation")
+disp("tan(alpha1) = rm*tan(alpha1m)/r")
+disp("Given")
+disp("Tip radius/Root radius = rt/rr = 1.4")
+disp("Therefore mean radius/root radius = 1.2")
+alpha1m = 25
+alpha1r = atan(1.2*tan(alpha1m*%pi/180))*180/%pi
+alpha1t = atan(tan(alpha1r*%pi/180)/1.4)*180/%pi
+disp("Velocity in m/s")
+disp("Cw2 = rm*x*Cw2m/rr = rm*Ca2/(rr*tan(alpha2m))")
+Cw2 = 1.2*250/tan(alpha1m*%pi/180)
+disp("Power developed in kW")
+W = m*U*Cw2/1000
-- 
cgit