summaryrefslogtreecommitdiff
path: root/2168/CH5/EX5.8/Chapter5_example8.sce
diff options
context:
space:
mode:
Diffstat (limited to '2168/CH5/EX5.8/Chapter5_example8.sce')
-rwxr-xr-x2168/CH5/EX5.8/Chapter5_example8.sce43
1 files changed, 43 insertions, 0 deletions
diff --git a/2168/CH5/EX5.8/Chapter5_example8.sce b/2168/CH5/EX5.8/Chapter5_example8.sce
new file mode 100755
index 000000000..6a6656729
--- /dev/null
+++ b/2168/CH5/EX5.8/Chapter5_example8.sce
@@ -0,0 +1,43 @@
+clc
+clear
+//Input data
+r=14//Compression ratio
+c=5//Cut off takes place at 5% of the stroke
+p1=1//Pressure at the beginning of compression in kg/cm^2. In texbook, it is given wrong as 10
+T1=90+273//Temperature at the beginning of compression in K
+p3=50//Maximum pressure in kg/cm^2
+R=29.3//Characteristic gas constant in mkg/kg/degree C
+cv=[0.171,0.00003]//Cv=0.171+0.00003T where Cv is Specific heat at constant volume and T is the temperature in K
+g1=1.4//Ratio of specific heats
+J=427//Mechanical equivalent of heat in kg.m/kcal
+
+//Calculations
+T2x=(T1*r^(g1-1))//Temperature in K
+function y=f1(T),y=cv(1)+(cv(2)*T),endfunction
+I1=intg(T1,T2x,f1)
+Cv=(1/(T2x-T1))*I1//Mean value of Cv in kJ/kg.K
+Cp=(Cv+(R/J))//Mean value of Cp in kJ/kg.K
+g=1.35//(Cp/Cv) value and rounded off to 2 decimal places for calculation purpose. Ratio of specific heats
+T2=(T1*r^(g-1))//Temperature in K
+I2=intg(T1,T2,f1)
+CV=(1/(T2-T1))*I2//Maen value of Cv in kJ/kg.K
+CP=(Cv+(R/J))//Mean value of Cp in kJ/kg.K
+g2=1.36//(Cp/Cv) value and rounded off to 2 decimal places for calculation purpose.Ratio of specific heats
+T2a=(T1*r^(g2-1))//Temperature in K
+p2=(p1*r*(T2a/T1))//Pressure in kg/cm^2
+T3=(T2a*(p3/p2))//Temperature in K
+T4=(((r-1)*(c/100))+1)*T3//Temperature in K
+g3=1.3//Assuming gamma as 1.3 for process 4-5
+T5=(T4/(r/(((r-1)*(c/100))+1))^(g3-1))//Temperature in K
+cV=cv(1)+(cv(2)/2)*(T5+T4)//Mean value of Cv in kJ/kg.K
+cP=cV+(R/J)//Mean value of Cp in kJ/kg.K
+g4=(cP/cV)//Ratio of specific heats
+T5a=(T4/(r/(((r-1)*(c/100))+1))^(g4-1))
+I3=intg(T2a,T3,f1)
+function y=f2(T),y=cv(1)+(R/J)+(cv(2)*T),endfunction
+I4=intg(T3,T4,f2)//Textbook answer is wrong
+q=I3+I4//Heat supplied per kg of working substance in kcal/kg
+
+//Output
+printf('(a) Temperatures at all the points of the cycle are: \n T1 = %i K \n T2 = %3.0f K \n T3 = %3.0f K \n T4 = %3.0f K \n T5 = %i K \n\n (b) heat supplied per kg of the working substance is %3.1f kcal/kg',T1,T2a,T3,T4,T5a,q)
+//Textbook answer is wrong