Application Center - Maplesoft

App Preview:

Thermodynamic Calculations of Cascade Vapor Compression Refrigeration Cycle

You can switch back to the summary page by clicking here.

Learn about Maple
Download Application




Thermodynamic Calculations of  Cascade Vapor Compression Refrigeration Cycle

Introduction

 

In the event that a high COP of a refrigeration cycle is of greater importance compared to other factors, it is possible to significantly increase the COP of a basic cycle through the use of a multistage vapor compression cycle. This is especially true when the pressure ratio between the heat rejection and heat absorption pressures is large 5 or more.
Multistaging involves one or more intermediate pressures between the heat rejection and heat absorption pressures, and a series of compressors operating between successive pressure intervals.
Cascade cycle is simply a chain of single stage vapor compression cycles operating in series, such that the condenser of a lower temperature cycle provides the heat input to the evaporator of the higher temperature cycle.

Creation functions on properties and processes of working fluids

 
> 

restart

> 

with(ThermophysicalData):

 

Pressure in the saturated region as a function of temperature, quality and working fluid

> 

PSTXwf := proc (T, X, wf) options operator, arrow; Property(P, temperature = T, Q = X, wf) end proc:

 

Temperature in the saturated region as a function of pressure, quality and working fluid

> 

TSPXwf := proc (p, X, wf) options operator, arrow; Property(T, pressure = p, Q = X, wf) end proc:

 

Specific enthalpy in the saturated region as a function of temperature, quality and working fluid

> 

HSTXwf := proc (T, X, wf) options operator, arrow; Property(enthalpy, temperature = T, Q = X, wf) end proc:

 

Specific enthalpy in the saturated region as a function of pressure, quality and working fluid

> 

HSPXwf := proc (p, X, wf) options operator, arrow; Property(enthalpy, pressure = p, Q = X, wf) end proc:

 

Specific entropy in the saturated region as a function of temperature, quality and working fluid

> 

SSTXwf := proc (T, X, wf) options operator, arrow; Property(entropy, temperature = T, Q = X, wf) end proc:

 

Specific entropy in the saturated region as a function of pressure, quality and working fluid

> 

SSPXwf := proc (p, X, wf) options operator, arrow; Property(entropy, pressure = p, Q = X, wf) end proc:

 

Density in the saturated region as a function of temperature, quality and working fluid

> 

DSTXwf := proc (T, X, wf) options operator, arrow; Property(D, temperature = T, Q = X, wf) end proc:

 

Density in the saturated region as a function of pressure, quality and working fluid

> 

DSPXwf := proc (p, X, wf) options operator, arrow; Property(D, pressure = p, Q = X, wf) end proc:

 

Temperature as a function of pressure, specific entropy and working fluid

> 

TPSwf := proc (p, s, wf) options operator, arrow; Property(T, pressure = p, entropy = s, wf) end proc:

 

Temperature as a function of pressure, specific enthalpy and working fluid

> 

TPHwf := proc (p, h, wf) options operator, arrow; Property(T, pressure = p, enthalpy = h, wf) end proc:

 

Specific enthalpy as a function of pressure, temperature and working fluid

> 

HPTwf := proc (p, T, wf) options operator, arrow; Property(H, pressure = p, temperature = T, wf) end proc:

 

Specific entropy as a function of pressure, temperature and working fluid

> 

SPTwf := proc (p, T, wf) options operator, arrow; Property(S, pressure = p, temperature = T, wf) end proc:

 

Density as a function of pressure, temperature and working fluid

> 

DPTwf := proc (p, T, wf) options operator, arrow; Property(D, pressure = p, temperature = T, wf) end proc:

 

Specific entropy as function of pressure and specific enthalpy and working fluid

> 

SPHwf := proc (p, h, wf) options operator, arrow; Property(S, pressure = p, enthalpy = h, wf) end proc:

Density as function of pressure and specific enthalpy and working fluid

> 

DPHwf := proc (p, h, wf) options operator, arrow; Property(D, pressure = p, enthalpy = h, wf) end proc:

 

> 

