Application Center - Maplesoft
  Submit your work
Application Center Applications Advanced Exergetic Analysis of Simple Vapor Compression Refrigeration Cycle Preview

App Preview:

Advanced Exergetic Analysis of Simple Vapor Compression Refrigeration Cycle

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

Learn about Maple
Download Application



Advanced Exergetic Analysis of Simple Vapor Compression Refrigeration Cycle

Introduction
 

Exergy analysis is without a doubt a powerful tool for developing, evaluating and improving thermal systems, particularly when this analysis is part of an exergoeconomic and exergoenvironmental analysis.

The exergetic analysis identifies the location, the magnitude and the causes of thermodynamic inefficiencies, assists in better understanding the operations in energy conversion processes, and enhances the creativity of engineers dealing with the improvement of such systems and processes.

Today it is not sufficient any more to just present the results of an exergy analysis without using them to improve the system being considered. A so called conventional exergetic analysis  cannot satisfactorily evaluate the mutual interdependencies among the system components, and does not consider the real potential for improving the system components [1].

An advanced exergetic analysis converts the most important limitations of the conventional analysis into strengths. An advanced analysis provides options to estimate (a) the real potential for improving the components (splitting the exergy destruction into unavoidable and avoidable parts); and (b) the mutual interdependencies among the exergy destructions within the components of a system (splitting the exergy destruction into endogenous and exogenous parts). Distinctions between avoidable and unavoidable exergy destruction on one side and endogenous and exogenous exergy destruction on the other side (with a further splitting of the exogenous exergy destruction) allow engineers to focus on the thermodynamic inefficiencies that can be avoided by simultaneously considering the interactions among system components. The avoidable endogenous and the avoidable exogenous exergy destructions provide the best guidance for improving the thermodynamic performance of energy conversion systems.

The drawbacks associated with the advanced exergetic analysis are the need for simulating some non-standard cycles (e.g., so-called hybrid cycles), the large number of required simulations and the associated engineering time.

Within this Maple-aided analysis, the potential for improving the thermodynamic efficiency of components and the overall vapor-compression refrigeration system as well as the interactions among components and their effect on the exergy destruction within each component can be estimated.

Reference:

1. 

Tsatsaronis G., 1999. Strengths and limitations of exergy analysis”. A. Bejan and E.Mamut, eds. “Thermodynamic Optimization of Complex Energy Systems”, Kluwer Academic Publishers, pp. 93-100.
2. 

Morosuk T., Tsatsaronis G. Strengths and limitations of advanced exergetic analyses. Proceedings of the ASME 2013 International Mechanical Engineering Congress & Exposition IMECE2013, November 15-21, 2013, San Diego, California
3. 

Morosuk T, Tsatsaronis G. Splitting the exergy destruction into endogenous and exogenous parts application to refrigeration machines. In: Frangopoulos C, Rakopoulos C, Tsatsaronis G, editors. Proceedings of the 19th international conference on efficiency, cost, optimization, simulation and environmental impact of energy systems, July 12–14; 2006, National Technical University of Athens, Greece, vol. 1; 2006. p. 165–172.
4. 

Morosuk T, Tsatsaronis G. The ‘''Cycle Method’'' used in the exergy analysis of refrigeration machines: from education to research. In: Frangopoulos C, Rakopoulos C, Tsatsaronis G, editors. Proceedings of the 19th international conference on efficiency, cost, optimization, simulation and environmental impact of energy systems, July 12–14; 2006, National Technical University of Athens, Greece, vol. 1; 2006. p. 157–163.
5. 

Kelly S, Tsatsaronis G, Morosuk T. Advanced exergetic analysis: approaches for splitting the exergy destruction into endogenous and exogenous parts. Energy Int J, 2009;34:384–391.
6. 

Morosuk T., Tsatsaronis G. Advanced Exergetic Evaluation of Refrigeration Machines Using Different Working Fluids // Energy –The International Journal. – 2009. – Vol. 34 – P. 2248-2258.
7. 

Morosuk T., Tsatsaronis G. New approach to the exergy analysis of absorption refrigeration machines // Energy – The International Journal. – 2008. – Vol. 33 – P. 890-907.

Fig. 1. Scheme of a Simple Vapor Compression Refrigeration 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

``

> 

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

 

Critical pressure as a function of working fluid

``

> 

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

Initial Data
 

Type of the machine (in case of heat pump set Type as 1, in case of refrigerator set Type as 2)

> 

Type := 2

2

 (3.1)

Capacity of the product

> 

Q[P] := Units:-Standard:-`*`(100, kW)

100*kW

 (3.2)

Temperature of heat source fluid at evaporator inlet

> 

T[EV, h, inlet] := Units:-Standard:-`*`(Units:-Standard:-`+`(-5, 273.15), Unit(''K''))

268.15*Units:-Unit(K)

 (3.3)

Pressure of heat source fluid at evaporator inlet

> 

p[EV, h, inlet] := Units:-Standard:-`*`(100000, Units:-Standard:-Unit(''Pa''))

100000*Units:-Unit(Pa)

 (3.4)

Temperature of heat source fluid at evaporator outlet

> 

T[EV, h, out] := Units:-Standard:-`*`(Units:-Standard:-`+`(-15, 273.15), Units:-Standard:-Unit(''K''))

258.15*Units:-Unit(K)

 (3.5)

Pressure of heat source fluid at evaporator outlet

> 

p[EV, h, out] := Units:-Standard:-`*`(100000, Units:-Standard:-Unit(''Pa''))

100000*Units:-Unit(Pa)

 (3.6)

 

Temperature of heat  consumption fluid at condenser inlet

> 

T[CD, c, inlet] := Units:-Standard:-`*`(Units:-Standard:-`+`(20, 273.15), Units:-Standard:-Unit(''K''))

293.15*Units:-Unit(K)

 (3.7)

Pressure of heat  consumption fluid at condenser inlet

> 

p[CD, c, inlet] := Units:-Standard:-`*`(100000, Units:-Standard:-Unit(''Pa''))

100000*Units:-Unit(Pa)

 (3.8)

 

Temperature of heat consumption fluid at condenser outlet

> 

T[CD, c, out] := Units:-Standard:-`*`(Units:-Standard:-`+`(30, 273.15), Units:-Standard:-Unit(''K''))

303.15*Units:-Unit(K)

 (3.9)

Pressure of heat consumption fluid at condenser outlet

> 

p[CD, c, out] := Units:-Standard:-`*`(100000, Units:-Standard:-Unit(''Pa''))

100000*Units:-Unit(Pa)

 (3.10)

Temperature difference at the cold end of the evaporator

> 

`&delta;T`[EV, c] := Units:-Standard:-`*`(10, Units:-Standard:-Unit(''K''))

10*Units:-Unit(K)

 (3.11)

Temperature difference at the cold end of the evaporator which is possible due to technological limitations

> 

`&delta;T`[EV, c, UN] := Units:-Standard:-`*`(.5, Units:-Standard:-Unit(''K''))

.5*Units:-Unit(K)

 (3.12)

 

Temperature difference at the hot end of the condenser

> 

`&delta;T`[CD, h] := Units:-Standard:-`*`(10, Units:-Standard:-Unit(''K''))

10*Units:-Unit(K)

 (3.13)

Temperature difference at the hot end of the condenser which is possible due to technological limitations

> 

`&delta;T`[CD, h, UN] := Units:-Standard:-`*`(.5, Units:-Standard:-Unit('K'))

.5*Units:-Unit(K)

 (3.14)

 

Isentropic efficiency of the compressor

> 

eta[comp] := .8

.8

 (3.15)

Max. isentropic efficiency of compressor which is possible due to technological limitations

> 

eta[comp, UN] := .95

.95

 (3.16)

 

Working fluid

> 

wf[cycle] := R717

R717

 (3.17)

Reference temperature for exergy analysis

> 

T[0] := Units:-Standard:-`*`(Units:-Standard:-`+`(20, 273.15), Units:-Standard:-Unit('K'))

293.15*Units:-Unit(K)

 (3.18)

Reference pressure for exergy analysis

> 

p[0] := Units:-Standard:-`*`(100000, Units:-Standard:-Unit('Pa'))

100000*Units:-Unit(Pa)

 (3.19)

Working fluid of the heat source

> 

wf[EV, h] := Air

Air

 (3.20)

Working fluid of the heat consumption

> 

wf[CD, c] := Air

Air

 (3.21)

Calculation
 

Specific enthalpy of working fluid at reference state

> 

h[0] := HPTwf(p[0], T[0], wf[cycle])

1679027.299*Units:-Unit(J/kg)

 (4.1)

Specific entropy of working fluid at reference state

> 

s[0] := SPTwf(p[0], T[0], wf[cycle])

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

 (4.2)

Specific enthalpy of heat source fluid at reference state

> 

h[0, EV, h] := HPTwf(p[0], T[0], wf[EV, h])

419408.0705*Units:-Unit(J/kg)

 (4.3)

Specific entropy of the heat source fluid at reference state

> 

s[0, EV, h] := SPTwf(p[0], T[0], wf[EV, h])

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

 (4.4)

Specific enthalpy of the heat consumption fluid at reference state

> 

h[0, CD, c] := HPTwf(p[0], T[0], wf[CD, c])

419408.0705*Units:-Unit(J/kg)

 (4.5)

Specific entropy of the heat consumption fluid at reference state

> 

s[0, CD, c] := SPTwf(p[0], T[0], wf[CD, c])

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

 (4.6)

Specific enthalpy of heat source fluid at evaporator inlet

> 

h[EV, h, inlet] := HPTwf(p[EV, h, inlet], T[EV, h, inlet], wf[EV, h])

394262.6750*Units:-Unit(J/kg)

 (4.7)

Specific enthalpy of heat source fluid at evaporator inlet

> 

s[EV, h, inlet] := SPTwf(p[EV, h, inlet], T[EV, h, inlet], wf[EV, h])

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

 (4.8)

Specific enthalpy of heat source fluid at evaporator outlet

> 

h[EV, h, out] := HPTwf(p[EV, h, out], T[EV, h, out], wf[EV, h])

384207.2149*Units:-Unit(J/kg)

 (4.9)

Specific enthalpy of heat source fluid at evaporator outlet

> 

s[EV, h, out] := SPTwf(p[EV, h, out], T[EV, h, out], wf[EV, h])

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

 (4.10)

Specific enthalpy of heat consumption fluid at condenser inlet

> 

h[CD, c, inlet] := HPTwf(p[CD, c, inlet], T[CD, c, inlet], wf[CD, c])

419408.0705*Units:-Unit(J/kg)

 (4.11)

Specific entropy of heat consumption fluid at condenser inlet

> 

s[CD, c, inlet] := SPTwf(p[CD, c, inlet], T[CD, c, inlet], wf[CD, c])

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

 (4.12)

Specific enthalpy of heat consumption fluid at condenser outlet

> 

h[CD, c, out] := HPTwf(p[CD, c, out], T[CD, c, out], wf[CD, c])

429470.9729*Units:-Unit(J/kg)

 (4.13)

Specific entropy of heat consumption fluid at condenser inlet

> 

s[CD, c, out] := SPTwf(p[CD, c, out], T[CD, c, out], wf[CD, c])

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

 (4.14)

Specific physical exergy of heat source fluid at evaporator inlet

> 

e[EV, h, inlet] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[EV, h, inlet], Units:-Standard:-`-`(h[0, EV, h])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[EV, h, inlet], Units:-Standard:-`-`(s[0, EV, h])))))

1137.24829*Units:-Unit(J/kg)

 (4.15)

Specific physical exergy of heat source fluid at evaporator outlet

> 

e[EV, h, out] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[EV, h, out], Units:-Standard:-`-`(h[0, EV, h])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[EV, h, out], Units:-Standard:-`-`(s[0, EV, h])))))

2284.95077*Units:-Unit(J/kg)

 (4.16)

Specific physical exergy of heat source fluid at condenser inlet

> 

e[CD, c, inlet] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[CD, c, inlet], Units:-Standard:-`-`(h[0, CD, c])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[CD, c, inlet], Units:-Standard:-`-`(s[0, CD, c])))))

-0.

 (4.17)

Specific physical exergy of heat source fluid at condenser outlet

> 

e[CD, c, out] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[CD, c, out], Units:-Standard:-`-`(h[0, CD, c])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[CD, c, out], Units:-Standard:-`-`(s[0, CD, c])))))

167.837563*Units:-Unit(J/kg)

 (4.18)

1. Real cycle

Temperature of the working fluid at the compressor outlet

> 

T[4, R] := Units:-Standard:-`+`(T[EV, h, out], Units:-Standard:-`-`(`&delta;T`[EV, c]))

248.15*Units:-Unit(K)

 (4.19)

Pressure of the working fluid at the compressor outlet

> 

p[2, R] := PSTXwf(Units:-Standard:-`+`(T[CD, c, out], `&delta;T`[CD, h]), 1, wf[cycle])

 

1555421.989*Units:-Unit(Pa)

 (4.20)

 

Pressure of  the working fluid at the condenser outlet

> 

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

 

1555421.989*Units:-Unit(Pa)

 (4.21)

 

Pressure of  the working fluid at the condenser outlet

> 

T[3, R] := TSPXwf(p[3, R], 0, wf[cycle])

 

313.1500000*Units:-Unit(K)

 (4.22)

 

 

Specific enthalpy of  the working fluid at the condenser outlet

> 

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

 

533791.2205*Units:-Unit(J/kg)

 (4.23)

 

Specific entropy of  the working fluid at the condenser outlet

> 

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

 

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

 (4.24)

 

Specific volume of  the working fluid at the condenser outlet

> 

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

 

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

 (4.25)

 

Specific exergy of  the working fluid at the condenser outlet

> 

e[3, R] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, R], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, R], Units:-Standard:-`-`(s[0])))))

 

299157.094*Units:-Unit(J/kg)

 (4.26)

 

Auxiliary values

> 

h[3, R, 1] := HPTwf(p[3, R], T[0], wf[cycle])

 

437261.1314*Units:-Unit(J/kg)

 (4.27)

 

> 

s[3, R, 1] := SPTwf(p[3, R], T[0], wf[cycle])

 

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

 (4.28)

 

Specific mechanical exergy of  the working fluid at the condenser outlet

> 

e[3, R, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[3, R, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[3, R, 1]))))

 

295986.156*Units:-Unit(J/kg)

 (4.29)

 

Specific thermal exergy of  the working fluid at the condenser outlet

> 

e[3, R, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, R], Units:-Standard:-`-`(h[3, R, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, R], Units:-Standard:-`-`(s[3, R, 1])))))

 

3170.93781*Units:-Unit(J/kg)

 (4.30)

 

 

Specific enthalpy of the working fluid at the evaporator intlet

> 

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

 

533791.2205*Units:-Unit(J/kg)

 (4.31)

 

Pressure of the working fluid at the evaporator intlet

 

> 

p[4, R] := Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`*`(PSTXwf(T[4, R], 0, wf[cycle]), Units:-Standard:-`+`(HSTXwf(T[4, R], 1, wf[cycle]), Units:-Standard:-`-`(h[4, R]))), Units:-Standard:-`*`(PSTXwf(T[4, R], 1, wf[cycle]), Units:-Standard:-`+`(h[4, R], Units:-Standard:-`-`(HSTXwf(T[4, R], 0, wf[cycle]))))), Units:-Standard:-`/`(Units:-Standard:-`+`(HSTXwf(T[4, R], 1, wf[cycle]), Units:-Standard:-`-`(HSTXwf(T[4, R], 0, wf[cycle])))))

 

151471.2936*Units:-Unit(Pa)

 (4.32)

 

Specific enthalpy of saturated liquid of  the working fluid in the evaporator

> 

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

229167.5145*Units:-Unit(J/kg)

 (4.33)

Specific enthalpy of saturated vapor of  the working fluid in the evaporator

> 

h[ss, 4, R] := HSPXwf(p[4, R], 1, wf[cycle])

1573806.240*Units:-Unit(J/kg)

 (4.34)

Quality of the working fluid at the evaporator intlet

> 

x[4, R] := Units:-Standard:-`*`(Units:-Standard:-`+`(h[4, R], Units:-Standard:-`-`(h[sl, 4, R])), Units:-Standard:-`/`(Units:-Standard:-`+`(h[ss, 4, R], Units:-Standard:-`-`(h[sl, 4, R]))))NULL

 

.2265468784

 (4.35)

 

Specific entropy of saturated liquid of  the working fluid in the evaporator

> 

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

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

 (4.36)

Specific entropy of saturated vapor of  the working fluid in the evaporator

> 

s[sv, 4, R] := SSPXwf(p[4, R], 1, wf[cycle])

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

 (4.37)

Specific entropy of  the working fluid at the evaporator inlet

> 

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

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

 (4.38)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

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

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

 (4.39)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

v[sv, 4, R] := Units:-Standard:-`/`(DSPXwf(p[4, R], 1, wf[cycle]))

.7716715235*Units:-Unit(m^3/kg)

 (4.40)

Specific volume of  the working fluid at the evaporator inlet

> 

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

.1759715514*Units:-Unit(m^3/kg)

 (4.41)

Specific exergy of  the working fluid at the evaporator inlet

> 

e[4, R] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, R], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, R], Units:-Standard:-`-`(s[0])))))

 

256035.355*Units:-Unit(J/kg)

 (4.42)

 

Auxiliary values

> 

h[4, R, 1] := HPTwf(p[4, R], T[0], wf[cycle])

 

1675673.399*Units:-Unit(J/kg)

 (4.43)

 

> 

s[4, R, 1] := SPTwf(p[4, R], T[0], wf[cycle])

 

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

 (4.44)

 

Specific mechanical exergy of  the working fluid at the evaporator inlet

> 

e[4, R, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[4, R, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[4, R, 1]))))

 

58604.35859*Units:-Unit(J/kg)

 (4.45)

 

Specific thermal exergy of  the working fluid at the evaporator inlet

> 

e[4, R, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, R], Units:-Standard:-`-`(h[4, R, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, R], Units:-Standard:-`-`(s[4, R, 1])))))

 

197430.996*Units:-Unit(J/kg)

 (4.46)

 

Pressure of  the working fluid at the compressor inlet

> 

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

151471.2936*Units:-Unit(Pa)

 (4.47)

Temperature of  the working fluid at the compressor inlet

> 

T[1, R] := TSPXwf(p[1, R], 1, wf[cycle])

248.1500000*Units:-Unit(K)

 (4.48)

Specific enthalpy of  the working fluid at the compressor inlet

> 

h[1, R] := HSPXwf(p[1, R], 1, wf[cycle])

1573806.240*Units:-Unit(J/kg)

 (4.49)

Specific entropy of  the working fluid at the compressor inlet

> 

s[1, R] := SSPXwf(p[1, R], 1, wf[cycle])

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

 (4.50)

Specific volume of  the working fluid at the compressor inlet

> 

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

 

.7716715235*Units:-Unit(m^3/kg)

 (4.51)

 

Specific exergy of  the working fluid at the compressor inlet

> 

e[1, R] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[1, R], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[1, R], Units:-Standard:-`-`(s[0])))))

 

67437.0235*Units:-Unit(J/kg)

 (4.52)

 

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

> 

h[2, R] := HcompPTPEFFwf(p[1, R], T[1, R], p[2, R], eta[comp], wf[cycle])

2022472.774*Units:-Unit(J/kg)

 (4.53)

Enhtalpy change in the compressor 1 after actual compression

> 

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

 

448666.534*Units:-Unit(J/kg)

 (4.54)

 

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

> 

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

 

456.9585107*Units:-Unit(K)

 (4.55)

 

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

> 

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

 

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

 (4.56)

 

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

> 

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

.1386479073*Units:-Unit(m^3/kg)

 (4.57)

Specific exergy of  the working fluid at the compressor inlet

> 

e[2, R] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[2, R], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[2, R], Units:-Standard:-`-`(s[0])))))

 

456185.3925*Units:-Unit(J/kg)

 (4.58)

 

 

Specific heat rejection in the condenser

> 

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

1488681.554*Units:-Unit(J/kg)

 (4.59)

Specific heat addition in the evaporator

> 

q[inlet, R] := Units:-Standard:-`+`(h[1, R], Units:-Standard:-`-`(h[4, R]))

1040015.020*Units:-Unit(J/kg)

 (4.60)

Coefficent of performance of a refrigerator

> 

COP[R, Ref] := Units:-Standard:-`*`(q[inlet, R], Units:-Standard:-`/`(l[comp, R]))

2.318013360

 (4.61)

Coefficent of performance of a heat pump

> 

COP[R, HP] := Units:-Standard:-`*`(q[out, R], Units:-Standard:-`/`(l[comp, R]))

3.318013361

 (4.62)

Mass rate of the working fluid

> 

m[wf, cycle, R] := proc (Type) if Type = 1 then Q[P]/q[out, R] else Q[P]/q[inlet, R] end if end proc:
> 

m[wf, cycle, R] := m[wf, cycle, R](Type)

0.9615245749e-4*kW*Units:-Unit(kg/J)

 (4.63)

Total heat addition in the evaporator

> 

Q[inlet, R] := Units:-Standard:-`*`(m[wf, cycle, R], q[inlet, R])

 

100.0000000*kW

 (4.64)

 

Total heat rejection in the condenser

> 

Q[out, R] := Units:-Standard:-`*`(m[wf, cycle, R], q[out, R])

 

143.1403898*kW

 (4.65)

 

Mass rate of the heat source fluid

> 

m[EV, h, R] := Units:-Standard:-`*`(Q[inlet, R], Units:-Standard:-`/`(Units:-Standard:-`+`(h[EV, h, inlet], Units:-Standard:-`-`(h[EV, h, out]))))

 

0.9944845786e-2*kW*Units:-Unit(kg/J)

 (4.66)

 

Mass rate of the heat consumption fluid

> 

m[CD, c, R] := Units:-Standard:-`*`(Q[out, R], Units:-Standard:-`/`(Units:-Standard:-`+`(h[CD, c, out], Units:-Standard:-`-`(h[CD, c, inlet]))))

 

0.1422456306e-1*kW*Units:-Unit(kg/J)

 (4.67)

 

Total exergy used in the compressor

> 

E[5, R] := Units:-Standard:-`*`(m[wf, cycle, R], l[comp, R])

 

43.14038984*kW

 (4.68)

 

Total physical exergy at the compressor intlet

> 

E[1, R] := Units:-Standard:-`*`(m[wf, cycle, R], e[1, R])

 

6.484235535*kW

 (4.69)

 

Total physical exergy at the compressor outlet

> 

E[2, R] := Units:-Standard:-`*`(m[wf, cycle, R], e[2, R])

 

43.86334656*kW

 (4.70)

 

Total physical exergy at the condenser outlet

> 

E[3, R] := Units:-Standard:-`*`(m[wf, cycle, R], e[3, R])

 

28.76468976*kW

 (4.71)

 

Total mechanical exergy at the condenser outlet

> 

E[3, R, M] := Units:-Standard:-`*`(m[wf, cycle, R], e[3, R, M])

 

28.45979628*kW

 (4.72)

 

Total thermal exergy at the condenser outlet

> 

E[3, R, T] := Units:-Standard:-`*`(m[wf, cycle, R], e[3, R, T])

 

.3048934630*kW

 (4.73)

 

Total physical exergy at the evaporator inlet

> 

E[4, R] := Units:-Standard:-`*`(m[wf, cycle, R], e[4, R])

 

24.61842859*kW

 (4.74)

 

Total mechanical exergy at the evaporator inlet

> 

E[4, R, M] := Units:-Standard:-`*`(m[wf, cycle, R], e[4, R, M])

 

5.634953098*kW

 (4.75)

 

Total thermal exergy at the evaporator inlet

> 

E[4, R, T] := Units:-Standard:-`*`(m[wf, cycle, R], e[4, R, T])

 

18.98347545*kW

 (4.76)

 

Total physical exergy of heat consumption fluid at the evaporator intlet 

> 

E[EV, h, inlet] := Units:-Standard:-`*`(m[EV, h, R], e[EV, h, inlet])

 

11.30975886*kW

 (4.77)

 

Total physical exergy of heat consumption fluid at the evaporator outlet

> 

E[EV, h, out] := Units:-Standard:-`*`(m[EV, h, R], e[EV, h, out])

 

22.72348304*kW

 (4.78)

 

Total physical exergy of heat consumption fluid at the condenser intlet 

> 

E[CD, c, inlet] := Units:-Standard:-`*`(m[CD, c, R], e[CD, c, inlet])

 

-0.

 (4.79)

 

Total physical exergy of heat consumption fluid at the condenser outlet

> 

E[CD, c, out] := Units:-Standard:-`*`(m[CD, c, R], e[CD, c, out])

 

2.387415999*kW

 (4.80)

 

Exergy of fuel in the compressor

> 

E[F, CM, R] := E[5, R]

 

43.14038984*kW

 (4.81)

 

Exergy of product in the compressor

> 

E[P, CM, R] := Units:-Standard:-`+`(E[2, R], Units:-Standard:-`-`(E[1, R]))

 

37.37911102*kW

 (4.82)

 

Exergy destruction in the compressor

> 

E[D, CM, R] := Units:-Standard:-`+`(E[F, CM, R], Units:-Standard:-`-`(E[P, CM, R]))

 

5.76127882*kW

 (4.83)

 

Exergetic efficiency in the compressor

> 

varepsilon[CM, R] := Units:-Standard:-`*`(E[P, CM, R], Units:-Standard:-`/`(E[F, CM, R]))

 

.8664527873

 (4.84)

 

Exergy of fuel in the condenser

> 

E[F, CD, R] := Units:-Standard:-`+`(E[2, R], Units:-Standard:-`-`(E[3, R]))

 

15.09865680*kW

 (4.85)

 

Exergy of product in the condenser

> 

E[P, CD, R] := Units:-Standard:-`+`(E[CD, c, out], Units:-Standard:-`-`(E[CD, c, inlet]))

 

2.387415999*kW

 (4.86)

 

Exergy destruction in the condenser

> 

E[D, CD, R] := Units:-Standard:-`+`(E[F, CD, R], Units:-Standard:-`-`(E[P, CD, R]))

 

12.71124080*kW

 (4.87)

 

Exergetic efficiency in the condenser

> 

varepsilon[CD, R] := Units:-Standard:-`*`(E[P, CD, R], Units:-Standard:-`/`(E[F, CD, R]))

 

.1581210853

 (4.88)

 

Exergy of fuel in the throttling valve

> 

E[F, TV, R] := Units:-Standard:-`+`(E[3, R, M], Units:-Standard:-`-`(E[4, R, M]))

 

22.82484318*kW

 (4.89)

 

Exergy of product in the throttling valve

> 

E[P, TV, R] := Units:-Standard:-`+`(E[4, R, T], Units:-Standard:-`-`(E[3, R, T]))

 

18.67858199*kW

 (4.90)

 

Exergy destruction in the throttling valve

> 

E[D, TV, R] := Units:-Standard:-`+`(E[F, TV, R], Units:-Standard:-`-`(E[P, TV, R]))

 

4.14626119*kW

 (4.91)

 

Exergetic efficiency in the throttling valve

> 

varepsilon[TV, R] := Units:-Standard:-`*`(E[P, TV, R], Units:-Standard:-`/`(E[F, TV, R]))

 

.8183443734

 (4.92)

 

Exergy of fuel in the evaporator

> 

