Next Article in Journal
Dynamic Banking Systemic Risk Accumulation under Multiple-Risk Exposures
Next Article in Special Issue
Investigation of Fluid Characteristic and Performance of an Ejector by a Wet Steam Model
Previous Article in Journal
Initial Solution Generation and Diversified Variable Picking in Local Search for (Weighted) Partial MaxSAT
Previous Article in Special Issue
Design and Investigation of a Dynamic Auto-Adjusting Ejector for the MED-TVC Desalination System Driven by Solar Energy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Back Pressure on Performances and Key Geometries of the Second Stage in a Highly Coupled Two-Stage Ejector

1
School of Civil Engineering and Architecture, Southwest University of Science and Technology, Mianyang 621010, China
2
School of Control Science and Engineering, Shandong University, Jinan 250061, China
3
Department of Mechanical, Aerospace and Civil Engineering, University of Manchester, Manchester M13 9PL, UK
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(12), 1847; https://doi.org/10.3390/e24121847
Submission received: 8 November 2022 / Revised: 12 December 2022 / Accepted: 16 December 2022 / Published: 18 December 2022
(This article belongs to the Special Issue Entropy and Exergy Analysis in Ejector-Based Systems)

Abstract

:
In this paper, for a highly coupled two-stage ejector-based cooling cycle, the optimization of primary nozzle length and angle of the second-stage ejector under varied primary nozzle diameters of the second stage was conducted first. Next, the evaluation for the influence of variable back pressure on ER of the two-stage ejector was performed. Last, the identification of the effect of the variable back pressure on the key geometries of the two-stage ejector was carried out. The results revealed that: (1) with the increase of the nozzle diameter at the second stage, the ER of both stages decreased with the increases of the length and angle of the converging section of the second-stage primary nozzle; (2) the pressure lift ratio range of the second-stage ejector in the critical mode gradually increased with the increase of the nozzle diameter of the second-stage; (3) when the pressure lift ratio increased from 102% to 106%, the peak ER of the second-stage decreased, and the influence of the area ratio and nozzle exit position of the second-stage ejector on its ER was reduced; (4) with the increase of nozzle diameter of the second-stage, the influence of area ratio and nozzle exit position of the second-stage on the second-stage performance decreased; and (5) the optimal AR of the second stage decreased but the optimal nozzle exit position of the second stage kept constant with the pressure lift ratio of the two-stage ejector.

1. Introduction

Nowadays, we are facing global warming and resource shortage [1]. Therefore, energy-saving and environmentally friendly refrigeration technologies have become a topic of widespread interest to refrigeration researchers and practitioners [2]. Ejector-based refrigeration systems have some advantages, such as no moving parts, waste heat driving and low operating costs [3,4]. Undoubtedly, ejector-based refrigeration systems are a promising industry [5,6].
In recent years, scholars have engaged in the research and development of two-stage ejector-based refrigeration systems. Kong et al. [7] investigated a supersonic two-stage ejector-diffuser system with numerical methods, and the system had four times better performance than a single ejector system. In addition, Kong et al. [8] predicted the flow phenomenon inside the two-stage ejector-diffuser. The results showed that the entrainment effects of the system greatly increased. Liu et al. [9] proposed a modified transcritical CO2 ejector enhanced two-stage compression cycle. The results showed that the heating coefficient of performance (COP) of the cycle, and so on, outperformed others. Wang et al. [10] presented a gas-fired air-to-water ejector heat pump. The system performance was improved with a high entrainment ratio. Liu et al. [11] proposed a novel two-stage compression transcritical CO2 refrigeration system with an ejector. The results indicated that the performance of the novel system was better than those of conventional systems. Yan et al. [12] presented a dual-ejector refrigeration system, and the area ratio (AR) had the most significant influences on the performance of the two-stage ejector. Ierin et al. [13] optimized a hybrid two-stage CO2 ejector-based cooling system, and the efficiency of the system increased by up to 32.7%. Ghorbani et al. [14] investigated a two-stage ejector cooling system, and the consumed power of the system decreased by 12.37%. Cao et al. [15] proposed a two-stage evaporation cycle and their numerical results disclosed that the COP of the cycle was improved, and the exergy was reduced. Sun et al. [16] claimed that the influence of phase transition in the ejector contributed to the ejector optimization. Chen et al. [17] studied the effect of the second-stage geometrical factors on the system performance and the optimized length to diameter ratio was 5. Similarly, Yadva et al. [18] numerically analyzed the performance of a two-stage ejector. Yang et al. [19] evaluated the exergy destruction characteristics inside a transcritical CO2 two-stage refrigeration system and showed that the system exergetic performance can be improved by enhancing the efficiency of the ejector. Yang et al. [20] also found that the gas cooler temperature had the greatest influence on ejector performance. Surendran et al. [21] explored a novel transcritical ejector regenerative refrigeration system and identified the system performance. Asfahan et al. [22] presented a system with two-stage ejectors and investigated them numerically. The results showed that the system had good performance when using the two-stage ejectors. Ding et al. [23] performed numerical studies using computational fluid dynamics (CFD) to predict two-stage ejector performance for subzero applications. Manjili et al. [24] used Engineering Equation Solver (EES) software to investigate a two-stage transcritical CO2 refrigeration cycle. It was found that the COP of new cycle was improved from 20% to 80% compared to the conventional cycle. Xue et al. [25] proposed and studied a two-stage vacuum ejector by comparing seven different ejector models. The results showed the two-stage ejector could provide superior suction pressure. Exposito-Carrillo et al. [26] optimized a two-stage CO2 refrigeration system, and the COP improved up to 13%. Wang et al. [27] developed a CFD model to investigate the performance of a proposed two-stage ejector. Viscito et al. [28] proposed a seasonal performance analysis of a hybrid ejector cooling system, they claimed that system required three or four ejectors for any reference climate, and they obtained an increase of the seasonal energy efficiency ratio up to 107%. Lillo et al. [29] presented a thermo-economic analysis of a waste heat recovery hybrid ejector cycle with a cooling load of 20 kW, the thermo-economic performance of this cycle has evident advantages over other waste heat driven system. Li et al. [30] carried out a numerical analysis of the influence of nozzle geometries on steam condensation and irreversibility in the ejector nozzle. The results indicated that the condensation of the steam makes a large amount of irreversible energy. Wen et al. [31] presented a two-stage ejector-based refrigeration system and optimized the two-stage ejector (TSE) geometries and system performance. Yan et al. [32] proposed another type of a highly coupled TSE-based system as shown in Figure 1, and they also optimized the key geometries such as area ratio (AR) and nozzle exit position (NXP) of the two-stage ejector, as illustrated in Figure 2, in which the mixture coming from the outlet of the first stage enters the primary nozzle of the second stage and entrains its secondary flow refrigerant.
However, no studies in the literature have mentioned the optimization of the key geometries, such as AR and NXP, of the highly coupled second-stage ejector under different back pressures and varied primary nozzles of the second stage (PNTD2). To bridge the gap, and based on our previous studies [33], further works in this study included:
  • CFD modelling and model validation of the highly coupled TSE;
  • Optimization of primary nozzle geometry of the second-stage ejector under varied primary nozzle of the second stage;
  • Evaluation the influence of variable back pressure on entrainment ratio (ER) of the TSE;
  • Identification of the effect of the variable back pressure on the key geometries of the second-stage ejector.