HcompPTPEFFwf := proc (p1, T1, p2, eff, wf) local h1, s1, s2_is, T2_is, h2_is; if TSPXwf(p1, 1, wf) < T1 then h1 := HPTwf(p1, T1, wf); s1 := SPTwf(p1, T1, wf); s2_is := s1; T2_is := TPSwf(p2, s2_is, wf); h2_is := HPTwf(p2, T2_is, wf) else h1 := HSTXwf(T1, 1, wf); s1 := SSTXwf(T1, 1, wf); s2_is := s1; T2_is := TPSwf(p2, s2_is, wf); h2_is := HPTwf(p2, T2_is, wf) end if; return h1+(h2_is-h1)/eff end proc:

 

 

> 

ScompPTPEFFwf := proc (p1, T1, p2, eff, wf) local h2, T2; h2 := HcompPTPEFFwf(p1, T1, p2, eff, wf); T2 := TPHwf(p2, h2, wf); return SPTwf(p2, T2, wf) end proc:

 

> 

TcompPTPEFFwf := proc (p1, T1, p2, eff, wf) local h2; h2 := HcompPTPEFFwf(p1, T1, p2, eff, wf); return TPHwf(p2, h2, wf) end proc:

 

> 

DcompPTPEFFwf := proc (p1, T1, p2, eff, wf) local h2, T2; h2 := HcompPTPEFFwf(p1, T1, p2, eff, wf); T2 := TPHwf(p2, h2, wf); return DPTwf(p2, T2, wf) end proc:

 

 

Critical temperature as a function of working fluid

NULL

> 

TCRwf := proc (wf) options operator, arrow; Property(T_critical, water, useunits) end proc:

 

Critical pressure as a function of working fluid

NULL

> 

PCRwf := proc (wf) options operator, arrow; Property(pcritical, water, useunits) end proc:

 

Input Data

 

Temperature of heat source

> 

T__in := Units:-Standard:-`*`(Units:-Standard:-`+`(-10, 273.15), Unit('K'))

263.15*Units:-Unit(K)

(3.1)

Temperature of heat consumption

> 

T__out := Units:-Standard:-`*`(Units:-Standard:-`+`(70, 273.15), Units:-Standard:-Unit('K'))

343.15*Units:-Unit(K)

(3.2)

Temperature difference in the evaporator

> 

`&delta;T__EV` := Units:-Standard:-`*`(2, Units:-Standard:-Unit('K'))

2*Units:-Unit(K)

(3.3)

Temperature difference in the condenser

> 

`&delta;T__CD` := Units:-Standard:-`*`(2, Units:-Standard:-Unit('K'))

2*Units:-Unit(K)

(3.4)

Temperature difference between working fluids in low pressure and high pressure circuits

> 

`&delta;T` := Units:-Standard:-`*`(2, Units:-Standard:-Unit('K'))

2*Units:-Unit(K)

(3.5)

Isentropic efficiency of the compressor

> 

eta[comp] := .8

.8

(3.6)

Pressure increase in compressor 1

> 

pi__1 := 3

3

(3.7)

Working fluid in low pressure circuit

> 

wf[L] := R410a

R410a

(3.8)

Working fluid in high pressure circuit

> 

wf[H] := R134a

R134a

(3.9)

 

 

Calculations

 

Temperature of  the working fluid at the evaporator outlet

> 

T[1, L] := Units:-Standard:-`+`(T__in, Units:-Standard:-`-`(`&delta;T__EV`))

261.15*Units:-Unit(K)

(4.1)

Pressure of  the working fluid at the evaporator outlet

> 

p[1, L] := PSTXwf(T[1, L], 1, wf[L])

534123.7895*Units:-Unit(Pa)

(4.2)

Specific enthalpy of  the working fluid at the evaporator outlet

> 

h[1, L] := HSTXwf(T[1, L], 1, wf[L])

417301.4457*Units:-Unit(J/kg)

(4.3)

Specific entropy of  the working fluid at the evaporator outlet

> 

`s__1,L` := SSTXwf(T[1, L], 1, wf[L])

1834.540688*Units:-Unit(J/(kg*K))

(4.4)

Specific volume of  the working fluid at the evaporator outlet

> 

`v__1,L` := Units:-Standard:-`/`(DSTXwf(T[1, L], 1, wf[L]))

0.4881666858e-1*Units:-Unit(m^3/kg)

(4.5)

 

Pressure of  the working fluid at the compressor 1 outlet

> 

p[2, L] := Units:-Standard:-`*`(p[1, L], pi__1)

1602371.368*Units:-Unit(Pa)

(4.6)

Pressure of  the working fluid at the condenser outlet 

> 

p[3, L] := p[2, L]

 