E[F, EV, R] := Units:-Standard:-`+`(E[4, R], Units:-Standard:-`-`(E[1, R]))

 

18.13419306*kW

 (4.93)

 

Exergy of product in the evaporator

> 

E[P, EV, R] := Units:-Standard:-`+`(E[EV, h, out], Units:-Standard:-`-`(E[EV, h, inlet]))

 

11.41372418*kW

 (4.94)

 

Exergy destruction in the evaporator

> 

E[D, EV, R] := Units:-Standard:-`+`(E[F, EV, R], Units:-Standard:-`-`(E[P, EV, R]))

 

6.72046888*kW

 (4.95)

 

Exergetic efficiency in the evaporator

> 

varepsilon[EV, R] := Units:-Standard:-`*`(E[P, EV, R], Units:-Standard:-`/`(E[F, EV, R]))

 

.6294034779

 (4.96)

 

Exergy of fuel of the overall system

> 

E[F, TOT, R] := E[5, R]

43.14038984*kW

 (4.97)

Exergy of product of the overall system 

> 

E[P, TOT, R] := proc (Type) if Type = 1 then E[P, CD, R] else E[P, EV, R] end if end proc:

 

> 

E[P, TOT, R] := E[P, TOT, R](Type)

 

11.41372418*kW

 (4.98)

 

Exergy destruction of the overall system

> 

E[D, TOT, R] := Units:-Standard:-`+`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CM, R], E[D, CD, R]), E[D, TV, R]), E[D, EV, R])

29.33924969*kW

 (4.99)

Exergetic efficiency of the overall system

> 

varepsilon[TOT, R] := Units:-Standard:-`*`(E[P, TOT, R], Units:-Standard:-`/`(E[F, TOT, R]))

.2645716514

 (4.100)

Exergy destruction ratio in the compressor related to the exergy of total fuel

> 

y[CM, R] := Units:-Standard:-`*`(E[D, CM, R], Units:-Standard:-`/`(E[F, TOT, R]))

 

.1335472128

 (4.101)

 

Exergy destruction ratio in the condenser related to the exergy of total fuel

> 

y[CD, R] := Units:-Standard:-`*`(E[D, CD, R], Units:-Standard:-`/`(E[F, TOT, R]))

.2946482600

 (4.102)

Exergy destruction ratio in the throttling valve related to the exergy of total fuel

> 

y[TV, R] := Units:-Standard:-`*`(E[D, TV, R], Units:-Standard:-`/`(E[F, TOT, R]))

0.9611088832e-1

 (4.103)

Exergy destruction ratio in the evaporator related to the exergy of total fuel

> 

y[EV, R] := Units:-Standard:-`*`(E[D, EV, R], Units:-Standard:-`/`(E[F, TOT, R]))

.1557813665

 (4.104)

Exergy destruction ratio in the system overall related to the exergy of total fuel

> 

y[TOT, R] := Units:-Standard:-`*`(E[D, TOT, R], Units:-Standard:-`/`(E[F, TOT, R]))

 

.6800877275

 (4.105)

 

Exergy destruction ratio in the compressor related to the total exergy destruction

> 

y[CM, R, 1] := Units:-Standard:-`*`(E[D, CM, R], Units:-Standard:-`/`(E[D, TOT, R]))

 

.1963676264

 (4.106)

 

Exergy destruction ratio in the condenser related to the total exergy destruction

> 

y[CD, R, 1] := Units:-Standard:-`*`(E[D, CD, R], Units:-Standard:-`/`(E[D, TOT, R]))

.4332503706

 (4.107)

Exergy destruction ratio in the throttling valve related to the total exergy destruction

> 

y[TV, R, 1] := Units:-Standard:-`*`(E[D, TV, R], Units:-Standard:-`/`(E[D, TOT, R]))

.1413213097

 (4.108)

Exergy destruction ratio in the evaporator related to the total exergy destruction

> 

y[EV, R, 1] := Units:-Standard:-`*`(E[D, EV, R], Units:-Standard:-`/`(E[D, TOT, R]))

.2290606935

 (4.109)

Exergy destruction ratio in the system overall related to the total exergy destruction

> 

y[TOT, R, 1] := Units:-Standard:-`*`(E[D, TOT, R], Units:-Standard:-`/`(E[D, TOT, R]))

 

1.000000000

 (4.110)

 

2. Unavoidable cycle

Temperature of the working fluid at the compressor outlet

> 

T[4, UN] := Units:-Standard:-`+`(T[EV, h, out], Units:-Standard:-`-`(`&delta;T`[EV, c, UN]))

257.65*Units:-Unit(K)

 (4.111)

Pressure of the working fluid at the compressor outlet

> 

p[2, UN] := PSTXwf(Units:-Standard:-`+`(T[CD, c, out], `&delta;T`[CD, h, UN]), 1, wf[cycle])

 

1184647.760*Units:-Unit(Pa)

 (4.112)

 

Pressure of  the working fluid at the condenser outlet

> 

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

 

1184647.760*Units:-Unit(Pa)

 (4.113)

 

Pressure of  the working fluid at the condenser outlet

> 

T[3, UN] := TSPXwf(p[3, UN], 0, wf[cycle])

 

303.6500000*Units:-Unit(K)

 (4.114)

 

 

Specific enthalpy of  the working fluid at the condenser outlet

> 

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

 

487332.5593*Units:-Unit(J/kg)

 (4.115)

 

Specific entropy of  the working fluid at the condenser outlet

> 

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

 

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

 (4.116)

 

Specific volume of  the working fluid at the condenser outlet

> 

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

 

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

 (4.117)

 

Specific exergy of  the working fluid at the condenser outlet

> 

e[3, UN] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, UN], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, UN], Units:-Standard:-`-`(s[0])))))

 

296260.646*Units:-Unit(J/kg)

 (4.118)

 

Auxiliary values

> 

h[3, UN, 1] := HPTwf(p[3, UN], T[0], wf[cycle])

 

437088.3990*Units:-Unit(J/kg)

 (4.119)

 

> 

s[3, UN, 1] := SPTwf(p[3, UN], T[0], wf[cycle])

 

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

 (4.120)

 

Specific mechanical exergy of  the working fluid at the condenser outlet

> 

e[3, UN, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[3, UN, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[3, UN, 1]))))

 

295378.982*Units:-Unit(J/kg)

 (4.121)

 

Specific thermal exergy of  the working fluid at the condenser outlet

> 

e[3, UN, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, UN], Units:-Standard:-`-`(h[3, UN, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, UN], Units:-Standard:-`-`(s[3, UN, 1])))))

 

881.66486*Units:-Unit(J/kg)

 (4.122)

 

 

Specific enthalpy of the working fluid at the evaporator intlet

> 

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

 

487332.5593*Units:-Unit(J/kg)

 (4.123)

 

Pressure of the working fluid at the evaporator intlet

 

> 

p[4, UN] := Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`*`(PSTXwf(T[4, UN], 0, wf[cycle]), Units:-Standard:-`+`(HSTXwf(T[4, UN], 1, wf[cycle]), Units:-Standard:-`-`(h[4, UN]))), Units:-Standard:-`*`(PSTXwf(T[4, UN], 1, wf[cycle]), Units:-Standard:-`+`(h[4, UN], Units:-Standard:-`-`(HSTXwf(T[4, UN], 0, wf[cycle]))))), Units:-Standard:-`/`(Units:-Standard:-`+`(HSTXwf(T[4, UN], 1, wf[cycle]), Units:-Standard:-`-`(HSTXwf(T[4, UN], 0, wf[cycle])))))

 

231194.2952*Units:-Unit(Pa)

 (4.124)

 

Temperature of the working fluid at the compressor outlet

> 

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

257.6500000*Units:-Unit(K)

 (4.125)

Specific enthalpy of saturated liquid of  the working fluid in the evaporator

> 

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

272097.9009*Units:-Unit(J/kg)

 (4.126)

Specific enthalpy of saturated vapor of  the working fluid in the evaporator

> 

h[sv, 4, UN] := HSPXwf(p[4, UN], 1, wf[cycle])

1586868.880*Units:-Unit(J/kg)

 (4.127)

Quality of the working fluid at the evaporator intlet

> 

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

 

.1637050572

 (4.128)

 

Specific entropy of saturated liquid of  the working fluid in the evaporator

> 

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

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

 (4.129)

Specific entropy of saturated vapor of  the working fluid in the evaporator

> 

s[sv, 4, UN] := SSPXwf(p[4, UN], 1, wf[cycle])

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

 (4.130)

Specific entropy of  the working fluid at the evaporator inlet

> 

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

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

 (4.131)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

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

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

 (4.132)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

v[sv, 4, UN] := Units:-Standard:-`/`(DSPXwf(p[4, UN], 1, wf[cycle]))

.5189588087*Units:-Unit(m^3/kg)

 (4.133)

Specific volume of  the working fluid at the evaporator inlet

> 

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

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

 (4.134)

Specific exergy of  the working fluid at the evaporator inlet

> 

e[4, UN] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, UN], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, UN], Units:-Standard:-`-`(s[0])))))

 

274922.645*Units:-Unit(J/kg)

 (4.135)

 

Auxiliary values

> 

h[4, UN, 1] := HPTwf(p[4, UN], T[0], wf[cycle])

 

1670374.872*Units:-Unit(J/kg)

 (4.136)

 

> 

s[4, UN, 1] := SPTwf(p[4, UN], T[0], wf[cycle])

 

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

 (4.137)

 

Specific mechanical exergy of  the working fluid at the evaporator inlet

> 

e[4, UN, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[4, UN, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[4, UN, 1]))))

 

117839.6126*Units:-Unit(J/kg)

 (4.138)

 

Specific thermal exergy of  the working fluid at the evaporator inlet

> 

e[4, UN, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, UN], Units:-Standard:-`-`(h[4, UN, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, UN], Units:-Standard:-`-`(s[4, UN, 1])))))

 

157083.032*Units:-Unit(J/kg)

 (4.139)

 

Pressure of  the working fluid at the compressor inlet

> 

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

231194.2952*Units:-Unit(Pa)

 (4.140)

Temperature of  the working fluid at the compressor inlet

> 

T[1, UN] := TSPXwf(p[1, UN], 1, wf[cycle])

257.6500000*Units:-Unit(K)

 (4.141)

Specific enthalpy of  the working fluid at the compressor inlet

> 

h[1, UN] := HSPXwf(p[1, UN], 1, wf[cycle])

1586868.880*Units:-Unit(J/kg)

 (4.142)

Specific entropy of  the working fluid at the compressor inlet

> 

s[1, UN] := SSPXwf(p[1, UN], 1, wf[cycle])

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

 (4.143)

Specific volume of  the working fluid at the compressor inlet

> 

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

 

.5189588090*Units:-Unit(m^3/kg)

 (4.144)

 

Specific exergy of  the working fluid at the compressor inlet

> 

e[1, UN] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[1, UN], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[1, UN], Units:-Standard:-`-`(s[0])))))

 

123424.3364*Units:-Unit(J/kg)

 (4.145)

 

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

> 

h[2, UN] := HcompPTPEFFwf(p[1, UN], T[1, UN], p[2, UN], eta[comp, UN], wf[cycle])

1836728.847*Units:-Unit(J/kg)

 (4.146)

Enhtalpy change in the compressor 1 after actual compression

> 

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

 

249859.967*Units:-Unit(J/kg)

 (4.147)

 

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

> 

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

 

379.5440924*Units:-Unit(K)

 (4.148)

 

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

> 

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

 

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

 (4.149)

 

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

> 

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

.1484683294*Units:-Unit(m^3/kg)

 (4.150)

Specific exergy of  the working fluid at the compressor inlet

> 

e[2, UN] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[2, UN], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[2, UN], Units:-Standard:-`-`(s[0])))))

 

363571.1677*Units:-Unit(J/kg)

 (4.151)

 

 

Specific heat rejection in the condenser

> 

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

1349396.288*Units:-Unit(J/kg)

 (4.152)

Specific heat addition in the evaporator

> 

q[inlet, UN] := Units:-Standard:-`+`(h[1, UN], Units:-Standard:-`-`(h[4, UN]))

1099536.321*Units:-Unit(J/kg)

 (4.153)

Coefficent of performance of a refrigerator

> 

COP[UN, Ref] := Units:-Standard:-`*`(q[inlet, UN], Units:-Standard:-`/`(l[comp, UN]))

4.400610207

 (4.154)

Coefficent of performance of a heat pump

> 

COP[UN, HP] := Units:-Standard:-`*`(q[out, UN], Units:-Standard:-`/`(l[comp, UN]))

5.400610207

 (4.155)

Mass rate of the working fluid

> 

m[wf, cycle, UN] := proc (Type) if Type = 1 then Q[P]/q[out, UN] else Q[P]/q[inlet, UN] end if end proc:

 

> 

m[wf, cycle, UN] := m[wf, cycle, UN](Type)

 

0.9094742765e-4*kW*Units:-Unit(kg/J)

 (4.156)

 

Total heat addition in the evaporator

> 

Q[inlet, UN] := Units:-Standard:-`*`(m[wf, cycle, UN], q[inlet, UN])

 

100.0000000*kW

 (4.157)

 

Total heat rejection in the condenser

> 

Q[out, UN] := Units:-Standard:-`*`(m[wf, cycle, UN], q[out, UN])

 

122.7241213*kW

 (4.158)

 

Mass rate of the heat source fluid

> 

m[EV, h, UN] := Units:-Standard:-`*`(Q[inlet, UN], Units:-Standard:-`/`(Units:-Standard:-`+`(h[EV, h, inlet], Units:-Standard:-`-`(h[EV, h, out]))))

 

0.9944845786e-2*kW*Units:-Unit(kg/J)

 (4.159)

 

Mass rate of the heat consumption fluid

> 

m[CD, c, UN] := Units:-Standard:-`*`(Q[out, UN], Units:-Standard:-`/`(Units:-Standard:-`+`(h[CD, c, out], Units:-Standard:-`-`(h[CD, c, inlet]))))

 

0.1219569826e-1*kW*Units:-Unit(kg/J)

 (4.160)

 

Total exergy used in the compressor

> 

E[5, UN] := Units:-Standard:-`*`(m[wf, cycle, UN], l[comp, UN])

 

22.72412127*kW

 (4.161)

 

Total physical exergy at the compressor intlet

> 

E[1, UN] := Units:-Standard:-`*`(m[wf, cycle, UN], e[1, UN])

 

11.22512590*kW

 (4.162)

 

Total physical exergy at the compressor outlet

> 

E[2, UN] := Units:-Standard:-`*`(m[wf, cycle, UN], e[2, UN])

 

33.06586247*kW

 (4.163)

 

Total physical exergy at the condenser outlet

> 

E[3, UN] := Units:-Standard:-`*`(m[wf, cycle, UN], e[3, UN])

 

26.94414367*kW

 (4.164)

 

Total mechanical exergy at the condenser outlet

> 

E[3, UN, M] := Units:-Standard:-`*`(m[wf, cycle, UN], e[3, UN, M])

 

26.86395859*kW

 (4.165)

 

Total thermal exergy at the condenser outlet

> 

E[3, UN, T] := Units:-Standard:-`*`(m[wf, cycle, UN], e[3, UN, T])

 

0.8018515107e-1*kW

 (4.166)

 

Total physical exergy at the evaporator inlet

> 

E[4, UN] := Units:-Standard:-`*`(m[wf, cycle, UN], e[4, UN])

 

25.00350737*kW

 (4.167)

 

Total mechanical exergy at the evaporator inlet

> 

``E[4, UN, M] := Units:-Standard:-`*`(m[wf, cycle, UN], e[4, UN, M])

 

10.71720964*kW

 (4.168)

 

Total thermal exergy at the evaporator inlet

> 

E[4, UN, T] := Units:-Standard:-`*`(m[wf, cycle, UN], e[4, UN, T])

 

14.28629769*kW

 (4.169)

 

Total physical exergy of heat consumption fluid at the evaporator intlet 

> 

E[EV, h, inlet] := Units:-Standard:-`*`(m[EV, h, UN], e[EV, h, inlet])

 

11.30975886*kW

 (4.170)

 

Total physical exergy of heat consumption fluid at the evaporator outlet

> 

E[EV, h, out] := Units:-Standard:-`*`(m[EV, h, UN], e[EV, h, out])

 

22.72348304*kW

 (4.171)

 

Total physical exergy of heat consumption fluid at the condenser intlet 

> 

E[CD, c, inlet] := Units:-Standard:-`*`(m[CD, c, UN], e[CD, c, inlet])

 

-0.

 (4.172)

 

Total physical exergy of heat consumption fluid at the condenser outlet

> 

E[CD, c, out] := Units:-Standard:-`*`(m[CD, c, UN], e[CD, c, out])

 

2.046896275*kW

 (4.173)

 

Exergy of fuel in the compressor

> 

E[F, CM, UN] := E[5, UN]

 

22.72412127*kW

 (4.174)

 

Exergy of product in the compressor

> 

E[P, CM, UN] := Units:-Standard:-`+`(E[2, UN], Units:-Standard:-`-`(E[1, UN]))

 

21.84073657*kW

 (4.175)

 

Exergy destruction in the compressor

> 

E[D, CM, UN] := Units:-Standard:-`+`(E[F, CM, UN], Units:-Standard:-`-`(E[P, CM, UN]))

 

.88338470*kW

 (4.176)

 

Exergetic efficiency in the compressor

> 

varepsilon[CM, UN] := Units:-Standard:-`*`(E[P, CM, UN], Units:-Standard:-`/`(E[F, CM, UN]))

 

.9611256828

 (4.177)

 

Exergy of fuel in the condenser

> 

E[F, CD, UN] := Units:-Standard:-`+`(E[2, UN], Units:-Standard:-`-`(E[3, UN]))

 

6.12171880*kW

 (4.178)

 

Exergy of product in the condenser

> 

E[P, CD, UN] := Units:-Standard:-`+`(E[CD, c, out], Units:-Standard:-`-`(E[CD, c, inlet]))

 

2.046896275*kW

 (4.179)

 

Exergy destruction in the condenser

> 

E[D, CD, UN] := Units:-Standard:-`+`(E[F, CD, UN], Units:-Standard:-`-`(E[P, CD, UN]))

 

4.074822525*kW

 (4.180)

 

Exergetic efficiency in the condenser

> 

varepsilon[CD, UN] := Units:-Standard:-`*`(E[P, CD, UN], Units:-Standard:-`/`(E[F, CD, UN]))

 

.3343662689

 (4.181)

 

Exergy of fuel in the throttling valve

> 

E[F, TV, UN] := Units:-Standard:-`+`(E[3, UN, M], Units:-Standard:-`-`(E[4, UN, M]))

 

16.14674895*kW

 (4.182)

 

Exergy of product in the throttling valve

> 

E[P, TV, UN] := Units:-Standard:-`+`(E[4, UN, T], Units:-Standard:-`-`(E[3, UN, T]))

 

14.20611254*kW

 (4.183)

 

Exergy destruction in the throttling valve

> 

E[D, TV, UN] := Units:-Standard:-`+`(E[F, TV, UN], Units:-Standard:-`-`(E[P, TV, UN]))

 

1.94063641*kW

 (4.184)

 

Exergetic efficiency in the throttling valve

> 

varepsilon[TV, UN] := Units:-Standard:-`*`(E[P, TV, UN], Units:-Standard:-`/`(E[F, TV, UN]))

 

.8798125606

 (4.185)

 

Exergy of fuel in the evaporator

> 

E[F, EV, UN] := Units:-Standard:-`+`(E[4, UN], Units:-Standard:-`-`(E[1, UN]))

 

13.77838147*kW

 (4.186)

 

Exergy of product in the evaporator

> 

E[P, EV, UN] := Units:-Standard:-`+`(E[EV, h, out], Units:-Standard:-`-`(E[EV, h, inlet]))

 

11.41372418*kW

 (4.187)

 

Exergy destruction in the evaporator

> 

E[D, EV, UN] := Units:-Standard:-`+`(E[F, EV, UN], Units:-Standard:-`-`(E[P, EV, UN]))

 

2.36465729*kW

 (4.188)

 

Exergetic efficiency in the evaporator

> 

varepsilon[EV, UN] := Units:-Standard:-`*`(E[P, EV, UN], Units:-Standard:-`/`(E[F, EV, UN]))

 

.8283791681

 (4.189)

 

Exergy of fuel of the overall system

> 

E[F, TOT, UN] := E[5, UN]

22.72412127*kW

 (4.190)

Exergy of product of the overall system 

> 

E[P, TOT, UN] := proc (Type) if Type = 1 then E[P, CD, UN] else E[P, EV, UN] end if end proc:

 

> 

E[P, TOT, UN] := E[P, TOT, UN](Type)

 

11.41372418*kW

 (4.191)

 

Exergy destruction of the overall system

> 

E[D, TOT, UN] := Units:-Standard:-`+`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CM, UN], E[D, CD, UN]), E[D, TV, UN]), E[D, EV, UN])

9.263500925*kW

 (4.192)

Exergetic efficiency of the overall system

> 

varepsilon[TOT, UN] := Units:-Standard:-`*`(E[P, TOT, UN], Units:-Standard:-`/`(E[F, TOT, UN]))

.5022735113

 (4.193)

Exergy destruction ratio in the compressor related to the exergy of total fuel

> 

y[CM, UN] := Units:-Standard:-`*`(E[D, CM, UN], Units:-Standard:-`/`(E[F, TOT, UN]))

 

0.3887431728e-1

 (4.194)

 

Exergy destruction ratio in the condenser related to the exergy of total fuel

> 

y[CD, UN] := Units:-Standard:-`*`(E[D, CD, UN], Units:-Standard:-`/`(E[F, TOT, UN]))

.1793170560

 (4.195)

Exergy destruction ratio in the throttling valve related to the exergy of total fuel

> 

y[TV, UN] := Units:-Standard:-`*`(E[D, TV, UN], Units:-Standard:-`/`(E[F, TOT, UN]))

0.8539984394e-1

 (4.196)

Exergy destruction ratio in the evaporator related to the exergy of total fuel

> 

y[EV, UN] := Units:-Standard:-`*`(E[D, EV, UN], Units:-Standard:-`/`(E[F, TOT, UN]))

.1040593501

 (4.197)

Exergy destruction ratio in the system overall related to the exergy of total fuel

> 

y[TOT, UN] := Units:-Standard:-`*`(E[D, TOT, UN], Units:-Standard:-`/`(E[F, TOT, UN]))

 

.4076505672

 (4.198)

 

Exergy destruction ratio in the compressor related to the total exergy destruction

> 

y[CM, UN, 1] := Units:-Standard:-`*`(E[D, CM, UN], Units:-Standard:-`/`(E[D, TOT, UN]))

 

0.9536186237e-1

 (4.199)

 

Exergy destruction ratio in the condenser related to the total exergy destruction

> 

y[CD, UN, 1] := Units:-Standard:-`*`(E[D, CD, UN], Units:-Standard:-`/`(E[D, TOT, UN]))

.4398793242

 (4.200)

Exergy destruction ratio in the throttling valve related to the total exergy destruction

> 

y[TV, UN, 1] := Units:-Standard:-`*`(E[D, TV, UN], Units:-Standard:-`/`(E[D, TOT, UN]))

.2094927637

 (4.201)

Exergy destruction ratio in the evaporator related to the total exergy destruction

> 

y[EV, UN, 1] := Units:-Standard:-`*`(E[D, EV, UN], Units:-Standard:-`/`(E[D, TOT, UN]))

.2552660501

 (4.202)

Exergy destruction ratio in the system overall related to the total exergy destruction

> 

y[TOT, UN, 1] := Units:-Standard:-`*`(E[D, TOT, UN], Units:-Standard:-`/`(E[D, TOT, UN]))

 

1.000000000

 (4.203)

 

3. Theoretical cycle

Temperature of the working fluid at the compressor outlet

> 

T[4, T] := Units:-Standard:-`+`(T[EV, h, out], 0)

258.15*Units:-Unit(K)

 (4.204)

Pressure of the working fluid at the compressor outlet

> 

p[2, T] := PSTXwf(Units:-Standard:-`+`(T[CD, c, out], 0), 1, wf[cycle])

 

1167197.054*Units:-Unit(Pa)

 (4.205)

 

Pressure of  the working fluid at the condenser outlet

> 

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

 

1167197.054*Units:-Unit(Pa)

 (4.206)

 

Pressure of  the working fluid at the condenser outlet

> 

T[3, T] := TSPXwf(p[3, T], 0, wf[cycle])

 

303.1500000*Units:-Unit(K)

 (4.207)

 

 

Specific enthalpy of  the working fluid at the condenser outlet

> 

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

 

484911.4924*Units:-Unit(J/kg)

 (4.208)

 

Specific entropy of  the working fluid at the condenser outlet

> 

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

 

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

 (4.209)

 

Specific volume of  the working fluid at the condenser outlet

> 

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

 

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

 (4.210)

 

Specific exergy of  the working fluid at the condenser outlet

> 

e[3, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, T], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, T], Units:-Standard:-`-`(s[0])))))

 

296150.505*Units:-Unit(J/kg)

 (4.211)

 

Auxiliary values

> 

h[3, T, 1] := HPTwf(p[3, T], T[0], wf[cycle])

 

437080.3036*Units:-Unit(J/kg)

 (4.212)

 

> 

s[3, T, 1] := SPTwf(p[3, T], T[0], wf[cycle])

 

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

 (4.213)

 

Specific mechanical exergy of  the working fluid at the condenser outlet

> 

e[3, T, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[3, T, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[3, T, 1]))))

 

295350.397*Units:-Unit(J/kg)

 (4.214)

 

Specific thermal exergy of  the working fluid at the condenser outlet

> 

e[3, T, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, T], Units:-Standard:-`-`(h[3, T, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, T], Units:-Standard:-`-`(s[3, T, 1])))))

 

800.10886*Units:-Unit(J/kg)

 (4.215)

 

 

Specific entropy of the working fluid at the evaporator intlet

> 

s[4, T] := s[3, T]

 

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

 (4.216)

 

Pressure of the working fluid at the evaporator intlet

 

> 

p[4, T] := Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`*`(PSTXwf(T[4, T], 0, wf[cycle]), Units:-Standard:-`+`(SSTXwf(T[4, T], 1, wf[cycle]), Units:-Standard:-`-`(s[4, T]))), Units:-Standard:-`*`(PSTXwf(T[4, T], 1, wf[cycle]), Units:-Standard:-`+`(s[4, T], Units:-Standard:-`-`(SSTXwf(T[4, R], 0, wf[cycle]))))), Units:-Standard:-`/`(Units:-Standard:-`+`(SSTXwf(T[4, T], 1, wf[cycle]), Units:-Standard:-`-`(SSTXwf(T[4, T], 0, wf[cycle])))))

 

244434.7345*Units:-Unit(Pa)

 (4.217)

 

Temperature of the working fluid at the compressor outlet

> 

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

ThermophysicalData:-Property(T, pressure = 244434.7345*Units:-Unit(Pa), enthalpy = h[4, T], R717)

 (4.218)

Specific enthalpy of saturated liquid of  the working fluid in the evaporator

> 

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

278068.0318*Units:-Unit(J/kg)

 (4.219)

Specific enthalpy of saturated vapor of  the working fluid in the evaporator

> 

h[sv, 4, T] := HSPXwf(p[4, T], 1, wf[cycle])

1588577.451*Units:-Unit(J/kg)

 (4.220)