2. System Description and Numerical Method

2.1. System and Initial TSE Geometries

The schematic of the highly coupled TSE-based refrigeration system is shown in Figure 1, and the initial geometrical parameters of the TSE are presented in Figure 3.

2.2. CFD Modelling

The flow inside the TSE is calculated by using governing equations [33,34], and Gambit 2.4 and Ansys 19.0 [35] are used in this simulation. Grids with 103,000 quadrilateral elements are created as shown in Figure 4.
R134a is the working fluid with parameters from NIST [36], the RNG k-ε turbulence model was selected in this study. The standard wall function was chosen, and the range of the first grid cell is in the region of 30 < y+ < 300. The residual convergence limit for each equation is below 10−5, except that for energy equation is set to less than 10−6. In addition, to ensure that the refrigerant liquid completely evaporates into refrigerant gas, three inlet streams are set as 10 K superheat. The primary fluid inlet and the secondary fluid inlet are both set as the pressure inlet, while the outlet is set as the pressure outlet [32], and the boundary conditions of the TSE are illustrated in Table 1.
Three levels of grids (71,000, 103,000 and 138,000) are used to validate the grid independence as illustrated in Figure 5. Since the three grid levels are quite close with each other, the medium one is finally used in the following simulation.

3. Model Validation

3.1. Experimental Setup

The experimental setup is presented in Figure 6, in the setup, Evaporator 1 is simulated as an air conditioner, and its evaporating temperature is set as 7 °C; Evaporator 2 is simulated as a refrigerator, and its evaporating temperature is arranged as −5 °C; whist, Evaporator 3 is simulated as a freezer, and its evaporating temperature is specified as −30 °C. The ambient temperature is valued at 36 °C. Based on the thermodynamic calculation, the individual required cooling loads for three evaporators are 1566.2 W, 609.4 W and 997.2 W, respectively. Other details can refer to our previous study [32]. The range and accuracy of sensors are presented in Table 2.

3.2. Validation of the CFD Model

Fifteen CFD simulation results, as illustrated in Table 3, were validated by the experimental data. The average and maximum discrepancy of ER1 were 7.2% and 11.9%, and those for ER2 are 5.9% and 10.6%, respectively; therefore, the models can be used in the following simulations.

4. Results and Discussion

4.1. Optimization of Nozzle Geometry of the Second-Stage Ejector

4.1.1. Optimization of Nozzle Length and Angle of the Second-Stage Ejector When PNTD2 Is 4.1 mm

Figure 7 shows the results of ER1 and ER2 with the second-stage ejector nozzle length (LC2) when PNTD2 is 4.1 mm. When LC2 changes from 5 mm to 40 mm, ER1 first rises from 0.595 to the maximum value of 0.646 (LC2 = 25 mm), and then slowly decreases to 0.636 (LC2 = 40 mm). ER2 rises slowly from 2.096 to 2.153 (LC2 = 20 mm) and then decreases rapidly until it reaches a minimum of 2.022. When LC2 changes from 5 mm to 40 mm, the maximum deviations of ER1 and ER2 reach 8.571% and 6.479%, respectively, which means that the change of LC2 has an impact on the performance of both the first stage and the second-stage, but the impact on the first-stage is more obvious. At the same time, it can be seen that when LC2 = 15–25 mm, the values of ER1 and ER2 are relatively large.
Figure 8 displays the results of ER1 and ER2 with the second-stage ejector nozzle angle (AC2) when PNTD2 is 4.1 mm. With the increase of AC2, ER1 first increases and then decreases, and the maximum value of ER1 is 0.653 (AC2 = 16°). As AC2 increases from 6° to 10°, ER2 increases from 2.146 to 2.171, and ER2 gradually decreases to 2.014 as AC2 continues to rise to 22°. Compared with the initial values of ER1 and ER2 (0.632 and 2.153), the maximum values of ER1 and ER2 increase by 0.021 and 0.019, respectively, and the maximum deviations of ER1 and ER2 are 11.054% and 7.795%, respectively. This means that AC2 has an impact on the performance of both stages, but ER1 is more sensitive to the changes of AC2. Moreover, it can be seen that AC2 has a greater impact on the two-stage performance than LC2. In addition, the AC2 range for which ER1 achieves large values is 10°–22°, while the AC2 range for which ER2 achieves large values is 6°–14°, indicating that the AC2 range for which both ER1 and ER2 obtain large values is 10°–14°.

