1 files changed, 24 insertions, 0 deletions
diff --git a/2063/CH1/EX1.15/1_15.sce b/2063/CH1/EX1.15/1_15.sce
new file mode 100755
index 000000000..d883e1753
--- /dev/null
+++ b/2063/CH1/EX1.15/1_15.sce
@@ -0,0 +1,24 @@
+clc
+clear
+//Input data
+r=18;//Compression ratio of diesel cycle
+Q=2000;//Heat added in kJ/kg
+T1=300;//Lowest temperature in the cycle in K
+p1=1;//Lowest pressure in the cycle in bar
+Cp=1;//Specific heat of air at constant pressure in kJ/kg K
+Cv=0.714;//Specific heat of air at constant volume in kJ/kg K
+
+//Calculations
+r1=Cp/Cv;//Isentropic ratio
+v1=((1-Cv)*T1)/(p1*10^5);//Initial volume at point 1 in the graph in m^3/kg
+v2=v1/r;//Volume at point 2 in m^3/kg
+p2=p1*(v1/v2)^(r1);//Pressure at point 2 in bar
+T2=T1*(v1/v2)^(r1-1);//Temperature at point 2 in K
+T3=(Q/Cp)+T2;//Temperature at point 3 in K
+v3=v2*(T3/T2);//Volume at point 3 in K
+v4=v1;//Since Constant volume heat rejection in m^3/kg
+T4=T3/(v4/v3)^(r1-1);//Temperature at point 4 in K for isentropic expansion
+p4=p1*(T4/T1);//Pressure at point 4 in bar
+
+//Output
+printf('(a)Pressure at point 1 in the cycle is %3.0f bar\n (b)Pressure at point 2 & 3 is %3.1f bar\n (c)Pressure at point 4 is %3.2f bar\n (d)Temperature at point 1 is %3.0f K\n (e)Temperature at point 2 is %3.0f K\n (f)Temperature at point 3 is %3.0f K\n (g)Temperature at point 4 is %3.0f K',p1,p2,p4,T1,T2,T3,T4)