Specific entropy of saturated liquid of  the working fluid in the evaporator

> 

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

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

 (4.221)

Specific entropy of saturated vapor of  the working fluid in the evaporator

> 

s[sv, 4, T] := SSPXwf(p[4, T], 1, wf[cycle])

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

 (4.222)

Quality of the working fluid at the evaporator intlet

> 

x[4, T] := Units:-Standard:-`*`(Units:-Standard:-`+`(s[4, T], Units:-Standard:-`-`(s[sl, 4, T])), Units:-Standard:-`/`(Units:-Standard:-`+`(s[sv, 4, T], Units:-Standard:-`-`(s[sl, 4, T]))))NULL

 

.1445814939

 (4.223)

 

 

Specific entropy of  the working fluid at the evaporator inlet

> 

h[4, T] := Units:-Standard:-`+`(h[sl, 4, T], Units:-Standard:-`*`(x[4, T], Units:-Standard:-`+`(h[sv, 4, T], Units:-Standard:-`-`(h[sl, 4, T]))))

467543.4414*Units:-Unit(J/kg)

 (4.224)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

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

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

 (4.225)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

v[sv, 4, T] := Units:-Standard:-`/`(DSPXwf(p[4, T], 1, wf[cycle]))

.4924834385*Units:-Unit(m^3/kg)

 (4.226)

Specific volume of  the working fluid at the evaporator inlet

> 

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

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

 (4.227)

Specific exergy of  the working fluid at the evaporator inlet

> 

e[4, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, T], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, T], Units:-Standard:-`-`(s[0])))))

 

278782.454*Units:-Unit(J/kg)

 (4.228)

 

Auxiliary values

> 

h[4, T, 1] := HPTwf(p[4, T], T[0], wf[cycle])

 

1669482.240*Units:-Unit(J/kg)

 (4.229)

 

> 

s[4, T, 1] := SPTwf(p[4, T], T[0], wf[cycle])

 

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

 (4.230)

 

Specific mechanical exergy of  the working fluid at the evaporator inlet

> 

e[4, T, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[4, T, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[4, T, 1]))))

 

125595.3239*Units:-Unit(J/kg)

 (4.231)

 

Specific thermal exergy of  the working fluid at the evaporator inlet

> 

e[4, T, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, T], Units:-Standard:-`-`(h[4, T, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, T], Units:-Standard:-`-`(s[4, T, 1])))))

 

153187.130*Units:-Unit(J/kg)

 (4.232)

 

Pressure of  the working fluid at the compressor inlet

> 

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

244434.7345*Units:-Unit(Pa)

 (4.233)

Temperature of  the working fluid at the compressor inlet

> 

T[1, T] := TSPXwf(p[1, T], 1, wf[cycle])

258.9630774*Units:-Unit(K)

 (4.234)

Specific enthalpy of  the working fluid at the compressor inlet

> 

h[1, T] := HSPXwf(p[1, T], 1, wf[cycle])

1588577.451*Units:-Unit(J/kg)

 (4.235)

Specific entropy of  the working fluid at the compressor inlet

> 

s[1, T] := SSPXwf(p[1, T], 1, wf[cycle])

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

 (4.236)

Specific volume of  the working fluid at the compressor inlet

> 

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

 

.4924834376*Units:-Unit(m^3/kg)

 (4.237)

 

Specific exergy of  the working fluid at the compressor inlet

> 

e[1, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[1, T], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[1, T], Units:-Standard:-`-`(s[0])))))

 

130789.5301*Units:-Unit(J/kg)

 (4.238)

 

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

> 

h[2, T] := HcompPTPEFFwf(p[1, T], T[1, T], p[2, T], 1, wf[cycle])

1814495.110*Units:-Unit(J/kg)

 (4.239)

Enhtalpy change in the compressor 1 after actual compression

> 

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

 

225917.659*Units:-Unit(J/kg)

 (4.240)

 

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

> 

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

 

370.4730620*Units:-Unit(K)

 (4.241)

 

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

> 

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

 

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

 (4.242)

 

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

> 

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

.1464462235*Units:-Unit(m^3/kg)

 (4.243)

Specific exergy of  the working fluid at the compressor inlet

> 

e[2, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[2, T], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[2, T], Units:-Standard:-`-`(s[0])))))

 

356707.1896*Units:-Unit(J/kg)

 (4.244)

 

 

Specific heat rejection in the condenser

> 

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

1329583.618*Units:-Unit(J/kg)

 (4.245)

Specific heat addition in the evaporator

> 

q[inlet, T] := Units:-Standard:-`+`(h[1, T], Units:-Standard:-`-`(h[4, T]))

1121034.010*Units:-Unit(J/kg)

 (4.246)

Coefficent of performance of a refrigerator

> 

COP[T, Ref] := Units:-Standard:-`*`(q[inlet, T], Units:-Standard:-`/`(l[comp, T]))

4.962135386

 (4.247)

Coefficent of performance of a heat pump

> 

COP[T, HP] := Units:-Standard:-`*`(q[out, T], Units:-Standard:-`/`(l[comp, T]))

5.885257593

 (4.248)

Mass rate of the working fluid

> 

m[wf, cycle, T] := proc (Type) if Type = 1 then Q[P]/q[out, T] else Q[P]/q[inlet, T] end if end proc:

 

> 

m[wf, cycle, T] := m[wf, cycle, T](Type)

 

0.8920335967e-4*kW*Units:-Unit(kg/J)

 (4.249)

 

Total heat addition in the evaporator

> 

Q[inlet, T] := Units:-Standard:-`*`(m[wf, cycle, T], q[inlet, T])

 

100.0000000*kW

 (4.250)

 

Total heat rejection in the condenser

> 

Q[out, T] := Units:-Standard:-`*`(m[wf, cycle, R], q[out, T])

 

127.8427323*kW

 (4.251)

 

Mass rate of the heat source fluid

> 

m[EV, h, T] := Units:-Standard:-`*`(Q[inlet, T], Units:-Standard:-`/`(Units:-Standard:-`+`(h[EV, h, inlet], Units:-Standard:-`-`(h[EV, h, out]))))

 

0.9944845786e-2*kW*Units:-Unit(kg/J)

 (4.252)

 

Mass rate of the heat consumption fluid

> 

m[CD, c, T] := Units:-Standard:-`*`(Q[out, T], Units:-Standard:-`/`(Units:-Standard:-`+`(h[CD, c, out], Units:-Standard:-`-`(h[CD, c, inlet]))))

 

0.1270435976e-1*kW*Units:-Unit(kg/J)

 (4.253)

 

Total exergy used in the compressor

> 

E[5, T] := Units:-Standard:-`*`(m[wf, cycle, T], l[comp, T])

 

20.15261419*kW

 (4.254)

 

Total physical exergy at the compressor intlet

> 

E[1, T] := Units:-Standard:-`*`(m[wf, cycle, T], e[1, T])

 

11.66686549*kW

 (4.255)

 

Total physical exergy at the compressor outlet

> 

E[2, T] := Units:-Standard:-`*`(m[wf, cycle, T], e[2, T])

 

31.81947973*kW

 (4.256)

 

Total physical exergy at the condenser outlet

> 

E[3, T] := Units:-Standard:-`*`(m[wf, cycle, T], e[3, T])

 

26.41762001*kW

 (4.257)

 

Total mechanical exergy at the condenser outlet

> 

E[3, T, M] := Units:-Standard:-`*`(m[wf, cycle, T], e[3, T, M])

 

26.34624769*kW

 (4.258)

 

Total thermal exergy at the condenser outlet

> 

E[3, T, T] := Units:-Standard:-`*`(m[wf, cycle, T], e[3, T, T])

 

0.7137239841e-1*kW

 (4.259)

 

Total physical exergy at the evaporator inlet

> 

E[4, T] := Units:-Standard:-`*`(m[wf, cycle, T], e[4, T])

 

24.86833151*kW

 (4.260)

 

Total mechanical exergy at the evaporator inlet

> 

E[4, T, M] := Units:-Standard:-`*`(m[wf, cycle, T], e[4, T, M])

 

11.20352485*kW

 (4.261)

 

Total thermal exergy at the evaporator inlet

> 

E[4, T, T] := Units:-Standard:-`*`(m[wf, cycle, T], e[4, T, T])

 

13.66480665*kW

 (4.262)

 

Total physical exergy of heat consumption fluid at the evaporator intlet 

> 

E[EV, h, inlet] := Units:-Standard:-`*`(m[EV, h, T], e[EV, h, inlet])

 

11.30975886*kW

 (4.263)

 

Total physical exergy of heat consumption fluid at the evaporator outlet

> 

E[EV, h, out] := Units:-Standard:-`*`(m[EV, h, T], e[EV, h, out])

 

22.72348304*kW

 (4.264)

 

Total physical exergy of heat consumption fluid at the condenser intlet 

> 

E[CD, c, inlet] := Units:-Standard:-`*`(m[CD, c, T], e[CD, c, inlet])

 

-0.

 (4.265)

 

Total physical exergy of heat consumption fluid at the condenser outlet

> 

E[CD, c, out] := Units:-Standard:-`*`(m[CD, c, T], e[CD, c, out])

 

2.132268782*kW

 (4.266)

 

Exergy of fuel in the compressor

> 

E[F, CM, T] := E[5, T]

 

20.15261419*kW

 (4.267)

 

Exergy of product in the compressor

> 

E[P, CM, T] := Units:-Standard:-`+`(E[2, T], Units:-Standard:-`-`(E[1, T]))

 

20.15261424*kW

 (4.268)

 

Exergy destruction in the compressor

> 

E[D, CM, T] := Units:-Standard:-`+`(E[F, CM, T], Units:-Standard:-`-`(E[P, CM, T]))

 

-0.5e-7*kW

 (4.269)

 

Exergetic efficiency in the compressor

> 

varepsilon[CM, T] := Units:-Standard:-`*`(E[P, CM, T], Units:-Standard:-`/`(E[F, CM, T]))

 

1.000000002

 (4.270)

 

Exergy of fuel in the condenser

> 

E[F, CD, T] := Units:-Standard:-`+`(E[2, T], Units:-Standard:-`-`(E[3, T]))

 

5.40185972*kW

 (4.271)

 

Exergy of product in the condenser

> 

E[P, CD, T] := Units:-Standard:-`+`(E[CD, c, out], Units:-Standard:-`-`(E[CD, c, inlet]))

 

2.132268782*kW

 (4.272)

 

Exergy destruction in the condenser

> 

E[D, CD, T] := Units:-Standard:-`+`(E[F, CD, T], Units:-Standard:-`-`(E[P, CD, T]))

 

3.269590938*kW

 (4.273)

 

Exergetic efficiency in the condenser

> 

varepsilon[CD, T] := Units:-Standard:-`*`(E[P, CD, T], Units:-Standard:-`/`(E[F, CD, T]))

 

.3947286476

 (4.274)

 

Exergy of fuel in the throttling valve

> 

E[F, TV, T] := Units:-Standard:-`+`(E[3, T, M], Units:-Standard:-`-`(E[4, T, M]))

 

15.14272284*kW

 (4.275)

 

Exergy of product in the throttling valve

> 

E[P, TV, T] := Units:-Standard:-`+`(E[3, T, M], Units:-Standard:-`-`(E[4, T, M]))

 

15.14272284*kW

 (4.276)

 

Exergy destruction in the throttling valve

> 

E[D, TV, T] := Units:-Standard:-`+`(E[F, TV, T], Units:-Standard:-`-`(E[P, TV, T]))

 

0.

 (4.277)

 

Exergetic efficiency in the throttling valve

> 

varepsilon[TV, T] := Units:-Standard:-`*`(E[P, TV, T], Units:-Standard:-`/`(E[F, TV, T]))

 

1.000000000

 (4.278)

 

Exergy of fuel in the evaporator

> 

E[F, EV, T] := Units:-Standard:-`+`(E[4, T], Units:-Standard:-`-`(E[1, T]))

 

13.20146602*kW

 (4.279)

 

Exergy of product in the evaporator

> 

E[P, EV, T] := Units:-Standard:-`+`(E[EV, h, out], Units:-Standard:-`-`(E[EV, h, inlet]))

 

11.41372418*kW

 (4.280)

 

Exergy destruction in the evaporator

> 

E[D, EV, T] := Units:-Standard:-`+`(E[F, EV, T], Units:-Standard:-`-`(E[P, EV, T]))

 

1.78774184*kW

 (4.281)

 

Exergetic efficiency in the evaporator

> 

varepsilon[EV, T] := Units:-Standard:-`*`(E[P, EV, T], Units:-Standard:-`/`(E[F, EV, T]))

 

.8645800521

 (4.282)

 

Exergy of fuel of the overall system

> 

E[F, TOT, T] := E[5, T]

20.15261419*kW

 (4.283)

Exergy of product of the overall system 

> 

E[P, TOT, T] := proc (Type) if Type = 1 then E[P, CD, T] else E[P, EV, T] end if end proc:

 

> 

E[P, TOT, T] := E[P, TOT, T](Type)

 

11.41372418*kW

 (4.284)

 

Exergy destruction of the overall system

> 

E[D, TOT, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CM, T], E[D, CD, T]), E[D, TV, T]), E[D, EV, T])

5.057332728*kW

 (4.285)

Exergetic efficiency of the overall system

> 

varepsilon[TOT, T] := Units:-Standard:-`*`(E[P, TOT, T], Units:-Standard:-`/`(E[F, TOT, T]))

.5663644464

 (4.286)

Exergy destruction ratio in the compressor related to the exergy of total fuel

> 

y[CM, T] := Units:-Standard:-`*`(E[D, CM, T], Units:-Standard:-`/`(E[F, TOT, T]))

 

-0.2481067693e-8

 (4.287)

 

Exergy destruction ratio in the condenser related to the exergy of total fuel

> 

y[CD, T] := Units:-Standard:-`*`(E[D, CD, T], Units:-Standard:-`/`(E[F, TOT, T]))

.1622415289

 (4.288)

Exergy destruction ratio in the throttling valve related to the exergy of total fuel

> 

y[TV, T] := Units:-Standard:-`*`(E[D, TV, T], Units:-Standard:-`/`(E[F, TOT, T]))

0.

 (4.289)

Exergy destruction ratio in the evaporator related to the exergy of total fuel

> 

y[EV, T] := Units:-Standard:-`*`(E[D, EV, T], Units:-Standard:-`/`(E[F, TOT, T]))

0.8871017045e-1

 (4.290)

Exergy destruction ratio in the system overall related to the exergy of total fuel

> 

y[TOT, T] := Units:-Standard:-`*`(E[D, TOT, T], Units:-Standard:-`/`(E[F, TOT, T]))

 

.2509516969

 (4.291)

 

Exergy destruction ratio in the compressor related to the total exergy destruction

> 

y[CM, T, 1] := Units:-Standard:-`*`(E[D, CM, T], Units:-Standard:-`/`(E[D, TOT, T]))

 

-0.9886634455e-8

 (4.292)

 

Exergy destruction ratio in the condenser related to the total exergy destruction

> 

y[CD, T, 1] := Units:-Standard:-`*`(E[D, CD, T], Units:-Standard:-`/`(E[D, TOT, T]))

.6465050084

 (4.293)

Exergy destruction ratio in the throttling valve related to the total exergy destruction

> 

y[TV, T, 1] := Units:-Standard:-`*`(E[D, TV, T], Units:-Standard:-`/`(E[D, TOT, T]))

0.

 (4.294)

Exergy destruction ratio in the evaporator related to the total exergy destruction

> 

y[EV, T, 1] := Units:-Standard:-`*`(E[D, EV, T], Units:-Standard:-`/`(E[D, TOT, T]))

.3534950014

 (4.295)

Exergy destruction ratio in the system overall related to the total exergy destruction

> 

y[TOT, T, 1] := Units:-Standard:-`*`(E[D, TOT, T], Units:-Standard:-`/`(E[D, TOT, T]))

 

1.000000000

 (4.296)

 

 

4. Hybrid cycle with irreversible component - Compressor

Temperature of the working fluid at the compressor outlet

> 

T[4, EN, CM] := Units:-Standard:-`+`(T[EV, h, out], 0)

258.15*Units:-Unit(K)

 (4.297)

Pressure of the working fluid at the compressor outlet

> 

p[2, `EN,CM`] := PSTXwf(Units:-Standard:-`+`(T[CD, c, out], 0), 1, wf[cycle])

 

1167197.054*Units:-Unit(Pa)

 (4.298)

 

Pressure of  the working fluid at the condenser outlet

> 

p[3, `EN,CM`] := p[2, `EN,CM`]

 

1167197.054*Units:-Unit(Pa)

 (4.299)

 

Pressure of  the working fluid at the condenser outlet

> 

T[3, `EN,CM`] := TSPXwf(p[3, `EN,CM`], 0, wf[cycle])

 

303.1500000*Units:-Unit(K)

 (4.300)

 

 

Specific enthalpy of  the working fluid at the condenser outlet

> 

h[3, `EN,CM`] := HSTXwf(T[3, `EN,CM`], 0, wf[cycle])

 

484911.4924*Units:-Unit(J/kg)

 (4.301)

 

Specific entropy of  the working fluid at the condenser outlet

> 

s[3, `EN,CM`] := SSTXwf(T[3, `EN,CM`], 0, wf[cycle])

 

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

 (4.302)

 

Specific volume of  the working fluid at the condenser outlet

> 

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

 

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

 (4.303)

 

Specific exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,CM`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, `EN,CM`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, `EN,CM`], Units:-Standard:-`-`(s[0])))))

 

296150.505*Units:-Unit(J/kg)

 (4.304)

 

Auxiliary values

> 

h[3, `EN,CM`, 1] := HPTwf(p[3, `EN,CM`], T[0], wf[cycle])

 

437080.3036*Units:-Unit(J/kg)

 (4.305)

 

> 

s[3, `EN,CM`, 1] := SPTwf(p[3, `EN,CM`], T[0], wf[cycle])

 

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

 (4.306)

 

Specific mechanical exergy of  the working fluid at the condenser outlet

> 

e[3, EN, CM, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[3, `EN,CM`, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[3, `EN,CM`, 1]))))

 

295350.397*Units:-Unit(J/kg)

 (4.307)

 

Specific thermal exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,CM`, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, `EN,CM`], Units:-Standard:-`-`(h[3, `EN,CM`, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, `EN,CM`], Units:-Standard:-`-`(s[3, `EN,CM`, 1])))))

 

800.10886*Units:-Unit(J/kg)

 (4.308)

 

Specific entropy of the working fluid at the evaporator intlet

> 

s[4, EN, CM] := s[3, `EN,CM`]

 

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

 (4.309)

 

Pressure of the working fluid at the evaporator intlet

 

> 

p[4, `EN,CM`] := Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`*`(PSTXwf(T[4, EN, CM], 0, wf[cycle]), Units:-Standard:-`+`(SSTXwf(T[4, EN, CM], 1, wf[cycle]), Units:-Standard:-`-`(s[4, EN, CM]))), Units:-Standard:-`*`(PSTXwf(T[4, EN, CM], 1, wf[cycle]), Units:-Standard:-`+`(s[4, EN, CM], Units:-Standard:-`-`(SSTXwf(T[4, EN, CM], 0, wf[cycle]))))), Units:-Standard:-`/`(Units:-Standard:-`+`(SSTXwf(T[4, EN, CM], 1, wf[cycle]), Units:-Standard:-`-`(SSTXwf(T[4, EN, CM], 0, wf[cycle])))))

 

236167.1501*Units:-Unit(Pa)

 (4.310)

 

Specific enthalpy of saturated liquid of  the working fluid in the evaporator

> 

h[sl, 4, `EN,CM`] := HSPXwf(p[4, `EN,CM`], 0, wf[cycle])

274370.1868*Units:-Unit(J/kg)

 (4.311)

Specific enthalpy of saturated vapor of  the working fluid in the evaporator

> 

h[sv, 4, `EN,CM`] := HSPXwf(p[4, `EN,CM`], 1, wf[cycle])

1587522.359*Units:-Unit(J/kg)

 (4.312)

Specific entropy of saturated liquid of  the working fluid in the evaporator

> 

s[sl, 4, `EN,CM`] := SSPXwf(p[4, `EN,CM`], 0, wf[cycle])

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

 (4.313)

Specific entropy of saturated vapor of  the working fluid in the evaporator

> 

s[sv, 4, `EN,CM`] := SSPXwf(p[4, `EN,CM`], 1, wf[cycle])

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

 (4.314)

Quality of the working fluid at the evaporator intlet

> 

x[4, `EN,CM`] := Units:-Standard:-`*`(Units:-Standard:-`+`(s[4, EN, CM], Units:-Standard:-`-`(s[sl, 4, `EN,CM`])), Units:-Standard:-`/`(Units:-Standard:-`+`(s[sv, 4, `EN,CM`], Units:-Standard:-`-`(s[sl, 4, `EN,CM`]))))NULL

 

.1466395135

 (4.315)

 

 

Specific entropy of  the working fluid at the evaporator inlet

> 

h[4, `EN,CM`] := Units:-Standard:-`+`(h[sl, 4, `EN,CM`], Units:-Standard:-`*`(x[4, `EN,CM`], Units:-Standard:-`+`(h[sv, 4, `EN,CM`], Units:-Standard:-`-`(h[sl, 4, `EN,CM`]))))

466930.1825*Units:-Unit(J/kg)

 (4.316)

Specific entropy of saturated vapor of  the working fluid in the evaporator

 

> 

s[sv, 4, `EN,CM`] := SSPXwf(p[4, `EN,CM`], 1, wf[cycle])

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

 (4.317)

Specific entropy of  the working fluid at the evaporator inlet

> 

s[4, `EN,CM`] := Units:-Standard:-`+`(s[sl, 4, `EN,CM`], Units:-Standard:-`*`(x[4, `EN,CM`], Units:-Standard:-`+`(s[sv, 4, `EN,CM`], Units:-Standard:-`-`(s[sl, 4, `EN,CM`]))))

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

 (4.318)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

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

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

 (4.319)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

v[sv, 4, `EN,CM`] := Units:-Standard:-`/`(DSPXwf(p[4, `EN,CM`], 1, wf[cycle]))

.5086794088*Units:-Unit(m^3/kg)

 (4.320)

Specific volume of  the working fluid at the evaporator inlet

> 

v[4, `EN,CM`] := Units:-Standard:-`+`(v[sl, 4, `EN,CM`], Units:-Standard:-`*`(x[4, `EN,CM`], Units:-Standard:-`+`(v[sv, 4, `EN,CM`], Units:-Standard:-`-`(v[sl, 4, `EN,CM`]))))

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

 (4.321)

Specific exergy of  the working fluid at the evaporator inlet

> 

e[4, `EN,CM`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, `EN,CM`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, `EN,CM`], Units:-Standard:-`-`(s[0])))))

 

278169.196*Units:-Unit(J/kg)

 (4.322)

 

Auxiliary values

> 

h[4, `EN,CM`, 1] := HPTwf(p[4, `EN,CM`], T[0], wf[cycle])

 

1670040.049*Units:-Unit(J/kg)

 (4.323)

 

> 

s[4, `EN,CM`, 1] := SPTwf(p[4, `EN,CM`], T[0], wf[cycle])

 

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

 (4.324)

 

Specific mechanical exergy of  the working fluid at the evaporator inlet

> 

e[4, `EN,CM`, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[4, `EN,CM`, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[4, `EN,CM`, 1]))))

 

120804.8345*Units:-Unit(J/kg)

 (4.325)

 

Specific thermal exergy of  the working fluid at the evaporator inlet

> 

e[4, `EN,CM`, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, `EN,CM`], Units:-Standard:-`-`(h[4, `EN,CM`, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, `EN,CM`], Units:-Standard:-`-`(s[4, `EN,CM`, 1])))))

 

157364.361*Units:-Unit(J/kg)

 (4.326)

 

Pressure of  the working fluid at the compressor inlet

> 

p[1, `EN,CM`] := p[4, `EN,CM`]

236167.1501*Units:-Unit(Pa)

 (4.327)

Temperature of  the working fluid at the compressor inlet

> 

T[1, `EN,CM`] := TSPXwf(p[1, `EN,CM`], 1, wf[cycle])

258.1500000*Units:-Unit(K)

 (4.328)

Specific enthalpy of  the working fluid at the compressor inlet

> 

h[1, `EN,CM`] := HSPXwf(p[1, `EN,CM`], 1, wf[cycle])

1587522.359*Units:-Unit(J/kg)

 (4.329)

Specific entropy of  the working fluid at the compressor inlet

> 

s[1, `EN,CM`] := SSPXwf(p[1, `EN,CM`], 1, wf[cycle])

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

 (4.330)

Specific volume of  the working fluid at the compressor inlet

> 

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

 

.5086794088*Units:-Unit(m^3/kg)

 (4.331)

 

Specific exergy of  the working fluid at the compressor inlet

> 

e[1, `EN,CM`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[1, `EN,CM`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[1, `EN,CM`], Units:-Standard:-`-`(s[0])))))

 

126239.2086*Units:-Unit(J/kg)

 (4.332)

 

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

> 

h[2, `EN,CM`] := HcompPTPEFFwf(p[1, `EN,CM`], T[1, `EN,CM`], p[2, `EN,CM`], eta[comp], wf[cycle])

1876772.826*Units:-Unit(J/kg)

 (4.333)

Enhtalpy change in the compressor 1 after actual compression

> 

l[comp, `EN,CM`] := Units:-Standard:-`+`(h[2, `EN,CM`], Units:-Standard:-`-`(h[1, `EN,CM`]))

 

289250.467*Units:-Unit(J/kg)

 (4.334)

 

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

> 

T[2, `EN,CM`] := TPHwf(p[2, `EN,CM`], h[2, `EN,CM`], wf[cycle])

 

395.4188136*Units:-Unit(K)

 (4.335)

 

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

> 

s[2, `EN,CM`] := SPTwf(p[2, `EN,CM`], T[2, `EN,CM`], wf[cycle])

 

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

 (4.336)

 

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

> 

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

.1583179873*Units:-Unit(m^3/kg)

 (4.337)

Specific exergy of  the working fluid at the compressor inlet

> 

e[2, `EN,CM`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[2, `EN,CM`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[2, `EN,CM`], Units:-Standard:-`-`(s[0])))))

 

371289.4760*Units:-Unit(J/kg)

 (4.338)

 

 

Specific heat rejection in the condenser

> 

q[out, `EN,CM`] := Units:-Standard:-`+`(h[2, `EN,CM`], Units:-Standard:-`-`(h[3, `EN,CM`]))

1391861.334*Units:-Unit(J/kg)

 (4.339)

Specific heat addition in the evaporator

> 

q[inlet, `EN,CM`] := Units:-Standard:-`+`(h[1, `EN,CM`], Units:-Standard:-`-`(h[4, `EN,CM`]))

1120592.176*Units:-Unit(J/kg)

 (4.340)

Coefficent of performance of a refrigerator

> 

COP[`EN,CM`, Ref] := Units:-Standard:-`*`(q[inlet, `EN,CM`], Units:-Standard:-`/`(l[comp, `EN,CM`]))

3.874124000

 (4.341)

Coefficent of performance of a heat pump

> 

COP[`EN,CM`, HP] := Units:-Standard:-`*`(q[out, `EN,CM`], Units:-Standard:-`/`(l[comp, `EN,CM`]))

4.811958814

 (4.342)

Mass rate of the working fluid

> 

