Refrigeration cycle with Freon 12
Housam Binous and Ahmed Bellagi

Thermophysical Properties of freon 12 from: https://webbook.nist.gov/cgi/fluid.cgi?ID=C75718&Action=Page
saturation properties - pressure increment = 10 psia ​
Temperature (F), Pressure (psia), Enthalpy (Btu/lbm), Entropy (Btu/lbm*R), Enthalpy (Btu/lbm), Entropy (Btu/lbm*R)
In[]:=
datasat={{-109.30,1.0000,-14.111,-0.036683,65.404,0.19026},{-33.428,11.000,1.3769,0.0032475,74.023,0.17368},{-5.8412,21.000,7.2436,0.016537,77.144,0.17056},{12.827,31.000,11.301,0.025253,79.215,0.16899},{27.326,41.000,14.506,0.031890,80.789,0.16800},{39.348,51.000,17.201,0.037315,82.064,0.16730},{49.707,61.000,19.554,0.041938,83.139,0.16677},{58.862,71.000,21.658,0.045989,84.066,0.16634},{67.101,81.000,23.572,0.049609,84.881,0.16600},{74.615,91.000,25.336,0.052892,85.606,0.16570},{81.539,101.00,26.978,0.055904,86.257,0.16543},{87.974,111.00,28.519,0.058691,86.846,0.16520},{93.995,121.00,29.974,0.061292,87.382,0.16498},{99.659,131.00,31.355,0.063732,87.873,0.16478},{105.01,141.00,32.673,0.066035,88.323,0.16459},{110.10,151.00,33.935,0.068217,88.737,0.16440},{114.94,161.00,35.147,0.070294,89.118,0.16422},{119.56,171.00,36.316,0.072278,89.470,0.16405},{123.99,181.00,37.445,0.074178,89.795,0.16387},{128.25,191.00,38.539,0.076003,90.096,0.16370},{132.34,201.00,39.600,0.077760,90.373,0.16352},{136.29,211.00,40.633,0.079457,90.628,0.16335},{140.10,221.00,41.639,0.081097,90.863,0.16317},{143.78,231.00,42.621,0.082687,91.079,0.16299},{147.36,241.00,43.581,0.084231,91.277,0.16280},{150.82,251.00,44.521,0.085731,91.457,0.16261},{154.18,261.00,45.442,0.087193,91.620,0.16242},{157.45,271.00,46.345,0.088619,91.766,0.16222},{160.63,281.00,47.234,0.090011,91.897,0.16201},{163.73,291.00,48.107,0.091373,92.013,0.16180},{166.75,301.00,48.967,0.092706,92.113,0.16158},{169.70,311.00,49.816,0.094013,92.198,0.16136},{172.57,321.00,50.653,0.095297,92.269,0.16112},{175.39,331.00,51.480,0.096558,92.325,0.16088},{178.13,341.00,52.297,0.097799,92.366,0.16062},{180.82,351.00,53.107,0.099021,92.392,0.16036},{183.45,361.00,53.909,0.10023,92.403,0.16008},{186.02,371.00,54.704,0.10142,92.399,0.15979},{188.55,381.00,55.494,0.10259,92.380,0.15950},{191.02,391.00,56.279,0.10376,92.345,0.15918},{193.44,401.00,57.060,0.10491,92.293,0.15885},{195.82,411.00,57.838,0.10605,92.224,0.15851},{198.15,421.00,58.614,0.10719,92.138,0.15815},{200.44,431.00,59.388,0.10832,92.032,0.15777},{202.69,441.00,60.163,0.10944,91.907,0.15737},{204.90,451.00,60.938,0.11056,91.761,0.15694},{207.07,461.00,61.716,0.11168,91.592,0.15649},{209.20,471.00,62.497,0.11280,91.398,0.15601},{211.29,481.00,63.284,0.11393,91.177,0.15550},{213.35,491.00,64.077,0.11506,90.926,0.15495},{215.38,501.00,64.881,0.11620,90.641,0.15436},{217.36,511.00,65.697,0.11736,90.317,0.15373},{219.32,521.00,66.530,0.11854,89.948,0.15303},{221.24,531.00,67.385,0.11974,89.526,0.15226},{223.13,541.00,68.270,0.12099,89.038,0.15141},{224.98,551.00,69.196,0.12229,88.467,0.15044},{226.80,561.00,70.184,0.12368,87.785,0.14932},{228.59,571.00,71.268,0.12520,86.942,0.14797},{230.34,581.00,72.524,0.12696,85.838,0.14626},{232.05,591.00,74.191,0.12932,84.177,0.14375}};
In[]:=
tblhl=Table[{datasat[[i]][[2]],datasat[[i]][[3]]},{i,1,Length[datasat]}];​​tblsl=Table[{datasat[[i]][[2]],datasat[[i]][[4]]},{i,1,Length[datasat]}];
In[]:=
fhl=Interpolation[tblhl];fsl=Interpolation[tblsl];
In[]:=
tblhv=Table[{datasat[[i]][[2]],datasat[[i]][[5]]},{i,1,Length[datasat]}];​​tblsv=Table[{datasat[[i]][[2]],datasat[[i]][[6]]},{i,1,Length[datasat]}];
In[]:=
fhv=Interpolation[tblhv];fsv=Interpolation[tblsv];
isobaric properties - pressure = 84.755 psia ​
{Temperature (F), Pressure (psia), Enthalpy (Btu/lbm), Entropy (Btu/lbm*R), State)}
In[]:=
dataisobaric={{10.000,84.755,10.732,0.023809,"liq"},{15.000,84.755,11.824,0.026122,"liq"},{20.000,84.755,12.922,0.028423,"liq"},{25.000,84.755,14.025,0.030712,"liq"},{30.000,84.755,15.135,0.032989,"liq"},{35.000,84.755,16.250,0.035256,"liq"},{40.000,84.755,17.372,0.037512,"liq"},{45.000,84.755,18.501,0.039759,"liq"},{50.000,84.755,19.636,0.041997,"liq"},{55.000,84.755,20.778,0.044228,"liq"},{60.000,84.755,21.928,0.046451,"liq"},{65.000,84.755,23.085,0.048667,"liq"},{69.999,84.755,24.250,0.050877,"liq"},{69.999,84.755,85.163,0.16588,"vap"},{70.000,84.755,85.163,0.16588,"vap"},{75.000,84.755,85.983,0.16742,"vap"},{80.000,84.755,86.800,0.16894,"vap"},{85.000,84.755,87.615,0.17044,"vap"},{90.000,84.755,88.427,0.17193,"vap"},{95.000,84.755,89.238,0.17340,"vap"},{100.00,84.755,90.048,0.17485,"vap"},{105.00,84.755,90.857,0.17629,"vap"},{110.00,84.755,91.666,0.17772,"vap"},{115.00,84.755,92.474,0.17913,"vap"},{120.00,84.755,93.283,0.18053,"vap"}};
In[]:=
pos=Position[dataisobaric,69.9990]
Out[]=
{{13,1},{14,1}}
In[]:=
tblhsv=Table[{dataisobaric[[i]][[1]],dataisobaric[[i]][[3]]},{i,14,Length[dataisobaric]}];​​tblssv=Table[{dataisobaric[[i]][[1]],dataisobaric[[i]][[4]]},{i,14,Length[dataisobaric]}];​​tblshsv=Table[{dataisobaric[[i]][[4]],dataisobaric[[i]][[3]]},{i,15,Length[dataisobaric]}];
In[]:=
fshsv=Interpolation[tblshsv];
In[]:=
tblhl1=Table[{datasat[[i]][[3]],datasat[[i]][[2]]},{i,1,Length[datasat]}];​​tblsl1=Table[{datasat[[i]][[4]],datasat[[i]][[1]]},{i,1,Length[datasat]}];
In[]:=
tblhv1=Table[{datasat[[i]][[5]],datasat[[i]][[2]]},{i,1,Length[datasat]}];​​tblsv1=Table[{datasat[[i]][[6]],datasat[[i]][[1]]},{i,1,Length[datasat]}];
In[]:=
fhsv=Interpolation[tblhsv];fssv=Interpolation[tblssv];
In[]:=
plt1=ListLogPlot[{tblhv1,tblhl1},JoinedTrue,FrameTrue,AspectRatio1,GridLinesAutomatic,FrameLabel{Style["Enthalpy (Btu/lbm)",14],Style["Pressure (psia)",14]}]
Out[]=
In[]:=
plt2=ListPlot[{tblsv1,tblsl1},JoinedTrue,FrameTrue,AspectRatio1,GridLinesAutomatic,FrameLabel{Style["Entropy (Btu/lbm*°F)",14],Style["Temperature (°F)",14]}]
Out[]=
We propose to study the ordinary vapor compression cycle for refrigerant FREON 12 GIVEN: TC = -10°F and TH = 70°F for full explanations see: https://www.youtube.com/watch?v=qD04XetRVzk&t=44s​​
​​First, let us determine the unknown variables for points {1, 2, 3, 4} of the Table below:​​
point
T °F
P psia
h Btu/lbm
s Btu/lbm*F
state
2
T2
P2
H2
S2
sat. vap.
3
T3
P3
H3
S3
super heated vap.
4
T4
P4
H4
S4
sat. liq.
1
T1
P1
H1
S1
quality
The OVC Cycle drawn on the H-P diagram
The OVC Cycle drawn on the T-S diagram
Let us determine the COP of the refrigeration cycle if compressor is only 70% efficient