1 files changed, 35 insertions, 0 deletions
diff --git a/3685/CH14/EX14.3/Ex14_3.sce b/3685/CH14/EX14.3/Ex14_3.sce
new file mode 100644
index 000000000..5715fc8c5
--- /dev/null
+++ b/3685/CH14/EX14.3/Ex14_3.sce
@@ -0,0 +1,35 @@
+clc
+h1 = 183.19 // Enthalpy at state 1 in kJ/kg
+h2 = 209.41 // Enthalpy at state 2 in kJ/kg
+h3 = 74.59 // Enthalpy at state 3 in kJ/kg
+h4 = h3 // Isenthalpic process
+T1 = 40 // Evaporator temperature in degree Celsius 
+T2 = -10 // Condenser temperature in degree Celsius
+W = 5 // Plant capacity in tonnes of refrigeration
+w = (W*14000/3600)/(h1-h4) // Refrigerant flow rate
+v1 = 0.077 // Specific volume of vapor in m^3/kg
+VFR = w*v1 // volume flow rate
+T = 48 // Compressor discharge temperature in degree Celsius
+P2 = 9.6066 // Pressure after compression
+P1 = 2.1912 // Pressure before compression
+rp = P2/P1 // Pressure ratio
+Q1 = w*(h2-h3) // Heat rejected in condenser
+hf = 26.87 // Enthalpy of fluid in kJ/kg
+hfg = 156.31// Latent heat of vaporization in kJ/kg
+x4 = (h4-hf)/hfg // quality of refrigerant
+COP_v = (h1-h4)/(h2-h1) // Actual coefficient of performance of cycle
+PI = w*(h2-h1) // Power input
+COP = (T2+273)/((T1+273)-(T2+273)) // Ideal coefficient of performance
+r = COP_v/COP
+printf("\n Example 14.3\n")
+printf("\n Refrigerant flow rate is %f kg/s",w)
+printf("\n Volume flow rate is %f m^3/s",VFR)
+printf("\n Compressor discharge temperature is %d degree Celsius ",T)
+printf("\n Pressure ratio is %f ",rp)
+printf("\n Heat rejected to the condenser is %f kW",Q1)
+printf("\n Flash gas percentage is %f percent",x4*100)
+printf("\n COP is %f kW",COP_v)
+printf("\n Power required to drive the compressor is %f kW",PI)
+printf("\n Ratio of COP of cycle with  Carnot refrigerator is %f",r)
+//The answers vary due to round off error
+