1602371.368*Units:-Unit(Pa)

(4.7)

 

Temperature of  the working fluid at the condenser outlet

> 

`T__3,L` := TSPXwf(p[3, L], 0, wf[L])

296.8903325*Units:-Unit(K)

(4.8)

Specific enthalpy of  the working fluid at the condenser outlet

> 

h[3, L] := HSTXwf(`T__3,L`, 0, wf[L])

237719.8158*Units:-Unit(J/kg)

(4.9)

Specific enthalpy of  the working fluid at the condenser outlet

> 

`s__3,L` := SSTXwf(`T__3,L`, 0, wf[L])

1129.808030*Units:-Unit(J/(kg*K))

(4.10)

Specific volume of  the working fluid at the condenser outlet

> 

`v__3,L` := Units:-Standard:-`/`(DSTXwf(`T__3,L`, 0, wf[L]))

0.9388227496e-3*Units:-Unit(m^3/kg)

(4.11)

Pressure of  the working fluid at the evaporator intlet

> 

p[4, L] := p[1, L]

534123.7895*Units:-Unit(Pa)

(4.12)

Specific enthalpy of  the working fluid at the evaporator intlet

> 

h[4, L] := h[3, L]

237719.8158*Units:-Unit(J/kg)

(4.13)

Temperature of  the working fluid at the evaporator intlet

> 

`T__4,L` := TPHwf(p[4, L], h[4, L], wf[L])

261.0766419*Units:-Unit(K)

(4.14)

Specific enthalpy of saturated liquid of  the working fluid at the evaporator inlet

> 

h[sl, 4, L] := HSPXwf(p[4, L], 0, wf[L])

181961.1572*Units:-Unit(J/kg)

(4.15)

Quality of the working fluid at the evaporator intlet

> 

x[4, L] := Units:-Standard:-`*`(Units:-Standard:-`+`(h[4, L], Units:-Standard:-`-`(h[sl, 4, L])), Units:-Standard:-`/`(Units:-Standard:-`+`(h[1, L], Units:-Standard:-`-`(h[sl, 4, L]))))

 

.2369278076

(4.16)

 

Specific entropy of saturated liquid of  the working fluid at the evaporator inlet

> 

`s__sl,4,L` := SSPXwf(p[4, L], 0, wf[L])

933.2933459*Units:-Unit(J/(kg*K))

(4.17)

Specific volume of saturated liquid of  the working fluid at the evaporator inlet

> 

`v__sl,4,L` := Units:-Standard:-`/`(DSPXwf(p[4, L], 0, wf[L]))

0.8220990714e-3*Units:-Unit(m^3/kg)

(4.18)

 

Specific entropy of  the working fluid at the evaporator inlet

> 

s[4, L] := Units:-Standard:-`+`(`s__sl,4,L`, Units:-Standard:-`*`(x[4, L], Units:-Standard:-`+`(`s__1,L`, Units:-Standard:-`-`(`s__sl,4,L`))))

1146.823903*Units:-Unit(J/(kg*K))

(4.19)

Specific volume of  the working fluid at the evaporator inlet

> 

v[4, L] := Units:-Standard:-`+`(`v__sl,4,L`, Units:-Standard:-`*`(x[4, L], Units:-Standard:-`+`(`v__1,L`, Units:-Standard:-`-`(`v__sl,4,L`))))

0.1219334720e-1*Units:-Unit(m^3/kg)

(4.20)

 

Specific entropy of  the working fluid at the compressor 1 outlet  after isentropic compression

> 

s[2, L, is] := `s__1,L`

 

1834.540688*Units:-Unit(J/(kg*K))

(4.21)

 

Temperature of  the working fluid at the compressor 1 outlet  after isentropic compression

> 

T[2, L, is] := TPSwf(p[2, L], s[2, L, is], wf[L])

 

313.7563374*Units:-Unit(K)

(4.22)

 

Specific enthalpy of  the working fluid at the compressor 1 outlet  after isentropic compression

> 

h[2, L, is] := HPTwf(p[2, L], T[2, L, is], wf[L])

 

447569.9024*Units:-Unit(J/kg)

(4.23)

 

Enhtalpy change in the compressor  after isentropic compression

> 

l[c1, is] := Units:-Standard:-`+`(h[2, L, is], Units:-Standard:-`-`(h[1, L]))

 

30268.4567*Units:-Unit(J/kg)