m[wf, cycle, `EN,CM`] := proc (Type) if Type = 1 then Q[P]/q[out, `EN,CM`] else Q[P]/q[inlet, `EN,CM`] end if end proc:

 

> 

m[wf, cycle, `EN,CM`] := m[wf, cycle, `EN,CM`](Type)

 

0.8923853132e-4*kW*Units:-Unit(kg/J)

 (4.343)

 

Total heat addition in the evaporator

> 

Q[inlet, `EN,CM`] := Units:-Standard:-`*`(m[wf, cycle, `EN,CM`], q[inlet, `EN,CM`])

 

99.99999999*kW

 (4.344)

 

Total heat rejection in the condenser

> 

Q[out, `EN,CM`] := Units:-Standard:-`*`(m[wf, cycle, R], q[out, `EN,CM`])

 

133.8308877*kW

 (4.345)

 

Mass rate of the heat source fluid

> 

m[EV, h, `EN,CM`] := Units:-Standard:-`*`(Q[inlet, `EN,CM`], Units:-Standard:-`/`(Units:-Standard:-`+`(h[EV, h, inlet], Units:-Standard:-`-`(h[EV, h, out]))))

 

0.9944845785e-2*kW*Units:-Unit(kg/J)

 (4.346)

 

Mass rate of the heat consumption fluid

> 

m[CD, c, `EN,CM`] := Units:-Standard:-`*`(Q[out, `EN,CM`], Units:-Standard:-`/`(Units:-Standard:-`+`(h[CD, c, out], Units:-Standard:-`-`(h[CD, c, inlet]))))

 

0.1329943215e-1*kW*Units:-Unit(kg/J)

 (4.347)

 

Total exergy used in the compressor

> 

E[5, `EN,CM`] := Units:-Standard:-`*`(m[wf, cycle, `EN,CM`], l[comp, `EN,CM`])

 

25.81228686*kW

 (4.348)

 

Total physical exergy at the compressor intlet

> 

E[1, `EN,CM`] := Units:-Standard:-`*`(m[wf, cycle, `EN,CM`], e[1, `EN,CM`])

 

11.26540157*kW

 (4.349)

 

Total physical exergy at the compressor outlet

> 

E[2, `EN,CM`] := Units:-Standard:-`*`(m[wf, cycle, `EN,CM`], e[2, `EN,CM`])

 

33.13332753*kW

 (4.350)

 

Total physical exergy at the condenser outlet

> 

E[3, `EN,CM`] := Units:-Standard:-`*`(m[wf, cycle, `EN,CM`], e[3, `EN,CM`])

 

26.42803612*kW

 (4.351)

 

Total mechanical exergy at the condenser outlet

> 

E[3, EN, CM, M] := Units:-Standard:-`*`(m[wf, cycle, `EN,CM`], e[3, EN, CM, M])

 

26.35663565*kW

 (4.352)

 

Total thermal exergy at the condenser outlet

> 

E[3, `EN,CM`, T] := Units:-Standard:-`*`(m[wf, cycle, `EN,CM`], e[3, `EN,CM`, T])

 

0.7140053956e-1*kW

 (4.353)

 

Total physical exergy at the evaporator inlet

> 

E[4, EN, CM] := Units:-Standard:-`*`(m[wf, cycle, `EN,CM`], e[4, `EN,CM`])

 

24.82341051*kW

 (4.354)

 

Total mechanical exergy at the evaporator inlet

> 

E[4, EN, CM, M] := Units:-Standard:-`*`(m[wf, cycle, `EN,CM`], e[4, `EN,CM`, M])

 

10.78044601*kW

 (4.355)

 

Total thermal exergy at the evaporator inlet

> 

E[4, `EN,CM`, T] := Units:-Standard:-`*`(m[wf, cycle, `EN,CM`], e[4, `EN,CM`, T])

 

14.04296446*kW

 (4.356)

 

Total physical exergy of heat consumption fluid at the evaporator intlet 

> 

E[EV, h, inlet] := Units:-Standard:-`*`(m[EV, h, `EN,CM`], e[EV, h, inlet])

 

11.30975886*kW

 (4.357)

 

Total physical exergy of heat consumption fluid at the evaporator outlet

> 

E[EV, h, out] := Units:-Standard:-`*`(m[EV, h, `EN,CM`], e[EV, h, out])

 

22.72348303*kW

 (4.358)

 

Total physical exergy of heat consumption fluid at the condenser intlet 

> 

E[CD, c, inlet] := Units:-Standard:-`*`(m[CD, c, `EN,CM`], e[CD, c, inlet])

 

-0.

 (4.359)

 

Total physical exergy of heat consumption fluid at the condenser outlet

> 

E[CD, c, out] := Units:-Standard:-`*`(m[CD, c, `EN,CM`], e[CD, c, out])

 

2.232144281*kW

 (4.360)

 

Exergy of fuel in the compressor

> 

E[F, CM, `EN,CM`] := E[5, `EN,CM`]

 

25.81228686*kW

 (4.361)

 

Exergy of product in the compressor

> 

E[P, CM, `EN,CM`] := Units:-Standard:-`+`(E[2, `EN,CM`], Units:-Standard:-`-`(E[1, `EN,CM`]))

 

21.86792596*kW

 (4.362)

 

Exergy destruction in the compressor

> 

E[D, CM, EN, CM] := Units:-Standard:-`+`(E[F, CM, `EN,CM`], Units:-Standard:-`-`(E[P, CM, `EN,CM`]))

 

3.94436090*kW

 (4.363)

 

Exergetic efficiency in the compressor

> 

varepsilon[CM, `EN,CM`] := Units:-Standard:-`*`(E[P, CM, `EN,CM`], Units:-Standard:-`/`(E[F, CM, `EN,CM`]))

 

.8471905677

 (4.364)

 

Exergy of fuel in the condenser

> 

E[F, CD, `EN,CM`] := Units:-Standard:-`+`(E[2, `EN,CM`], Units:-Standard:-`-`(E[3, `EN,CM`]))

 

6.70529141*kW

 (4.365)

 

Exergy of product in the condenser

> 

E[P, CD, `EN,CM`] := Units:-Standard:-`+`(E[CD, c, out], Units:-Standard:-`-`(E[CD, c, inlet]))

 

2.232144281*kW

 (4.366)

 

Exergy destruction in the condenser

> 

E[D, CD, `EN,CM`] := Units:-Standard:-`+`(E[F, CD, `EN,CM`], Units:-Standard:-`-`(E[P, CD, `EN,CM`]))

 

4.473147129*kW

 (4.367)

 

Exergy of fuel in the throttling valve

> 

E[F, TV, `EN,CM`] := Units:-Standard:-`+`(E[3, EN, CM, M], Units:-Standard:-`-`(E[4, EN, CM, M]))

 

15.57618964*kW

 (4.368)

 

Exergy of product in the throttling valve

> 

E[P, TV, `EN,CM`] := Units:-Standard:-`+`(E[3, EN, CM, M], Units:-Standard:-`-`(E[4, EN, CM, M]))

 

15.57618964*kW

 (4.369)

 

Exergy destruction in the throttling valve

> 

E[D, TV, `EN,CM`] := Units:-Standard:-`+`(E[F, TV, `EN,CM`], Units:-Standard:-`-`(E[P, TV, `EN,CM`]))

 

0.

 (4.370)

 

Exergetic efficiency in the throttling valve

> 

varepsilon[TV, `EN,CM`] := Units:-Standard:-`*`(E[P, TV, `EN,CM`], Units:-Standard:-`/`(E[F, TV, `EN,CM`]))

 

1.000000000

 (4.371)

 

Exergy of fuel in the evaporator

> 

E[F, EV, `EN,CM`] := Units:-Standard:-`+`(E[4, EN, CM], Units:-Standard:-`-`(E[1, `EN,CM`]))

 

13.55800894*kW

 (4.372)

 

Exergy of product in the evaporator

> 

E[P, EV, `EN,CM`] := Units:-Standard:-`+`(E[EV, h, out], Units:-Standard:-`-`(E[EV, h, inlet]))

 

11.41372417*kW

 (4.373)

 

Exergy destruction in the evaporator

> 

E[D, EV, `EN,CM`] := Units:-Standard:-`+`(E[F, EV, `EN,CM`], Units:-Standard:-`-`(E[P, EV, `EN,CM`]))

 

2.14428477*kW

 (4.374)

 

Exergetic efficiency in the evaporator

> 

varepsilon[EV, `EN,CM`] := Units:-Standard:-`*`(E[P, EV, `EN,CM`], Units:-Standard:-`/`(E[F, EV, `EN,CM`]))

 

.8418436822

 (4.375)

 

Exergy of fuel of the overall system

> 

E[F, TOT, `EN,CM`] := E[5, `EN,CM`]

25.81228686*kW

 (4.376)

Exergy of product of the overall system 

> 

E[P, TOT, `EN,CM`] := proc (Type) if Type = 1 then E[P, CD, `EN,CM`] else E[P, EV, `EN,CM`] end if end proc:

 

> 

E[P, TOT, `EN,CM`] := E[P, TOT, `EN,CM`](Type)

 

11.41372417*kW

 (4.377)

 

Exergy destruction of the overall system

> 

E[D, TOT, `EN,CM`] := Units:-Standard:-`+`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CM, EN, CM], E[D, CD, `EN,CM`]), E[D, TV, `EN,CM`]), E[D, EV, `EN,CM`])

10.56179280*kW

 (4.378)

Exergetic efficiency of the overall system

> 

varepsilon[TOT, `EN,CM`] := Units:-Standard:-`*`(E[P, TOT, `EN,CM`], Units:-Standard:-`/`(E[F, TOT, `EN,CM`]))

.4421818272

 (4.379)

Exergy destruction ratio in the compressor related to the exergy of total fuel

> 

y[CM, `EN,CM`] := Units:-Standard:-`*`(E[D, CM, EN, CM], Units:-Standard:-`/`(E[F, TOT, `EN,CM`]))

 

.1528094322

 (4.380)

 

Exergy destruction ratio in the condenser related to the exergy of total fuel

> 

y[CD, `EN,CM`] := Units:-Standard:-`*`(E[D, CD, `EN,CM`], Units:-Standard:-`/`(E[F, TOT, `EN,CM`]))

.1732952664

 (4.381)

Exergy destruction ratio in the throttling valve related to the exergy of total fuel

> 

y[TV, `EN,CM`] := Units:-Standard:-`*`(E[D, TV, `EN,CM`], Units:-Standard:-`/`(E[F, TOT, `EN,CM`]))

0.

 (4.382)

Exergy destruction ratio in the evaporator related to the exergy of total fuel

> 

y[EV, `EN,CM`] := Units:-Standard:-`*`(E[D, EV, `EN,CM`], Units:-Standard:-`/`(E[F, TOT, `EN,CM`]))

0.8307225088e-1

 (4.383)

Exergy destruction ratio in the system overall related to the exergy of total fuel

> 

y[TOT, `EN,CM`] := Units:-Standard:-`*`(E[D, TOT, `EN,CM`], Units:-Standard:-`/`(E[F, TOT, `EN,CM`]))

 

.4091769496

 (4.384)

 

Exergy destruction ratio in the compressor related to the total exergy destruction

> 

y[CM, `EN,CM`, 1] := Units:-Standard:-`*`(E[D, CM, EN, CM], Units:-Standard:-`/`(E[D, TOT, `EN,CM`]))

 

.3734556220

 (4.385)

 

Exergy destruction ratio in the condenser related to the total exergy destruction

> 

y[CD, `EN,CM`, 1] := Units:-Standard:-`*`(E[D, CD, `EN,CM`], Units:-Standard:-`/`(E[D, TOT, `EN,CM`]))

.4235215757

 (4.386)

Exergy destruction ratio in the throttling valve related to the total exergy destruction

> 

y[TV, `EN,CM`, 1] := Units:-Standard:-`*`(E[D, TV, `EN,CM`], Units:-Standard:-`/`(E[D, TOT, `EN,CM`]))

0.

 (4.387)

Exergy destruction ratio in the evaporator related to the total exergy destruction

> 

y[EV, `EN,CM`, 1] := Units:-Standard:-`*`(E[D, EV, `EN,CM`], Units:-Standard:-`/`(E[D, TOT, `EN,CM`]))

.2030228021

 (4.388)

Exergy destruction ratio in the system overall related to the total exergy destruction

> 

y[TOT, `EN,CM`, 1] := Units:-Standard:-`*`(E[D, TOT, `EN,CM`], Units:-Standard:-`/`(E[D, TOT, `EN,CM`]))

 

1.000000000

 (4.389)

 

5. Hybrid cycle with irreversible component - Condenser

Temperature of the working fluid at the compressor outlet

> 

T[4, EN, CD] := Units:-Standard:-`+`(T[EV, h, out], 0)

258.15*Units:-Unit(K)

 (4.390)

Pressure of the working fluid at the compressor outlet

> 

p[2, EN, CD] := PSTXwf(Units:-Standard:-`+`(T[CD, c, out], `&delta;T`[CD, h]), 1, wf[cycle])

 

1555421.989*Units:-Unit(Pa)

 (4.391)

 

Pressure of  the working fluid at the condenser outlet

> 

p[3, EN, CD] := p[2, EN, CD]

 

1555421.989*Units:-Unit(Pa)

 (4.392)

 

Pressure of  the working fluid at the condenser outlet

> 

T[3, `EN,CD`] := TSPXwf(p[3, EN, CD], 0, wf[cycle])

 

313.1500000*Units:-Unit(K)

 (4.393)

 

 

Specific enthalpy of  the working fluid at the condenser outlet

> 

h[3, `EN,CD`] := HSTXwf(T[3, `EN,CD`], 0, wf[cycle])

 

533791.2205*Units:-Unit(J/kg)

 (4.394)

 

Specific entropy of  the working fluid at the condenser outlet

> 

s[3, `EN,CD`] := SSTXwf(T[3, `EN,CD`], 0, wf[cycle])

 

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

 (4.395)

 

Specific volume of  the working fluid at the condenser outlet

> 

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

 

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

 (4.396)

 

Specific exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,CD`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, `EN,CD`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, `EN,CD`], Units:-Standard:-`-`(s[0])))))

 

299157.094*Units:-Unit(J/kg)

 (4.397)

 

Auxiliary values

> 

h[3, `EN,CD`, 1] := HPTwf(p[3, EN, CD], T[0], wf[cycle])

 

437261.1314*Units:-Unit(J/kg)

 (4.398)

 

> 

s[3, `EN,CD`, 1] := SPTwf(p[3, EN, CD], T[0], wf[cycle])

 

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

 (4.399)

 

Specific mechanical exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,CD`, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[3, `EN,CD`, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[3, `EN,CD`, 1]))))

 

295986.156*Units:-Unit(J/kg)

 (4.400)

 

Specific thermal exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,CD`, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, `EN,CD`], Units:-Standard:-`-`(h[3, `EN,CD`, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, `EN,CD`], Units:-Standard:-`-`(s[3, `EN,CD`, 1])))))

 

3170.93781*Units:-Unit(J/kg)

 (4.401)

 

Specific entropy of the working fluid at the evaporator intlet

> 

s[4, EN, CD] := s[3, `EN,CD`]

 

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

 (4.402)

 

Pressure of the working fluid at the evaporator intlet

 

> 

p[4, EN, CD] := Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`*`(PSTXwf(T[4, EN, CD], 0, wf[cycle]), Units:-Standard:-`+`(SSTXwf(T[4, EN, CD], 1, wf[cycle]), Units:-Standard:-`-`(s[4, EN, CD]))), Units:-Standard:-`*`(PSTXwf(T[4, EN, CD], 1, wf[cycle]), Units:-Standard:-`+`(s[4, EN, CD], Units:-Standard:-`-`(SSTXwf(T[4, EN, CD], 0, wf[cycle]))))), Units:-Standard:-`/`(Units:-Standard:-`+`(SSTXwf(T[4, EN, CD], 1, wf[cycle]), Units:-Standard:-`-`(SSTXwf(T[4, EN, CD], 0, wf[cycle])))))

 

236167.1499*Units:-Unit(Pa)

 (4.403)

 

Specific enthalpy of saturated liquid of  the working fluid in the evaporator

> 

h[sl, 4, `EN,CD`] := HSPXwf(p[4, EN, CD], 0, wf[cycle])

274370.1868*Units:-Unit(J/kg)

 (4.404)

Specific enthalpy of saturated vapor of  the working fluid in the evaporator

> 

h[sv, 4, `EN,CD`] := HSPXwf(p[4, EN, CD], 1, wf[cycle])

1587522.359*Units:-Unit(J/kg)

 (4.405)

Specific entropy of saturated liquid of  the working fluid in the evaporator

> 

s[sl, 4, EN, CD] := SSPXwf(p[4, EN, CD], 0, wf[cycle])

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

 (4.406)

Specific entropy of saturated vapor of  the working fluid in the evaporator

> 

s[sv, 4, EN, CD] := SSPXwf(p[4, EN, CD], 1, wf[cycle])

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

 (4.407)

Quality of the working fluid at the evaporator intlet

> 

x[4, `EN,CD`] := Units:-Standard:-`*`(Units:-Standard:-`+`(s[4, EN, CD], Units:-Standard:-`-`(s[sl, 4, EN, CD])), Units:-Standard:-`/`(Units:-Standard:-`+`(s[sv, 4, EN, CD], Units:-Standard:-`-`(s[sl, 4, EN, CD]))))NULL

 

.1774022976

 (4.408)

 

 

Specific entropy of  the working fluid at the evaporator inlet

> 

h[4, `EN,CD`] := Units:-Standard:-`+`(h[sl, 4, `EN,CD`], Units:-Standard:-`*`(x[4, `EN,CD`], Units:-Standard:-`+`(h[sv, 4, `EN,CD`], Units:-Standard:-`-`(h[sl, 4, `EN,CD`]))))

507326.3992*Units:-Unit(J/kg)

 (4.409)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

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

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

 (4.410)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

v[sv, 4, `EN,CD`] := Units:-Standard:-`/`(DSPXwf(p[4, EN, CD], 1, wf[cycle]))

.5086794091*Units:-Unit(m^3/kg)

 (4.411)

Specific volume of  the working fluid at the evaporator inlet

> 

v[4, `EN,CD`] := Units:-Standard:-`+`(v[sl, 4, `EN,CD`], Units:-Standard:-`*`(x[4, `EN,CD`], Units:-Standard:-`+`(v[sv, 4, `EN,CD`], Units:-Standard:-`-`(v[sl, 4, `EN,CD`]))))

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

 (4.412)

Specific exergy of  the working fluid at the evaporator inlet

> 

e[4, EN, CD] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, `EN,CD`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, EN, CD], Units:-Standard:-`-`(s[0])))))

 

272692.272*Units:-Unit(J/kg)

 (4.413)

 

Auxiliary values

> 

h[4, `EN,CD`, 1] := HPTwf(p[4, EN, CD], T[0], wf[cycle])

 

1670040.049*Units:-Unit(J/kg)

 (4.414)

 

> 

s[4, `EN,CD`, 1] := SPTwf(p[4, EN, CD], T[0], wf[cycle])

 

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

 (4.415)

 

Specific mechanical exergy of  the working fluid at the evaporator inlet

> 

e[4, `EN,CD`, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[4, `EN,CD`, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[4, `EN,CD`, 1]))))

 

120804.8342*Units:-Unit(J/kg)

 (4.416)

 

Specific thermal exergy of  the working fluid at the evaporator inlet

> 

e[4, `EN,CD`, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, `EN,CD`], Units:-Standard:-`-`(h[4, `EN,CD`, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, EN, CD], Units:-Standard:-`-`(s[4, `EN,CD`, 1])))))

 

151887.438*Units:-Unit(J/kg)

 (4.417)

 

Pressure of  the working fluid at the compressor inlet

> 

p[1, `EN,CD`] := p[4, EN, CD]

236167.1499*Units:-Unit(Pa)

 (4.418)

Temperature of  the working fluid at the compressor inlet

> 

T[1, `EN,CD`] := TSPXwf(p[1, `EN,CD`], 1, wf[cycle])

258.1500000*Units:-Unit(K)

 (4.419)

Specific enthalpy of  the working fluid at the compressor inlet

> 

h[1, `EN,CD`] := HSPXwf(p[1, `EN,CD`], 1, wf[cycle])

1587522.359*Units:-Unit(J/kg)

 (4.420)

Specific entropy of  the working fluid at the compressor inlet

> 

s[1, `EN,CD`] := SSPXwf(p[1, `EN,CD`], 1, wf[cycle])

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

 (4.421)

Specific volume of  the working fluid at the compressor inlet

> 

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

 

.5086794088*Units:-Unit(m^3/kg)

 (4.422)

 

Specific exergy of  the working fluid at the compressor inlet

> 

e[1, `EN,CD`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[1, `EN,CD`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[1, `EN,CD`], Units:-Standard:-`-`(s[0])))))

 

126239.2086*Units:-Unit(J/kg)

 (4.423)

 

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

> 

h[2, `EN,CD`] := HcompPTPEFFwf(p[1, `EN,CD`], T[1, `EN,CD`], p[2, EN, CD], 1, wf[cycle])

1869829.772*Units:-Unit(J/kg)

 (4.424)

Enhtalpy change in the compressor 1 after actual compression

> 

l[comp, `EN,CD`] := Units:-Standard:-`+`(h[2, `EN,CD`], Units:-Standard:-`-`(h[1, `EN,CD`]))

 

282307.413*Units:-Unit(J/kg)

 (4.425)

 

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

> 

T[2, `EN,CD`] := TPHwf(p[2, EN, CD], h[2, `EN,CD`], wf[cycle])

 

396.8264918*Units:-Unit(K)

 (4.426)

 

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

> 

s[2, `EN,CD`] := SPTwf(p[2, EN, CD], T[2, `EN,CD`], wf[cycle])

 

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

 (4.427)

 

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

> 

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

.1174758710*Units:-Unit(m^3/kg)

 (4.428)

Specific exergy of  the working fluid at the compressor inlet

> 

e[2, `EN,CD`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[2, `EN,CD`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[2, `EN,CD`], Units:-Standard:-`-`(s[0])))))

 

408546.6219*Units:-Unit(J/kg)

 (4.429)

 

 

Specific heat rejection in the condenser

> 

q[out, `EN,CD`] := Units:-Standard:-`+`(h[2, `EN,CD`], Units:-Standard:-`-`(h[3, `EN,CD`]))

1336038.552*Units:-Unit(J/kg)

 (4.430)

Specific heat addition in the evaporator

> 

q[inlet, `EN,CD`] := Units:-Standard:-`+`(h[1, `EN,CD`], Units:-Standard:-`-`(h[4, `EN,CD`]))

1080195.960*Units:-Unit(J/kg)

 (4.431)

Coefficent of performance of a refrigerator

> 

COP[`EN,CD`, Ref] := Units:-Standard:-`*`(q[inlet, `EN,CD`], Units:-Standard:-`/`(l[comp, `EN,CD`]))

3.826311001

 (4.432)

Coefficent of performance of a heat pump

> 

COP[`EN,CD`, HP] := Units:-Standard:-`*`(q[out, `EN,CD`], Units:-Standard:-`/`(l[comp, `EN,CD`]))

4.732566311

 (4.433)

Mass rate of the working fluid

> 

m[wf, cycle, `EN,CD`] := proc (Type) if Type = 1 then Q[P]/q[out, `EN,CD`] else Q[P]/q[inlet, `EN,CD`] end if end proc:

 

> 

m[wf, cycle, `EN,CD`] := m[wf, cycle, `EN,CD`](Type)

 

0.9257579523e-4*kW*Units:-Unit(kg/J)

 (4.434)

 

Total heat addition in the evaporator

> 

Q[inlet, `EN,CD`] := Units:-Standard:-`*`(m[wf, cycle, `EN,CD`], q[inlet, `EN,CD`])

 

100.0000000*kW

 (4.435)

 

Total heat rejection in the condenser

> 

Q[out, `EN,CD`] := Units:-Standard:-`*`(m[wf, cycle, R], q[out, `EN,CD`])

 

128.4633901*kW

 (4.436)

 

Mass rate of the heat source fluid

> 

m[EV, h, `EN,CD`] := Units:-Standard:-`*`(Q[inlet, `EN,CD`], Units:-Standard:-`/`(Units:-Standard:-`+`(h[EV, h, inlet], Units:-Standard:-`-`(h[EV, h, out]))))

 

0.9944845786e-2*kW*Units:-Unit(kg/J)

 (4.437)

 

Mass rate of the heat consumption fluid

> 

m[CD, c, `EN,CD`] := Units:-Standard:-`*`(Q[out, `EN,CD`], Units:-Standard:-`/`(Units:-Standard:-`+`(h[CD, c, out], Units:-Standard:-`-`(h[CD, c, inlet]))))

 

0.1276603757e-1*kW*Units:-Unit(kg/J)

 (4.438)

 

Total exergy used in the compressor

> 

E[5, `EN,CD`] := Units:-Standard:-`*`(m[wf, cycle, `EN,CD`], l[comp, `EN,CD`])

 

26.13483326*kW

 (4.439)

 

Total physical exergy at the compressor intlet

> 

E[1, EN, CD] := Units:-Standard:-`*`(m[wf, cycle, `EN,CD`], e[1, `EN,CD`])

 

11.68669513*kW

 (4.440)

 

Total physical exergy at the compressor outlet

> 

E[2, `EN,CD`] := Units:-Standard:-`*`(m[wf, cycle, `EN,CD`], e[2, `EN,CD`])

 

37.82152841*kW

 (4.441)

 

Total physical exergy at the condenser outlet

> 

E[3, `EN,CD`] := Units:-Standard:-`*`(m[wf, cycle, `EN,CD`], e[3, `EN,CD`])

 

27.69470588*kW

 (4.442)

 

Total mechanical exergy at the condenser outlet

> 

E[3, `EN,CD`, M] := Units:-Standard:-`*`(m[wf, cycle, `EN,CD`], e[3, `EN,CD`, M])

 

27.40115377*kW

 (4.443)

 

Total thermal exergy at the condenser outlet

> 

E[3, `EN,CD`, T] := Units:-Standard:-`*`(m[wf, cycle, `EN,CD`], e[3, `EN,CD`, T])

 

.2935520894*kW

 (4.444)

 

Total physical exergy at the evaporator inlet

> 

E[4, EN, CD] := Units:-Standard:-`*`(m[wf, cycle, `EN,CD`], e[4, EN, CD])

 

25.24470393*kW

 (4.445)

 

Total mechanical exergy at the evaporator inlet

> 

E[4, EN, CD, M] := Units:-Standard:-`*`(m[wf, cycle, `EN,CD`], e[4, `EN,CD`, M])

 

11.18360359*kW

 (4.446)

 

Total thermal exergy at the evaporator inlet

> 

E[4, `EN,CD`, T] := Units:-Standard:-`*`(m[wf, cycle, `EN,CD`], e[4, `EN,CD`, T])

 

14.06110036*kW

 (4.447)

 

Total physical exergy of heat consumption fluid at the evaporator intlet 

> 

E[EV, h, inlet] := Units:-Standard:-`*`(m[EV, h, `EN,CD`], e[EV, h, inlet])

 

11.30975886*kW

 (4.448)

 

Total physical exergy of heat consumption fluid at the evaporator outlet

> 

E[EV, h, out] := Units:-Standard:-`*`(m[EV, h, `EN,CD`], e[EV, h, out])

 

22.72348304*kW

 (4.449)

 

Total physical exergy of heat consumption fluid at the condenser intlet 

> 

E[CD, c, inlet] := Units:-Standard:-`*`(m[CD, c, `EN,CD`], e[CD, c, inlet])

 

-0.

 (4.450)

 

Total physical exergy of heat consumption fluid at the condenser outlet

> 

E[CD, c, out] := Units:-Standard:-`*`(m[CD, c, `EN,CD`], e[CD, c, out])

 

2.142620635*kW

 (4.451)

 

Exergy of fuel in the compressor

> 

E[F, CM, `EN,CD`] := E[5, `EN,CD`]

 

26.13483326*kW

 (4.452)

 

Exergy of product in the compressor

> 

E[P, CM, `EN,CD`] := Units:-Standard:-`+`(E[2, `EN,CD`], Units:-Standard:-`-`(E[1, EN, CD]))

 

26.13483328*kW

 (4.453)

 

Exergy destruction in the compressor

> 

E[D, CM, `EN,CD`] := Units:-Standard:-`+`(E[F, CM, `EN,CD`], Units:-Standard:-`-`(E[P, CM, `EN,CD`]))

 

-0.2e-7*kW

 (4.454)

 

Exergetic efficiency in the compressor

> 

varepsilon[CM, `EN,CD`] := Units:-Standard:-`*`(E[P, CM, `EN,CD`], Units:-Standard:-`/`(E[F, CM, `EN,CD`]))

 

1.000000001

 (4.455)

 

Exergy of fuel in the condenser

> 

E[F, CD, `EN,CD`] := Units:-Standard:-`+`(E[2, `EN,CD`], Units:-Standard:-`-`(E[3, `EN,CD`]))

 

10.12682253*kW

 (4.456)

 

Exergy of product in the condenser

> 

E[P, CD, `EN,CD`] := Units:-Standard:-`+`(E[CD, c, out], Units:-Standard:-`-`(E[CD, c, inlet]))

 

2.142620635*kW

 (4.457)

 

Exergy destruction in the condenser

> 

E[D, CD, EN, CD] := Units:-Standard:-`+`(E[F, CD, `EN,CD`], Units:-Standard:-`-`(E[P, CD, `EN,CD`]))

 

7.984201895*kW

 (4.458)

 

Exergy of fuel in the throttling valve

> 

E[F, TV, `EN,CD`] := Units:-Standard:-`+`(E[3, `EN,CD`, M], Units:-Standard:-`-`(E[4, EN, CD, M]))

 

16.21755018*kW

 (4.459)

 

Exergy of product in the throttling valve

> 

E[P, TV, `EN,CD`] := Units:-Standard:-`+`(E[3, `EN,CD`, M], Units:-Standard:-`-`(E[4, EN, CD, M]))

 

16.21755018*kW

 (4.460)

 

Exergy destruction in the throttling valve

> 

E[D, TV, `EN,CD`] := Units:-Standard:-`+`(E[F, TV, `EN,CD`], Units:-Standard:-`-`(E[P, TV, `EN,CD`]))

 

0.

 (4.461)

 

Exergetic efficiency in the throttling valve       

> 

varepsilon[TV, `EN,CD`] := Units:-Standard:-`*`(E[P, TV, `EN,CD`], Units:-Standard:-`/`(E[F, TV, `EN,CD`]))

 

1.000000000

 (4.462)

 

Exergy of fuel in the evaporator

> 

E[F, EV, `EN,CD`] := Units:-Standard:-`+`(E[4, EN, CD], Units:-Standard:-`-`(E[1, EN, CD]))

 

13.55800880*kW

 (4.463)

 

Exergy of product in the evaporator

> 

E[P, EV, `EN,CD`] := Units:-Standard:-`+`(E[EV, h, out], Units:-Standard:-`-`(E[EV, h, inlet]))

 

11.41372418*kW

 (4.464)

 

Exergy destruction in the evaporator

> 

E[D, EV, `EN,CD`] := Units:-Standard:-`+`(E[F, EV, `EN,CD`], Units:-Standard:-`-`(E[P, EV, `EN,CD`]))

 

2.14428462*kW

 (4.465)

 

Exergetic efficiency in the evaporator

> 

varepsilon[EV, `EN,CD`] := Units:-Standard:-`*`(E[P, EV, `EN,CD`], Units:-Standard:-`/`(E[F, EV, `EN,CD`]))

 

.8418436916

 (4.466)

 

Exergy of fuel of the overall system

> 

E[F, TOT, `EN,CD`] := E[5, `EN,CD`]

26.13483326*kW

 (4.467)

Exergy of product of the overall system 

> 

E[P, TOT, `EN,CD`] := proc (Type) if Type = 1 then E[P, CD, `EN,CD`] else E[P, EV, `EN,CD`] end if end proc:

 

> 

E[P, TOT, `EN,CD`] := E[P, TOT, `EN,CD`](Type)

 

11.41372418*kW

 (4.468)

 

Exergy destruction of the overall system

> 

E[D, TOT, `EN,CD`] := Units:-Standard:-`+`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CM, `EN,CD`], E[D, CD, EN, CD]), E[D, TV, `EN,CD`]), E[D, EV, `EN,CD`])

10.12848650*kW

 (4.469)

Exergetic efficiency of the overall system

> 

varepsilon[TOT, `EN,CD`] := Units:-Standard:-`*`(E[P, TOT, `EN,CD`], Units:-Standard:-`/`(E[F, TOT, `EN,CD`]))

.4367245839

 (4.470)

Exergy destruction ratio in the compressor related to the exergy of total fuel

> 

y[CM, `EN,CD`] := Units:-Standard:-`*`(E[D, CM, `EN,CD`], Units:-Standard:-`/`(E[F, TOT, `EN,CD`]))

 

-0.7652622002e-9

 (4.471)

 

Exergy destruction ratio in the condenser related to the exergy of total fuel

> 

y[CD, `EN,CD`] := Units:-Standard:-`*`(E[D, CD, EN, CD], Units:-Standard:-`/`(E[F, TOT, `EN,CD`]))

.3055003955

 (4.472)

Exergy destruction ratio in the throttling valve related to the exergy of total fuel

> 

y[TV, `EN,CD`] := Units:-Standard:-`*`(E[D, TV, `EN,CD`], Units:-Standard:-`/`(E[F, TOT, `EN,CD`]))

0.

 (4.473)

Exergy destruction ratio in the evaporator related to the exergy of total fuel

> 

y[EV, `EN,CD`] := Units:-Standard:-`*`(E[D, EV, `EN,CD`], Units:-Standard:-`/`(E[F, TOT, `EN,CD`]))

0.8204699831e-1

 (4.474)

Exergy destruction ratio in the system overall related to the exergy of total fuel

> 

y[TOT, `EN,CD`] := Units:-Standard:-`*`(E[D, TOT, `EN,CD`], Units:-Standard:-`/`(E[F, TOT, `EN,CD`]))

 

.3875473932

 (4.475)

 

Exergy destruction ratio in the compressor related to the total exergy destruction

> 

y[CM, `EN,CD`, 1] := Units:-Standard:-`*`(E[D, CM, `EN,CD`], Units:-Standard:-`/`(E[D, TOT, `EN,CD`]))

 

-0.1974628687e-8

 (4.476)

 

Exergy destruction ratio in the condenser related to the total exergy destruction

> 

y[CD, `EN,CD`, 1] := Units:-Standard:-`*`(E[D, CD, EN, CD], Units:-Standard:-`/`(E[D, TOT, `EN,CD`]))

.7882917053

 (4.477)

Exergy destruction ratio in the throttling valve related to the total exergy destruction

> 

y[TV, `EN,CD`, 1] := Units:-Standard:-`*`(E[D, TV, `EN,CD`], Units:-Standard:-`/`(E[D, TOT, `EN,CD`]))

0.

 (4.478)

Exergy destruction ratio in the evaporator related to the total exergy destruction

> 

y[EV, `EN,CD`, 1] := Units:-Standard:-`*`(E[D, EV, `EN,CD`], Units:-Standard:-`/`(E[D, TOT, `EN,CD`]))

.2117082962

 (4.479)

Exergy destruction ratio in the system overall related to the total exergy destruction

> 

y[TOT, `EN,CD`, 1] := Units:-Standard:-`*`(E[D, TOT, `EN,CD`], Units:-Standard:-`/`(E[D, TOT, `EN,CD`]))

 

1.000000000

 (4.480)

 

6. Hybrid cycle with irreversible component - Throttling Valve

Temperature of the working fluid at the compressor outlet

> 

T[4, EN, TV] := Units:-Standard:-`+`(T[EV, h, out], 0)

258.15*Units:-Unit(K)

 (4.481)

Pressure of the working fluid at the compressor outlet

> 

p[2, EN, TV] := PSTXwf(Units:-Standard:-`+`(T[CD, c, out], 0), 1, wf[cycle])

 

1167197.054*Units:-Unit(Pa)

 (4.482)

 

Pressure of  the working fluid at the condenser outlet

> 

p[3, EN, TV] := p[2, EN, TV]

 

1167197.054*Units:-Unit(Pa)

 (4.483)

 

Pressure of  the working fluid at the condenser outlet

> 

T[3, `EN,TV`] := TSPXwf(p[3, EN, TV], 0, wf[cycle])

 

303.1500000*Units:-Unit(K)

 (4.484)

 

 

Specific enthalpy of  the working fluid at the condenser outlet

> 

h[3, `EN,TV`] := HSTXwf(T[3, `EN,TV`], 0, wf[cycle])

 

484911.4924*Units:-Unit(J/kg)

 (4.485)

 

Specific entropy of  the working fluid at the condenser outlet

> 

s[3, `EN,TV`] := SSTXwf(T[3, `EN,TV`], 0, wf[cycle])

 

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

 (4.486)

 

Specific volume of  the working fluid at the condenser outlet

> 

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

 

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

 (4.487)

 

Specific exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,TV`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, `EN,TV`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, `EN,TV`], Units:-Standard:-`-`(s[0])))))

 

