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/CH18/EX18.7/Ex18_7.sce | |
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/CH18/EX18.7/Ex18_7.sce')
-rwxr-xr-x | 2223/CH18/EX18.7/Ex18_7.sce | 25 |
1 files changed, 25 insertions, 0 deletions
diff --git a/2223/CH18/EX18.7/Ex18_7.sce b/2223/CH18/EX18.7/Ex18_7.sce new file mode 100755 index 000000000..cbf36c7ff --- /dev/null +++ b/2223/CH18/EX18.7/Ex18_7.sce @@ -0,0 +1,25 @@ +// scilab Code Exa 18.7 RHF of a three stage turbine
+
+p1=1.0; // Initial Pressure in bar
+gamma=1.4;
+T1=1500; // Initial Temperature in K
+s=3; // number of stages
+opr=11; // Overall Pressure Ratio
+pr=opr^(1/s); // equal Pressure Ratio in each stage
+n_T=0.88; // Overall Efficiency
+delTa=T1*(1-opr^(-((gamma-1)/gamma)))*n_T;
+T2=T1-delTa;
+n_p=(log(T1/T2))/(((gamma-1)/gamma)*(log(opr))); // polytropic or small stage efficiency
+cp=1.005; // Specific Heat at Constant Pressure in kJ/(kgK)
+n_st=(1-pr^(n_p*(-((gamma-1)/gamma))))/(1-pr^(-((gamma-1)/gamma))); // stage efficiency
+T(1)=T1;
+for i=1:3
+ delT(i)=T(i)*(1-pr^(n_p*(-((gamma-1)/gamma))));
+ delw_s(i)=delT(i)*cp/n_st;
+ T(i+1)=T(i)-delT(i);
+end
+w_a=cp*delTa;
+w_s=w_a/n_T;
+RHF=(delw_s(1)+delw_s(2)+delw_s(3))/w_s;
+disp(RHF,"the reheat factor is")
+
|