Scribd
Upload
28 views6 pages

EES Sample Program

Uploaded by

Abdurrafie H. Abdulcader
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
28 views6 pages

EES Sample Program

Uploaded by

Abdurrafie H. Abdulcader
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 6

"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

You might also like