296150.505*Units:-Unit(J/kg)

 (4.488)

 

Auxiliary values

> 

h[3, `EN,TV`, 1] := HPTwf(p[3, EN, TV], T[0], wf[cycle])

 

437080.3036*Units:-Unit(J/kg)

 (4.489)

 

> 

s[3, `EN,TV`, 1] := SPTwf(p[3, EN, TV], T[0], wf[cycle])

 

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

 (4.490)

 

Specific mechanical exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,TV`, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[3, `EN,TV`, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[3, `EN,TV`, 1]))))

 

295350.397*Units:-Unit(J/kg)

 (4.491)

 

Specific thermal exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,TV`, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, `EN,TV`], Units:-Standard:-`-`(h[3, `EN,TV`, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, `EN,TV`], Units:-Standard:-`-`(s[3, `EN,TV`, 1])))))

 

800.10886*Units:-Unit(J/kg)

 (4.492)

 

Specific enthalpy of the working fluid at the evaporator intlet

> 

h[4, EN, TV] := h[3, `EN,TV`]

 

484911.4924*Units:-Unit(J/kg)

 (4.493)

 

Pressure of the working fluid at the evaporator intlet

 

> 

p[4, EN, TV] := Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`*`(PSTXwf(T[4, EN, TV], 0, wf[cycle]), Units:-Standard:-`+`(HSTXwf(T[4, EN, TV], 1, wf[cycle]), Units:-Standard:-`-`(h[4, EN, TV]))), Units:-Standard:-`*`(PSTXwf(T[4, EN, TV], 1, wf[cycle]), Units:-Standard:-`+`(h[4, EN, TV], Units:-Standard:-`-`(HSTXwf(T[4, EN, TV], 0, wf[cycle]))))), Units:-Standard:-`/`(Units:-Standard:-`+`(HSTXwf(T[4, EN, TV], 1, wf[cycle]), Units:-Standard:-`-`(HSTXwf(T[4, EN, TV], 0, wf[cycle])))))

 

236167.1501*Units:-Unit(Pa)

 (4.494)

 

Specific enthalpy of saturated liquid of  the working fluid in the evaporator

> 

h[sl, 4, EN, TV] := HSPXwf(p[4, EN, TV], 0, wf[cycle])

274370.1868*Units:-Unit(J/kg)

 (4.495)

Specific enthalpy of saturated vapor of  the working fluid in the evaporator

> 

h[sv, 4, EN, TV] := HSPXwf(p[4, EN, TV], 1, wf[cycle])

1587522.359*Units:-Unit(J/kg)

 (4.496)

Quality of the working fluid at the evaporator intlet

> 

x[4, `EN,TV`] := Units:-Standard:-`*`(Units:-Standard:-`+`(h[4, EN, TV], Units:-Standard:-`-`(h[sl, 4, EN, TV])), Units:-Standard:-`/`(Units:-Standard:-`+`(h[sv, 4, EN, TV], Units:-Standard:-`-`(h[sl, 4, EN, TV]))))``

 

.1603327551

 (4.497)

 

Specific entropy of saturated liquid of  the working fluid in the evaporator

> 

s[sl, 4, `EN,TV`] := SSPXwf(p[4, EN, TV], 0, wf[cycle])

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

 (4.498)

Specific entropy of saturated vapor of  the working fluid in the evaporator

> 

s[sv, 4, `EN,TV`] := SSPXwf(p[4, EN, TV], 1, wf[cycle])

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

 (4.499)

Specific entropy of  the working fluid at the evaporator inlet

> 

s[4, `EN,TV`] := Units:-Standard:-`+`(s[sl, 4, `EN,TV`], Units:-Standard:-`*`(x[4, `EN,TV`], Units:-Standard:-`+`(s[sv, 4, `EN,TV`], Units:-Standard:-`-`(s[sl, 4, `EN,TV`]))))

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

 (4.500)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

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

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

 (4.501)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

v[sv, 4, `EN,TV`] := Units:-Standard:-`/`(DSPXwf(p[4, EN, TV], 1, wf[cycle]))

.5086794088*Units:-Unit(m^3/kg)

 (4.502)

Specific volume of  the working fluid at the evaporator inlet

> 

v[4, `EN,TV`] := Units:-Standard:-`+`(v[sl, 4, `EN,TV`], Units:-Standard:-`*`(x[4, `EN,TV`], Units:-Standard:-`+`(v[sv, 4, `EN,TV`], Units:-Standard:-`-`(v[sl, 4, `EN,TV`]))))

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

 (4.503)

Specific exergy of  the working fluid at the evaporator inlet

> 

e[4, EN, TV] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, EN, TV], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, `EN,TV`], Units:-Standard:-`-`(s[0])))))

 

275731.287*Units:-Unit(J/kg)

 (4.504)

 

Auxiliary values

> 

h[4, `EN,TV`, 1] := HPTwf(p[4, EN, TV], T[0], wf[cycle])

 

1670040.049*Units:-Unit(J/kg)

 (4.505)

 

> 

s[4, `EN,TV`, 1] := SPTwf(p[4, EN, TV], T[0], wf[cycle])

 

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

 (4.506)

 

Specific mechanical exergy of  the working fluid at the evaporator inlet

> 

e[4, `EN,TV`, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[4, `EN,TV`, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[4, `EN,TV`, 1]))))

 

120804.8345*Units:-Unit(J/kg)

 (4.507)

 

Specific thermal exergy of  the working fluid at the evaporator inlet

> 

e[4, `EN,TV`, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, EN, TV], Units:-Standard:-`-`(h[4, `EN,TV`, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, `EN,TV`], Units:-Standard:-`-`(s[4, `EN,TV`, 1])))))

 

154926.453*Units:-Unit(J/kg)

 (4.508)

 

Pressure of  the working fluid at the compressor inlet

> 

p[1, `EN,TV`] := p[4, EN, TV]

236167.1501*Units:-Unit(Pa)

 (4.509)

Temperature of  the working fluid at the compressor inlet

> 

T[1, `EN,TV`] := TSPXwf(p[1, `EN,TV`], 1, wf[cycle])

258.1500000*Units:-Unit(K)

 (4.510)

Specific enthalpy of  the working fluid at the compressor inlet

> 

h[1, EN, TV] := HSPXwf(p[1, `EN,TV`], 1, wf[cycle])

1587522.359*Units:-Unit(J/kg)

 (4.511)

Specific entropy of  the working fluid at the compressor inlet

> 

s[1, `EN,TV`] := SSPXwf(p[1, `EN,TV`], 1, wf[cycle])

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

 (4.512)

Specific volume of  the working fluid at the compressor inlet

> 

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

 

.5086794088*Units:-Unit(m^3/kg)

 (4.513)

 

Specific exergy of  the working fluid at the compressor inlet

> 

e[1, `EN,TV`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[1, EN, TV], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[1, `EN,TV`], Units:-Standard:-`-`(s[0])))))

 

126239.2086*Units:-Unit(J/kg)

 (4.514)

 

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

> 

h[2, `EN,TV`] := HcompPTPEFFwf(p[1, `EN,TV`], T[1, `EN,TV`], p[2, EN, TV], 1, wf[cycle])

1818922.733*Units:-Unit(J/kg)

 (4.515)

Enhtalpy change in the compressor 1 after actual compression

> 

l[comp, `EN,TV`] := Units:-Standard:-`+`(h[2, `EN,TV`], Units:-Standard:-`-`(h[1, EN, TV]))

 

231400.374*Units:-Unit(J/kg)

 (4.516)

 

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

> 

T[2, `EN,TV`] := TPHwf(p[2, EN, TV], h[2, `EN,TV`], wf[cycle])

 

372.2312166*Units:-Unit(K)

 (4.517)

 

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

> 

s[2, `EN,TV`] := SPTwf(p[2, EN, TV], T[2, `EN,TV`], wf[cycle])

 

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

 (4.518)

 

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

> 

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

.1472972458*Units:-Unit(m^3/kg)

 (4.519)

Specific exergy of  the working fluid at the compressor inlet

> 

e[2, `EN,TV`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[2, `EN,TV`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[2, `EN,TV`], Units:-Standard:-`-`(s[0])))))

 

357639.5831*Units:-Unit(J/kg)

 (4.520)

 

 

Specific heat rejection in the condenser

> 

q[out, `EN,TV`] := Units:-Standard:-`+`(h[2, `EN,TV`], Units:-Standard:-`-`(h[3, `EN,TV`]))

1334011.241*Units:-Unit(J/kg)

 (4.521)

Specific heat addition in the evaporator

> 

q[inlet, `EN,TV`] := Units:-Standard:-`+`(h[1, EN, TV], Units:-Standard:-`-`(h[4, EN, TV]))

1102610.867*Units:-Unit(J/kg)

 (4.522)

Coefficent of performance of a refrigerator

> 

COP[`EN,TV`, Ref] := Units:-Standard:-`*`(q[inlet, `EN,TV`], Units:-Standard:-`/`(l[comp, `EN,TV`]))

4.764948509

 (4.523)

Coefficent of performance of a heat pump

> 

COP[`EN,TV`, HP] := Units:-Standard:-`*`(q[out, `EN,TV`], Units:-Standard:-`/`(l[comp, `EN,TV`]))

5.764948509

 (4.524)

Mass rate of the working fluid

> 

m[wf, cycle, `EN,TV`] := proc (Type) if Type = 1 then Q[P]/q[out, `EN,TV`] else Q[P]/q[inlet, `EN,TV`] end if end proc:

 

> 

m[wf, cycle, `EN,TV`] := m[wf, cycle, `EN,TV`](Type)

 

0.9069382771e-4*kW*Units:-Unit(kg/J)

 (4.525)

 

Total heat addition in the evaporator

> 

Q[inlet, `EN,TV`] := Units:-Standard:-`*`(m[wf, cycle, `EN,TV`], q[inlet, `EN,TV`])

 

100.0000000*kW

 (4.526)

 

Total heat rejection in the condenser

> 

Q[out, `EN,TV`] := Units:-Standard:-`*`(m[wf, cycle, R], q[out, `EN,TV`])

 

128.2684591*kW

 (4.527)

 

Mass rate of the heat source fluid

> 

m[EV, h, `EN,TV`] := Units:-Standard:-`*`(Q[inlet, `EN,TV`], Units:-Standard:-`/`(Units:-Standard:-`+`(h[EV, h, inlet], Units:-Standard:-`-`(h[EV, h, out]))))

 

0.9944845786e-2*kW*Units:-Unit(kg/J)

 (4.528)

 

Mass rate of the heat consumption fluid

> 

m[CD, c, `EN,TV`] := Units:-Standard:-`*`(Q[out, `EN,TV`], Units:-Standard:-`/`(Units:-Standard:-`+`(h[CD, c, out], Units:-Standard:-`-`(h[CD, c, inlet]))))

 

0.1274666632e-1*kW*Units:-Unit(kg/J)

 (4.529)

 

Total exergy used in the compressor

> 

E[5, `EN,TV`] := Units:-Standard:-`*`(m[wf, cycle, `EN,TV`], l[comp, `EN,TV`])

 

20.98658565*kW

 (4.530)

 

Total physical exergy at the compressor intlet

> 

E[1, EN, TV] := Units:-Standard:-`*`(m[wf, cycle, `EN,TV`], e[1, `EN,TV`])

 

11.44911704*kW

 (4.531)

 

Total physical exergy at the compressor outlet

> 

E[2, `EN,TV`] := Units:-Standard:-`*`(m[wf, cycle, `EN,TV`], e[2, `EN,TV`])

 

32.43570273*kW

 (4.532)

 

Total physical exergy at the condenser outlet

> 

E[3, `EN,TV`] := Units:-Standard:-`*`(m[wf, cycle, `EN,TV`], e[3, `EN,TV`])

 

26.85902288*kW

 (4.533)

 

Total mechanical exergy at the condenser outlet

> 

E[3, `EN,TV`, M] := Units:-Standard:-`*`(m[wf, cycle, `EN,TV`], e[3, `EN,TV`, M])

 

26.78645802*kW

 (4.534)

 

Total thermal exergy at the condenser outlet

> 

E[3, `EN,TV`, T] := Units:-Standard:-`*`(m[wf, cycle, `EN,TV`], e[3, `EN,TV`, T])

 

0.7256493510e-1*kW

 (4.535)

 

Total physical exergy at the evaporator inlet

> 

E[4, EN, TV] := Units:-Standard:-`*`(m[wf, cycle, `EN,TV`], e[4, EN, TV])

 

25.00712584*kW

 (4.536)

 

Total mechanical exergy at the evaporator inlet

> 

E[4, `EN,TV`, M] := Units:-Standard:-`*`(m[wf, cycle, `EN,TV`], e[4, `EN,TV`, M])

 

10.95625285*kW

 (4.537)

 

Total thermal exergy at the evaporator inlet

> 

E[4, `EN,TV`, T] := Units:-Standard:-`*`(m[wf, cycle, `EN,TV`], e[4, `EN,TV`, T])

 

14.05087304*kW

 (4.538)

 

Total physical exergy of heat consumption fluid at the evaporator intlet 

> 

E[EV, h, inlet] := Units:-Standard:-`*`(m[EV, h, `EN,TV`], e[EV, h, inlet])

 

11.30975886*kW

 (4.539)

 

Total physical exergy of heat consumption fluid at the evaporator outlet

> 

E[EV, h, out] := Units:-Standard:-`*`(m[EV, h, `EN,TV`], e[EV, h, out])

 

22.72348304*kW

 (4.540)

 

Total physical exergy of heat consumption fluid at the condenser intlet 

> 

E[CD, c, inlet] := Units:-Standard:-`*`(m[CD, c, `EN,TV`], e[CD, c, inlet])

 

-0.

 (4.541)

 

Total physical exergy of heat consumption fluid at the condenser outlet

> 

E[CD, c, out] := Units:-Standard:-`*`(m[CD, c, `EN,TV`], e[CD, c, out])

 

2.139369412*kW

 (4.542)

 

Exergy of fuel in the compressor

> 

E[F, CM, `EN,TV`] := E[5, `EN,TV`]

 

20.98658565*kW

 (4.543)

 

Exergy of product in the compressor

> 

E[P, CM, `EN,TV`] := Units:-Standard:-`+`(E[2, `EN,TV`], Units:-Standard:-`-`(E[1, EN, TV]))

 

20.98658569*kW

 (4.544)

 

Exergy destruction in the compressor

> 

E[D, CM, `EN,TV`] := Units:-Standard:-`+`(E[F, CM, `EN,TV`], Units:-Standard:-`-`(E[P, CM, `EN,TV`]))

 

-0.4e-7*kW

 (4.545)

 

Exergetic efficiency in the compressor

> 

varepsilon[CM, `EN,TV`] := Units:-Standard:-`*`(E[P, CM, `EN,TV`], Units:-Standard:-`/`(E[F, CM, `EN,TV`]))

 

1.000000002

 (4.546)

 

Exergy of fuel in the condenser

> 

E[F, CD, `EN,TV`] := Units:-Standard:-`+`(E[2, `EN,TV`], Units:-Standard:-`-`(E[3, `EN,TV`]))

 

5.57667985*kW

 (4.547)

 

Exergy of product in the condenser

> 

E[P, CD, `EN,TV`] := Units:-Standard:-`+`(E[CD, c, out], Units:-Standard:-`-`(E[CD, c, inlet]))

 

2.139369412*kW

 (4.548)

 

Exergy destruction in the condenser

> 

E[D, CD, `EN,TV`] := Units:-Standard:-`+`(E[F, CD, `EN,TV`], Units:-Standard:-`-`(E[P, CD, `EN,TV`]))

 

3.437310438*kW

 (4.549)

 

Exergy of fuel in the throttling valve

> 

E[F, TV, `EN,TV`] := Units:-Standard:-`+`(E[3, `EN,TV`, M], Units:-Standard:-`-`(E[4, `EN,TV`, M]))

 

15.83020517*kW

 (4.550)

 

Exergy of product in the throttling valve

> 

E[P, TV, `EN,TV`] := Units:-Standard:-`+`(E[4, `EN,TV`, T], Units:-Standard:-`-`(E[3, `EN,TV`, T]))

 

13.97830810*kW

 (4.551)

 

Exergy destruction in the throttling valve

> 

E[D, TV, EN, TV] := Units:-Standard:-`+`(E[F, TV, `EN,TV`], Units:-Standard:-`-`(E[P, TV, `EN,TV`]))

 

1.85189707*kW

 (4.552)

 

Exergetic efficiency in the throttling valve 

> 

varepsilon[TV, `EN,TV`] := Units:-Standard:-`*`(E[P, TV, `EN,TV`], Units:-Standard:-`/`(E[F, TV, `EN,TV`]))

 

.8830149672

 (4.553)

 

Exergy of fuel in the evaporator

> 

E[F, EV, `EN,TV`] := Units:-Standard:-`+`(E[4, EN, TV], Units:-Standard:-`-`(E[1, EN, TV]))

 

13.55800880*kW

 (4.554)

 

Exergy of product in the evaporator

> 

E[P, EV, `EN,TV`] := Units:-Standard:-`+`(E[EV, h, out], Units:-Standard:-`-`(E[EV, h, inlet]))

 

11.41372418*kW

 (4.555)

 

Exergy destruction in the evaporator

> 

E[D, EV, EN, TV] := Units:-Standard:-`+`(E[F, EV, `EN,TV`], Units:-Standard:-`-`(E[P, EV, `EN,TV`]))

 

2.14428462*kW

 (4.556)

 

Exergetic efficiency in the evaporator

