diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /2223/CH8/EX8.5 | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '2223/CH8/EX8.5')
-rwxr-xr-x | 2223/CH8/EX8.5/Ex8_5.sav | bin | 0 -> 56640 bytes | |||
-rwxr-xr-x | 2223/CH8/EX8.5/Ex8_5.sce | 46 |
2 files changed, 46 insertions, 0 deletions
diff --git a/2223/CH8/EX8.5/Ex8_5.sav b/2223/CH8/EX8.5/Ex8_5.sav Binary files differnew file mode 100755 index 000000000..67fda5ef3 --- /dev/null +++ b/2223/CH8/EX8.5/Ex8_5.sav diff --git a/2223/CH8/EX8.5/Ex8_5.sce b/2223/CH8/EX8.5/Ex8_5.sce new file mode 100755 index 000000000..2dd135a13 --- /dev/null +++ b/2223/CH8/EX8.5/Ex8_5.sce @@ -0,0 +1,46 @@ +// scilab Code Exa 8.5 compressor type radial cascade tunnel
+
+M=0.7; // Mach Number
+pr=0.721; // pr=pt/p0 From isentropic gas tables
+t_opt=0.911; // t_opt=Tt/T0
+pa=1.013; // Atmospheric Pressure in bar
+Ta=306; // in K
+n_c=0.65; // efficiency
+R=288;
+gamma=1.4;
+alpha=30;
+dm=45/100; // mean diameter of the impeller blade in m
+b=10/100; // blade width in m
+cp_a=1.008; // Specific Heat of air at Constant Pressure in kJ/(kgK)
+
+// part(a) pressure ratio of the compressor
+pr_c=1/pr;
+disp(pr_c,"(a)pressure ratio of the compressor is")
+
+// part(b) stagnation pressure in the settling chamber
+p02=pa*pr_c;
+disp("bar",p02,"(b)stagnation pressure in the settling chamber is")
+
+// part(c)test section conditions(static pressure, temperature and velocity)
+n=(gamma-1)/gamma;
+T02s=Ta*(pr_c^((gamma-1)/gamma));
+T02=Ta+((T02s-Ta)/n_c);
+T_t=t_opt*T02;
+p_t=pr*p02;
+c_t=M*sqrt(gamma*R*T_t);
+disp("(c)test section conditions are given by: ")
+disp("bar",p_t,"static pressure of air in the test section is")
+disp("K",T_t,"static temperature of air in the test section is")
+disp("m/s",c_t,"velocity of air in the test section is")
+
+// part(d) determining mass flow rate
+c_r=c_t*sind(alpha);
+ro_t=p_t*1e5/(R*T_t); // density of air in kg/m3
+A_t=%pi*dm*b;
+m=ro_t*A_t*c_r;
+disp("kg/s",m,"(d) mass flow rate of compressor is")
+
+// part (e) determining the power required to drive the air compressor
+delh_s=cp_a*(T02-Ta);
+P=m*delh_s;
+disp("kW",P,"(e)Power required to drive the air compressor is")
|