九州大学学術情報リポジトリ

Kyushu University Institutional Repository

Thermodynamic Feasibility Evaluation of a Novel

Lowtemperature Ejector-based Trans-critical
R744 Refrigeration Cycle
Mojtaba Purjam
Department of Energy and Environmental Engineering, Interdisciplinary Graduate School of
Engineering Sciences, Kyushu University

Kyaw Thu
Department of Energy and Environmental Engineering, Interdisciplinary Graduate School of
Engineering Sciences, Kyushu University

Miyazaki, Takahiko
Department of Energy and Environmental Engineering, Interdisciplinary Graduate School of
Engineering Sciences, Kyushu University

https://doi.org/10.5109/4372280

出版情報：Evergreen. 8 (1), pp.204-212, 2021-03. Transdisciplinary Research and Education

Center for Green Technologies, Kyushu University
バージョン：
権利関係：
 EVERGREEN Joint Journal of Novel Carbon Resource Sciences & Green Asia Strategy, Vol. 08, Issue 01, pp204-212, March 2021

 
Thermodynamic Feasibility Evaluation of a Novel Low-
temperature Ejector-based Trans-critical R744 Refrigeration
Cycle
Mojtaba Purjam1, *, Kyaw Thu1,2, Takahiko Miyazaki1,2
1
Department of Energy and Environmental Engineering, Interdisciplinary Graduate School of Engineering
Sciences, Kyushu University, 6-1 Kasuga-koen, Kasuga-shi, Fukuoka, 816-8580, Japan
2
Thermal Science and Engineering Division, International Institute of Carbon-Neutral Energy Research
(I2CNER), Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395, Japan

*Author to whom correspondence should be addressed:

E-mail: [email protected]

(Received October 30, 2020; Revised March 25, 2021; accepted March 26, 2021).

Abstract: In conventional ejector refrigeration cycles, the ejector mostly operates at fixed
pressure points in the secondary stream and outlet. Studying the sensitivity of the system to the
pressure alternation of the mentioned points has not been investigated. To cover this gap, an extra
turboexpander, compressor, and gas cooler were introduced to the conventional ejector layout. It
was found that pressurizing the secondary entrance of the ejector could be beneficial to COP if it is
not warmer than the gas coolers. Also using the expansion process right after the ejector has the
potential to increase the efficiency of the system if the discharge pressure of the ejector reaches
supercritical or superheat condition.

Keywords: Ejector; Trans-critical CO2; Low-Temperature Refrigeration; Turboexpander

1. Introduction and empirical research over an R744 ejector cycle and

1),2) announced an 8% boost in COP compared to a
Lorentzen and Petteresen considered carbon dioxide
simple throttling valve system. Zhu et al.19) studied a
as an unharmful and environmentally friendly refrigerant
transcritical R744 ejector heat pump-based water
that has several positive characteristics such as low
heating cycle in the laboratory and proved 10.3% of
toxicity, inexplosive, high capacity to mass ratio,
increased performance over the single-stage cycle. In a
lubricants-compatible, cheap, and unique heat transfer
transcritical cycle, gas cooler pressure has a significant
properties. In comparison to conventional heat pump
effect on COP. There is an optimum pressure that
cycles, the compressors of trans-critical CO2 cycles must
maximizes the COP of the system20),21),22). Elbel17)
provide higher lift3). Because of this, their COP is usually
reviewed and studied few unconventional layouts for
smaller than other ones with organofluorine-based
transcritical R744 cycles in a nutshell.
refrigerants. In order to enhance the COP, researchers
One of the most energy-consuming parts of the
tried different layouts and components such as internal
refrigeration process in the industries is the low-
heat exchangers4), cascade cycles5),6), turboexpanders7),8),
temperature refrigeration23). These systems have
multistage compression9),10),11), vortex tube expansions12),
widespread application, energy-intensive, and usually
subcooling boost13),14).
work non-stop at high capacities. In evaporation
Along with the described components and methods, the
temperatures below -35℃ the COPs of most refrigeration
ejectors’ implementation approach is favored. Li and
cycles drop under 1. Until now there is not a significant
Groll15) and Deng et al.16) suppositionally investigated an
specific investigation for using CO2 cycles in this
ejector in a transcritical CO2 cycle. Li and
temperature domain. So, proposing and studying new
Groll15) operated a dual-phase ejector for a
cycles for these temperatures will fill some gaps in
transcritical HVAC system and reported a performance
transcritical CO2 heat pumping knowledge. From a
gain of up to 16%. It was surfaced that the highest COP in
thermodynamic point of view, if a cycle can present
the ejector cycle was 22% better than other simple vapor
reasonable performance in low temperatures, it has a good
compression cycles16). Elbel17) reported that by using an
chance to be highly efficient in conventional cooling
ejector COP and cooling capacity enhanced up to 7% and
temperatures. Investigation on the reaction of the cycle to
8% respectively. Banasiak et al.18) conducted a simulation

