initial commit / add all books

1 files changed, 42 insertions, 0 deletions

diff --git a/3407/CH6/EX6.1/Ex6_1.sce b/3407/CH6/EX6.1/Ex6_1.sce
new file mode 100644
index 000000000..dd23b1d81
--- /dev/null
+++ b/3407/CH6/EX6.1/Ex6_1.sce
@@ -0,0 +1,42 @@
+clear;

+clc;

+funcprot(0);

+

+//given data

+dt = 1.0;//tip diameter in m

+dh = 0.9;//hub diameter in m

+alpha1 = 30;//in deg

+beta1 = 60;//in deg

+alpha2 = 60;//in deg

+beta2 = 30;//in deg

+N = 6000;//rotational speed in rev/min

+rhog = 1.5;//gas density in kg/m^3

+Rt = 0.5;//degree of reaction at the tip

+

+//Calculations

+omega = 2*%pi*N/60;

+Ut = omega*0.5*dt;

+Uh = omega*0.5*dh;

+cx = Ut/(tan(alpha1*%pi/180) + tan(beta1*%pi/180));

+mdot = %pi*((0.5*dt)^2 - (0.5*dh)^2)*rhog*cx;

+Wcdot = mdot*Ut*cx*(tan(alpha2*%pi/180)- tan(alpha1*%pi/180));

+ctheta1t = cx*tan(alpha1*%pi/180);

+ctheta1h = ctheta1t*(dt/dh);

+ctheta2t = cx*tan(alpha2*%pi/180);

+ctheta2h = ctheta2t*(dt/dh);

+alpha1_ = (180/%pi)*atan(ctheta1h/cx);

+beta1_ = (180/%pi)*atan((Uh/cx) - tan(alpha1_*%pi/180));

+alpha2_ = (180/%pi)*atan(ctheta2h/cx);

+beta2_ = (180/%pi)*atan((Uh/cx) - tan(alpha2_*%pi/180));

+k = Rt*(0.5*dt)^2;

+Rh = 1 - (k/(0.5*dh)^2);

+

+//Results

+printf('(i)The axial velocity, cx = %d m/s',cx);

+printf('\n (ii)The mass flow rate = %.1f kg/s',mdot);

+printf('\n (iii)The power absorbed by the stage = %.1f MW',Wcdot/(10^6));

+printf('\n (iv)The flow angles at the hub are:\n alpha1 = %.2f deg,\n beta1 = %.2f deg,\n alpha2 = %.1f deg, and\n beta2 = %.2f deg.',alpha1_,beta1_,alpha2_,beta2_);

+printf('\n (v)The reaction ratio of the stage at the hub, R = %.3f.',Rh);

+

+

+//there are small errors in the answers given in textbook