diff options
Diffstat (limited to '1808/CH3/EX3.25')
-rw-r--r-- | 1808/CH3/EX3.25/Chapter3_Exampl25.sce | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/1808/CH3/EX3.25/Chapter3_Exampl25.sce b/1808/CH3/EX3.25/Chapter3_Exampl25.sce new file mode 100644 index 000000000..ba73a018d --- /dev/null +++ b/1808/CH3/EX3.25/Chapter3_Exampl25.sce @@ -0,0 +1,32 @@ +clc
+clear
+//INPUT DATA
+Rc=20;//compression ratio
+p1=101.325;//Pressure in kN/m^2
+t1=300;//temperature in K
+p32=2;//pressure ratio of heating process
+v43=1.5;//volume ratio of heating process
+g=1.4;//constant
+cp=1.005;//specific pressure
+cv=0.718;//specific volume
+R=0.287;//gas constant
+
+//CALCULATIONS
+t2=t1*(Rc^(g-1));//Temperature in K
+p2=p1*(Rc^g);//pressure in kN/m^2
+p3=p32*p2;//pressure in kN/m^2
+t3=p32*t2;//Temperature in K
+t4=v43*t3;//Temperature in K
+p5=p3/((Rc/v43)^g);//pressure in kN/m^2
+t5=t4/((Rc/v43)^(g-1));//Temperature in K
+nd=(((cv*(t3-t2))+(cp*(t4-t3))-(cv*(t5-t1)))/((cv*(t3-t2)+(cp*(t4-t3)))))*100;//The thermal efficiency in percentage
+x=((cv*(t3-t2)+(cp*(t4-t3))));//workdone
+y=(cv*(t5-t1));//workdone
+v1=R*t1/p1;//specific volume in m^3/kg
+pm=(x-y)/(v1*(1-(1/Rc)));//mean effective pressure in kN/m^2
+
+//OUTPUT
+printf('(a)The thermal efficiency is %3.2f percentage \n (b)The mean effective pressure is %3.1f kN/m^2',nd,pm)
+
+
+
|