4.1.2. Optimization of Nozzle Length and Angle of Second-Stage Ejector When PNTD2 Is 4.7 mm

Figure 9 reveals the results of ER1 and ER2 with LC2 when PNTD2 is 4.7 mm. It can be seen that when LC2 changes from 5 mm to 45 mm, ER1 first rises from 1.031 to the maximum value of 1.098 (LC2 = 35 mm), and then rapidly decreases to 1.074 (LC2 = 45 mm). However, ER2 slowly increases to the maximum value of 1.861 (LC2 = 15 mm), follows by a rapid decline in ER2 until it decreases to the lowest value of 1.799. When LC2 changes from 5 mm to 45 mm, the maximum deviations of ER1 and ER2 reach 6.499% and 3.446%, respectively, which reflects that the performance of the first-stage and the second-stage are affected by the change of LC2, but the first-stage can be affected more obviously. Furthermore, it can be seen that when LC2 = 15–25 mm, the values of ER1 and ER2 are relatively large.
Figure 10 indicates the results of ER1 and ER2 with AC2 when PNTD2 is 4.7 mm. As shown in the figure, with the increase of AC2, ER1 increases first and then decreases, and the highest value of ER1 is 1.091 (AC2 = 14°). With the increase of AC2 from 6° to 8°, ER2 increases from 1.849 to 1.856; and with the increase of AC2 to 22°, ER2 decreases to 1.811. The maximum deviations of ER1 and ER2 are 5.106% and 2.485%, respectively, which means that AC2 has an impact on the performance of both stages, but ER1 is more sensitive to the change of AC2. In addition, the range of AC2 where both ER1 and ER2 are at large values is still within 10°–14°.

4.1.3. Optimization of Nozzle Length and Angle of Second-Stage Ejector When PNTD2 Is 5.3 mm

Figure 11 presents the results of ER1 and ER2 with LC2 when PNTD2 is 5.3 mm. As illustrated in the figure, when LC2 changes from 5 mm to 40 mm, ER1 increases from 1.502 to the maximum value 1.585. ER2 decreases from 1.645 to a minimum of 1.598. Compared with the initial values of ER1 and ER2 (1.573 and 1.627), the maximum values of ER1 and ER2 increase by 0.012 and 0.018. In addition, when LC2 changes from 5 mm to 40 mm, the maximum deviations of ER1 and ER2 are 5.526% and 2.941%, respectively, which means that the performance of both stages is affected by the change of LC2, and the performance of the first-stage is slightly more affected. In contrast, when LC2 is 15–25 mm, the values of ER1 and ER2 are at large values.
Figure 12 shows the changes of ER1 and ER2 with AC2 when PNTD2 is 5.3 mm. It can be seen that with the increase of AC2, ER1 first increases and then decreases, and the maximum value of ER1 is 1.573 (AC2 = 14°). As AC2 increases from 6° to 10°, ER2 increases from 1.629 to 1.632, when AC2 continues to increase to 18°, ER2 gradually decreases to 1.621. The maximum deviations of ER1 and ER2 are 3.897% and 1.527%, respectively, which indicates that AC2 has an impact on the performance of both stages, but ER1 is slightly more sensitive to AC2. In addition, the range of AC2 is still within 10°–14° when both ER1 and ER2 are at large values.
In conclusion, both LC2 and AC2 have certain effects on the performance of the TSE. With the increase of PNTD2, the influence of LC2 and AC2 on the performance of the two stages is gradually weakened. With a comprehensive consideration, LC2 = 25 mm and AC2 = 12° are selected as the optimized geometries of the second-stage ejector converging nozzle to carry out the following study on the influence of variable back pressure on the ER of TSE.

4.2. Influence of Variable Back Pressure on ERs of the TSE

