initial commit / add all books

6 files changed, 179 insertions, 0 deletions

diff --git a/2175/CH11/EX11.1/11_1.sce b/2175/CH11/EX11.1/11_1.sce
new file mode 100755
index 000000000..6006e92d6
--- /dev/null
+++ b/2175/CH11/EX11.1/11_1.sce
@@ -0,0 +1,40 @@
+clc;

+Cai=900;

+Cb=300;

+alpha=20*%pi/180;

+Cri=(Cai^2+Cb^2-2*Cb*Cai*cos(alpha))^0.5;

+b=asin(Cai*sin(alpha)/Cri);

+Beta=180*b/%pi

+disp("the blade inlet angle is:");

+disp("degree",Beta)

+

+//part II

+k=0.7;

+Cre=k*Cri

+AD=Cri*cos(b);

+AE=Cre*cos(b);

+

+Cw=AD+AE;

+disp("driving force on wheel is:");

+m=1;

+Df=m*Cw

+disp("N per kg/s",Df);

+

+

+//part III

+Cfi=Cri*sin(b);

+Cfe=Cre*sin(b);

+Cf=Cfi-Cfe;

+At=m*Cf;

+disp("axial thrust is:");

+disp("N per kg/s",At)

+

+//part IV

+Dp=Cb*Cw;

+disp("diagram power per unit mass flow rate:");

+disp("kW",Dp/1000);

+

+//part V

+De=Cb*Cw/(Cai^2);

+disp("Diagram efficiency is");

+disp("%",De*100);

diff --git a/2175/CH11/EX11.10/11_10.sce b/2175/CH11/EX11.10/11_10.sce
new file mode 100755
index 000000000..7fbab3d7a
--- /dev/null
+++ b/2175/CH11/EX11.10/11_10.sce
@@ -0,0 +1,17 @@
+clc;

+T=20+273;

+y=1.4;

+Ti=T*4^([y-1]/y)

+ir=Ti-T;

+actual_r=ir/0.8;

+cp=1.005;

+P=cp*actual_r;

+Cai=150;

+Cbi=15000*%pi*250/(60*10^3);

+Cwi=Cai*sin(25*%pi/180);

+Cbe=15000*%pi*590/(60*10^3);

+Cwe=Cbe;

+P=178.9*10^3;

+C_we=(P+Cbi*Cwi)/(Cbe);

+Sf=C_we/Cwe;

+disp(Sf,"Slip factor is:");

diff --git a/2175/CH11/EX11.2/11_2.sce b/2175/CH11/EX11.2/11_2.sce
new file mode 100755
index 000000000..f83cc0195
--- /dev/null
+++ b/2175/CH11/EX11.2/11_2.sce
@@ -0,0 +1,50 @@
+clc;

+k=0.9;

+Cri1=486;//m/s

+Cri2=187.5;//m/s

+Caei=327;//m/s

+Cre1=k*Cri1;

+Cre2=k*Cri2;

+Cai2=k*Caei;

+//from velocity diagram;

+disp("inlet blade angle firls row of moving blades");

+Bi1=20;

+disp("degree",Bi1)

+

+disp("inlet blade angle fixed blades");

+alpha=20;

+disp("degree",alpha)

+

+disp("inlet blade angle second row of moving blades");

+Bi2=34.5;

+disp("degree",Bi2);

+

+//part II

+m=1;

+Cw1=874;

+Cw2=292.5;

+disp("N",m*Cw1,"driving force on first row:");

+disp("N",m*Cw2,"driving force on second row:");

+

+Cfi1=167;

+Cfe1=135;

+Cfi2=106;

+Cfe2=97;

+At1=m*(Cfi1-Cfe1);

+At2=m*(Cfi2-Cfe2);

+disp("N per kg/s",(At1+At2),"Total axial thrust:");

+

+//part III

+T_df=Cw1+Cw2

+disp("N per kg/s",T_df,"total driving force");

+bv=120

+P=T_df*bv/10^3;

+Cai1=600;

+E=m*Cai1^2/(2*10^3);

+De=P/E;

+disp("%",De*100,"diagramefficiency is");

+

+//partIV

+alpha_i=16*%pi/180;

+M=cos(alpha_i)^2;

+disp("%",M*100,"Maximum diagram efficienty is:");

diff --git a/2175/CH11/EX11.3/11_3.sce b/2175/CH11/EX11.3/11_3.sce
new file mode 100755
index 000000000..3dc3626f9
--- /dev/null
+++ b/2175/CH11/EX11.3/11_3.sce
@@ -0,0 +1,28 @@
+clc;

+Cai=600;

+alpha_i=16*%pi/180;

+l=25/1000;

+m=5;

+vi=0.375;

+n=m*vi/(Cai*sin(alpha_i)*l);

+disp("m",n,"length of nozzle arc is:");

+

+//part II

+p=0.025;

+Beta_1=18*%pi/180;

+Cre=437;

+t=0.0005;

+l1=m*vi*p/n/(p*sin(Beta_1)-t)/Cre;

+bhm=l1;

+

+Beta_2=21*%pi/180;

+Crf=294;

+lf=m*vi*p/[n*(p*sin(Beta_2)-t)*Crf];

+bhf=lf

+

+Beta_3=35*%pi/180;

+Crf2=169;

+l2=m*vi*p/n/(p*sin(Beta_3)-t)/Crf2;

+

+disp("Blade height at exit of first row, fixed and second row is respectively");

+disp("mm",l2*1000,"mm",bhf*1000,"mm",bhm*1000);

diff --git a/2175/CH11/EX11.4/11_4.sce b/2175/CH11/EX11.4/11_4.sce
new file mode 100755
index 000000000..570ac9a2c
--- /dev/null
+++ b/2175/CH11/EX11.4/11_4.sce
@@ -0,0 +1,31 @@
+clc;

+Cai=90;

+alpha=20*%pi/180;

+Cf=Cai*sin(alpha)

+

+Cb=4*Cf/3;

+

+v=0.6686;//m^3/kg

+m=9000/3600;

+A=m*v/Cf

+h=0.04;

+r=A/(2*%pi*h)

+N=Cb/(A/h)

+disp("rev/s",N,"Wheel speed is:")

+

+//partII

+Cw=2*Cai*cos(alpha)-Cb;

+DP=m*Cb*Cw;

+disp("kW",DP/1000,"diagram powar is:");

+

+//part III

+R=Cb*Cw

+Cri=[(Cai^2)+(Cb^2)-(2*Cai*Cb*cos(alpha))]^0.5

+Ei=Cai^2-(Cri^2/2)

+DE=R/Ei

+disp("%",DE*100,"diagram efficiency is:");

+

+//part IV

+Ed=(Cai^2-Cri^2)/2;

+Td=2*Ed;

+disp("kJ/kg",Td/1000,"total enthalpy drop per stage:")

diff --git a/2175/CH11/EX11.5/11_5.sce b/2175/CH11/EX11.5/11_5.sce
new file mode 100755
index 000000000..3dfa7cdcb
--- /dev/null
+++ b/2175/CH11/EX11.5/11_5.sce
@@ -0,0 +1,13 @@
+clc;

+Cw=115;//m/s

+Cb=200;//m/s

+wf=0.86;

+P=(Cw*Cb*wf)/1000;

+CP=12*P;

+T=20+273;

+y=1.4;

+ET=T*6^[(y-1)/y];

+cp=1.005;

+sp=cp*(ET-T);

+Ce=sp/CP;

+disp("%",Ce*100,"compressor isintropic efficiency is:");