(4.24)

 

Enhtalpy change in the compressor 1 after actual compression

> 

l[c1] := Units:-Standard:-`*`(Units:-Standard:-`+`(h[2, L, is], Units:-Standard:-`-`(h[1, L])), Units:-Standard:-`/`(eta[comp]))

 

37835.57088*Units:-Unit(J/kg)

(4.25)

 

Specific enthalpy of  the working fluid at the compressor 1 outlet  after actual compression

> 

h[2, L] := Units:-Standard:-`+`(h[1, L], l[c1])

 

455137.0166*Units:-Unit(J/kg)

(4.26)

 

 

Specific enthalpy of  the working fluid at the compressor 1 outlet  after actual compression

> 

T[2, L] := TPHwf(p[2, L], h[2, L], wf[L])

 

320.3318365*Units:-Unit(K)

(4.27)

 

 

Specific entropy of  the working fluid at the compressor 1 outlet  after actual compression 

> 

s[2, L] := SPTwf(p[2, L], T[2, L], wf[L])

 

1858.410794*Units:-Unit(J/(kg*K))

(4.28)

 

 

Specific volume of  the working fluid at the compressor 1 outlet  after actual compression 

> 

v[2, L] := Units:-Standard:-`/`(DPTwf(p[2, L], T[2, L], wf[L]))

 

0.1878649972e-1*Units:-Unit(m^3/kg)

(4.29)

 

 

Temperature of  the working fluid at the compressor 2 intlet

> 

T[1, H] := Units:-Standard:-`+`(TSPXwf(p[2, L], 1, wf[L]), Units:-Standard:-`-`(`&delta;T`))

 

295.0058134*Units:-Unit(K)

(4.30)

 

Pressure of  the working fluid at the compressor 2 intlet

> 

p[1, H] := PSTXwf(T[1, H], 1, wf[H])

 

605226.3010*Units:-Unit(Pa)

(4.31)

 

 

 

Specific enthalpy of  the working fluid at the compressor 2 intlet 

> 

h[1, H] := HSTXwf(T[1, H], 1, wf[H])

 

410719.0452*Units:-Unit(J/kg)

(4.32)

 

Specific entropy of  the working fluid at the compressor 2 intlet

> 

s[1, H] := SSTXwf(T[1, H], 1, wf[H])

 

1717.354341*Units:-Unit(J/(kg*K))

(4.33)

 

Specific volume of  the working fluid at the compressor 2 intlet

> 

v[1, H] := Units:-Standard:-`/`(DSTXwf(T[1, H], 1, wf[H]))

 

0.3400315959e-1*Units:-Unit(m^3/kg)

(4.34)

 

Temperature of the working fluid at the condenser outlet

> 

T[3, H] := Units:-Standard:-`+`(T__out, `&delta;T__CD`)

345.15*Units:-Unit(K)

(4.35)

Pressure of  the working fluid at the condenser outlet

> 

p[3, H] := PSTXwf(T[3, H], 0, wf[H])

 

2213229.244*Units:-Unit(Pa)

(4.36)

 

Specific enthalpy of  the working fluid at the condenser outlet

> 

h[3, H] := HSTXwf(T[3, H], 0, wf[H])

 

307777.3420*Units:-Unit(J/kg)

(4.37)

 

Specific entropy of  the working fluid at the condenser outlet

> 

s[3, H] := SSTXwf(T[3, H], 0, wf[H])

 

1343.048574*Units:-Unit(J/(kg*K))

(4.38)

 

Specific volume of  the working fluid at the condenser outlet

> 

v[3, H] := Units:-Standard:-`/`(DSTXwf(T[3, H], 0, wf[H]))

 

0.1016510305e-2*Units:-Unit(m^3/kg)

(4.39)

 

 

Pressure of  the working fluid at the compressor 2 outlet

> 

p[2, H] := p[3, H]

 

2213229.244*Units:-Unit(Pa)

(4.40)

Pressure increase in the compressor 2  

> 

pi__2 := Units:-Standard:-`*`(p[2, H], Units:-Standard:-`/`(p[1, H]))

3.656862302

(4.41)

Pressure of  the working fluid at the expansion valve 2 outlet

> 

p[4, H] := p[1, H]

605226.3010*Units:-Unit(Pa)

(4.42)

 

Specific enthalpy of  the working fluid at the expansion valve 2 outlet

> 

h[4, H] := h[3, H]

