1 files changed, 30 insertions, 0 deletions
diff --git a/1808/CH3/EX3.18/Chapter3_Exampl18.sce b/1808/CH3/EX3.18/Chapter3_Exampl18.sce
new file mode 100644
index 000000000..e91c3483e
--- /dev/null
+++ b/1808/CH3/EX3.18/Chapter3_Exampl18.sce
@@ -0,0 +1,30 @@
+clc
+clear
+//INPUT DATA
+p1=101.325;//Pressure in kN/m^2
+t1=303;//Temperature in K
+g=1.4;//constant
+cp=1.005;//specific pressure in kJ/kgK
+cv=0.718;//specific volume in kJ/kgK
+R=0.287;//gas constant in kJ/kgK
+r=2.5;//cut off ratio
+v2=0.1;//clearance volume in m^3/kg
+
+
+//CALCULATIONS
+v1=R*t1/p1;//volume at state 1 in m^3/kg
+p2=p1*((v1/v2)^(g));//pressure at state 2 in kN/m^2
+t2=t1*((v1/v2));//temperature at state 2 in K
+v3=r*v2;//volume at state 3 in m^3/kg
+t3=t2*(v3/v2);//Temperature at state 3 in K
+t4=t3*((v3/v1)^(g-1));//Temperature at state 4 in K
+p4=p2*((v3/v1)^g);//Pressure at state 4 in kN/m^2
+c=(v2/(v1-v2))*100;//percentage clearance
+Qs=cp*(t3-t2);//heat added in kJ/kg
+Qr=cv*(t4-t1);//heat rejected in kJ/kg
+nd=((Qs-Qr)/Qs)*100;//Thermal efficiency in percentage
+pm=((Qs-Qr)/(v1-v2));//Mean effective pressure in kN/m^2
+
+//OUTPUT
+printf('(a)pressure at state 2 is %3.2f kN/m^2 \n  temperature at state 2 is %3.2f K \n  volume at state 3 is %3.2f m^3/kg\n  Temperature at state 3 is %3.2f K \n  Temperature at state 4 is %3.2f K \n Pressure at state 4 is %3.2f kN/m^2 \n (b)percentage clearance is %3.2f percentage \n(c)Thermal efficiency is %3.2f percentage \n (d)Mean effective pressure is %3.2f kN/m^2',p2,t2,v3,t3,t4,p4,c,nd,pm)
+