> 

varepsilon[EV, `EN,TV`] := Units:-Standard:-`*`(E[P, EV, `EN,TV`], Units:-Standard:-`/`(E[F, EV, `EN,TV`]))

 

.8418436916

 (4.557)

 

Exergy of fuel of the overall system

> 

E[F, TOT, `EN,TV`] := E[5, `EN,TV`]

20.98658565*kW

 (4.558)

Exergy of product of the overall system 

> 

E[P, TOT, `EN,TV`] := proc (Type) if Type = 1 then E[P, CD, `EN,TV`] else E[P, EV, `EN,TV`] end if end proc:

 

> 

E[P, TOT, `EN,TV`] := E[P, TOT, `EN,TV`](Type)

 

11.41372418*kW

 (4.559)

 

Exergy destruction of the overall system

> 

E[D, TOT, `EN,TV`] := Units:-Standard:-`+`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CM, `EN,TV`], E[D, CD, `EN,TV`]), E[D, TV, EN, TV]), E[D, EV, EN, TV])

7.433492088*kW

 (4.560)

Exergetic efficiency of the overall system

> 

varepsilon[TOT, `EN,TV`] := Units:-Standard:-`*`(E[P, TOT, `EN,TV`], Units:-Standard:-`/`(E[F, TOT, `EN,TV`]))

.5438580801

 (4.561)

Exergy destruction ratio in the compressor related to the exergy of total fuel

> 

y[CM, `EN,TV`] := Units:-Standard:-`*`(E[D, CM, `EN,TV`], Units:-Standard:-`/`(E[F, TOT, `EN,TV`]))

 

-0.1905979404e-8

 (4.562)

 

Exergy destruction ratio in the condenser related to the exergy of total fuel

> 

y[CD, `EN,TV`] := Units:-Standard:-`*`(E[D, CD, `EN,TV`], Units:-Standard:-`/`(E[F, TOT, `EN,TV`]))

.1637860725

 (4.563)

Exergy destruction ratio in the throttling valve related to the exergy of total fuel

> 

y[TV, `EN,TV`] := Units:-Standard:-`*`(E[D, TV, EN, TV], Units:-Standard:-`/`(E[F, TOT, `EN,TV`]))

0.8824194183e-1

 (4.564)

Exergy destruction ratio in the evaporator related to the exergy of total fuel

> 

y[EV, `EN,TV`] := Units:-Standard:-`*`(E[D, EV, EN, TV], Units:-Standard:-`/`(E[F, TOT, `EN,TV`]))

.1021740580

 (4.565)

Exergy destruction ratio in the system overall related to the exergy of total fuel

> 

y[TOT, `EN,TV`] := Units:-Standard:-`*`(E[D, TOT, `EN,TV`], Units:-Standard:-`/`(E[F, TOT, `EN,TV`]))

 

.3542020704

 (4.566)

 

Exergy destruction ratio in the compressor related to the total exergy destruction

> 

y[CM, `EN,TV`, 1] := Units:-Standard:-`*`(E[D, CM, `EN,TV`], Units:-Standard:-`/`(E[D, TOT, `EN,TV`]))

 

-0.5381050996e-8

 (4.567)

 

Exergy destruction ratio in the condenser related to the total exergy destruction

> 

y[CD, `EN,TV`, 1] := Units:-Standard:-`*`(E[D, CD, `EN,TV`], Units:-Standard:-`/`(E[D, TOT, `EN,TV`]))

.4624085689

 (4.568)

Exergy destruction ratio in the throttling valve related to the total exergy destruction

> 

y[TV, `EN,TV`, 1] := Units:-Standard:-`*`(E[D, TV, EN, TV], Units:-Standard:-`/`(E[D, TOT, `EN,TV`]))

.2491288143

 (4.569)

Exergy destruction ratio in the evaporator related to the total exergy destruction

> 

y[EV, `EN,TV`, 1] := Units:-Standard:-`*`(E[D, EV, EN, TV], Units:-Standard:-`/`(E[D, TOT, `EN,TV`]))

.2884626223

 (4.570)

Exergy destruction ratio in the system overall related to the total exergy destruction

> 

y[TOT, `EN,TV`, 1] := Units:-Standard:-`*`(E[D, TOT, `EN,TV`], Units:-Standard:-`/`(E[D, TOT, `EN,TV`]))

 

1.000000000

 (4.571)

 

7. Hybrid cycle with irreversible component - Evaporator

Temperature of the working fluid at the compressor outlet

> 

T[4, EN, EV] := Units:-Standard:-`+`(T[EV, h, out], Units:-Standard:-`-`(`&delta;T`[EV, c]))

248.15*Units:-Unit(K)

 (4.572)

Pressure of the working fluid at the compressor outlet

> 

p[2, EN, EV] := PSTXwf(Units:-Standard:-`+`(T[CD, c, out], 0), 1, wf[cycle])

 

1167197.054*Units:-Unit(Pa)

 (4.573)

 

Pressure of  the working fluid at the condenser outlet

> 

p[3, EN, EV] := p[2, EN, EV]

 

1167197.054*Units:-Unit(Pa)

 (4.574)

 

Pressure of  the working fluid at the condenser outlet

> 

T[3, `EN,EV`] := TSPXwf(p[3, EN, EV], 0, wf[cycle])

 

303.1500000*Units:-Unit(K)

 (4.575)

 

 

Specific enthalpy of  the working fluid at the condenser outlet

> 

h[3, `EN,EV`] := HSTXwf(T[3, `EN,EV`], 0, wf[cycle])

 

484911.4924*Units:-Unit(J/kg)

 (4.576)

 

Specific entropy of  the working fluid at the condenser outlet

> 

s[3, `EN,EV`] := SSTXwf(T[3, `EN,EV`], 0, wf[cycle])

 

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

 (4.577)

 

Specific volume of  the working fluid at the condenser outlet

> 

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

 

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

 (4.578)

 

Specific exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,EV`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, `EN,EV`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, `EN,EV`], Units:-Standard:-`-`(s[0])))))

 

296150.505*Units:-Unit(J/kg)

 (4.579)

 

Auxiliary values

> 

h[3, `EN,EV`, 1] := HPTwf(p[3, EN, EV], T[0], wf[cycle])

 

437080.3036*Units:-Unit(J/kg)

 (4.580)

 

> 

s[3, `EN,EV`, 1] := SPTwf(p[3, EN, EV], T[0], wf[cycle])

 

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

 (4.581)

 

Specific mechanical exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,EV`, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[3, `EN,EV`, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[3, `EN,EV`, 1]))))

 

295350.397*Units:-Unit(J/kg)

 (4.582)

 

Specific thermal exergy of  the working fluid at the condenser outlet

> 

e[3, `EN,EV`, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[3, `EN,EV`], Units:-Standard:-`-`(h[3, `EN,EV`, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[3, `EN,EV`], Units:-Standard:-`-`(s[3, `EN,EV`, 1])))))

 

800.10886*Units:-Unit(J/kg)

 (4.583)

 

Specific entropy of the working fluid at the evaporator intlet

> 

s[4, EN, EV] := s[3, `EN,EV`]

 

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

 (4.584)

 

Pressure of the working fluid at the evaporator intlet

 

> 

p[4, EN, EV] := Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`*`(PSTXwf(T[4, EN, EV], 0, wf[cycle]), Units:-Standard:-`+`(SSTXwf(T[4, EN, EV], 1, wf[cycle]), Units:-Standard:-`-`(s[4, EN, EV]))), Units:-Standard:-`*`(PSTXwf(T[4, EN, EV], 1, wf[cycle]), Units:-Standard:-`+`(s[4, EN, EV], Units:-Standard:-`-`(SSTXwf(T[4, EN, EV], 0, wf[cycle]))))), Units:-Standard:-`/`(Units:-Standard:-`+`(SSTXwf(T[4, EN, EV], 1, wf[cycle]), Units:-Standard:-`-`(SSTXwf(T[4, EN, EV], 0, wf[cycle])))))

 

151471.2936*Units:-Unit(Pa)

 (4.585)

 

Specific enthalpy of saturated liquid of  the working fluid in the evaporator

> 

h[sl, 4, `EN,EV`] := HSPXwf(p[4, EN, EV], 0, wf[cycle])

229167.5145*Units:-Unit(J/kg)

 (4.586)

Specific enthalpy of saturated vapor of  the working fluid in the evaporator

> 

h[sv, 4, `EN,EV`] := HSPXwf(p[4, EN, EV], 1, wf[cycle])

1573806.240*Units:-Unit(J/kg)

 (4.587)

Specific entropy of saturated liquid of  the working fluid in the evaporator

> 

s[sl, 4, EN, EV] := SSPXwf(p[4, EN, EV], 0, wf[cycle])

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

 (4.588)

Specific entropy of saturated vapor of  the working fluid in the evaporator

> 

s[sv, 4, EN, EV] := SSPXwf(p[4, EN, EV], 1, wf[cycle])

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

 (4.589)

Quality of the working fluid at the evaporator intlet

> 

x[4, `EN,EV`] := Units:-Standard:-`*`(Units:-Standard:-`+`(s[4, EN, EV], Units:-Standard:-`-`(s[sl, 4, EN, EV])), Units:-Standard:-`/`(Units:-Standard:-`+`(s[sv, 4, EN, EV], Units:-Standard:-`-`(s[sl, 4, EN, EV]))))NULL

 

.1705216311

 (4.590)

 

 

Specific enthalpy of  the working fluid at the evaporator inlet

> 

h[4, `EN,EV`] := Units:-Standard:-`+`(h[sl, 4, `EN,EV`], Units:-Standard:-`*`(x[4, `EN,EV`], Units:-Standard:-`+`(h[sv, 4, `EN,EV`], Units:-Standard:-`-`(h[sl, 4, `EN,EV`]))))

458457.5033*Units:-Unit(J/kg)

 (4.591)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

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

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

 (4.592)

Specific volume of saturated liquid of  the working fluid in the evaporator

> 

v[sv, 4, `EN,EV`] := Units:-Standard:-`/`(DSPXwf(p[4, EN, EV], 1, wf[cycle]))

.7716715235*Units:-Unit(m^3/kg)

 (4.593)

Specific volume of  the working fluid at the evaporator inlet

> 

v[4, `EN,EV`] := Units:-Standard:-`+`(v[sl, 4, `EN,EV`], Units:-Standard:-`*`(x[4, `EN,EV`], Units:-Standard:-`+`(v[sv, 4, `EN,EV`], Units:-Standard:-`-`(v[sl, 4, `EN,EV`]))))

.1328218926*Units:-Unit(m^3/kg)

 (4.594)

Specific exergy of  the working fluid at the evaporator inlet

> 

e[4, `EN,EV`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, `EN,EV`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, EN, EV], Units:-Standard:-`-`(s[0])))))

 

269696.516*Units:-Unit(J/kg)

 (4.595)

 

Auxiliary values

> 

h[4, `EN,EV`, 1] := HPTwf(p[4, EN, EV], T[0], wf[cycle])

 

1675673.399*Units:-Unit(J/kg)

 (4.596)

 

> 

s[4, `EN,EV`, 1] := SPTwf(p[4, EN, EV], T[0], wf[cycle])

 

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

 (4.597)

 

Specific mechanical exergy of  the working fluid at the evaporator inlet

> 

e[4, `EN,EV`, M] := Units:-Standard:-`+`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[0], Units:-Standard:-`-`(s[4, `EN,EV`, 1]))), Units:-Standard:-`-`(Units:-Standard:-`+`(h[0], Units:-Standard:-`-`(h[4, `EN,EV`, 1]))))

 

58604.35859*Units:-Unit(J/kg)

 (4.598)

 

Specific thermal exergy of  the working fluid at the evaporator inlet

> 

e[4, `EN,EV`, T] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[4, `EN,EV`], Units:-Standard:-`-`(h[4, `EN,EV`, 1])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[4, EN, CD], Units:-Standard:-`-`(s[4, `EN,EV`, 1])))))

 

165219.018*Units:-Unit(J/kg)

 (4.599)

 

Pressure of  the working fluid at the compressor inlet

> 

p[1, EN, EV] := p[4, EN, EV]

151471.2936*Units:-Unit(Pa)

 (4.600)

Temperature of  the working fluid at the compressor inlet

> 

T[1, EN, EV] := TSPXwf(p[1, EN, EV], 1, wf[cycle])

248.1500000*Units:-Unit(K)

 (4.601)

Specific enthalpy of  the working fluid at the compressor inlet

> 

h[1, `EN,EV`] := HSPXwf(p[1, EN, EV], 1, wf[cycle])

1573806.240*Units:-Unit(J/kg)

 (4.602)

Specific entropy of  the working fluid at the compressor inlet

> 

s[1, `EN,EV`] := SSPXwf(p[1, EN, EV], 1, wf[cycle])

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

 (4.603)

Specific volume of  the working fluid at the compressor inlet

> 

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

 

.7716715235*Units:-Unit(m^3/kg)

 (4.604)

 

Specific exergy of  the working fluid at the compressor inlet

> 

e[1, `EN,EV`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[1, `EN,EV`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[1, `EN,EV`], Units:-Standard:-`-`(s[0])))))

 

67437.0235*Units:-Unit(J/kg)

 (4.605)

 

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

> 

h[2, `EN,EV`] := HcompPTPEFFwf(p[1, EN, EV], T[1, EN, EV], p[2, EN, EV], 1, wf[cycle])

1877968.466*Units:-Unit(J/kg)

 (4.606)

Enhtalpy change in the compressor 1 after actual compression

> 

l[comp, `EN,EV`] := Units:-Standard:-`+`(h[2, `EN,EV`], Units:-Standard:-`-`(h[1, `EN,EV`]))

 

304162.226*Units:-Unit(J/kg)

 (4.607)

 

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

> 

T[2, EN, EV] := TPHwf(p[2, EN, EV], h[2, `EN,EV`], wf[cycle])

 

395.9009657*Units:-Unit(K)

 (4.608)

 

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

> 

s[2, `EN,EV`] := SPTwf(p[2, EN, EV], T[2, EN, EV], wf[cycle])

 

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

 (4.609)

 

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

> 

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

.1585436940*Units:-Unit(m^3/kg)

 (4.610)

Specific exergy of  the working fluid at the compressor inlet

> 

e[2, `EN,EV`] := Units:-Standard:-`+`(Units:-Standard:-`+`(h[2, `EN,EV`], Units:-Standard:-`-`(h[0])), Units:-Standard:-`-`(Units:-Standard:-`*`(T[0], Units:-Standard:-`+`(s[2, `EN,EV`], Units:-Standard:-`-`(s[0])))))

 

371599.2495*Units:-Unit(J/kg)

 (4.611)

 

Specific heat rejection in the condenser

> 

q[out, `EN,EV`] := Units:-Standard:-`+`(h[2, `EN,EV`], Units:-Standard:-`-`(h[3, `EN,EV`]))

1393056.974*Units:-Unit(J/kg)

 (4.612)

Specific heat addition in the evaporator

> 

q[inlet, `EN,EV`] := Units:-Standard:-`+`(h[1, `EN,EV`], Units:-Standard:-`-`(h[4, `EN,EV`]))

1115348.737*Units:-Unit(J/kg)

 (4.613)

Coefficent of performance of a refrigerator

> 

COP[`EN,EV`, Ref] := Units:-Standard:-`*`(q[inlet, `EN,EV`], Units:-Standard:-`/`(l[comp, `EN,EV`]))

3.666953493

 (4.614)

Coefficent of performance of a heat pump

> 

COP[`EN,EV`, HP] := Units:-Standard:-`*`(q[out, `EN,EV`], Units:-Standard:-`/`(l[comp, `EN,EV`]))

4.579980205

 (4.615)

Mass rate of the working fluid

> 

m[wf, cycle, `EN,EV`] := proc (Type) if Type = 1 then Q[P]/q[out, `EN,EV`] else Q[P]/q[inlet, `EN,EV`] end if end proc:

 

> 

m[wf, cycle, `EN,EV`] := m[wf, cycle, `EN,EV`](Type)

 

0.8965805643e-4*kW*Units:-Unit(kg/J)

 (4.616)

 

Total heat addition in the evaporator

> 

Q[inlet, `EN,EV`] := Units:-Standard:-`*`(m[wf, cycle, `EN,EV`], q[inlet, `EN,EV`])

 

100.0000000*kW

 (4.617)

 

Total heat rejection in the condenser

> 

Q[out, `EN,EV`] := Units:-Standard:-`*`(m[wf, cycle, R], q[out, `EN,EV`])

 

133.9458515*kW

 (4.618)

 

Mass rate of the heat source fluid

> 

m[EV, h, `EN,EV`] := Units:-Standard:-`*`(Q[inlet, `EN,EV`], Units:-Standard:-`/`(Units:-Standard:-`+`(h[EV, h, inlet], Units:-Standard:-`-`(h[EV, h, out]))))

 

0.9944845786e-2*kW*Units:-Unit(kg/J)

 (4.619)

 

Mass rate of the heat consumption fluid

> 

m[CD, c, `EN,EV`] := Units:-Standard:-`*`(Q[out, `EN,EV`], Units:-Standard:-`/`(Units:-Standard:-`+`(h[CD, c, out], Units:-Standard:-`-`(h[CD, c, inlet]))))

 

0.1331085667e-1*kW*Units:-Unit(kg/J)

 (4.620)

 

Total exergy used in the compressor

> 

E[5, EN, EV] := Units:-Standard:-`*`(m[wf, cycle, `EN,EV`], l[comp, `EN,EV`])

 

27.27059402*kW

 (4.621)

 

Total physical exergy at the compressor intlet

> 

E[1, `EN,EV`] := Units:-Standard:-`*`(m[wf, cycle, `EN,EV`], e[1, `EN,EV`])

 

6.046272458*kW

 (4.622)

 

Total physical exergy at the compressor outlet

> 

E[2, `EN,EV`] := Units:-Standard:-`*`(m[wf, cycle, `EN,EV`], e[2, `EN,EV`])

 

33.31686648*kW

 (4.623)

 

Total physical exergy at the condenser outlet

> 

E[3, `EN,EV`] := Units:-Standard:-`*`(m[wf, cycle, `EN,EV`], e[3, `EN,EV`])

 

26.55227869*kW

 (4.624)

 

Total mechanical exergy at the condenser outlet

> 

E[3, `EN,EV`, M] := Units:-Standard:-`*`(m[wf, cycle, `EN,EV`], e[3, `EN,EV`, M])

 

26.48054256*kW

 (4.625)

 

Total thermal exergy at the condenser outlet

> 

E[3, `EN,EV`, T] := Units:-Standard:-`*`(m[wf, cycle, `EN,EV`], e[3, `EN,EV`, T])

 

0.7173620532e-1*kW

 (4.626)

 

Total physical exergy at the evaporator inlet

> 

E[4, EN, EV] := Units:-Standard:-`*`(m[wf, cycle, `EN,EV`], e[4, `EN,EV`])

 

24.18046545*kW

 (4.627)

 

Total mechanical exergy at the evaporator inlet

> 

E[4, EN, EV, M] := Units:-Standard:-`*`(m[wf, cycle, `EN,EV`], e[4, `EN,EV`, M])

 

5.254352890*kW

 (4.628)

 

Total thermal exergy at the evaporator inlet

> 

E[4, `EN,EV`, T] := Units:-Standard:-`*`(m[wf, cycle, `EN,EV`], e[4, `EN,EV`, T])

 

14.81321604*kW

 (4.629)

 

Total physical exergy of heat consumption fluid at the evaporator intlet 

> 

E[EV, h, inlet] := Units:-Standard:-`*`(m[EV, h, `EN,EV`], e[EV, h, inlet])

 

11.30975886*kW

 (4.630)

 

Total physical exergy of heat consumption fluid at the evaporator outlet

> 

E[EV, h, out] := Units:-Standard:-`*`(m[EV, h, `EN,EV`], e[EV, h, out])

 

22.72348304*kW

 (4.631)

 

Total physical exergy of heat consumption fluid at the condenser intlet 

> 

E[CD, c, inlet] := Units:-Standard:-`*`(m[CD, c, `EN,EV`], e[CD, c, inlet])

 

-0.

 (4.632)

 

Total physical exergy of heat consumption fluid at the condenser outlet

> 

E[CD, c, out] := Units:-Standard:-`*`(m[CD, c, `EN,EV`], e[CD, c, out])

 

2.234061745*kW

 (4.633)

 

Exergy of fuel in the compressor

> 

E[F, CM, `EN,EV`] := E[5, EN, EV]

 

27.27059402*kW

 (4.634)

 

Exergy of product in the compressor

> 

E[P, CM, `EN,EV`] := Units:-Standard:-`+`(E[2, `EN,EV`], Units:-Standard:-`-`(E[1, `EN,EV`]))

 

27.27059402*kW

 (4.635)

 

Exergy destruction in the compressor

> 

E[D, CM, `EN,EV`] := Units:-Standard:-`+`(E[F, CM, `EN,EV`], Units:-Standard:-`-`(E[P, CM, `EN,EV`]))

 

0.

 (4.636)

 

Exergetic efficiency in the compressor

> 

varepsilon[CM, `EN,EV`] := Units:-Standard:-`*`(E[P, CM, `EN,EV`], Units:-Standard:-`/`(E[F, CM, `EN,EV`]))

 

1.000000000

 (4.637)

 

Exergy of fuel in the condenser

> 

E[F, CD, `EN,EV`] := Units:-Standard:-`+`(E[2, `EN,EV`], Units:-Standard:-`-`(E[3, `EN,EV`]))

 

6.76458779*kW

 (4.638)

 

Exergy of product in the condenser

> 

E[P, CD, `EN,EV`] := Units:-Standard:-`+`(E[CD, c, out], Units:-Standard:-`-`(E[CD, c, inlet]))

 

2.234061745*kW

 (4.639)

 

Exergy destruction in the condenser

> 

E[D, CD, `EN,EV`] := Units:-Standard:-`+`(E[F, CD, `EN,EV`], Units:-Standard:-`-`(E[P, CD, `EN,EV`]))

 

4.530526045*kW

 (4.640)

 

Exergy of fuel in the throttling valve

> 

E[F, TV, EN, EV] := Units:-Standard:-`+`(E[3, `EN,EV`, M], Units:-Standard:-`-`(E[4, EN, EV, M]))

 

21.22618967*kW

 (4.641)

 

Exergy of product in the throttling valve

> 

E[P, TV, EN, EV] := Units:-Standard:-`+`(E[3, `EN,EV`, M], Units:-Standard:-`-`(E[4, EN, EV, M]))

 

21.22618967*kW

 (4.642)

 

Exergy destruction in the throttling valve

> 

E[D, TV, EN, EV] := Units:-Standard:-`+`(E[F, TV, EN, EV], Units:-Standard:-`-`(E[P, TV, EN, EV]))

 

0.

 (4.643)

 

Exergetic efficiency in the throttling valve       

> 

varepsilon[TV, `EN,EV`] := Units:-Standard:-`*`(E[P, TV, EN, EV], Units:-Standard:-`/`(E[F, TV, EN, EV]))

 

1.000000000

 (4.644)

 

Exergy of fuel in the evaporator

> 

E[F, EV, `EN,EV`] := Units:-Standard:-`+`(E[4, EN, EV], Units:-Standard:-`-`(E[1, `EN,EV`]))

 

18.13419299*kW

 (4.645)

 

Exergy of product in the evaporator

> 

E[P, EV, `EN,EV`] := Units:-Standard:-`+`(E[EV, h, out], Units:-Standard:-`-`(E[EV, h, inlet]))

 

11.41372418*kW

 (4.646)

 

Exergy destruction in the evaporator

> 

E[D, EV, EN, EV] := Units:-Standard:-`+`(E[F, EV, `EN,EV`], Units:-Standard:-`-`(E[P, EV, `EN,EV`]))

 

6.72046881*kW

 (4.647)

 

Exergetic efficiency in the evaporator

> 

varepsilon[EV, `EN,EV`] := Units:-Standard:-`*`(E[P, EV, `EN,EV`], Units:-Standard:-`/`(E[F, EV, `EN,EV`]))

 

.6294034803

 (4.648)

 

Exergy of fuel of the overall system

> 

E[F, TOT, `EN,EV`] := E[5, EN, EV]

27.27059402*kW

 (4.649)

Exergy of product of the overall system 

> 

E[P, TOT, `EN,EV`] := proc (Type) if Type = 1 then E[P, CD, `EN,EV`] else E[P, EV, `EN,EV`] end if end proc:

 

> 

E[P, TOT, `EN,EV`] := E[P, TOT, `EN,EV`](Type)

 

11.41372418*kW

 (4.650)

 

Exergy destruction of the overall system

> 

E[D, TOT, `EN,EV`] := Units:-Standard:-`+`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CM, `EN,EV`], E[D, CD, `EN,EV`]), E[D, TV, EN, EV]), E[D, EV, EN, EV])

11.25099486*kW

 (4.651)

Exergetic efficiency of the overall system

> 

varepsilon[TOT, `EN,EV`] := Units:-Standard:-`*`(E[P, TOT, `EN,EV`], Units:-Standard:-`/`(E[F, TOT, `EN,EV`]))

.4185359575

 (4.652)

Exergy destruction ratio in the compressor related to the exergy of total fuel

> 

y[CM, `EN,EV`] := Units:-Standard:-`*`(E[D, CM, `EN,EV`], Units:-Standard:-`/`(E[F, TOT, `EN,EV`]))

 

0.

 (4.653)

 

Exergy destruction ratio in the condenser related to the exergy of total fuel

> 

y[CD, `EN,EV`] := Units:-Standard:-`*`(E[D, CD, `EN,EV`], Units:-Standard:-`/`(E[F, TOT, `EN,EV`]))

.1661322831

 (4.654)

Exergy destruction ratio in the throttling valve related to the exergy of total fuel

> 

y[TV, `EN,EV`] := Units:-Standard:-`*`(E[D, TV, EN, EV], Units:-Standard:-`/`(E[F, TOT, `EN,EV`]))

0.

 (4.655)

Exergy destruction ratio in the evaporator related to the exergy of total fuel

> 

y[EV, `EN,EV`] := Units:-Standard:-`*`(E[D, EV, EN, EV], Units:-Standard:-`/`(E[F, TOT, `EN,EV`]))

.2464364658

 (4.656)

Exergy destruction ratio in the system overall related to the exergy of total fuel

> 

y[TOT, `EN,EV`] := Units:-Standard:-`*`(E[D, TOT, `EN,EV`], Units:-Standard:-`/`(E[F, TOT, `EN,EV`]))

 

.4125687490

 (4.657)

 

Exergy destruction ratio in the compressor related to the total exergy destruction

> 

y[CM, `EN,EV`, 1] := Units:-Standard:-`*`(E[D, CM, `EN,EV`], Units:-Standard:-`/`(E[D, TOT, `EN,EV`]))

 

0.

 (4.658)

 

Exergy destruction ratio in the condenser related to the total exergy destruction

> 

y[CD, `EN,EV`, 1] := Units:-Standard:-`*`(E[D, CD, `EN,EV`], Units:-Standard:-`/`(E[D, TOT, `EN,EV`]))

.4026778166

 (4.659)

Exergy destruction ratio in the throttling valve related to the total exergy destruction

> 

y[TV, `EN,EV`, 1] := Units:-Standard:-`*`(E[D, TV, EN, EV], Units:-Standard:-`/`(E[D, TOT, `EN,EV`]))

