diff options
Diffstat (limited to '1808/CH5/EX5.3/Chapter5_Exampl3.sce')
-rw-r--r-- | 1808/CH5/EX5.3/Chapter5_Exampl3.sce | 37 |
1 files changed, 37 insertions, 0 deletions
diff --git a/1808/CH5/EX5.3/Chapter5_Exampl3.sce b/1808/CH5/EX5.3/Chapter5_Exampl3.sce new file mode 100644 index 000000000..0dbacbe00 --- /dev/null +++ b/1808/CH5/EX5.3/Chapter5_Exampl3.sce @@ -0,0 +1,37 @@ +clc
+clear
+//INPUT DATA
+t1=293;//temperature in K
+p1=1;//pressure in bar
+p2=8;//pressure in bar
+v1=80;//volume in m^3
+g=1.4;//for isentropic compression
+n=1.25;//Adiabatic compression
+
+//CALCULATIONS
+
+//Isothermal compression
+v2=p1*v1/p2;//volume in m^3
+wn=p1*10^2*v1*log(p2/p1);//net work done in kJ/min
+P=wn/60;//Power required in kW
+
+//Adiabatic compression
+v121=(p2/p1)^(1/g);//volume in m^3/min
+v21=v1/v121;//volume in m^3/min
+t211=t1*(p2/v121);//temperature at the end of the compression
+wn1=(g/(g-1))*p1*10^2*v1*(((p2/p1)^((g-1)/g))-1);//net work done in kJ
+P1=wn1/60;//Power in kW
+
+//Polytropic process
+v12=(p2/p1)^(1/n);//volume in m^3/min
+v22=v1/v12;//volume in m^3/min
+t2=t1*(p2/v12);//temperature in K
+wn2=(n/(n-1))*p1*10^2*v1*(((p2/p1)^((n-1)/n))-1);//net work done in kJ
+P2=wn2/60;//power required in kW
+Qc=((g-n)/(g-1))*P2;//Heat transferred in kW
+
+//OUTPUT
+printf('(a)Isothermal compression \n temperature at the end of the compression is 293K \n net work done is %3.2f kJ/min \n Power required is %3.3f kW \n ',wn,P)
+printf('(b)Adiabatic compression \n temperature at the end of the compression is %3.2f K \n net work done is %3.2f kJ/min \n Power required is %3.3f kW \n ',t211,wn1,P1)
+printf('(c)Polytropic compression \n temperature at the end of the compression is %3.2f K \n net work done is %3.2f kJ/min \n Power required is %3.3f kW \n Heat rejected is %3.2f kW\n ',t2,wn2,P2,Qc)
+
|