Boundary conditions of the TSE except the back pressure are kept constant, and the change of the back pressure of the TSE is expressed as the percentage of pressure lift, namely the change of PLR (the ratio of the back pressure to the secondary inlet pressure of the second-stage ejector). The initial PLR is 108%. When PNTD2 is 4.1 mm and PLR changes in the range of 102–118%, the influence of the changed PLR on ER1 and ER2 is shown in Figure 13. It can be seen that when PLR increases from 102% to 118%, ER1 gradually decreases, but its maximum value and minimum value are 0.644 and 0.643, respectively. Therefore, ER1 almost does not change. ER2 decreases almost linearly from 2.442 to 1.105, and the maximum deviation of ER2 is 121.0%. Therefore, changes in PLR have a significant impact on ER2.
When PNTD2 = 4.7 mm, the influence of changing PLR on ER1 and ER2 is shown in Figure 14. ER1 is almost unaffected by PLR, with the highest and lowest values of 1.092 and 1.090, respectively. In addition, when PLR increases from 102% to 104%, ER2 remains at 1.938; when PLR increases from 104% to 118%, ER2 decreases to the minimum value of 1.253, and the maximum deviation of ER2 is 54.7%. Therefore, the change of PLR has a relatively obvious impact on ER2. Moreover, it can be seen that compared with the ER value of PNTD2 = 4.1 mm, ER1 increases a lot, while ER2 decreases a little.
The effect of the changing PLR on ER1 and ER2 at PNTD2 = 5.3 mm is shown in Figure 15. With the increase of PLR, ER1 is still almost unaffected. However, compared with PNTD2 = 4.7 mm, ER1 increases largely. The maximum value of ER2 is lower than that of PNTD2 = 4.7 mm. When PLR increases from 102% to 106%, ER2 remains at 1.654 and it still shows a downward trend while PLR increases from 106% to 118%, and its maximum and minimum values are 1.654 and 1.213, respectively. The maximum deviation of ER2 is 36.4%. It is noted that at PNTD2 = 5.3 mm, the PLR affects ER2 to a smaller extent than at PNTD2 = 4.7 mm.
In summary, the change of PLR basically has no effect on the performance of the first stage and a significant effect on the performance of the second stage. Furthermore, the effect of PLR on the second-stage performance gradually diminishes with increasing PNTD2. It can also be seen that when PNTD2 = 4.1 mm, the ejector with a PLR of 102% is already in subcritical mode. When PNTD2 = 4.7 mm and PLR is 102–104%, the second-stage ejector is in critical mode, and when PLR is greater than 104%, it is in subcritical mode. When PNTD2 = 5.3 mm, the PLR is in the range of 102–106%, and the ejector is in critical mode. Therefore, the critical back pressure increases with the increase of PNTD2. The reason for this phenomenon is that, when the PNTD2 increases, which means the area ratio of the secondary stage reduces, normally the entrainment ratio increases with the area ratio when the back pressure keeps unchanged; as a result, the critical pressure increases with the increase of PNTD2. For different PNTD2, the next study will be carried out for the optimization of the key geometries of the second-stage ejector under different PLR, such as AR2 and NXP2, to identify the influence of back pressure and PNTD2 on the best AR2 and NXP2.

4.3. Effect of the Variable Back Pressure on the Key Geometries of the TSE

4.3.1. Optimized AR2

The effect of AR2 at PNTD2 = 4.1 mm on ER2 at different PLR is shown in Figure 16. It can be seen that the change trend of ER2 affected by AR2 under the three PLR is consistent, that is, ER2 first increases, and then decreases as AR2 increases. At PLR = 102%, ER2 reaches its maximum value of 3.354 at AR2 = 23.5, so the best value for AR2 is 23.5, which is 13.0 times more than the optimal value of AR2 at PLR = 108%, and ER2 increases by 37.3% compared to the optimal value of 2.442 at PLR = 108%. When PLR is 104%, ER2 reaches a maximum value of 2.682 at AR2 = 17.5, then the optimal value of AR2 is 17.5, which is 7.0 times greater than the optimal value of AR2 at PLR = 108%; the optimal value of ER2 increases by 13.1% compared to the optimal value of 2.372 for ER2 at PLR = 108%. At PLR = 106%, the optimal value of 2.480 for ER2 is obtained at AR2 = 14.5, so that the optimal value for AR2 is 14.5, which is 4.0 times greater than the optimal value for AR2 at PLR = 108%, and the optimal ER2 increases by 9.4% over the optimal value for ER2 (2.266) at PLR = 108%. In summary, with the PLR increasing from 102% to 108%, the maximum ER2 decreases from 3.354 to 2.266, and the corresponding optimal AR2 decreases from 23.5 to 10.5.
Figure 17 shows the effect of AR2 at PNTD2 = 4.7 mm on ER2 under the three PLR. When PLR is 102%, ER2 increases first and then decreases with the change of AR2, and its maximum value of 3.022 is obtained at AR2 = 20.1, which is 55.9% higher than the optimal value of ER2 (1.938) at PLR of 108%; and the optimal AR2 is 12.0 times larger than that at PLR of 108%. With a PLR of 104%, when AR2 changes, the optimal value for AR2 is 16.1, and its corresponding maximum ER2 of 2.598, and ER2 increases by 0.662 over the optimal value of ER2 (1.936) at a PLR of 108%. When the PLR is 106%, and when AR2 changes, ER2 rises first and then decreases; its optimal value of 2.283 is obtained at AR2 = 13.1; and the optimal value of ER2 is increased by 0.368 compared to the optimal value of ER2 of 1.915 (PLR of 108%). Similar to PNTD2 = 4.1 mm, the maximum value of ER2 decreases and the corresponding optimal AR2 decreases as the PLR increases from 102% to 108%.
Figure 18 shows the change of ER2 at PNTD2 = 5.3 mm with AR2 under the different PLR conditions. At PLR = 102%, as AR2 changes from 2.9 to 19.9, ER2 increases from 0.336 to a maximum of 2.663 at AR2 = 16.9, thus, the best value for AR2 is 16.9. Furthermore, the optimal value of ER2 is increased by 1.008 when compared to the optimal ER2 (1.655) at PLR = 108%. At PLR = 104%, ER2 achieves a maximum value of 2.299 at AR2 = 12.9. Therefore, the optimal AR2 is 12.9. In addition, the optimal value of ER2 is increased by 0.645 when compared to the best ER2 (1.654) at PLR = 108%. When PLR = 106%, ER2 shows a trend of increasing first and then decreasing with the increase of AR2. Its maximum value 2.050 is obtained at AR2 = 10.9, hence, the optimum of AR2 is 10.9. Compared with the optimal ER2 (1.653) when PLR is 108%, the maximum ER2 increases by 0.397. It can be seen that when PNTD2 is 5.3 mm and PLR changes from 102% to 108%, the maximum value of ER2 and the corresponding optimal AR2 also show a decreasing trend. This phenomenon can be probably explained as follows: when the PLR increases, the pressure difference between the outlet and the secondary flow inlet pressure increases, in order to maintain the increased pressure lift; thus, the ER2 usually drops and the suitable AR2 reduces accordingly.