0.

 (4.660)

Exergy destruction ratio in the evaporator related to the total exergy destruction

> 

y[EV, `EN,EV`, 1] := Units:-Standard:-`*`(E[D, EV, EN, EV], Units:-Standard:-`/`(E[D, TOT, `EN,EV`]))

.5973221829

 (4.661)

Exergy destruction ratio in the system overall related to the total exergy destruction

> 

y[TOT, `EN,EV`, 1] := Units:-Standard:-`*`(E[D, TOT, `EN,EV`], Units:-Standard:-`/`(E[D, TOT, `EN,EV`]))

 

1.000000000

 (4.662)

 

Endogenous exergy destruction of the overall system

> 

E[D, TOT, EN] := Units:-Standard:-`+`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CM, EN, CM], E[D, CD, EN, CD]), E[D, TV, EN, TV]), E[D, EV, EN, EV])

 

20.50092868*kW

 (4.663)

 

 

Defining function for computing endogenous unavoidable exergy destructions of the system componens using the proc commandNULL
 
> 

E[D, UN, EN] := proc (Compomemt) local `&eta;_comp`, `&delta;T_EV_c`, `&delta;T_CD_h`, T4, p2, p3, T3, h3, s3, e3, h3_0, s3_0, e3_M, e3_T, h4, p4, hsl_4, hsv_4, x4, ssl_4, ssv_4, s4, e4, h4_0, s4_0, e4_M, e4_T, p1, T1, h1, s1, e1, h2, l_comp, T2, s2, e2, q_out, q_inlet, m_w_f, Q_inlet, Q_out, m_EV, m_CD, E_EV_inlet, E_EV_out, E5, E1, E2, E3, E3_M, E3_T, E4, E4_M, E4_T, E_F, E_P, E_D; if Compomemt = CM then `&eta;_comp` := eta[comp, UN]; `&delta;T_EV_c` := 0; `&delta;T_CD_h` := 0 else if Compomemt = CD then `&eta;_comp` := 1; `&delta;T_EV_c` := 0; `&delta;T_CD_h` := `&delta;T`[CD, h, UN] else if Compomemt = TV then `&eta;_comp` := 1; `&delta;T_EV_c` := 0; `&delta;T_CD_h` := 0 else if Compomemt = EV then `&eta;_comp` := 1; `&delta;T_EV_c` := `&delta;T`[EV, c, UN]; `&delta;T_CD_h` := 0 end if end if end if end if; T4 := T[EV, h, out]-`&delta;T_EV_c`; p2 := PSTXwf(T[CD, c, out]+`&delta;T_CD_h`, 1, wf[cycle]); p3 := p2; T3 := TSPXwf(p3, 0, wf[cycle]); h3 := HSTXwf(T3, 0, wf[cycle]); s3 := SSTXwf(T3, 0, wf[cycle]); e3 := h3-h[0]-T[0]*(s3-s[0]); h3_0 := HPTwf(p3, T[0], wf[cycle]); s3_0 := SPTwf(p3, T[0], wf[cycle]); e3_M := T[0]*(s[0]-s3_0)-h[0]+h3_0; e3_T := h3-h3_0-T[0]*(s3-s3_0); if Compomemt = TV then h4 := h3; p4 := (PSTXwf(T4, 0, wf[cycle])*(HSTXwf(T4, 1, wf[cycle])-h4)+PSTXwf(T4, 1, wf[cycle])*(h4-HSTXwf(T4, 0, wf[cycle])))/(HSTXwf(T4, 1, wf[cycle])-HSTXwf(T4, 0, wf[cycle])); hsl_4 := HSPXwf(p4, 0, wf[cycle]); hsv_4 := HSPXwf(p4, 1, wf[cycle]); x4 := (h4-hsl_4)/(hsv_4-hsl_4); ssl_4 := SSPXwf(p4, 0, wf[cycle]); ssv_4 := SSPXwf(p4, 1, wf[cycle]); s4 := ssl_4+x4*(ssv_4-ssl_4); e4 := h4-h[0]-T[0]*(s4-s[0]); h4_0 := HPTwf(p4, T[0], wf[cycle]); s4_0 := SPTwf(p4, T[0], wf[cycle]); e4_M := T[0]*(s[0]-s4_0)-h[0]+h4_0; e4_T := h4-h4_0-T[0]*(s4-s4_0) else s4 := s3; p4 := (PSTXwf(T4, 0, wf[cycle])*(SSTXwf(T4, 1, wf[cycle])-s4)+PSTXwf(T4, 1, wf[cycle])*(s4-SSTXwf(T4, 0, wf[cycle])))/(SSTXwf(T4, 1, wf[cycle])-SSTXwf(T4, 0, wf[cycle])); ssl_4 := SSPXwf(p4, 0, wf[cycle]); ssv_4 := SSPXwf(p4, 1, wf[cycle]); x4 := (s4-ssl_4)/(ssv_4-ssl_4); hsl_4 := HSPXwf(p4, 0, wf[cycle]); hsv_4 := HSPXwf(p4, 1, wf[cycle]); h4 := hsl_4+x4*(hsv_4-hsl_4); e4 := h4-h[0]-T[0]*(s4-s[0]); h4_0 := HPTwf(p4, T[0], wf[cycle]); s4_0 := SPTwf(p4, T[0], wf[cycle]); e4_M := T[0]*(s[0]-s4_0)-h[0]+h4_0; e4_T := h4-h4_0-T[0]*(s4-s4_0) end if; p1 := p4; T1 := TSPXwf(p1, 1, wf[cycle]); h1 := HSPXwf(p1, 1, wf[cycle]); s1 := SSPXwf(p1, 1, wf[cycle]); e1 := h1-h[0]-T[0]*(s1-s[0]); h2 := HcompPTPEFFwf(p1, T1, p2, `&eta;_comp`, wf[cycle]); l_comp := h2-h1; T2 := TPHwf(p2, h2, wf[cycle]); s2 := SPTwf(p2, T2, wf[cycle]); e2 := h2-h[0]-T[0]*(s2-s[0]); q_out := h2-h3; q_inlet := h1-h4; if Type = 1 then m_w_f := Q[P]/q_out else m_w_f := Q[P]/q_inlet end if; Q_inlet := m_w_f*q_inlet; Q_out := m_w_f*q_out; m_EV := Q_inlet/(h[EV, h, inlet]-h[EV, h, out]); m_CD := Q_out/(h[CD, c, out]-h[CD, c, inlet]); E_EV_inlet := m_EV*e[EV, h, inlet]; E_EV_out := m_EV*e[EV, h, out]; E5 := l_comp*m_w_f; E1 := e1*m_w_f; E2 := e2*m_w_f; E3 := e3*m_w_f; E3_M := e3_M*m_w_f; E3_T := e3_T*m_w_f; E4 := e4*m_w_f; E4_M := e4_M*m_w_f; E4_T := e4_T*m_w_f; if Compomemt = CM then E_F := E5; E_P := E2-E1; E_D := E_F-E_P else if Compomemt = CD then E_F := E2-E3; E_P := E[CD, c, out]-E[CD, c, inlet]; E_D := E_F-E_P else if Compomemt = TV then E_F := E3_M-E4_M; E_P := E4_T-E3_T; E_D := E_F-E_P else if Compomemt = EV then E_F := E4-E1; E_P := E_EV_out-E_EV_inlet; E_D := E_F-E_P end if end if end if end if; return E_D end proc:

Endogenous unaviodable exergy destruction of the overall system

> 

E[D, TOT, UN, EN] := Units:-Standard:-`+`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, UN, EN](CM), E[D, UN, EN](CD)), E[D, UN, EN](TV)), E[D, UN, EN](EV))

 

8.548693595*kW

 (4.1.1)

 

 

Table 1: Conventional exergetic analysis

Component

Exergy of fuel

Exergy of product

Exergy destruction

Exergetic efficiency

Exergy destruction related to the total exergy fuel

Exergy destruction related to the total exergy destruction

Compressor,kW
> 

E[F, CM, R]

43.14038984*kW

 (4.664)
> 

E[P, CM, R]

37.37911102*kW

 (4.665)
> 

E[D, CM, R]

5.76127882*kW

 (4.666)

-

-

-

%
> 

Units:-Standard:-`*`(E[F, CM, R], Units:-Standard:-`/`(E[F, TOT, R]))

1.000000000

 (4.667)
> 

Units:-Standard:-`*`(E[P, CM, R], Units:-Standard:-`/`(E[P, TOT, R]))

3.274926784

 (4.668)
> 

Units:-Standard:-`*`(E[D, CM, R], Units:-Standard:-`/`(E[D, TOT, R]))

.1963676264

 (4.669)
> 

varepsilon[CM, R]

.8664527873

 (4.670)
> 

y[CM, R]

.1335472128

 (4.671)
> 

y[CM, R, 1]

.1963676264

 (4.672)

Condenser, kW
> 

E[F, CD, R]

15.09865680*kW

 (4.673)
> 

E[P, CD, R]

2.387415999*kW

 (4.674)
> 

E[D, CD, R]

12.71124080*kW

 (4.675)

-

-

-

%
> 

Units:-Standard:-`*`(E[F, CD, R], Units:-Standard:-`/`(E[F, TOT, R]))

.3499888818

 (4.676)
> 

Units:-Standard:-`*`(E[P, CD, R], Units:-Standard:-`/`(E[P, TOT, R]))

.2091706407

 (4.677)
> 

Units:-Standard:-`*`(E[D, CD, R], Units:-Standard:-`/`(E[D, TOT, R]))

.4332503706

 (4.678)
> 

varepsilon[CD, R]

.1581210853

 (4.679)
> 

y[CD, R]

.2946482600

 (4.680)
> 

y[CD, R, 1]

.4332503706

 (4.681)

Throttling valve, kW
> 

E[F, TV, R]

22.82484318*kW

 (4.682)
> 

E[P, TV, R]

18.67858199*kW

 (4.683)
> 

E[D, TV, R]

4.14626119*kW

 (4.684)

-

-

-

%
> 

Units:-Standard:-`*`(E[F, TV, R], Units:-Standard:-`/`(E[F, TOT, R]))

.5290829143

 (4.685)
> 

Units:-Standard:-`*`(E[P, TV, R], Units:-Standard:-`/`(E[P, TOT, R]))

1.636501960

 (4.686)
> 

Units:-Standard:-`*`(E[D, TV, R], Units:-Standard:-`/`(E[D, TOT, R]))

.1413213097

 (4.687)
> 

varepsilon[TV, R]

.8183443734

 (4.688)
> 

y[TV, R]

0.9611088832e-1

 (4.689)
> 

y[TV, R, 1]

.1413213097

 (4.690)

Evaporator, kW
> 

E[F, EV, R]

18.13419306*kW

 (4.691)
> 

E[P, EV, R]

11.41372418*kW

 (4.692)
> 

E[D, EV, R]

6.72046888*kW

 (4.693)

-

-

-

%
> 

Units:-Standard:-`*`(E[F, EV, R], Units:-Standard:-`/`(E[F, TOT, R]))

.4203530179

 (4.694)
> 

Units:-Standard:-`*`(E[P, EV, R], Units:-Standard:-`/`(E[P, TOT, R]))

1.000000000

 (4.695)
> 

Units:-Standard:-`*`(E[D, EV, R], Units:-Standard:-`/`(E[D, TOT, R]))

.2290606935

 (4.696)
> 

varepsilon[EV, R]

.6294034779

 (4.697)
> 

y[EV, R]

.1557813665

 (4.698)
> 

y[EV, R, 1]

.2290606935

 (4.699)

Overall system, kW
> 

E[F, TOT, R]

43.14038984*kW

 (4.700)
> 

E[P, TOT, R]

11.41372418*kW

 (4.701)
> 

E[D, TOT, R]

29.33924969*kW

 (4.702)

-

-

-

%
> 

Units:-Standard:-`*`(E[F, TOT, R], Units:-Standard:-`/`(E[F, TOT, R]))

1.000000000

 (4.703)
> 

Units:-Standard:-`*`(E[P, TOT, R], Units:-Standard:-`/`(E[P, TOT, R]))

1.000000000

 (4.704)
> 

Units:-Standard:-`*`(E[D, TOT, R], Units:-Standard:-`/`(E[D, TOT, R]))

1.000000000

 (4.705)
> 

varepsilon[TOT, R]

.2645716514

 (4.706)
> 

y[TOT, R]

.6800877275

 (4.707)
> 

y[TOT, R, 1]

1.000000000

 (4.708)

``

Table 2: Splitting the exergy destruction for components in advanced exergetic analysis

Component

Total in Real Cycle

Theoretical

Unavoidable

Avoidable

Endogenous

Exogenous

Compressor,kW
> 

E[D, CM, R]

5.76127882*kW

 (4.709)
> 

E[D, CM, T]

-0.5e-7*kW

 (4.710)
> 

E[D, CM, UN]

.88338470*kW

 (4.711)
> 

Units:-Standard:-`+`(E[D, CM, R], Units:-Standard:-`-`(E[D, CM, UN]))

4.87789412*kW

 (4.712)
> 

E[D, CM, EN, CM]

3.94436090*kW

 (4.713)
> 

Units:-Standard:-`+`(E[D, CM, R], Units:-Standard:-`-`(E[D, CM, EN, CM]))

1.81691792*kW

 (4.714)

%
> 

Units:-Standard:-`*`(E[D, CM, R], Units:-Standard:-`/`(E[D, CM, R]))

1.000000000

 (4.715)
> 

Units:-Standard:-`*`(E[D, CM, T], Units:-Standard:-`/`(E[D, CM, R]))

-0.8678628750e-8

 (4.716)
> 

Units:-Standard:-`*`(E[D, CM, UN], Units:-Standard:-`/`(E[D, CM, R]))

.1533313571

 (4.717)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, CM, R], Units:-Standard:-`-`(E[D, CM, UN])), Units:-Standard:-`/`(E[D, CM, R]))

.8466686430

 (4.718)
> 

Units:-Standard:-`*`(E[D, CM, EN, CM], Units:-Standard:-`/`(E[D, CM, R]))

.6846328781

 (4.719)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, CM, R], Units:-Standard:-`-`(E[D, CM, EN, CM])), Units:-Standard:-`/`(E[D, CM, R]))

.3153671219

 (4.720)

Condenser, kW
> 

E[D, CD, R]

12.71124080*kW

 (4.721)
> 

E[D, CD, T]

3.269590938*kW

 (4.722)
> 

E[D, CD, UN]

4.074822525*kW

 (4.723)
> 

Units:-Standard:-`+`(E[D, CD, R], Units:-Standard:-`-`(E[D, CD, UN]))

8.636418275*kW

 (4.724)
> 

E[D, CD, EN, CD]

7.984201895*kW

 (4.725)
> 

Units:-Standard:-`+`(E[D, CD, R], Units:-Standard:-`-`(E[D, CD, EN, CD]))

4.727038905*kW

 (4.726)

%
> 

Units:-Standard:-`*`(E[D, CD, R], Units:-Standard:-`/`(E[D, CD, R]))

1.000000000

 (4.727)
> 

Units:-Standard:-`*`(E[D, CD, T], Units:-Standard:-`/`(E[D, CD, R]))

.2572204389

 (4.728)
> 

Units:-Standard:-`*`(E[D, CD, UN], Units:-Standard:-`/`(E[D, CD, R]))

.3205684314

 (4.729)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, CD, R], Units:-Standard:-`-`(E[D, CD, UN])), Units:-Standard:-`/`(E[D, CD, R]))

.6794315685

 (4.730)
> 

Units:-Standard:-`*`(E[D, CD, EN, CD], Units:-Standard:-`/`(E[D, CD, R]))

.6281213628

 (4.731)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, CD, R], Units:-Standard:-`-`(E[D, CD, EN, CD])), Units:-Standard:-`/`(E[D, CD, R]))

.3718786371

 (4.732)

Throttling valve, kW
> 

E[D, TV, R]

4.14626119*kW

 (4.733)
> 

E[D, TV, T]

0.

 (4.734)
> 

E[D, TV, UN]

1.94063641*kW

 (4.735)
> 

Units:-Standard:-`+`(E[D, TV, R], Units:-Standard:-`-`(E[D, TV, UN]))

2.20562478*kW

 (4.736)
> 

E[D, TV, EN, TV]

1.85189707*kW

 (4.737)
> 

Units:-Standard:-`+`(E[D, TV, R], Units:-Standard:-`-`(E[D, TV, EN, TV]))

2.29436412*kW

 (4.738)

%
> 

Units:-Standard:-`*`(E[D, TV, R], Units:-Standard:-`/`(E[D, TV, R]))

1.000000000

 (4.739)
> 

Units:-Standard:-`*`(E[D, TV, T], Units:-Standard:-`/`(E[D, TV, R]))

0.

 (4.740)
> 

Units:-Standard:-`*`(E[D, TV, UN], Units:-Standard:-`/`(E[D, TV, R]))

.4680449015

 (4.741)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, TV, R], Units:-Standard:-`-`(E[D, TV, UN])), Units:-Standard:-`/`(E[D, TV, R]))

.5319550986

 (4.742)
> 

Units:-Standard:-`*`(E[D, TV, EN, TV], Units:-Standard:-`/`(E[D, TV, R]))

.4466426463

 (4.743)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, TV, R], Units:-Standard:-`-`(E[D, TV, EN, TV])), Units:-Standard:-`/`(E[D, TV, R]))

.5533573538

 (4.744)

Evaporator, kW
> 

E[D, EV, R]

6.72046888*kW

 (4.745)
> 

E[D, EV, T]

1.78774184*kW

 (4.746)
> 

E[D, EV, UN]

2.36465729*kW

 (4.747)
> 

Units:-Standard:-`+`(E[D, EV, R], Units:-Standard:-`-`(E[D, EV, UN]))

4.35581159*kW

 (4.748)
> 

E[D, EV, EN, EV]

6.72046881*kW

 (4.749)
> 

Units:-Standard:-`+`(E[D, EV, R], Units:-Standard:-`-`(E[D, EV, EN, EV]))

0.7e-7*kW

 (4.750)

%
> 

Units:-Standard:-`*`(E[D, EV, R], Units:-Standard:-`/`(E[D, EV, R]))

.9999999998

 (4.751)
> 

Units:-Standard:-`*`(E[D, EV, T], Units:-Standard:-`/`(E[D, EV, R]))

.2660144510

 (4.752)
> 

Units:-Standard:-`*`(E[D, EV, UN], Units:-Standard:-`/`(E[D, EV, R]))

.3518589747

 (4.753)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, EV, R], Units:-Standard:-`-`(E[D, EV, UN])), Units:-Standard:-`/`(E[D, EV, R]))

.6481410251

 (4.754)
> 

Units:-Standard:-`*`(E[D, EV, EN, EV], Units:-Standard:-`/`(E[D, EV, R]))

.9999999894

 (4.755)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, EV, R], Units:-Standard:-`-`(E[D, EV, EN, EV])), Units:-Standard:-`/`(E[D, EV, R]))

0.1041593990e-7

 (4.756)

Overall system, kW
> 

E[D, TOT, R]

29.33924969*kW

 (4.757)
> 

E[D, TOT, T]

5.057332728*kW

 (4.758)
> 

E[D, TOT, UN]

9.263500925*kW

 (4.759)
> 

Units:-Standard:-`+`(E[D, TOT, R], Units:-Standard:-`-`(E[D, TOT, UN]))

20.07574876*kW

 (4.760)
> 

E[D, TOT, EN]

20.50092868*kW

 (4.761)
> 

Units:-Standard:-`+`(E[D, TOT, R], Units:-Standard:-`-`(E[D, TOT, EN]))

8.83832101*kW

 (4.762)

%
> 

Units:-Standard:-`*`(E[D, TOT, R], Units:-Standard:-`/`(E[D, TOT, R]))

1.000000000

 (4.763)
> 

Units:-Standard:-`*`(E[D, TOT, T], Units:-Standard:-`/`(E[D, TOT, R]))

.1723743034

 (4.764)
> 

Units:-Standard:-`*`(E[D, TOT, UN], Units:-Standard:-`/`(E[D, TOT, R]))

.3157374856

 (4.765)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, TOT, R], Units:-Standard:-`-`(E[D, TOT, UN])), Units:-Standard:-`/`(E[D, TOT, R]))

.6842625144

 (4.766)
> 

Units:-Standard:-`*`(E[D, TOT, EN], Units:-Standard:-`/`(E[D, TOT, R]))

.6987543614

 (4.767)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, TOT, R], Units:-Standard:-`-`(E[D, TOT, EN])), Units:-Standard:-`/`(E[D, TOT, R]))

.3012456387

 (4.768)

 

Table 3: Splitting the exergy destruction for components in advanced exergetic analysis

Component

Total in Real Cycle

Unavoidable Endogenous

Unavoidable

 Exogenous

Avoidable

 Endogenous

Avoidable

 Exogenous

Compressor,kW
> 

E[D, CM, R]

5.76127882*kW

 (4.769)
> 

E[D, UN, EN](CM)

.85040775*kW

 (4.770)
> 

Units:-Standard:-`+`(E[D, CM, UN], Units:-Standard:-`-`(E[D, UN, EN](CM)))

0.3297695e-1*kW

 (4.771)
> 

Units:-Standard:-`+`(E[D, CM, EN, CM], Units:-Standard:-`-`(E[D, UN, EN](CM)))

3.09395315*kW

 (4.772)
> 

Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CM, R], Units:-Standard:-`-`(E[D, CM, UN])), Units:-Standard:-`-`(Units:-Standard:-`+`(E[D, CM, EN, CM], Units:-Standard:-`-`(E[D, UN, EN](CM)))))

1.78394097*kW

 (4.773)

%
> 

Units:-Standard:-`*`(E[D, CM, R], Units:-Standard:-`/`(E[D, CM, R]))

1.000000000

 (4.774)
> 

Units:-Standard:-`*`(E[D, UN, EN](CD), Units:-Standard:-`/`(E[D, CM, R]))

.6043330872

 (4.775)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, CM, UN], Units:-Standard:-`-`(E[D, UN, EN](CM))), Units:-Standard:-`/`(E[D, CM, R]))

0.5723894127e-2

 (4.776)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, CM, EN, CM], Units:-Standard:-`-`(E[D, UN, EN](CM))), Units:-Standard:-`/`(E[D, CM, R]))

.5370254152

 (4.777)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CM, R], Units:-Standard:-`-`(E[D, CM, UN])), Units:-Standard:-`-`(Units:-Standard:-`+`(E[D, CM, EN, CM], Units:-Standard:-`-`(E[D, UN, EN](CM))))), Units:-Standard:-`/`(E[D, CM, R]))

.3096432278

 (4.778)

Condenser, kW
> 

E[D, CD, R]

12.71124080*kW

 (4.779)
> 

E[D, UN, EN](CD)

3.481731415*kW

 (4.780)
> 

Units:-Standard:-`+`(E[D, CD, UN], Units:-Standard:-`-`(E[D, UN, EN](CD)))

.593091110*kW

 (4.781)
> 

Units:-Standard:-`+`(E[D, CD, EN, CD], Units:-Standard:-`-`(E[D, UN, EN](CD)))

4.502470480*kW

 (4.782)
> 

Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CD, R], Units:-Standard:-`-`(E[D, CD, UN])), Units:-Standard:-`-`(Units:-Standard:-`+`(E[D, CD, EN, CD], Units:-Standard:-`-`(E[D, UN, EN](CD)))))

4.133947795*kW

 (4.783)

%
> 

Units:-Standard:-`*`(E[D, CD, R], Units:-Standard:-`/`(E[D, CD, R]))

1.000000000

 (4.784)
> 

Units:-Standard:-`*`(E[D, UN, EN](CD), Units:-Standard:-`/`(E[D, CD, R]))

.2739096418

 (4.785)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, CD, UN], Units:-Standard:-`-`(E[D, UN, EN](CD))), Units:-Standard:-`/`(E[D, CD, R]))

0.4665878960e-1

 (4.786)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, CD, EN, CD], Units:-Standard:-`-`(E[D, UN, EN](CD))), Units:-Standard:-`/`(E[D, CD, R]))

.3542117210

 (4.787)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, CD, R], Units:-Standard:-`-`(E[D, CD, UN])), Units:-Standard:-`-`(Units:-Standard:-`+`(E[D, CD, EN, CD], Units:-Standard:-`-`(E[D, UN, EN](CD))))), Units:-Standard:-`/`(E[D, CD, R]))

.3252198475

 (4.788)

Throttling valve, kW
> 

E[D, TV, R]

4.14626119*kW

 (4.789)
> 

E[D, UN, EN](TV)

1.85189699*kW

 (4.790)
> 

Units:-Standard:-`+`(E[D, TV, UN], Units:-Standard:-`-`(E[D, UN, EN](TV)))

0.8873942e-1*kW

 (4.791)
> 

Units:-Standard:-`+`(E[D, TV, EN, TV], Units:-Standard:-`-`(E[D, UN, EN](TV)))

0.8e-7*kW

 (4.792)
> 

Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, TV, R], Units:-Standard:-`-`(E[D, TV, UN])), Units:-Standard:-`-`(Units:-Standard:-`+`(E[D, TV, EN, TV], Units:-Standard:-`-`(E[D, UN, EN](TV)))))

2.20562470*kW

 (4.793)

%
> 

Units:-Standard:-`*`(E[D, TV, R], Units:-Standard:-`/`(E[D, TV, R]))

1.000000000

 (4.794)
> 

Units:-Standard:-`*`(E[D, UN, EN](TV), Units:-Standard:-`/`(E[D, TV, R]))

.4466426270

 (4.795)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, TV, UN], Units:-Standard:-`-`(E[D, UN, EN](TV))), Units:-Standard:-`/`(E[D, TV, R]))

0.2140227447e-1

 (4.796)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, TV, EN, TV], Units:-Standard:-`-`(E[D, UN, EN](TV))), Units:-Standard:-`/`(E[D, TV, R]))

0.1929449119e-7

 (4.797)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, TV, R], Units:-Standard:-`-`(E[D, TV, UN])), Units:-Standard:-`-`(Units:-Standard:-`+`(E[D, TV, EN, TV], Units:-Standard:-`-`(E[D, UN, EN](TV))))), Units:-Standard:-`/`(E[D, TV, R]))

.5319550793

 (4.798)

Evaporator, kW
> 

E[D, EV, R]

6.72046888*kW

 (4.799)
> 

E[D, UN, EN](EV)

2.36465744*kW

 (4.800)
> 

Units:-Standard:-`+`(E[D, EV, UN], Units:-Standard:-`-`(E[D, UN, EN](EV)))

-0.15e-6*kW

 (4.801)
> 

Units:-Standard:-`+`(E[D, EV, EN, EV], Units:-Standard:-`-`(E[D, UN, EN](EV)))

4.35581137*kW

 (4.802)
> 

Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, EV, R], Units:-Standard:-`-`(E[D, EV, UN])), Units:-Standard:-`-`(Units:-Standard:-`+`(E[D, EV, EN, EV], Units:-Standard:-`-`(E[D, UN, EN](EV)))))

0.22e-6*kW

 (4.803)

%
> 

Units:-Standard:-`*`(E[D, EV, R], Units:-Standard:-`/`(E[D, EV, R]))

.9999999998

 (4.804)
> 

Units:-Standard:-`*`(E[D, UN, EN](EV), Units:-Standard:-`/`(E[D, EV, R]))

