"Single Stage Vapor Compression Refrigeration System"
t_evap = -10
t_cond = 40
p_evap = pressure(R22, t = t_evap, x = 0)
p_cond = pressure(R22, t = t_cond, x = 0)
DSC = 11
Q_ref = 1
"State 1: Condenser Exit, Saturated Liquid"
h_1 = enthalpy(R22, t = t_cond, x = 0)
h[1] = h_1
v_1 = volume(R22, t = t_cond, x = 0)
v[1] = v_1
s_1 = entropy(R22, t = t_cond, x = 0)
s[1] = s_1
t[1] = t_cond
p[1] = p_cond
"State 2: Expansion Valve Exit, Saturated Mixture"
h_2 = h_1
h[2] = h_2
v_2 = volume(R22, p = p_evap, h = h_2)
v[2 ] = v_2
s_2 = entropy(R22, p = p_evap, h = h_2)
s[2] = s_2
t[2] = t_evap
p[2] = p_evap
x_2 = quality(R22,p = p_evap, h = h_2)
"State 3: Evaporator Exit, Saturated Vapor"
h_3 = enthalpy(R22, t = t_evap, x = 1)
h[3] = h_3
v_3 = volume(R22, t = t_evap, x = 1)
v[3] = v_3
s_3 = entropy(R22, t = t_evap, x = 1)
s[3] = s_3
t[3] = t_evap
p[3] = p_evap
"State 4: Compressor Exit, Superheated Vapor"
s_4 = s_3
h_4 = enthalpy(R22, p = p_cond, s = s_4)
h[4] = h_4
v_4 = volume(R22, p = p_cond, s = s_4)
v[4] = v_4
t_4 = temperature(R22, p = p_cond, s = s_4)
p[4] = p_cond
Q_ref = m_ref*(h_3 - h_2)
W_comp = m_ref*(h_4 - h_3)
COP = Q_ref/W_comp
"MAIN SYSTEM"
t_evap_main = -15
t_cond_main = 40
p_evap_main = pressure(R22, t = t_evap_main, x = 0)
p_cond_main = pressure(R22, t = t_cond_main, x = 0)
DSC = 10
m_ref = 0.01719
W_comp_actual= 1.28
"Subsystem"
t_evap_subsystem= 10
t_cond_subsystem= 40
p_evap_subsystem= pressure(R22, t = t_evap_subsytem, x = 0)
p_cond_subsystem= pressure(R22, t = t_cond_subsystem, x = 0)
"State 1: Evaporator Exit, Saturated vapor"
h_1 = enthalpy(R22, t = t_evap_main, x = 1)
h[1] = h_1
v_1 = volume(R22, t = t_evap_main, x = 1)
v[1] = v_1
s_1 = entropy(R22, t = t_evap_main, x = 1)
s[1] = s_1
t[1] = t_evap_main
p[1] = p_evap_main
"State 2: Compressor Exit, Superheated Vapor"
s_2 = s_1
h_2 = enthalpy(R22, p = p_cond_main, s = s_2)
h[2] = h_2
v_2 = volume(R22, p = p_cond_main, s = s_2)
v[2] = v_2
t_2 = temperature(R22, p = p_cond_main, s = s_2)
p[2] = p_cond_main
"State 3: Condenser Exit, Saturated liquid "
h_3 = enthalpy(R22, t = t_cond_main, x = 0)
h[3] = h_3
v_3 = volume(R22, t = t_cond_main, x = 0)
v[3] = v_3
s_3 = entropy(R22, t = t_cond_main, x = 0)
s[3] = s_3
t[3] = t_cond_main
p[3] = p_cond_main
"State 4: Subcooler exit, Saturated liquid"
h_4 = enthalpy(R22, t = t_cond_main-DSC, x = 0)
h[4] = h_4
v_4 = volume(R22, t = t_cond_main-DSC, x = 0)
v[4] = v_4
s_4 = entropy(R22, t = t_cond_main-DSC, x = 0)
s[4] = s_4
t[4] = t_cond_main
p[4] = p_cond_main
"State 5: Expansion Valve Exit, Saturated mixture"
h_5 = h_4
h[5] = h_5
v_5 = volume(R22, p = p_evap_main , h = h_5)
v[5] = v_5
s_5 = entropy(R22, p = p_evap_main , h = h_5)
s[5] = s_5
t[5] = t_evap_main
p[5] = p_evap_main
x_5 = quality(R22,p = p_evap_main, h = h_5)
"State a: Evaporator Exit, Saturated Vapor"
h_6= enthalpy(R22, t = t_evap_subsytem, x = 1)
h[6] = h_6
v_6 = volume(R22, t = t_evap_subsystem, x = 1)
v[6] = v_6
s_6 = entropy(R22, t = t_evap_subsystem, x = 1)
s[6] = s_6
t[6] = t_evap_subsystem
p[6] = p_evap_subsystem
"State b: Compressor Exit, Superheated Vapor"
s_7 = s_7
h_7 = enthalpy(R22, p = p_cond_subsystem, s = s_7)
h[7] = h_7
v_7 = volume(R22, p = p_cond_subsystem, s = s_7)
v[7] = v_7
t_7 = temperature(R22, p = p_cond_subsystem, s = s_7)
p[7] = p_cond_subsystem
"State c: Condenser Exit, Saturated liquid "
h_8= enthalpy(R22, t = t_cond_subsytem, x = 0)
h[8] = h_8
v_8 = volume(R22, t = t_cond_subsystem, x = 0)
v[8] = v_8
s_8 = entropy(R22, t = t_cond_subsystem, x = 0)
s[8] = s_8
t[8] = t_cond_subsystem
p[8] = p_cond_subsytem
"State d: Expansion Valve Exit, Saturated mixture"
h_9 = h_8
h[9] = h_9
v_9 = volume(R22, p = p_evap_subsystem , h = h_9)
v[9] = v_9
s_9 = entropy(R22, p = p_evap_subsystem , h = h_9)
s[9] = s_9
t[9] = t_evap_subsystem
p[9] = p_evap_subsystem
x_9 = quality(R22,p = p_evap_subsystem, h = h_9)
Q_ref_main= m_ref*(h_1 - h_4)
W_comp_main_ideal = m_ref*(h_2 - h_1)
eta_main= W_comp_main_ideal/W_comp_actual
Q_main_Subcooler=m_ref(h_3 - h_4)
Q_main_Subcooler= Q_evap_subsytem
Q_evap_subsytem=m_subcooler*(
COP = Q_ref_main/W_comp_main_ideal