4.3.2. Optimized NXP2

Figure 19 shows the influence of NXP2 on ER2 under different PLR when PNTD2 is 4.1 mm. With given PLR = 102%, as NXP2 increases from 18 mm to 30 mm, the change trend of ER2 rises first and then decreases. Its maximum value of 2.486 is obtained when NXP2 = 26 mm, and the maximum deviation of ER2 is 1.944%. With given PLR = 104%, as NXP2 increases from 18 mm to 26 mm, ER2 increases from 2.341 to 2.382 and then decreases to 2.369, so the maximum deviation of ER2 is 1.756%. With given PLR = 106%, ER2 also shows a trend of first increasing and then decreasing, and its peak value of 2.269 is obtained at NXP2 = 26 mm, and the maximum deviation of ER2 is 1.584%. Therefore, it can be seen that when PNTD2 = 4.1 mm, the optimal value of NXP2 does not change with the change of PLR, and thus the influence of NXP2 on ER2 is far less than that of AR2.
Figure 20 demonstrates the influence of NXP2 on ER2 under different PLR when PNTD2 is 4.7 mm. When PLR is 102%, as NXP2 changes from 18 mm to 24 mm, ER2 rises from 1.932 to 1.950, and then decreases to 1.929 when NXP2 is 30 mm, and the maximum deviation of ER2 is 1.057%. When PLR is 104%, the change trend of ER2 also increases first and then decreases, and its peak value of 1.936 appears at the value of 24 mm of NXP2, while the maximum deviation of ER2 is 0.868%. When PLR is 106%, the variation trend of ER2 is similar to the previous two, and its maximum value of 1.921 is obtained at NXP2 = 24 mm, and the maximum deviation of ER2 is 0.744%. Therefore, when PNTD2 is 4.7 mm, the optimal value of NXP2 is 24 mm, which does not change with the change of PLR.
Figure 21 shows the influence of NXP2 on ER2 under different PLR when PNTD2 is 5.3 mm. The changing trend of ER2 is the same under the three PLR conditions, that is, ER2 increases first and then decreases, and all the maximum value of ER2 appear at NXP2 = 22 mm. Therefore, when PNTD2 is 5.3 mm, the best value of NXP2 is 22 mm, that is, it is not affected by the change of PLR.

5. Conclusions

In this paper, the CFD simulation method was first used to optimize the nozzle geometry of the second-stage ejector under different PNTD2 with the operating conditions given in Table 1. Then, the effect of variable back pressure on the ejector performance was studied. Finally, three PLRs that place the ejector in critical or near critical mode were selected to study the influence of AR2 and NXP2 on the performance of the second-stage ejector. The main findings obtained are as follows:
(1)
When LC2 and AC2 change, the maximum values of ER1 and ER2 do not appear at the same length or angle. LC2 = 25 mm and AC2 = 12° are the relative optimal combination of values for the second-stage ejector nozzle;
(2)
With the increase of PNTD2, ER1 and ER2 decrease with the increase of LC2 and AC2, and the PLR range of the ejector in the critical mode gradually increases;
(3)
The change of PLR has no effect on the performance of the first-stage ejector, but has a significant effect on the performance of the second-stage ejector; with the increase of PNTD2, the influence of PLR on the performance of the second stage is gradually weakened;
(4)
When PNTD2 is 4.1 mm, the optimal value of AR2 decreases from 23.5 to 14.5 with the increase of PLR, and the peak value of ER2 decreases from 3.354 to 2.480. When PNTD2 is 4.7 mm, the optimal value of AR2 decreases from 20.1 to 13.1 with the increase of PLR, and the maximum value of ER2 decreases from 3.022 to 2.383. When PNTD2 is 5.3 mm, with the change of PLR from 102% to 106%, the optimal value of AR2 is from 16.9 to 10.9, and the peak value of ER2 is reduced from 3.354 to 2.382;
(5)
The optimal value of NXP2 is not affected by the change of PLR. When PNTD2 is 4.1 mm, 4.7 mm and 5.3 mm, the corresponding optimal value of NXP2 is 26 mm, 24 mm and 22 mm, respectively.

Author Contributions

Conceptualization, J.Y.; methodology, C.W.; software, J.Y. and Y.S.; validation, C.W.; and investigation, J.Y. and Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

The work was funded by the National Natural Science Foundation of China under Grant Number: 51806235.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Symbols
PNTD 1 primary nozzle throat diameter of first-stage ejector, mm
PNTD 2 primary nozzle throat diameter of second-stage ejector, mm
LC 1 convergent section length of primary nozzle, mm
LC 2 divergent section length of primary nozzle, mm
AC 1 convergent section angle of primary nozzle, mm
AC 2 divergent section angle of primary nozzle, mm
m ˙ mass flow rate, g s−1
Ppressure, kPa
Ttemperature, K or °C
Abbreviations
ARarea ratio
CFDComputational Fluid Dynamics
COPcoefficient of performance
EMERCejector-based multi-evaporator refrigeration cycle
ERentrainment ratio
ERCejector-based refrigeration cycle
NXPnozzle exit position
PRRpressure recovery ratio
PRVpressure regulating valve
TSEtwo-stage ejector
Subscripts
pprimary flow
ssecondary flow
ooutflow, optimum
1first-stage ejector
2second-stage ejector