.3518589970

 (4.805)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, EV, UN], Units:-Standard:-`-`(E[D, UN, EN](EV))), Units:-Standard:-`/`(E[D, EV, R]))

-0.2231987122e-7

 (4.806)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, EV, EN, EV], Units:-Standard:-`-`(E[D, UN, EN](EV))), Units:-Standard:-`/`(E[D, EV, R]))

.6481409924

 (4.807)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, EV, R], Units:-Standard:-`-`(E[D, EV, UN])), Units:-Standard:-`-`(Units:-Standard:-`+`(E[D, EV, EN, EV], Units:-Standard:-`-`(E[D, UN, EN](EV))))), Units:-Standard:-`/`(E[D, EV, R]))

0.3273581113e-7

 (4.808)

Overall system, kW
> 

E[D, TOT, R]

29.33924969*kW

 (4.809)
> 

E[D, TOT, UN, EN]

8.548693595*kW

 (4.810)
> 

Units:-Standard:-`+`(E[D, TOT, UN], Units:-Standard:-`-`(E[D, TOT, UN, EN]))

.714807330*kW

 (4.811)
> 

Units:-Standard:-`+`(E[D, TOT, EN], Units:-Standard:-`-`(E[D, TOT, UN, EN]))

11.95223508*kW

 (4.812)
> 

Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, TOT, R], Units:-Standard:-`-`(E[D, TOT, UN])), Units:-Standard:-`-`(Units:-Standard:-`+`(E[D, TOT, EN], Units:-Standard:-`-`(E[D, TOT, UN, EN]))))

8.12351368*kW

 (4.813)

%
> 

Units:-Standard:-`*`(E[D, TOT, R], Units:-Standard:-`/`(E[D, TOT, R]))

1.000000000

 (4.814)
> 

Units:-Standard:-`*`(E[D, TOT, UN, EN], Units:-Standard:-`/`(E[D, TOT, R]))

.2913739678

 (4.815)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, TOT, UN], Units:-Standard:-`-`(E[D, TOT, UN, EN])), Units:-Standard:-`/`(E[D, TOT, R]))

0.2436351773e-1

 (4.816)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(E[D, TOT, EN], Units:-Standard:-`-`(E[D, TOT, UN, EN])), Units:-Standard:-`/`(E[D, TOT, R]))

.4073803934

 (4.817)
> 

Units:-Standard:-`*`(Units:-Standard:-`+`(Units:-Standard:-`+`(E[D, TOT, R], Units:-Standard:-`-`(E[D, TOT, UN])), Units:-Standard:-`-`(Units:-Standard:-`+`(E[D, TOT, EN], Units:-Standard:-`-`(E[D, TOT, UN, EN])))), Units:-Standard:-`/`(E[D, TOT, R]))

.2768821209

 (4.818)

Functions for plotting Ts-diagram
 

Critical temperature

> 

T[cr] := TCRwf(wf[cycle]):

Critical pressure

> 

p[cr] := PCRwf(wf[cycle]):

 

> 

Ts[sl] := plot([Units:-Standard:-`*`(SSTXwf(t, 0, wf[cycle]), Units:-Standard:-`/`(1000)), t, t = Units:-Standard:-`*`(TSPXwf(p[4, R], 0, wf[cycle]), Units:-Standard:-`/`(Unit('K'))) .. Units:-Standard:-`*`(T[cr], Units:-Standard:-`/`(Unit('K')))], color = "Black", thickness = 1, linestyle = "longdash", legend = "Liquid on saturated line"):

 

> 

Ts[sv] := plot([Units:-Standard:-`*`(SSTXwf(t, 1, wf[cycle]), Units:-Standard:-`/`(1000)), t, t = Units:-Standard:-`*`(TSPXwf(p[4, R], 1, wf[cycle]), Units:-Standard:-`/`(Unit('K'))) .. Units:-Standard:-`*`(T[cr], Units:-Standard:-`/`(Unit('K')))], color = "Black", thickness = 1, linestyle = "solid", legend = "Vapor on saturated line"):

1. Real cycle

> 

Ts[1, Units:-Standard:-`+`(R, -2), R] := plot([[Units:-Standard:-`*`(ScompPTPEFFwf(Units:-Standard:-`*`(p[1, R], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, R], Units:-Standard:-`/`(Unit('K'))), p, eta[comp], wf[cycle]), Units:-Standard:-`/`(1000)), TcompPTPEFFwf(Units:-Standard:-`*`(p[1, R], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, R], Units:-Standard:-`/`(Unit('K'))), p, eta[comp], wf[cycle]), p = Units:-Standard:-`*`(p[1, R], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[2, R], Units:-Standard:-`/`(Unit('Pa')))]], color = "Red", thickness = 3, linestyle = "solid", legend = "Compression (real cycle)"):
> 

Ts[2, Units:-Standard:-`+`(R, Units:-Standard:-`-`(sv)), 2, R] := plot([Units:-Standard:-`*`(SPTwf(Units:-Standard:-`*`(p[2, R], Units:-Standard:-`/`(Unit('Pa'))), t, wf[cycle]), Units:-Standard:-`/`(1000)), t, t = Units:-Standard:-`*`(TSPXwf(p[2, R], 1, wf[cycle]), Units:-Standard:-`/`(Unit('K'))) .. Units:-Standard:-`*`(T[2, R], Units:-Standard:-`/`(Unit('K')))], color = "Magenta", thickness = 3, linestyle = "solid", legend = "Cooling in condenser (real cycle)"):

``

> 

Ts[sv, 2, Units:-Standard:-`+`(R, -3), R] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[2, R], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[2, R], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = 1 .. 0]], color = "Brown", thickness = 3, linestyle = "solid", legend = "Condensing (real cycle)"):

 

> 

Ts[3, Units:-Standard:-`+`(R, -4), R] := plot([[Units:-Standard:-`*`(SPHwf(p, Units:-Standard:-`*`(h[4, R], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`('kg'))))), wf[cycle]), Units:-Standard:-`/`(1000)), TPHwf(p, Units:-Standard:-`*`(h[4, R], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`('kg'))))), wf[cycle]), p = Units:-Standard:-`*`(p[3, R], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[4, R], Units:-Standard:-`/`(Unit('Pa')))]], color = "Blue", thickness = 3, linestyle = "solid", legend = "Expansion (real cycle)"):

 

> 

Ts[4, Units:-Standard:-`+`(R, -1), R] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[1, R], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[1, R], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = x[4, R] .. 1]], color = "DeepSkyBlue", thickness = 3, linestyle = "solid", legend = "Evaporation (real cycle)"):

2. Unavoidable cycle

> 

Ts[1, Units:-Standard:-`+`(UN, -2), UN] := plot([[Units:-Standard:-`*`(ScompPTPEFFwf(Units:-Standard:-`*`(p[1, UN], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, UN], Units:-Standard:-`/`(Unit('K'))), p, eta[comp, UN], wf[cycle]), Units:-Standard:-`/`(1000)), TcompPTPEFFwf(Units:-Standard:-`*`(p[1, UN], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, UN], Units:-Standard:-`/`(Unit('K'))), p, eta[comp, UN], wf[cycle]), p = Units:-Standard:-`*`(p[1, UN], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[2, UN], Units:-Standard:-`/`(Unit('Pa')))]], color = "Red", thickness = 3, linestyle = "dash", legend = "Compression (unavoidable cycle)"):
> 

Ts[2, Units:-Standard:-`+`(UN, Units:-Standard:-`-`(sv)), 2, UN] := plot([Units:-Standard:-`*`(SPTwf(Units:-Standard:-`*`(p[2, UN], Units:-Standard:-`/`(Unit('Pa'))), t, wf[cycle]), Units:-Standard:-`/`(1000)), t, t = Units:-Standard:-`*`(TSPXwf(p[2, UN], 1, wf[cycle]), Units:-Standard:-`/`(Unit('K'))) .. Units:-Standard:-`*`(T[2, UN], Units:-Standard:-`/`(Unit('K')))], color = "Magenta", thickness = 3, linestyle = "dash", legend = "Cooling in condenser (unavoidable cycle)"):

NULL

> 

Ts[sv, 2, Units:-Standard:-`+`(UN, -3), UN] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[2, UN], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[2, UN], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = 1 .. 0]], color = "Brown", thickness = 3, linestyle = "dash", legend = "Condensing (unavoidable cycle)"):

 

> 

Ts[3, Units:-Standard:-`+`(UN, -4), UN] := plot([[Units:-Standard:-`*`(SPHwf(p, Units:-Standard:-`*`(h[4, UN], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`('kg'))))), wf[cycle]), Units:-Standard:-`/`(1000)), TPHwf(p, Units:-Standard:-`*`(h[4, UN], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`('kg'))))), wf[cycle]), p = Units:-Standard:-`*`(p[3, UN], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[4, UN], Units:-Standard:-`/`(Unit('Pa')))]], color = "Blue", thickness = 3, linestyle = "dash", legend = "Expansion (unavoidable cycle)"):

 

> 

Ts[4, Units:-Standard:-`+`(UN, -1), UN] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[1, UN], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[1, UN], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = x[4, UN] .. 1]], color = "DeepSkyBlue", thickness = 3, linestyle = "dash", legend = "Evaporation (unavoidable cycle)"):

3. Theoretical cycle

> 

Ts[1, Units:-Standard:-`+`(T, -2), T] := plot([[Units:-Standard:-`*`(ScompPTPEFFwf(Units:-Standard:-`*`(p[1, T], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, T], Units:-Standard:-`/`(Unit('K'))), p, 1, wf[cycle]), Units:-Standard:-`/`(1000)), TcompPTPEFFwf(Units:-Standard:-`*`(p[1, T], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, T], Units:-Standard:-`/`(Unit('K'))), p, 1, wf[cycle]), p = Units:-Standard:-`*`(p[1, T], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[2, T], Units:-Standard:-`/`(Unit('Pa')))]], color = "Red", thickness = 3, linestyle = "dot", legend = "Compression (theoretical cycle)"):
> 

Ts[2, Units:-Standard:-`+`(T, Units:-Standard:-`-`(sv)), 2, T] := plot([Units:-Standard:-`*`(SPTwf(Units:-Standard:-`*`(p[2, T], Units:-Standard:-`/`(Unit('Pa'))), t, wf[cycle]), Units:-Standard:-`/`(1000)), t, t = Units:-Standard:-`*`(TSPXwf(p[2, T], 1, wf[cycle]), Units:-Standard:-`/`(Unit('K'))) .. Units:-Standard:-`*`(T[2, T], Units:-Standard:-`/`(Unit('K')))], color = "Magenta", thickness = 3, linestyle = "dot", legend = "Cooling in condenser (theoretical cycle)"):

NULL

> 

Ts[sv, 2, Units:-Standard:-`+`(T, -3), T] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[2, T], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[2, T], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = 1 .. 0]], color = "Brown", thickness = 3, linestyle = "dot", legend = "Condensing (theoretical cycle)"):

 

> 

Ts[3, Units:-Standard:-`+`(T, -4), T] := plot([[Units:-Standard:-`*`(Units:-Standard:-`*`(s[4, T], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`(Units:-Standard:-`*`('kg', 'K')))))), Units:-Standard:-`/`(1000)), TPSwf(p, Units:-Standard:-`*`(s[4, T], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`(Units:-Standard:-`*`('kg', 'K')))))), wf[cycle]), p = Units:-Standard:-`*`(p[3, T], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[4, T], Units:-Standard:-`/`(Unit('Pa')))]], color = "Blue", thickness = 3, linestyle = "dot", legend = "Expansion (theoretical cycle)"):

 

> 

Ts[4, Units:-Standard:-`+`(T, -1), T] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[1, T], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[1, T], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = x[4, T] .. 1]], color = "DeepSkyBlue", thickness = 3, linestyle = "dot", legend = "Evaporation (theoretical cycle)"):

4. Hybrid cycle with irreversible component - Compressor

> 

Ts[1, Units:-Standard:-`+`(`EN,CM`, -2), `EN,CM`] := plot([[Units:-Standard:-`*`(ScompPTPEFFwf(Units:-Standard:-`*`(p[1, `EN,CM`], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, `EN,CM`], Units:-Standard:-`/`(Unit('K'))), p, eta[comp], wf[cycle]), Units:-Standard:-`/`(1000)), TcompPTPEFFwf(Units:-Standard:-`*`(p[1, `EN,CM`], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, `EN,CM`], Units:-Standard:-`/`(Unit('K'))), p, eta[comp], wf[cycle]), p = Units:-Standard:-`*`(p[1, `EN,CM`], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[2, `EN,CM`], Units:-Standard:-`/`(Unit('Pa')))]], color = "Red", thickness = 3, linestyle = "dash", legend = "Compression (cycle with endogenues exergy destruction in compressor)"):
> 

Ts[2, Units:-Standard:-`+`(`EN,CM`, Units:-Standard:-`-`(sv)), 2, `EN,CM`] := plot([Units:-Standard:-`*`(SPTwf(Units:-Standard:-`*`(p[2, `EN,CM`], Units:-Standard:-`/`(Unit('Pa'))), t, wf[cycle]), Units:-Standard:-`/`(1000)), t, t = Units:-Standard:-`*`(TSPXwf(p[2, `EN,CM`], 1, wf[cycle]), Units:-Standard:-`/`(Unit('K'))) .. Units:-Standard:-`*`(T[2, `EN,CM`], Units:-Standard:-`/`(Unit('K')))], color = "Magenta", thickness = 3, linestyle = "dash", legend = "Cooling in condenser (cycle with endogenues exergy destruction in compressor)"):

NULL

> 

Ts[sv, 2, Units:-Standard:-`+`(`EN,CM`, -3), `EN,CM`] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[2, `EN,CM`], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[2, `EN,CM`], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = 1 .. 0]], color = "Brown", thickness = 3, linestyle = "dash", legend = "Condensing (cycle with endogenues exergy destruction in compressor)"):

 

> 

Ts[3, Units:-Standard:-`+`(`EN,CM`, -4), `EN,CM`] := plot([[Units:-Standard:-`*`(Units:-Standard:-`*`(s[4, `EN,CM`], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`(Units:-Standard:-`*`('kg', 'K')))))), Units:-Standard:-`/`(1000)), TPSwf(p, Units:-Standard:-`*`(s[4, `EN,CM`], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`(Units:-Standard:-`*`('kg', 'K')))))), wf[cycle]), p = Units:-Standard:-`*`(p[3, `EN,CM`], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[4, `EN,CM`], Units:-Standard:-`/`(Unit('Pa')))]], color = "Blue", thickness = 3, linestyle = "dash", legend = "Expansion (cycle with endogenues exergy destruction in compressor)"):

 

> 

Ts[4, Units:-Standard:-`+`(`EN,CM`, -1), `EN,CM`] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[1, `EN,CM`], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[1, `EN,CM`], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = x[4, `EN,CM`] .. 1]], color = "DeepSkyBlue", thickness = 3, linestyle = "dash", legend = "Evaporation (cycle with endogenues exergy destruction in compressor)"):

5. Hybrid cycle with irreversible component - Condenser

> 

Ts[1, EN, Units:-Standard:-`+`(CD, -2), EN, CD] := plot([[Units:-Standard:-`*`(ScompPTPEFFwf(Units:-Standard:-`*`(p[1, `EN,CD`], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, `EN,CD`], Units:-Standard:-`/`(Unit('K'))), p, 1, wf[cycle]), Units:-Standard:-`/`(1000)), TcompPTPEFFwf(Units:-Standard:-`*`(p[1, `EN,CD`], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, `EN,CD`], Units:-Standard:-`/`(Unit('K'))), p, 1, wf[cycle]), p = Units:-Standard:-`*`(p[1, `EN,CD`], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[2, EN, CD], Units:-Standard:-`/`(Unit('Pa')))]], color = "Red", thickness = 7, linestyle = "dot", legend = "Compression (cycle with endogenues exergy destruction in condenser)"):
> 

Ts[2, EN, Units:-Standard:-`+`(CD, Units:-Standard:-`-`(sv)), 2, EN, CD] := plot([Units:-Standard:-`*`(SPTwf(Units:-Standard:-`*`(p[2, EN, CD], Units:-Standard:-`/`(Unit('Pa'))), t, wf[cycle]), Units:-Standard:-`/`(1000)), t, t = Units:-Standard:-`*`(TSPXwf(p[2, EN, CD], 1, wf[cycle]), Units:-Standard:-`/`(Unit('K'))) .. Units:-Standard:-`*`(T[2, `EN,CD`], Units:-Standard:-`/`(Unit('K')))], color = "Magenta", thickness = 7, linestyle = "dot", legend = "Cooling in condenser (cycle with endogenues exergy destruction in condenser)"):

NULL

> 

Ts[sv, 2, EN, Units:-Standard:-`+`(CD, -3), EN, CD] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[2, EN, CD], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[2, EN, CD], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = 1 .. 0]], color = "Brown", thickness = 7, linestyle = "dot", legend = "Condensing (cycle with endogenues exergy destruction in condenser)"):

 

> 

Ts[3, EN, Units:-Standard:-`+`(CD, -4), EN, CD] := plot([[Units:-Standard:-`*`(Units:-Standard:-`*`(s[4, EN, CD], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`(Units:-Standard:-`*`('kg', 'K')))))), Units:-Standard:-`/`(1000)), TPSwf(p, Units:-Standard:-`*`(s[4, EN, CD], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`(Units:-Standard:-`*`('kg', 'K')))))), wf[cycle]), p = Units:-Standard:-`*`(p[3, EN, CD], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[4, EN, CD], Units:-Standard:-`/`(Unit('Pa')))]], color = "Blue", thickness = 7, linestyle = "dot", legend = "Expansion (cycle with endogenues exergy destruction in condenser)"):

 

> 

Ts[4, EN, Units:-Standard:-`+`(CD, -1), EN, CD] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[1, `EN,CD`], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[1, `EN,CD`], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = x[4, `EN,CD`] .. 1]], color = "DeepSkyBlue", thickness = 7, linestyle = "dot", legend = "Evaporation (cycle with endogenues exergy destruction in condenser)"):

6. Hybrid cycle with irreversible component - Throttling Valve

> 

Ts[1, EN, Units:-Standard:-`+`(TV, -2), EN, TV] := plot([[Units:-Standard:-`*`(ScompPTPEFFwf(Units:-Standard:-`*`(p[1, `EN,TV`], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, `EN,TV`], Units:-Standard:-`/`(Unit('K'))), p, 1, wf[cycle]), Units:-Standard:-`/`(1000)), TcompPTPEFFwf(Units:-Standard:-`*`(p[1, `EN,TV`], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, `EN,TV`], Units:-Standard:-`/`(Unit('K'))), p, 1, wf[cycle]), p = Units:-Standard:-`*`(p[1, `EN,TV`], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[2, EN, TV], Units:-Standard:-`/`(Unit('Pa')))]], color = "Red", thickness = 1, linestyle = "solid", symbol = "asterisk", legend = "Compression (cycle with endogenues exergy destruction in throttling valve)"):
> 

Ts[2, EN, Units:-Standard:-`+`(TV, Units:-Standard:-`-`(sv)), 2, EN, TV] := plot([Units:-Standard:-`*`(SPTwf(Units:-Standard:-`*`(p[2, EN, TV], Units:-Standard:-`/`(Unit('Pa'))), t, wf[cycle]), Units:-Standard:-`/`(1000)), t, t = Units:-Standard:-`*`(TSPXwf(p[2, EN, TV], 1, wf[cycle]), Units:-Standard:-`/`(Unit('K'))) .. Units:-Standard:-`*`(T[2, `EN,TV`], Units:-Standard:-`/`(Unit('K')))], color = "Magenta", thickness = 3, linestyle = "dash", legend = "Cooling in condenser (cycle with endogenues exergy destruction in throttling valve)"):

NULL

> 

Ts[sv, 2, EN, Units:-Standard:-`+`(TV, -3), EN, TV] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[2, EN, TV], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[2, EN, TV], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = 1 .. 0]], color = "Brown", thickness = 3, linestyle = "dash", legend = "Condensing (cycle with endogenues exergy destruction in throttling valve)"):

 

> 

Ts[3, EN, Units:-Standard:-`+`(TV, -4), EN, TV] := plot([[Units:-Standard:-`*`(SPHwf(p, Units:-Standard:-`*`(h[4, EN, TV], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`('kg'))))), wf[cycle]), Units:-Standard:-`/`(1000)), TPHwf(p, Units:-Standard:-`*`(h[4, EN, TV], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`('kg'))))), wf[cycle]), p = Units:-Standard:-`*`(p[3, EN, TV], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[4, EN, TV], Units:-Standard:-`/`(Unit('Pa')))]], color = "Blue", thickness = 3, linestyle = "dash", legend = "Expansion (cycle with endogenues exergy destruction in throttling valve)"):

 

> 

Ts[4, EN, Units:-Standard:-`+`(TV, -1), EN, TV] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[1, `EN,TV`], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[1, `EN,TV`], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = x[4, `EN,TV`] .. 1]], color = "DeepSkyBlue", thickness = 3, linestyle = "dash", legend = "Evaporation (cycle with endogenues exergy destruction in throttling valve)"):

 

  NULL

7. Hybrid cycle with irreversible component - Evaporator

> 

Ts[1, EN, Units:-Standard:-`+`(EV, -2), EN, EV] := plot([[Units:-Standard:-`*`(ScompPTPEFFwf(Units:-Standard:-`*`(p[1, EN, EV], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, EN, EV], Units:-Standard:-`/`(Unit('K'))), p, 1, wf[cycle]), Units:-Standard:-`/`(1000)), TcompPTPEFFwf(Units:-Standard:-`*`(p[1, EN, EV], Units:-Standard:-`/`(Unit('Pa'))), Units:-Standard:-`*`(T[1, EN, EV], Units:-Standard:-`/`(Unit('K'))), p, 1, wf[cycle]), p = Units:-Standard:-`*`(p[1, EN, EV], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[2, EN, EV], Units:-Standard:-`/`(Unit('Pa')))]], color = "Red", thickness = 1, linestyle = "dot", legend = "Compression (cycle with endogenues exergy destruction in evaporator)"):
> 

Ts[2, EN, Units:-Standard:-`+`(EV, Units:-Standard:-`-`(sv)), 2, EN, EV] := plot([Units:-Standard:-`*`(SPTwf(Units:-Standard:-`*`(p[2, EN, EV], Units:-Standard:-`/`(Unit('Pa'))), t, wf[cycle]), Units:-Standard:-`/`(1000)), t, t = Units:-Standard:-`*`(TSPXwf(p[2, EN, EV], 1, wf[cycle]), Units:-Standard:-`/`(Unit('K'))) .. Units:-Standard:-`*`(T[2, EN, EV], Units:-Standard:-`/`(Unit('K')))], color = "Magenta", thickness = 1, linestyle = "dot", legend = "Cooling in condenser (cycle with endogenues exergy destruction in evaporator)"):

NULL

> 

Ts[sv, 2, EN, Units:-Standard:-`+`(EV, -3), EN, EV] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[2, EN, EV], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[2, EN, EV], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = 1 .. 0]], color = "Brown", thickness = 1, linestyle = "dot", legend = "Condensing (cycle with endogenues exergy destruction in evaporator)"):

 

> 

Ts[3, EN, Units:-Standard:-`+`(EV, -4), EN, EV] := plot([[Units:-Standard:-`*`(Units:-Standard:-`*`(s[4, EN, EV], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`(Units:-Standard:-`*`('kg', 'K')))))), Units:-Standard:-`/`(1000)), TPSwf(p, Units:-Standard:-`*`(s[4, EN, EV], Units:-Standard:-`/`(Unit(Units:-Standard:-`*`('J', Units:-Standard:-`/`(Units:-Standard:-`*`('kg', 'K')))))), wf[cycle]), p = Units:-Standard:-`*`(p[3, EN, EV], Units:-Standard:-`/`(Unit('Pa'))) .. Units:-Standard:-`*`(p[4, EN, EV], Units:-Standard:-`/`(Unit('Pa')))]], color = "Blue", thickness = 1, linestyle = "dot", legend = "Expansion (cycle with endogenues exergy destruction in evaporator)"):

 

> 

Ts[4, EN, Units:-Standard:-`+`(EV, -1), EN, EV] := plot([[Units:-Standard:-`*`(SSPXwf(Units:-Standard:-`*`(p[1, EN, EV], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), Units:-Standard:-`/`(1000)), TSPXwf(Units:-Standard:-`*`(p[1, EN, EV], Units:-Standard:-`/`(Unit('Pa'))), X, wf[cycle]), X = x[4, `EN,EV`] .. 1]], color = "DeepSkyBlue", thickness = 1, linestyle = "dot", legend = "Evaporation (cycle with endogenues exergy destruction in evaporator)"):

Ploting the Refrigeration Cycle on a Ts-diagram
 
 
> 

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

 
> 

plots:-display(Ts[sl], Ts[sv], Ts[1, Units:-Standard:-`+`(R, -2), R], Ts[2, Units:-Standard:-`+`(R, Units:-Standard:-`-`(sv)), 2, R], Ts[sv, 2, Units:-Standard:-`+`(R, -3), R], Ts[3, Units:-Standard:-`+`(R, -4), R], Ts[4, Units:-Standard:-`+`(R, -1), R], Ts[1, Units:-Standard:-`+`(`EN,CM`, -2), `EN,CM`], Ts[2, Units:-Standard:-`+`(`EN,CM`, Units:-Standard:-`-`(sv)), 2, `EN,CM`], Ts[sv, 2, Units:-Standard:-`+`(`EN,CM`, -3), `EN,CM`], Ts[3, Units:-Standard:-`+`(`EN,CM`, -4), `EN,CM`], Ts[4, Units:-Standard:-`+`(`EN,CM`, -1), `EN,CM`], Ts[1, EN, Units:-Standard:-`+`(CD, -2), EN, CD], Ts[2, EN, Units:-Standard:-`+`(CD, Units:-Standard:-`-`(sv)), 2, EN, CD], Ts[sv, 2, EN, Units:-Standard:-`+`(CD, -3), EN, CD], Ts[3, EN, Units:-Standard:-`+`(CD, -4), EN, CD], Ts[4, EN, Units:-Standard:-`+`(CD, -1), EN, CD], Ts[1, EN, Units:-Standard:-`+`(TV, -2), EN, TV], Ts[2, EN, Units:-Standard:-`+`(TV, Units:-Standard:-`-`(sv)), 2, EN, TV], Ts[sv, 2, EN, Units:-Standard:-`+`(TV, -3), EN, TV], Ts[3, EN, Units:-Standard:-`+`(TV, -4), EN, TV], Ts[4, EN, Units:-Standard:-`+`(TV, -1), EN, TV], Ts[1, EN, Units:-Standard:-`+`(EV, -2), EN, EV], Ts[2, EN, Units:-Standard:-`+`(EV, Units:-Standard:-`-`(sv)), 2, EN, EV], Ts[sv, 2, EN, Units:-Standard:-`+`(EV, -3), EN, EV], Ts[3, EN, Units:-Standard:-`+`(EV, -4), EN, EV], Ts[4, EN, Units:-Standard:-`+`(EV, -1), EN, EV], labels = [Units:-Standard:-`*`(s, [Units:-Standard:-`*`(kJ, Units:-Standard:-`/`(Units:-Standard:-`*`(kg, K)))]), Units:-Standard:-`*`(T, [K])], size = [800, 900], gridlines)