- 204 -
 EVERGREEN Joint Journal of Novel Carbon Resource Sciences & Green Asia Strategy, Vol. 08, Issue 01, pp204-212, March 2021

 
changing the pressure of 2 inlets and 1 outlet of the ejector pressure manipulation gave the research team the chance
is another target of this study. to investigate the feasibility of simultaneous application
This paper proposed 3 extra components to the of ejector and turboexpanders in this layout.
conventional CO2 ejector refrigeration cycle and tried to In the future, full parametric 1st law, 2nd law24) and
study the feasibility of the new layout by study the thermo-economics analysis will be done and
pressure of points before and after the ejector. It is the first optimization25) algorithm will be considered in the codes.
step for the parametric analysis.
The coexistence of ejectors and expansion devices in a 2. The cycle’s mathematical model
refrigeration cycle has never been reported. It’s mainly
2.1. The schematics and governing equations
since the pressure difference between the flash chamber
and evaporator is not significant. However, there is not Fig.1 presents the schematics of the purposed ejector
any mathematical model that proves this. Magnificent cycle for the deep-freezing cycle.
pressure lifts in the transcritical R744 cycle and the ejector

Fig. 1: Schematics of the simulated cycle

Fig. 2: Program’s flowchart in Phyton

- 205 -
 Thermodynamic Feasibility Evaluation of a Novel Low-temperature Ejector-based Trans-critical R744 Refrigeration Cycle

 
2.2.1 Compressors and expansion devices
During compression, the fluid’s pressure is increased
by compressors. During expansion, the fluid’s pressure is
reduced by the turboexpanders. Both processes are
adiabatic. However, compressors are the work consumers
and turboexpanders are the work producers.
The governing equations for compression process:

ℎ 𝑃 ,𝑠 ℎ
ℎ ℎ (3)
𝜂

𝑤 𝑥 ℎ ℎ 1 𝑥 ℎ ℎ (4)
Fig. 3: T-s diagram of the cycle according to Table 1 (1st
scenario) The governing equations for the expansion process:

ℎ ℎ 𝜂 ℎ ℎ 𝑃 ,𝑠 (5)

𝑤 ℎ ℎ 1 𝑥 ℎ ℎ (6)

2.2.2 Evaporator
The specific cooling capacity of the cycle:

𝑞 1 𝑥 ℎ ℎ (7)

The mass flow rate of refrigerant in the evaporator per

Fig. 4: P-h diagram of the cycle according to Table 1 (1st every refrigeration ton:
scenario)  

3516.9
2.2 1st law governing equations 𝑚 (8)
𝑞
Capacity-independent thermodynamic simulation is the
target of this study. Though, the machinery and friction
loss is neglected. 2.2.3 Ejector’s governing equations
Fig. 3 and Fig. 4 show the T-s and P-h diagram of the
In the ejector, the primary flow’s pressure-drop after the
cycle according to Table 1.
primary nozzle creates a vacuum which in turn sucks in
The first law of thermodynamics for a steady-state
the secondary flow. After that both flows are mixed in the
steady flow control volume with the negligible potential
mixing section and finally are discharged out of the ejector
energy of the flow is:
by the diffuser. All of the mentioned internal parts of
ejector were assumed to be adiabatic except the mixing
𝑉
𝑄 𝑊 𝑚 ℎ section which is also isobaric.
2
(1) 2.2.3.1 Nozzles’ modeling
𝑉
𝑚 ℎ Nozzles have an isentropic efficiency. With current
2
assumptions, the nozzles’ thermodynamic oriented
simulation can be used for the modeling of supersonic and
Mass independent form of the above equation: subsonic velocities.