307777.3420*Units:-Unit(J/kg)

(4.43)

Temperature of  the working fluid at the expansion valve 2 outlet

> 

T[4, H] := TPHwf(p[4, H], h[4, H], wf[H])

295.0058134*Units:-Unit(K)

(4.44)

Specific enthalpy of saturated liquid of  the working fluid at the expansion valve 2 outlet

> 

h[sl, 4, H] := HSPXwf(p[4, H], 0, wf[H])

230083.8526*Units:-Unit(J/kg)

(4.45)

Quality of the working fluid at the expansion valve 2 outlet

> 

x[4, H] := Units:-Standard:-`*`(Units:-Standard:-`+`(h[4, H], Units:-Standard:-`-`(h[sl, 4, H])), Units:-Standard:-`/`(Units:-Standard:-`+`(h[1, H], Units:-Standard:-`-`(h[sl, 4, H]))))

 

.4301126944

(4.46)

 

Specific entropy of saturated liquid of  the working fluid at the expansion valve 2 outlet

> 

s[sl, 4, H] := SSPXwf(p[4, H], 0, wf[H])

1105.043721*Units:-Unit(J/(kg*K))

(4.47)

Specific volume of saturated liquid of  the working fluid at the expansion valve 2 outlet

> 

v[sl, 4, H] := Units:-Standard:-`/`(DSPXwf(p[4, H], 0, wf[H]))

0.8206876245e-3*Units:-Unit(m^3/kg)

(4.48)

 

Specific entropy of  the working fluid at the expansion valve 2 outlet

> 

s[4, H] := Units:-Standard:-`+`(s[sl, 4, H], Units:-Standard:-`*`(x[4, H], Units:-Standard:-`+`(s[1, H], Units:-Standard:-`-`(s[sl, 4, H]))))

1368.406292*Units:-Unit(J/(kg*K))

(4.49)

Specific volume of  the working fluid at the expansion valve 2 outlet

> 

v[4, H] := Units:-Standard:-`+`(v[sl, 4, H], Units:-Standard:-`*`(x[4, H], Units:-Standard:-`+`(v[1, H], Units:-Standard:-`-`(v[sl, 4, H]))))

0.1509289005e-1*Units:-Unit(m^3/kg)

(4.50)

Specific entropy of  the working fluid at the compressor 2 outlet  after isentropic compression

> 

s[2, H, is] := s[1, H]

 

1717.354341*Units:-Unit(J/(kg*K))

(4.51)

 

Temperature of  the working fluid at the compressor 2 outlet  after isentropic compression

> 

T[2, H, is] := TPSwf(p[2, H], s[2, H, is], wf[H])

 

350.3634493*Units:-Unit(K)

(4.52)

 

Specific enthalpy of  the working fluid at the compressor 2 outlet  after isentropic compression

> 

h[2, H, is] := HPTwf(p[2, H], T[2, H, is], wf[H])

 

437028.8950*Units:-Unit(J/kg)

(4.53)

 

Enhtalpy change in the compressor 2 after isentropic compression

> 

l[c2, is] := Units:-Standard:-`+`(h[2, H, is], Units:-Standard:-`-`(h[1, H]))

 

26309.8498*Units:-Unit(J/kg)

(4.54)

 

Enhtalpy change in the compressor 2 after actual compression

> 

l[c2] := Units:-Standard:-`*`(Units:-Standard:-`+`(h[2, H, is], Units:-Standard:-`-`(h[1, H])), Units:-Standard:-`/`(eta[comp]))

 

32887.31225*Units:-Unit(J/kg)

(4.55)

 

Specific enthalpy of  the working fluid at the compressor 2 outlet  after actual compression

> 

h[2, H] := Units:-Standard:-`+`(h[1, H], l[c2])

 

443606.3574*Units:-Unit(J/kg)

(4.56)

 

Specific enthalpy of  the working fluid at the compressor 2 outlet  after actual compression

> 

T[2, H] := TPHwf(p[2, H], h[2, H], wf[H])

 

354.9381885*Units:-Unit(K)

(4.57)

 

Specific entropy of  the working fluid at the compressor 2 outlet  after actual compression 

> 

s[2, H] := SPTwf(p[2, H], T[2, H], wf[H])

 

1736.007363*Units:-Unit(J/(kg*K))

(4.58)

 

Specific volume of  the working fluid at the compressor 2 outlet  after actual compression 

