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)

 

``