𝑉 𝑉 ℎ ℎ 𝜂 ℎ ℎ 𝑃 ,𝑠 (9)
𝑞 𝑤 ℎ ℎ (2)
2 2

.
𝑉 2 ℎ ℎ (10)
The thermophysical properties of R744 were extracted
from CoolProp libraries based on non-ideal (non-real gas)
conditions26). The simulation environment was mostly
Python. The same equations were applied to the secondary
nozzle.

- 206 -
 EVERGREEN Joint Journal of Novel Carbon Resource Sciences & Green Asia Strategy, Vol. 08, Issue 01, pp204-212, March 2021

 
2.2.3.2 Speed of sound in the 2-phase fluid Conservation of energy:
Niguen et al.27) addressed a mathematical model to
calculate the speed of sound in the 2-phase condition. 𝑉 𝑉
ℎ ℎ (19)
2 2
.
1 𝛼 𝛼𝜌
𝐶 1 𝛼 Thermodynamic state equation for density:
𝐶 𝜌 𝐶
. (11)
𝛼 1 𝛼 𝜌 𝜌 𝜌 𝑃 ,ℎ (20)
𝛼
𝐶 𝜌𝐶
Density’s initial guessing value.
The α is the void fraction of sub-critical dual-phase
fluid. 𝛾 1 𝑀
𝜌 𝜌 (21)
𝑥𝜌 𝛾 1 𝑀 2
𝛼 (12)
𝑥𝜌 1 𝑥 𝜌
In which γ for single-phase:
2.2.3.3 Mixing section of ejector
𝐶 𝑃 ,ℎ
In the current study, the entrainment ratio of the ejector 𝛾 (22)
is decided by the vapor quality before the flash chamber. 𝐶 𝑃 ,ℎ

𝑚 1 𝑥 For dual phase:

𝑅 (13)
𝑚 𝑥
𝐶 𝑃 ,𝑥 1
𝛾 𝑥
From the conservation of momentum in an isobaric 𝐶 𝑃 ,𝑥 1
(23)
medium. 𝐶 𝑃 ,𝑥 0
1 𝑥
𝐶 𝑃 ,𝑥 0
𝑉 𝑅 𝑉 1 𝑅 𝑉 . (14)

2.3.3.5 Diffuser
Conservation of energy:
𝑉
𝑉 𝑉 ℎ ℎ (24)
ℎ 𝑅 ℎ 2
2 2
𝑉 ,
1 𝑅 ℎ ,
(15) ℎ , ℎ 𝜂 ℎ ℎ (25)
2
𝑉
1 𝑅 ℎ 𝑃 𝑃 ℎ ,𝑠 (26)
2 ,

𝑉 𝑇 𝑇 𝑃 ,ℎ (27)
𝜂 (16)
𝑉 ,
  2.2.4 Coefficient of performance

𝑉 and ℎ are the results in real condition. 𝑞

𝐶𝑂𝑃 (28)
𝑤 𝑤
2.2.3.4 Normal Shockwave inside the constant area
section
𝑇
If the flow velocity after the mixing section exceeds the COPCarnot (29)
speed of the sound, a normal shock will happen in the 𝑇Sink 𝑇
constant area section. It starts with the conservation of
mass28):
The sub-critical 2-phase fluid enters the flash chamber.
The saturated vapor is separated from the saturated liquid
𝜌 𝑉 𝜌 𝑉 (17)
by the flash chamber in an adiabatic isobaric process.
Evaporator and gas coolers are isobaric.
Conservation of momentum: The fixed values in this paper are efficiencies and
source and sink temperatures. The efficiencies for nozzles,
𝑃 𝜌 𝑉 𝑃 𝜌 𝑉 𝑉 (18) mixing section, diffuser, compressors, and
turboexpanders are 0.9, 0.85, 0.85, 0.75, 0.63 respectively.

- 207 -
 Thermodynamic Feasibility Evaluation of a Novel Low-temperature Ejector-based Trans-critical R744 Refrigeration Cycle

 