> 

v[2, H] := Units:-Standard:-`/`(DPTwf(p[2, H], T[2, H], wf[H]))

 

0.9158027894e-2*Units:-Unit(m^3/kg)

(4.59)

 

Ratio of mass flow rates between the high pressure circuit and that of the low pressure circuit

> 

phi := Units:-Standard:-`*`(Units:-Standard:-`+`(h[2, L], Units:-Standard:-`-`(h[3, L])), Units:-Standard:-`/`(Units:-Standard:-`+`(h[1, H], Units:-Standard:-`-`(h[4, H]))))

 

2.112042001

(4.60)

Functions for plotting Tq-diagram in the condencer/evaporator

> 

plots:-display, "expecting plot structure but received: %1", QT[H]

 

 

 

 

 

NULL

Heat rejection in the condenser referred to 1 kg of refrigerant in the low pressure circuit

> 

q[out] := Units:-Standard:-`*`(Units:-Standard:-`+`(h[2, H], Units:-Standard:-`-`(h[3, H])), phi)

286876.5855*Units:-Unit(J/kg)

(4.61)

Heat addition in the evaporator referred to 1 kg of refrigerant in the low pressure circuit

> 

`#msub(mi("q",fontstyle = "normal"),mo("in"))` := Units:-Standard:-`+`(h[1, L], Units:-Standard:-`-`(h[4, L]))

179581.6299*Units:-Unit(J/kg)

(4.62)

The total work of compressors referred to 1 kg of the working fluid in the low pressure circuit

> 

l[c] := Units:-Standard:-`+`(l[c1], Units:-Standard:-`*`(l[c2], phi))

107294.9556*Units:-Unit(J/kg)

(4.63)

Coefficent of performance of a refrigerator

> 

for X1 from x[4, H] by 0.1e-1 to 1 do if Units:-Standard:-`<`(Units:-Standard:-`+`(Tcdevh(X1, Units:-Standard:-`*`(p[4, H], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[3, H], Units:-Standard:-`/`(Unit('K'))), Units:-Standard:-`*`(h[2, L], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`('kg'))))), Units:-Standard:-`*`(h[1, H], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`('kg'))))), Units:-Standard:-`*`(p[2, L], Units:-Standard:-`/`(Unit('Pa')))), Units:-Standard:-`-`(Tcdevc(X1, Units:-Standard:-`*`(p[4, H], Units:-Standard:-`/`(Unit('Pa')))))), 0) then error "The temperature of hot stream is lower the temperature of cold stream in the intermediate condenser/evaporator" else COP[R] := Units:-Standard:-`*`(`#msub(mi("q",fontstyle = "normal"),mo("in"))`, Units:-Standard:-`/`(l[c])) end if end do; -1; COP[R]

 

1.673719225

(4.64)

 

Coefficent of performance of a heat pump

> 

for X1 from x[4, H] by 0.1e-1 to 1 do if Units:-Standard:-`<`(Units:-Standard:-`+`(Tcdevh(X1, Units:-Standard:-`*`(p[4, H], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[3, H], Units:-Standard:-`/`(Unit('K'))), Units:-Standard:-`*`(h[2, L], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`('kg'))))), Units:-Standard:-`*`(h[1, H], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`('kg'))))), Units:-Standard:-`*`(p[2, L], Units:-Standard:-`/`(Unit('Pa')))), Units:-Standard:-`-`(Tcdevc(X1, Units:-Standard:-`*`(p[4, H], Units:-Standard:-`/`(Unit('Pa')))))), 0) then error "The temperature of hot stream is lower the temperature of cold stream in the intermediate condenser/evaporator" else COP[HP] := Units:-Standard:-`*`(q[out], Units:-Standard:-`/`(l[c])) end if end do; -1; COP[HP]

2.673719225

(4.65)

Functions for plotting Ts-diagram

 

Ploting the Refrigeration Cycle on a Ts-diagram

 