References

  1. Hong, B.; Li, X.; Li, Y.; Chen, S.; Tan, Y.; Fan, D.; Song, S.; Zhu, B.; Gong, J. An improved hydraulic model of gathering pipeline network integrating pressure-exchange ejector. Energy 2022, 260, 125101. [Google Scholar] [CrossRef]
  2. Zhao, Y.; Yu, J. Thermodynamic analysis of a modified vapor-injection heat pump cycle using an ejector. Int. J. Refrig. 2023, 145, 137–147. [Google Scholar] [CrossRef]
  3. Zhang, Y.; He, C.; Sun, L. Optimization of an ejector to mitigate cavitation phenomena with coupled CFD/BP neural network and particle swarm optimization algorithm. Prog. Nucl. Energy 2022, 153, 104412. [Google Scholar] [CrossRef]
  4. Wilhelmsen, O.; Aasen, A.; Banasiak, K.; Herlyng, H.; Hafner, A. One-dimensional mathematical modeling of two-phase ejectors: Extension to mixtures and mapping of the local exergy destruction. Appl. Therm. Eng. 2022, 217, 119228. [Google Scholar] [CrossRef]
  5. Kulkarni, A.; Yadav, S.; Kumar, A. A comprehensive study on solar ejector cooling system: A review. Mater. Today Proc. 2022, 69, 463–467. [Google Scholar] [CrossRef]
  6. Kumar, V.; Yadav, S.; Ram, S. A comprehensive studies on constant area mixing (CAM) and constant pressure mixing (CPM) Ejectors: A review. Mater. Today Proc. 2022, 69, 513–518. [Google Scholar] [CrossRef]
  7. Kong, F.; Kim, H. Optimization study of a two-stage ejector–diffuser system. Int. J. Heat Mass Transf. 2016, 101, 1151–1162. [Google Scholar] [CrossRef]
  8. Kong, F.; Kim, H. Analytical and computational studies on the performance of a two-stage ejector–diffuser system. Int. J. Heat Mass Transf. 2015, 85, 71–87. [Google Scholar] [CrossRef]
  9. Liu, Y.; Liu, J.; Yu, J. Theoretical analysis on a novel two-stage compression transcritical CO2 dual-evaporator refrigeration cycle with an ejector. Int. J. Refrig. 2020, 119, 268–275. [Google Scholar] [CrossRef]
  10. Wang, P.; Ma, H.; Spitzenberger, J.; Abu-Heiba, A.; Nawaz, K. Thermal performance of an absorption-assisted two-stage ejector air-to-water heat pump. Energy Convers. Manag. 2021, 230, 113761. [Google Scholar] [CrossRef]
  11. Liu, J.; Zhou, L.; Lyu, N.; Lin, Z.; Zhang, S.; Zhang, X. Analysis of a modified transcritical CO2 two-stage ejector-compression cycle for domestic hot water production. Energy Convers. Manag. 2022, 269, 116094. [Google Scholar] [CrossRef]
  12. Yan, J.; Wen, N.; Wang, L.; Li, X.; Liu, Z.; Li, S. Optimization on ejector key geometries of a two-stage ejector-based mul-ti-evaporator refrigeration system. Energ. Convers. Manag. 2018, 175, 142–150. [Google Scholar] [CrossRef]
  13. Ierin, V.; Chen, G.; Volovyk, O.; Shestopalov, K. Hybrid two–stage CO2 transcritical mechanical compression–ejector cooling cycle: Thermodynamic analysis and optimization. Int. J. Refrig. 2021, 132, 45–55. [Google Scholar] [CrossRef]
  14. Ghorbani, B.; Ebrahimi, A.; Moradi, M.; Ziabasharhagh, M. Continuous production of cryogenic energy at low-temperature using two-stage ejector cooling system, Kalina power cycle, cold energy storage unit, and photovoltaic system. Energy Convers. Manag. 2021, 227, 113541. [Google Scholar] [CrossRef]
  15. Cao, X.; Liang, X.; Shao, L.; Zhang, C. Performance analysis of an ejector-assisted two-stage evaporation single-stage va-por-compression cycle. Appl. Therm. Eng. 2022, 205, 118005. [Google Scholar] [CrossRef]
  16. Sun, W.; Liu, C.; Zhang, H.; Sun, W.; Xue, H.; Jia, L. Numerical analysis of two-stage vacuum ejector performance considering the influence of phase transition and non-condensable gases. Appl. Therm. Eng. 2022, 213, 118730. [Google Scholar] [CrossRef]
  17. Chen, W.; Xue, K.; Wang, Y.; Chong, D.; Yan, J. Numerical assessment on the performance of two-stage ejector to boost the low-pressure natural gas. J. Nat. Gas Sci. Eng. 2016, 34, 575–584. [Google Scholar] [CrossRef]
  18. Yadav, S.; Pandey, K.; Kumar, V.; Gupta, R. Computational analysis of a supersonic two-stage ejector. Mater. Today Proc. 2021, 38, 2700–2705. [Google Scholar] [CrossRef]
  19. Yang, D.; Li, Y.; Xie, J.; Wang, J. Exergy destruction characteristics of a transcritical carbon dioxide two-stage compres-sion/ejector refrigeration system for low-temperature cold storage. Energy Rep. 2022, 8, 8546–8562. [Google Scholar] [CrossRef]
  20. Yang, D.; Jie, Z.; Zhang, Q.; Li, Y.; Xie, J. Evaluation of the ejector two-stage compression refrigeration cycle with work perfor-mance from energy, conventional exergy and advanced exergy perspectives. Energy Rep. 2022, 8, 12944–12957. [Google Scholar] [CrossRef]
  21. Surendran, A.; Seshadri, S. An ejector based Transcritical Regenerative Series Two-Stage Organic Rankine Cycle for du-al/multi-source heat recovery applications. Therm. Sci. Eng. Prog. 2022, 27, 101158. [Google Scholar] [CrossRef]
  22. Ali, E.S.; Asfahan, H.M.; Sultan, M.; Askalany, A.A. A novel ejectors integration with two-stages adsorption desalination: Away to scavenge the ambient energy. Sustain. Energy Technol. Assess. 2021, 48, 101658. [Google Scholar] [CrossRef]
  23. Ding, Z.; Wang, L.; Zhao, H.; Zhang, H.; Wang, C. Numerical study and design of a two-stage ejector for subzero refrigeration. Appl. Therm. Eng. 2016, 108, 436–448. [Google Scholar] [CrossRef]
  24. Eskandari Manjili, F.; Cheraghi, M. Performance of a new two-stage transcritical CO2 refrigeration cycle with two ejectors. Appl. Therm. Eng. 2019, 156, 402–409. [Google Scholar] [CrossRef]
  25. Xue, H.; Wang, L.; Jia, L.; Xie, C.; Lv, Q. Design and investigation of a two-stage vacuum ejector for MED-TVC system. Appl. Therm. Eng. 2020, 167, 114713. [Google Scholar] [CrossRef]
  26. Expósito-Carrillo, J.A.; La Flor, F.J.S.-D.; Perís-Pérez, B.; Salmerón-Lissén, J.M. Thermodynamic analysis of the optimal operating conditions for a two-stage CO2 refrigeration unit in warm climates with and without ejector. Appl. Therm. Eng. 2021, 185, 116284. [Google Scholar] [CrossRef]
  27. Wang, X.; Wang, L.; Song, Y.; Deng, J.; Zhan, Y. Optimal design of two-stage ejector for subzero refrigeration system on fishing vessel. Appl. Therm. Eng. 2021, 187, 116565. [Google Scholar] [CrossRef]
  28. Viscito, L.; Lillo, G.; Napoli, G.; Mauro, A.W. Waste Heat Driven Multi-Ejector Cooling Systems: Optimization of Design at Partial Load; Seasonal Performance and Cost Evaluation. Energies 2021, 14, 5663. [Google Scholar] [CrossRef]
  29. Lillo, G.; Mastrullo, R.; Mauro, A.W.; Trinchieri, R.; Viscito, L. Thermo-Economic Analysis of a Hybrid Ejector Refrigerating System Based on a Low Grade Heat Source. Energies 2020, 13, 562. [Google Scholar] [CrossRef] [Green Version]
  30. Li, H.; Wang, X.; Huang, H.; Ning, J.; Tu, J. A Numerical Analysis of the Influence of Nozzle Geometric Structure on Spontaneous Steam Condensation and Irreversibility in the Steam Ejector Nozzle. Appl. Sci. 2021, 11, 11954. [Google Scholar] [CrossRef]
  31. Wen, N.; Wang, L.; Yan, J.; Li, X.; Liu, Z.; Li, S.; Zou, G. Effects of operating conditions and cooling loads on two-stage ejector perfor-mances. Appl Therm. Eng. 2019, 150, 770–780. [Google Scholar] [CrossRef]
  32. Yan, J.; Cai, Q.; Wen, H. Optimization on key geometries of a highly coupled two-stage ejector. Appl. Therm. Eng. 2021, 197, 117362. [Google Scholar] [CrossRef]
  33. Allouche, Y.; Bouden, C.; Varga, S. A CFD analysis of the flow structure inside a steam ejector to identify the suitable experi-mental operating conditions for a solardriven refrigeration system. Int. J. Refrig. 2014, 39, 186–195. [Google Scholar] [CrossRef]
  34. He, S.; Li, Y.; Wang, R. Progress of mathematical modeling on ejectors. Renew. Sustain. Energy Rev. 2009, 13, 1760–1780. [Google Scholar] [CrossRef]
  35. Wen, H.; Yan, J.; Li, X. Influence of liquid volume fraction on ejector performance: A numerical study. Appl. Therm. Eng. 2021, 190, 116845. [Google Scholar] [CrossRef]
  36. NIST Chemistry WebBook. NIST Standard Reference Database Number 69. Available online: http://webbook.nist.gov/chemistry (accessed on 1 February 2010).