The sink temperature is 303K and the source temperature 3. Results and discussions
is 233K. The least temperature difference for gas coolers
Table 1 shows the properties of every point in the
and the evaporator are 10 and 5 degrees respectively.
schematics with some extra information. This is an
Fig. 2 presents the flow chart that calculates the COP of
estimated optimal operation point of this cycle for the
the cycle. The entrainment ratio is the most important part
fixed source and sink temperatures.
of the simulation and it is must be calculated before
calculation of the COP.
Table 1: The IIR referenced thermo-physical information the points in the schematics in the first scenario
State # T (K) P (MPa) h (kJ kg-1) s (kJ kg-1 K-1) Phase or x (–) Mach Number
1 313.0 9.750 316.3 1.368 Supercritical 0
aa 278.4 3.990 304.2 1.373 0.4256 0.87
2 313.0 4.200 477.4 1.975 Superheat 0
bb 309.3 3.990 475.3 1.979 Superheat 0.27
cc 278.4 3.990 355.8 1.558 0.6664 0.63
dd 278.4 3.990 355.8 1.558 0.6664 0.63
3 290.0 5.324 364.0 1.562 0.7042 0
4 290.0 5.324 364.0 1.562 0.7042 0
5 290.0 5.324 413.7 1.734 1 0
6 337.3 9.750 440.6 1.754 Supercritical 0
7 290.0 5.324 245.8 1.155 0 0
8 228.0 0.827 232.2 1.190 0.391 0
9 228.0 0.827 434.1 2.075 1 0
10 361.4 4.20 534.0 2.147 Superheat 0
COP 1.342 𝐶𝑂𝑃 3.329
𝑚 17.4 gr/s.TR R 0.42

Table 2: The IIR referenced thermo-physical information the points in the schematics in the second scenario
State # T (K) P (MPa) h (kJ kg-1) s (kJ kg-1 K-1) Phase or x (–) Mach Number
1 313.0 10.50 306.1 1.332 Supercritical 0
aa 298.3 6.460 300.3 1.334 0.2095 0.58
2 313.0 6.800 436.2 1.781 Superheat 0
bb 309.2 6.460 434.5 1.782 Superheat 0.27
cc 298.3 6.460 351.3 1.505 0.6403 0.49
dd 298.3 6.460 351.3 1.505 0.6403 0.49
3 306.1 7.660 355.2 1.507 Supercritical 0
4 289.1 5.200 350.2 1.517 0.623 0
5 289.1 5.200 415.2 1.742 1 0
6 344.7 10.50 447.3 1.765 Supercritical 0
7 289.1 5.200 242.7 1.145 0 0
8 228.0 0.827 229.6 1.179 0.3833 0
9 228.0 0.827 434.1 2.075 1 0
10 408.5 6.80 571.0 2.162 Superheat 0
COP 1.249 𝐶𝑂𝑃 3.329
𝑚 17.2 gr/s.TR R 0.605

- 208 -
 EVERGREEN Joint Journal of Novel Carbon Resource Sciences & Green Asia Strategy, Vol. 08, Issue 01, pp204-212, March 2021

 
The computed COP illustrates that the current study
seems to have a remarkable performance at Table 1
temperatures, even though the assumptions are almost
pessimistic. The results look reasonable and they predict
no normal shockwaves inside the ejector in these
operating conditions. Because of this, future empirical
investigations seem to near results.
For the sensitivity analysis, there are supposition and
constraints. The pressure after Exp1 cannot be larger than
the ejector’s outlet pressure. Also, the inlet of the flash
chamber must be in the 2-phase state.
Two scenarios for the cycle’s operation were assumed.
This assumption is based on the behavior of the ejector Fig. 5: P-h diagram of the cycle according to Table 2 (2nd
outlet to the pressure alternation. The first scenario for the scenario)
time that the purged flow of the ejector is in the 2-phase
condition and the second for the time that the purged Installing compressor after evaporator is beneficial to
pressure becomes a single-phase fluid. Depending on the the performance. It will significantly increase the COP up
purge pressure, it can be superheated vapor or to an optimal point. However, the necessity of a
supercritical fluid. The target of this study is the turboexpander depends on the scenario. If the ejector
investigation of the effects of pressure on the COP. The drives the cycle to the 2nd scenario, the role of the first
effect of the evaporation and gas cooler temperature is not turboexpander is vital. It not only increases the COP to a
covered here. In Fig. 2, the point that scenario identified certain amount but also prepares the fluid for the
is shown. separation in the flash chamber. On the other hand, if the
The Flash chamber is the separator of the flow of the discharge is in 2 phase condition, the first expander
two loops. If in the sensitivity analysis the inlet flow to the becomes unnecessary due to the fact that the optimal
flash chamber turns to single-phase the computation is operating pressure happens at the end of the curve.
stopped. These boundaries are shown with the words “1st
Scenario” or “2nd Scenario” in Fig. 7.
Fig. 5 and Fig. 6 are the P-h and T-s diagrams of the
second scenario according to Table 2.
Fig. 7 shows the COP based on three pressure variable
inputs. The same color was chosen for the same purge
pressure lines. An important factor in COP is the
entrainment ratio which depends on directly the quality of
vapor after turboexpander. The optimum maximum
pressure of the cycle is between 9 to 11 MPa.
The shortest continuous line in Fig. 7 shows the
situation without Com2 and gas cooler. The endpoint of
this curve is the situation that there is no first
Fig. 6: T-s diagram of the cycle according to Table 2 (2nd
turboexpander in the layout. It can be concluded that this
scenario)
point shows a simple ejector cycle. COP in this point is
almost the same as previous studies which is a validation
for the calculation of this study. Comparing the COP of
the conventional ejector refrigeration cycle with the
maximum performance of the current layout, it can be
concluded that this novel layout can considerably increase
the performance at low-temperature refrigeration.

- 209 -
 Thermodynamic Feasibility Evaluation of a Novel Low-temperature Ejector-based Trans-critical R744 Refrigeration Cycle

 
Nomenclature
C: Speed of sound (m s-1)
C : Isobaric specific heat capacity (kJ kg-1 K-1)
C : Isochoric specific heat capacity (kJ kg-1 K-1)
COP: Coefficient of performance (–)
Com: Compressor
Exp: Expansion device or turboexpander
Eva: Evaporator
FC: Flash chamber
h: Specific enthalpy (kJ kg-1)
m: Mass flow rate (g s-1)
M: Mach number (–)
q: Specific heat energy (kW kg-1)
Q: Heat energy (kW)
R: Entrainment ratio
P: pressure (MPa)
Fig. 7: Pressure-based sensitivity analysis for Ppri = 10.5 MPa s: Specific entropy (kJ kg-1 K-1)
T: Temperature (K)
Adding the compressor inevitably increases the TR: Ton of refrigeration
pressure of the flash chamber which further boosts the V: Speed of flow (m s-1)
work recovery of Exp2 and the cooling capacity of the w: Specific Work (kW kg-1)
evaporator. As it was mentioned in the introduction, W: Work (kW)
simultaneous utilization of the ejector and the x: Quality of vapor (–)
turboexpander has not been reported. However, in the
current cycle, replacing the throttling valve between the Greek symbols
flash chamber and the evaporator with turboexpander α: Void fraction (–)
seems thermodynamically interesting enough to be γ: Specific heats ratio (–)
economic. ρ: Density (kg m-3)
Worthy of mention, the elevation of pressure in the
flash chamber leads to shrinkage of the compression ratio Subscripts
in the second compressor which in turn decreases its 1 ~ 10: State points number
energy consumption and mechanical problems. aa, bb, cc, dd: State points inside the ejector
These results are only for one working condition. In the d: diffuser
future, several plans must be executed for this research. id: Ideal
Full 1st and 2nd law analysis of cycle, thermo-economics is: Isentropic
and exergoeconomic analysis, improving the cycle in: inlet
performance by making the entrainment ratio more l: Liquid
controllable, combining this cycle with other heat m: Mixing
pumping and heat engine cycles in parallel and cascade n: nozzle
formation and studying the performance of cycle as a heat out: outlet
pump must be performed in the future. pri: primary flow into ejector
purge: purged flow after the ejector
4. Conclusion sec: secondary flow into ejector
A novel recommendation for the R744 deep freezing v: Vapor
application was presented. In the mathematical model, a
simple ejector cycle was modified by adding a compressor References
and gas cooler between the secondary inlet flow of the 1) G. Lorentzen, “Revival of carbon dioxide as a
ejector the evaporator, and one turboexpander after the refrigerant,” Int. J. Refrig., 17 (5) 292–301 (1994).
ejector’s outlet. The coefficient of the performance of the doi:10.1016/0140-7007(94)90059-0.
cycle is remarkable (around 1.3) at -45℃. 2) J. Lorentzen, Gustav, Pettersen, “A new , efficient
After the pressure-based sensitivity analysis of the and environmentally benign system for car un
cycle, it was shown that utilizing these components nouveau syst6me de conditionnement d ’ air
especially the compressor is beneficial. The automobile efficace et sans nuisance sur l ’
turboexpander is unnecessary if the discharge of the environnement,” Int. J.. Refrig., 16 (July) 4–12
ejector is in the 2-phase condition. (1993).
3) J. Ko, N. Takata, K. Thu, and T. Miyazaki, “Dynamic

- 210 -
 EVERGREEN Joint Journal of Novel Carbon Resource Sciences & Green Asia Strategy, Vol. 08, Issue 01, pp204-212, March 2021

 
modeling and validation of a carbon dioxide heat (2013). doi:10.1002/er.1879.
pump system,” Evergreen, 7 (2) 172–194 (2020). 15) D. Li, and E.A. Groll, “Transcritical co2 refrigeration
doi:10.5109/4055215. cycle with ejector-expansion device,” Int. J. Refrig.,
4) C. Aprea, and A. Maiorino, “An experimental 28 (5) 766–773 (2005).
evaluation of the transcritical co2 refrigerator doi:10.1016/j.ijrefrig.2004.10.008.
performances using an internal heat exchanger,” Int. 16) J. qiang Deng, P. xue Jiang, T. Lu, and W. Lu,
J. Refrig., 31 (6) 1006–1011 (2008). “Particular characteristics of transcritical co 2
doi:10.1016/j.ijrefrig.2007.12.016. refrigeration cycle with an ejector,” Appl. Therm.
5) S. Bhattacharyya, S. Mukhopadhyay, A. Kumar, R.K. Eng., 27 (2–3) 381–388 (2007).
Khurana, and J. Sarkar, “Optimization of a co2-c3h8 doi:10.1016/j.applthermaleng.2006.07.016.
cascade system for refrigeration and heating,” Int. J. 17) S. Elbel, “Historical and present developments of
Refrig., 28 (8) 1284–1292 (2005). ejector refrigeration systems with emphasis on
doi:10.1016/j.ijrefrig.2005.08.010. transcritical carbon dioxide air-conditioning,” Int. J.
6) M. Purjam, and K. Goudarzi, “High efficiency sub- Refrig., 34 (7) 1545–1561 (2011).
critical carbon dioxide supplementary heat pump for doi:10.1016/j.ijrefrig.2010.11.011.
low temperature climates (energy and exergy 18) K. Banasiak, A. Hafner, and T. Andresen,
analysis),” Renew. Energy, (2019). “Experimental and numerical investigation of the
doi:10.1016/j.renene.2018.10.024. influence of the two-phase ejector geometry on the
7) D.M. Robinson, and E.A. Groll, “Efficiencies of performance of the r744 heat pump,” Int. J. Refrig.,
transcritical co2 cycles with and without an 35 (6) 1617–1625 (2012).
expansion turbine: rendement de cycles transcritiques doi:10.1016/j.ijrefrig.2012.04.012.
au co2 avec et sans turbine d’expansion,” Int. J. 19) Y. Zhu, Y. Huang, C. Li, F. Zhang, and P.X. Jiang,
Refrig., 21 (7) 577–589 (1998). doi:10.1016/S0140- “Experimental investigation on the performance of
7007(98)00024-3. transcritical co2 ejector-expansion heat pump water
8) J.S. Baek, E.A. Groll, and P.B. Lawless, “Piston- heater system,” Energy Convers. Manag., 167
cylinder work producing expansion device in a (March) 147–155 (2018).
transcritical carbon dioxide cycle . part i : doi:10.1016/j.enconman.2018.04.081.
experimental investigation ´ tente et re ´ cupe ´ ration 20) J. Sarkar, “Optimization of ejector-expansion
de travail par piston dans un cycle de ´ rimentale transcritical co2 heat pump cycle,” Energy, 33 (9)
transcritique au dioxyde de carbone . partie i :,” Int. J. 1399–1406 (2008).
Refrig., 28 141–151 (2005). doi:10.1016/j.energy.2008.04.007.
doi:10.1016/j.ijrefrig.2004.08.006. 21) C. Lucas, and J. Koehler, “Experimental investigation
9) Y. Hwang, P. D, A. Celik, R. Radermacher, and C. of the cop improvement of a refrigeration cycle by use
Park, “Performance of co 2 cycles with a two-stage of an ejector,” Int. J. Refrig., 35 (6) 1595–1603
compressor,” Refrig. Air Cond., 1–8 (2004). (2012). doi:10.1016/j.ijrefrig.2012.05.010.
10) H. Fritschi, F. Tillenkamp, R. Löhrer, and M. Brügger, 22) Y. He, J. Deng, F. Yang, and Z. Zhang, “An optimal
“Efficiency increase in carbon dioxide refrigeration multivariable controller for transcritical co2
technology with parallel compression,” Int. J. Low- refrigeration cycle with an adjustable ejector,”
Carbon Technol., 12 (2) 171–180 (2017). Energy Convers. Manag., 142 466–476 (2017).
doi:10.1093/ijlct/ctw002. doi:10.1016/j.enconman.2017.03.070.
11) Z. Zhang, L. Tong, and X. Wang, “Thermodynamic 23) M.A. Islam, A. Pal, K. Thu, and B.B. Saha, “Study on
analysis of double-stage compression transcritical performance and environmental impact of
co<inf>2</inf> refrigeration cycles with an supermarket refrigeration system in japan,”
expander,” Entropy, 17 (4) 2544–2555 (2015). Evergreen, 6 (2) 168–176 (2019).
doi:10.3390/e17042544. doi:10.5109/2321014.
12) J. Sarkar, “Cycle parameter optimization of vortex 24) M.I. Alhamid, N. Nasruddin, Budihardjo, E. Susanto,
tube expansion transcritical co2 system,” Int. J. T.F. Vickary, and M. Arif Budiyanto, “Refrigeration
Therm. Sci., 48 (9) 1823–1828 (2009). cycle exergy-based analysis of hydrocarbon (r600a)
doi:10.1016/j.ijthermalsci.2009.01.016. refrigerant for optimization of household refrigerator,”
13) R. Llopis, R. Cabello, D. Sánchez, and E. Torrella, Evergreen, 6 (1) 71–77 (2019). doi:10.5109/2321015.
“Energy improvements of co<inf>2</inf> 25) S.A. Shaedi, N. Mohd-Ghazali, J.T. Oh, R. Ahmad,
transcritical refrigeration cycles using dedicated and Y. Mohd-Yunos, “Entropy generation
mechanical subcooling,” Int. J. Refrig., 55 129–141 minimization of two-phase flow in a mini channel
(2015). doi:10.1016/j.ijrefrig.2015.03.016. with genetic algorithm,” Evergreen, 6 (1) 39–43
14) J. Sarkar, “Performance optimization of transcritical (2019). doi:10.5109/2321004.
co 2 refrigeration cycle with thermoelectric 26) R. Span, and W. Wagner, “A new equation of state
subcooler,” Int. J. Energy Res., 37 (2) 121–128 for carbon dioxide covering the fluid region from the

- 211 -
 Thermodynamic Feasibility Evaluation of a Novel Low-temperature Ejector-based Trans-critical R744 Refrigeration Cycle

 
triple‐point temperature to 1100 k at pressures up to
800 mpa,” J. Phys. Chem. Ref. Data, 25 (6) 1509–
1596 (1996). doi:10.1063/1.555991.
27) D.L. Nguyen, E.R.F. Winter, and M. Greiner, “Sonic
velocity in two-phase systems,” Int. J. Multiph. Flow,
7 (3) 311–320 (1981). doi:10.1016/0301-
9322(81)90024-0.
28) T. Oka, T. Handa, F. Akagi, S. Yamaguchi, T. Aoki,
K. Yamabe, and Y. Kihara, “Steady-state analysis of
supersonic mixing enhanced by a three-dimensional
cavity flow,” Evergreen, 4 (1) 44–51 (2017).
doi:10.5109/1808452.

- 212 -