plots:-display(Ts[sl, L], Ts[sv, L], Ts[1, Units:-Standard:-`+`(L, -2), L], Ts[2, Units:-Standard:-`+`(L, Units:-Standard:-`-`(sv)), 2, L], Ts[sv, 2, Units:-Standard:-`+`(L, -3), L], Ts[3, Units:-Standard:-`+`(L, -4), L], Ts[4, Units:-Standard:-`+`(L, -1), L], Ts[sl, H], Ts[sv, H], Ts[1, Units:-Standard:-`+`(H, -2), H], Ts[2, Units:-Standard:-`+`(H, Units:-Standard:-`-`(sv)), 2, H], Ts[sv, 2, Units:-Standard:-`+`(H, -3), H], Ts[3, Units:-Standard:-`+`(H, -4), H], Ts[4, Units:-Standard:-`+`(H, -1), H], labels = [Units:-Standard:-`*`(s, [Units:-Standard:-`*`(kJ, Units:-Standard:-`/`(Units:-Standard:-`*`(kg, K)))]), Units:-Standard:-`*`(T, [K])], size = [800, 500], gridlines)

 

Functions for plotting ph-diagram

 

Ploting the Refrigeration Cycle on a ph-diagram

 
> 

plots:-display(ph[sl, L], ph[sv, L], ph[1, Units:-Standard:-`+`(L, -2), L], ph[2, Units:-Standard:-`+`(L, Units:-Standard:-`-`(sv)), 2, L], ph[sv, 2, Units:-Standard:-`+`(L, -3), L], ph[3, Units:-Standard:-`+`(L, -4), L], ph[4, Units:-Standard:-`+`(L, -1), L], ph[sl, H], ph[sv, H], ph[1, Units:-Standard:-`+`(H, -2), H], ph[2, Units:-Standard:-`+`(H, Units:-Standard:-`-`(sv)), 2, H], ph[sv, 2, Units:-Standard:-`+`(H, -3), H], ph[3, Units:-Standard:-`+`(H, -4), H], ph[4, Units:-Standard:-`+`(H, -1), H], labels = [Units:-Standard:-`*`(h, [Units:-Standard:-`*`(kJ, Units:-Standard:-`/`(Units:-Standard:-`*`(kg, K)))]), Units:-Standard:-`*`(pressure, [MPa])], size = [800, 500], gridlines)

 

Functions for plotting hs-diagram

 

Ploting the Refrigeration Cycle on a hs-diagram

 
> 

plots:-display(hs[sl, L], hs[sv, L], hs[1, Units:-Standard:-`+`(L, -2), L], hs[2, Units:-Standard:-`+`(L, Units:-Standard:-`-`(sv)), 2, L], hs[sv, 2, Units:-Standard:-`+`(L, -3), L], hs[3, Units:-Standard:-`+`(L, -4), L], hs[4, Units:-Standard:-`+`(L, -1), L], hs[sl, H], hs[sv, H], hs[1, Units:-Standard:-`+`(H, -2), H], hs[2, Units:-Standard:-`+`(H, Units:-Standard:-`-`(sv)), 2, H], hs[sv, 2, Units:-Standard:-`+`(H, -3), H], hs[3, Units:-Standard:-`+`(H, -4), H], hs[4, Units:-Standard:-`+`(H, -1), H], labels = [Units:-Standard:-`*`(s, [Units:-Standard:-`*`(kJ, Units:-Standard:-`/`(Units:-Standard:-`*`(kg, K)))]), Units:-Standard:-`*`(h, [Units:-Standard:-`*`(kJ, Units:-Standard:-`/`(kg))])], size = [800, 500], gridlines)

 

Functions for plotting pv-diagram

 

Ploting the Refrigeration Cycle on a pv-diagram

 
> 

plots:-display(pv[sl, L], pv[sv, L], pv[1, Units:-Standard:-`+`(L, -2), L], pv[2, Units:-Standard:-`+`(L, Units:-Standard:-`-`(sv)), 2, L], pv[sv, 2, Units:-Standard:-`+`(L, -3), L], pv[3, Units:-Standard:-`+`(L, -4), L], pv[4, Units:-Standard:-`+`(L, -1), L], pv[sl, H], pv[sv, H], pv[1, Units:-Standard:-`+`(H, -2), H], pv[2, Units:-Standard:-`+`(H, Units:-Standard:-`-`(sv)), 2, H], pv[sv, 2, Units:-Standard:-`+`(H, -3), H], pv[3, Units:-Standard:-`+`(H, -4), H], pv[4, Units:-Standard:-`+`(H, -1), H], labels = [Units:-Standard:-`*`(v, [Units:-Standard:-`*`(m3, Units:-Standard:-`/`(kg))]), Units:-Standard:-`*`(p, [MPa])], size = [800, 500], gridlines)

 
> 

``