Figure 1. The schematic of the highly coupled TSE-based refrigeration system: (a) schematic; and (b) pressure-enthalpy diagram.
Figure 1. The schematic of the highly coupled TSE-based refrigeration system: (a) schematic; and (b) pressure-enthalpy diagram.
Entropy 24 01847 g001
Figure 2. The schematic of the highly coupled TSE.
Figure 2. The schematic of the highly coupled TSE.
Entropy 24 01847 g002
Figure 3. Initial geometrical parameters of the TSE.
Figure 3. Initial geometrical parameters of the TSE.
Entropy 24 01847 g003
Figure 4. 2D axisymmetric quadrilateral grids of the TSE.
Figure 4. 2D axisymmetric quadrilateral grids of the TSE.
Entropy 24 01847 g004
Figure 5. Axial static pressure with grid levels of the TSE.
Figure 5. Axial static pressure with grid levels of the TSE.
Entropy 24 01847 g005
Figure 6. Photograph of experimental setup.
Figure 6. Photograph of experimental setup.
Entropy 24 01847 g006
Figure 7. Changes of ER1 and ER2 with LC2 when PNTD2 is 4.1 mm.
Figure 7. Changes of ER1 and ER2 with LC2 when PNTD2 is 4.1 mm.
Entropy 24 01847 g007
Figure 8. Changes of ER1 and ER2 with AC2 when PNTD2 is 4.1 mm.
Figure 8. Changes of ER1 and ER2 with AC2 when PNTD2 is 4.1 mm.
Entropy 24 01847 g008
Figure 9. Changes of ER1 and ER2 with LC2 when PNTD2 is 4.7 mm.
Figure 9. Changes of ER1 and ER2 with LC2 when PNTD2 is 4.7 mm.
Entropy 24 01847 g009
Figure 10. Changes of ER1 and ER2 with AC2 when PNTD2 is 4.7 mm.
Figure 10. Changes of ER1 and ER2 with AC2 when PNTD2 is 4.7 mm.
Entropy 24 01847 g010
Figure 11. Changes of ER1 and ER2 with LC2 when PNTD2 is 5.3 mm.
Figure 11. Changes of ER1 and ER2 with LC2 when PNTD2 is 5.3 mm.
Entropy 24 01847 g011
Figure 12. Changes of ER1 and ER2 with AC2 when PNTD2 is 5.3 mm.
Figure 12. Changes of ER1 and ER2 with AC2 when PNTD2 is 5.3 mm.
Entropy 24 01847 g012
Figure 13. Effect of varied PLR on ER1 and ER2 (PNTD2 = 4.1 mm).
Figure 13. Effect of varied PLR on ER1 and ER2 (PNTD2 = 4.1 mm).
Entropy 24 01847 g013
Figure 14. Effect of varied PLR on ER1 and ER2 (PNTD2 = 4.7 mm).
Figure 14. Effect of varied PLR on ER1 and ER2 (PNTD2 = 4.7 mm).
Entropy 24 01847 g014
Figure 15. Effect of varied PLR on ER1 and ER2 (PNTD2 = 5.3 mm).
Figure 15. Effect of varied PLR on ER1 and ER2 (PNTD2 = 5.3 mm).
Entropy 24 01847 g015
Figure 16. The variation of ER2 with AR2 under different PLR (PNTD2 = 4.1 mm).
Figure 16. The variation of ER2 with AR2 under different PLR (PNTD2 = 4.1 mm).
Entropy 24 01847 g016
Figure 17. The variation of ER2 with AR2 under different PLR (PNTD2 = 4.7 mm).
Figure 17. The variation of ER2 with AR2 under different PLR (PNTD2 = 4.7 mm).
Entropy 24 01847 g017
Figure 18. The variation of ER2 with AR2 under different PLR (PNTD2 = 5.3 mm).
Figure 18. The variation of ER2 with AR2 under different PLR (PNTD2 = 5.3 mm).
Entropy 24 01847 g018
Figure 19. The variation of ER2 with NXP2 under different PLR (PNTD2 = 4.1 mm).
Figure 19. The variation of ER2 with NXP2 under different PLR (PNTD2 = 4.1 mm).
Entropy 24 01847 g019
Figure 20. The variation of ER2 with NXP2 under different PLR (PNTD2 = 4.7 mm).
Figure 20. The variation of ER2 with NXP2 under different PLR (PNTD2 = 4.7 mm).
Entropy 24 01847 g020
Figure 21. The variation of ER2 with NXP2 under different PLR (PNTD2 = 5.3 mm).
Figure 21. The variation of ER2 with NXP2 under different PLR (PNTD2 = 5.3 mm).
Entropy 24 01847 g021
Table 1. Boundary conditions of the TSE.
Table 1. Boundary conditions of the TSE.
Boundary ConditionsP (kPa)T (K)
First-stage primary flow374.6290
First-stage secondary flow243.4278
Second-stage secondary flow84.4253
Outlet flow91.15260
Table 2. Position, range and accuracy of sensors.
Table 2. Position, range and accuracy of sensors.
SensorsPositionUnitRangeAccuracy
TemperatureT1, T2, T3°C−40~40 ± 0.3
T4, T5, T6°C0~100 ± 0.3
PressureP1, P2, P3Bar−1~8 ± 0.5 %
P4, P5, P6Bar−1~16 ± 0.5 %
Volume flow rateFlowmeter 1 L / h 6~60 ±1.6%
Flowmeter 2 L / h 6~60±1.6%
Flowmeter 3 L / h 6~60±1.6%
Table 3. Operating conditions of the TSE for CFD model validation.
Table 3. Operating conditions of the TSE for CFD model validation.
First-Stage Primary FlowFirst-Stage Secondary FlowSecond-Stage Secondary FlowOutflow
PTPTPTP
(kPa)(K)(kPa)(K)(kPa)(K)(kPa)
Group 1414.6294243.427884.425391.15
394.6292
374.6290
354.6288
334.6286
Group 2374.6290283.428284.4253
263.4280
243.4278
223.4276
203.4274
Group 3374.6290243.427888.4254
86.4253.5
84.4253
82.4252.5
80.4252
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yan, J.; Shu, Y.; Wang, C. Effect of Back Pressure on Performances and Key Geometries of the Second Stage in a Highly Coupled Two-Stage Ejector. Entropy 2022, 24, 1847. https://doi.org/10.3390/e24121847

AMA Style

Yan J, Shu Y, Wang C. Effect of Back Pressure on Performances and Key Geometries of the Second Stage in a Highly Coupled Two-Stage Ejector. Entropy. 2022; 24(12):1847. https://doi.org/10.3390/e24121847

Chicago/Turabian Style

Yan, Jia, Yuetong Shu, and Chen Wang. 2022. "Effect of Back Pressure on Performances and Key Geometries of the Second Stage in a Highly Coupled Two-Stage Ejector" Entropy 24, no. 12: 1847. https://doi.org/10.3390/e24121847

APA Style

Yan, J., Shu, Y., & Wang, C. (2022). Effect of Back Pressure on Performances and Key Geometries of the Second Stage in a Highly Coupled Two-Stage Ejector. Entropy, 24(12), 1847. https://doi.org/10.3390/e24121847

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop