Numerical Analysis of Steam Ejector Performance with Non-Equilibrium Condensation for Refrigeration Applications

: In recent years, there has been great interest in developing cooling systems with humidity-and temperature-independent control capabilities that can operate efﬁciently at varying temperatures. This paper proposes a bi-loop double-evaporator ejection–compression cycle, which utilizes low-grade heat and is suitable for the construction industry. The proposed cycle involves the concurrent operation of a vapor compression cycle and an ejector refrigeration cycle that enables it to handle altered pressure levels and operate with varying compression ratios all the way to a common condenser pressure. Conventional computational ﬂuid dynamics (CFD) approaches often model steam as an ideal gas with single-phase ﬂow. In contrast, this research employs the wet steam model to optimize ejector geometry. The wet steam model takes into account non-equilibrium water vapor condensation, thus providing a more precise assessment of spontaneous condensation behavior and its impact on ejector performance. When compared to the conventional dry gas model, the use of the wet steam model dramatically decreases the entrainment ratio error from 16.24% for single-phase steam to 3.92% when compared to experimental data. This study concentrates on four critical attributes of wet steam, including Mach number, droplet nucleation rate, average droplet radius, and liquid mass fraction, to develop a strategy for enhancing ejector performance and efﬁciency. The study demonstrates that optimal area and primary nozzle diameter ratios for the steam ejector are 5 and 2.4, respectively. Increasing the area ratio mitigates condensation intensity, thereby reducing the liquid mass fraction in the diffuser. Overall, this paper provides valuable insights into improving and optimizing ejector performance, thus highlighting the importance of considering the behavior of spontaneous condensation in ejector design and modeling.


Introduction
More than 60% of energy generated by fossil fuels is wasted, with over 50% being low-grade heat at temperatures below 275 • C [1].As a result, researchers and engineers are increasingly turning their attention to alternative energy sources such as waste heat and geothermal energy to replace electricity in running refrigeration systems [2][3][4].This trend has become especially prevalent in the construction industry, where low-grade heat is typically utilized in ejector refrigeration cycles (ERC), desiccant systems, and absorption systems [5].While LiBr-water and ammonia-water absorption systems are currently available on the market (Coefficient of Performance (COP) of absorption refrigeration system changes within the range of 0.5 to 0.8 [6]), ejector cycles offer several advantages over these options.For instance, they can avoid issues related to internal corrosion and crystallization of the solution [7].Moreover, ERC is the only refrigeration system that can use nontoxic, nonflammable, and environmentally safe refrigerants (such as water) as its single-component working fluid, which happens to be both inexpensive and ecologically friendly [8].In regions where electricity is scarce or expensive, ejector refrigeration systems can provide an efficient and cost-effective alternative to conventional refrigeration technologies [9].This is because ejector refrigeration systems do not require any electrical power to operate.Instead, they rely on a combination of thermal energy input (such as hot water or steam) and the pressure differential created by the ejector to generate cooling.Scholars such as Braimakis [10], Tashtoush et al. [11], and Besagni et al. [12] have provided in-depth descriptions of the latest practical technologies in this field.
Despite numerous advantages, ERC face significant challenges in the market compared to vapor compression cycles (VCC).One major contributing factor is their typically lower COP which is often less than 0.5.Combining ERC and VCC into an ejection-compression refrigeration cycle, such as enhanced ERC and cascade ERC, is a promising solution to address this issue [13].For example, Sokolov [14] introduced a modified ERC that mechanically boosts the vapor to leave the evaporator, with the goal of increasing the ejector secondary pressure.Sun [15] built a cascade cycle that combines a compression subcycle and an ejection subcycle, with the intercooler serving as both an evaporator and a condenser for both subsystems.Xu et al. [16] proposed a hybrid ejection-compression refrigeration cycle for cooling in multi-story buildings, which aimed to reduce heat consumption and solar collector area.This system entrains the ejector secondary flow at an intermediate pressure from a separator, and the vapor exiting the ejector is mixed with the compressor discharge before entering the condenser.Results show superior electrical COP, with a maximal increase of 22.8%.Additionally, Ghaebi and Rostamzadeh [4] proposed a novel dual-evaporator ejection-compression refrigeration cycle for freezing and air-conditioning applications.Their analysis considered energy, exergy, exergoeconomic, and sustainability viewpoints, and computed an optimum range of COP (1.077-1.745)for various working fluids.Overall, these innovative designs represent promising developments in improving the efficiency and performance of ejection-compression refrigeration cycles.
The ejector geometry plays a critical role in determining the performance of an ejector refrigeration system.For certain operating conditions, there is an optimal or appropriate ejector geometry that provides the maximum entrainment ratio and critical condensation pressure.Several mathematical models for ejector geometrical parameters have been proposed based on the "one-dimensional flow ejector theory" [17].For example, Huang et al. [18] established a variable geometry ejector model to evaluate ejector performance under critical mode operation, assuming constant pressure mixing in an equal area section.Meanwhile, Chen et al. [19] proposed a fixed geometry ejector model that is suitable for both critical and sub-critical operations.Furthermore, detailed ejector geometric optimization can be performed using computational fluid dynamics (CFD) simulations.CFD can provide a detailed understanding of the flow behavior within the ejector and can help optimize the ejector geometry to achieve maximum performance for a given set of operating conditions.Besagni and Inzoli [20,21] provided a comprehensive review of ejector design, highlighting the importance of two key geometric factors (see Table 1): the area ratio of the ejector constant area chamber to the main throat area (AR) and the nozzle exit position (NXP).However, other geometrical parameters (such as the nozzle diameter and length, the diffuser angle, and the constant area chamber length) also play important roles in determining ejector performance.
Fixed geometry ejectors have limited applicability due to their sensitivity to operational conditions and lack of controllability.For instance, changes in external temperature necessitate the need for the primary and secondary flow rate control in ERC, while the cooling load and heat source power of solar ejector refrigeration cycle (SERC)fluctuate over time.Therefore, the variable geometry ejector (VGE) is proposed and analyzed [26].Movable spindles are a widely applied solution, allowing AR to be controlled by adjusting the front and rear positions of the spindles [27].Changing the position of the motive nozzle is another option.This technique implies that the NXP can remain unchanged if the ejector is designed and optimized properly.In contrast, the spindle position affects all the operating variables, such as cooling capacity, COP, generator pressure, generator heat, as well as mass flow rates.Table 2 briefly summarizes the literature of this study.Despite numerous theoretical studies, a gap still remains between their results and practical applications, primarily because of the use of dry gas assumptions and the neglect of the phase transition process.Ignoring the impact of the condensation shock on fluid density renders some conclusions highly inaccurate and some results even invalid.

Model Area Ratio and Nozzle Exit Position
Chen et al. [22] CFD, natural gas The NXP varies between 3.6 and 7.2 mm to obtain the maximum entrainment ratio.In addition, it is displayed between 1.2 and 7.2 mm for the optimal pressure ratio.

Pei et al. [24] CFD, Hydrogen gas
When the NXP exceeds the optimal range, the hydrogen entrainment ratio drops sharply in the entire operating range.The optimum AR ranges from 3 to 3.54.

Bai et al. [25]
Experiment, R23/R600a The main factors affecting the available freezing temperature and the cooling rate are the mixing chamber length and the nozzle throat, rather than the mixing chamber diameter and the position of the nozzle outlet.

R152a
R600a By adjusting the position of the spindle, an improvement in the entrainment ratio of up to 177% is achieved.

Pereira et al. [29]
Experimental Spindle R600a The COP is increased up to 85% with the variable geometry ejector.
Galindo et al. [30] Theoretical Spindle R1234yf There was an improvement in COP from 0.34 and 0.31 (in July and May, respectively) to 0.42 and 0.48.
Van Nguyen et al. [31] Experimental Spindle, NXP R600a T e = 8.5 The COP is increased by 24% compared to the fixed geometry ejector.
Besagni and Cristiani [27] Theoretical Spindle R290 Increasing the area ratio by 33% resulted in an average COP enhancement of 57.1%.
Tang et al. [32] conducted an experiment to visualize two-phase choking flow in ejectors.Their study revealed that an increased primary fluid pressure led to a reduction in the choking area, which, in turn, resulted in a decreased entrainment ratio as the choking area moved upstream.Several ejector CFD optimization studies, such as those conducted by Fu et al. [33], Dong et al. [34], Hou et al. [35], and Yan et al. [36], do not account for the effects of spontaneous steam condensation on ejector design.As supersonic ejectors have a significantly short residence time and high cooling rates, a considerable deviation from the equilibrium process occurs.Consequently, an accurate condensation model is essential for predicting ejector performance reliably.Ariafar et al. [37] compared the ideal-gas model with the condensation model and concluded that the latter resulted in a 10% higher entrainment ratio and a 7% higher critical backpressure.In a similar study, Wen et al. [38] found that the dry gas model leads to an unrealistic temperature distribution, resulting in an underestimation of entropy loss by 15% compared to non-equilibrium condensation models.Conversely, Ding et al. [39] found that the dry gas model led to an overestimation of entropy destruction by approximately 16%.Despite several studies in recent decades that focused on spontaneous condensation inside the nozzle of steam ejectors, these studies were limited to specific components such as the nozzle [40][41][42][43].As for the entire steam ejector, researchers have only paid attention to the effect of steam spontaneous condensation on internal flow and the optimization of numerical algorithms (Table 3).A thorough comprehension of the performance characteristics of steam ejector flows is essential for extending their industrial applications.Hence, the steam condensation effect inside the entire ejector needs to be carefully examined, especially with regard to optimizing the key geometric parameters while taking into account non-equilibrium condensation (NEC).The NEC process is of great significance to the ejector performance evaluation for ERC.Nevertheless, the optimization of the ejector structure through the NEC model has rarely been considered.

Reference Remarks
Zhang et al. [40] This study presented a modified condensation model to optimize steam ejector performance.The primary nozzle is optimized using the Multi-Objective Genetic Algorithm method, resulting in a 27.5% increase in entrainment ratio.
Li et al. [44] This study investigated the relationship between double choking characteristics and ejector performance.Results showed that the double choking mode occurs when the minimum distance between the sonic velocity line and the wall is 0.21 mm, and the second choking position is 4 mm downstream of the diffuser entrance.
Han et al. [45] This study analyzed the effect of phase change on the performance of a hydrogen recirculation ejector in proton exchange membrane fuel cell systems.The results showed that a higher entrainment ratio predicted by the two-phase flow model compared to the single-phase flow model.Droplet nucleation occurs at the junction of primary and secondary flow, with the nucleation rate increasing with primary flow pressure.
Dadpour et al. [46] This study investigated the effect of droplets injection at secondary flow on the performance of ejector refrigeration cycle.Results showed that a decrease in COP, and entrainment ratio with increasing wetness and number of droplets in the secondary flow.
Li et al. [47] The study examined how variable motive pressures affect the performance and shock waves in the system.Results indicated that increasing the motive pressure led to an improvement in the entrainment ratio and a reduction in the strength of the condensation shock wave, but it also resulted in an increase in the total pressure loss coefficient.
In the past, researchers have concentrated on enhancing the performance of ejectioncompression refrigeration cycles by incorporating economizers or other modifications.However, limited research has been carried out on utilizing dual-loop bi-evaporator ejection-compression refrigeration cycles to offer different cooling capacities for each loop.An air conditioning system that provides temperature-and humidity-independent control can remove moisture load through the low-temperature evaporator and heat load using the high-temperature evaporator.The application of the proposed cycle is not limited to either freezing or air-conditioning systems.Instead, it holds a broader relevance in the development of energy-efficient and environmentally friendly refrigeration systems and could spark the interest of both manufacturers and researchers in the field.As far as the authors know, no research has been conducted on ejector design and the application of dual-loop bi-evaporator ERC for office air conditioning, indicating a lack of research in investigating the dual-loop bi-evaporator air conditioning performance of the ejection-compression air conditioner.Moreover, research has primarily focused on the impact of spontaneous steam condensation on the internal flow of steam ejectors and on enhancing numerical algorithms for individual ejector components.However, comprehending the characteristics of steam ejector flow holistically is critical for extending their industrial applicability.Thus, thoroughly examining the steam condensation effect throughout the entire ejector and optimizing characteristic geometrical parameters while factoring in non-equilibrium condensation (NEC) is paramount.While NEC is critical for evaluating ejector performance in ERC, little attention has been paid towards optimizing the ejector structure via NEC model.
To sum up, the novelty of this research is as follows: • The phase change condensing model is established for ejector refrigeration systems considering phase transition processes with the condensation and evaporation of massive droplets in supersonic flows.

•
The optimization of pivotal ejector parameters and ejector performance is investigated considering the behavior of spontaneous condensation phenomenon.

•
A novel conceptual configuration of dual-loop bi-evaporator ejector refrigeration cycles with independent temperature and humidity control strategy based on climate adaptation is proposed (see Figure 1) to deal with the varying sensible and latent loads in different climates.
The main purposes of present study are multi-fold and can be summarized as follows: • Proposing a novel dual-loop bi-evaporator ejection-compression refrigeration cycle.

•
Designing the first set of ejector dimensions for a dual-loop bi-evaporator refrigeration system.

•
Developing a non-equilibrium condensation model for the steam ejector.

•
Validating the wet steam model by comparing it with experimental data and dry gas model.

•
Studying the effect of the ejector area ratio and primary nozzle diameter ratio on the ejector performance using the wet steam model for a certain operating condition.
Buildings 2023, 13, x FOR PEER REVIEW 5 of 28 development of energy-efficient and environmentally friendly refrigeration systems and could spark the interest of both manufacturers and researchers in the field.As far as the authors know, no research has been conducted on ejector design and the application of dual-loop bi-evaporator ERC for office air conditioning, indicating a lack of research in investigating the dual-loop bi-evaporator air conditioning performance of the ejectioncompression air conditioner.Moreover, research has primarily focused on the impact of spontaneous steam condensation on the internal flow of steam ejectors and on enhancing numerical algorithms for individual ejector components.However, comprehending the characteristics of steam ejector flow holistically is critical for extending their industrial applicability.Thus, thoroughly examining the steam condensation effect throughout the entire ejector and optimizing characteristic geometrical parameters while factoring in nonequilibrium condensation (NEC) is paramount.While NEC is critical for evaluating ejector performance in ERC, little attention has been paid towards optimizing the ejector structure via NEC model.
To sum up, the novelty of this research is as follows: • The phase change condensing model is established for ejector refrigeration systems considering phase transition processes with the condensation and evaporation of massive droplets in supersonic flows.

•
The optimization of pivotal ejector parameters and ejector performance is investigated considering the behavior of spontaneous condensation phenomenon.

•
A novel conceptual configuration of dual-loop bi-evaporator ejector refrigeration cycles with independent temperature and humidity control strategy based on climate adaptation is proposed (see Figure 1) to deal with the varying sensible and latent loads in different climates.
(a) The main purposes of present study are multi-fold and can be summarized as follows: • Proposing a novel dual-loop bi-evaporator ejection-compression refrigeration cycle.

•
Designing the first set of ejector dimensions for a dual-loop bi-evaporator refrigeration system.

•
Developing a non-equilibrium condensation model for the steam ejector.

•
Validating the wet steam model by comparing it with experimental data and dry gas model.

•
Studying the effect of the ejector area ratio and primary nozzle diameter ratio on the ejector performance using the wet steam model for a certain operating condition.

Description of Bi-Evaporator Ejection-Compression Refrigeration Cycle
Figure 1a illustrates the proposed bi-evaporator ejection-compression refrigeration cycle (BECRC) plant layout and corresponding P-h diagram.In such refrigeration systems, each loop operates at a different pressure level and is compressed to a common pressure (condenser pressure) with differing compression ratios.The VCC stage involves compressing the vapor that emerges from the low-pressure evaporator (state 12) to high pressure (state 13) using a compressor, and then mixing it with the outlet of the ejector (state 4).In the condenser, the mixed flow (state 5) is condensed into a saturated liquid by rejecting heat to the surroundings (state 6).The resulting liquid is then split into two streams, with one stream going to the low-pressure evaporator after throttling (state 11), completing the low-pressure loop.In the ERC stage, the remaining flow is divided into two streams.One of the streams passes through the expansion valve and is throttled to a lowerpressure, two-phase flow (state 8), while the other stream enters the generator after being pumped up to a higher pressure (state 7).The high-pressure vapor from the generator (state 1) entrains the low-pressure vapor from the high-temperature evaporator (state 7), and the resulting superheated flow leaves the ejector and enters the condenser (state 4), thereby completing the high-pressure loop.Such a dual-temperature system has potential for broad applications, where the radiant panel cooling system employs a high-temperature cold source (18/23 °C) and the air handling unit uses a low-temperature cold source (7/12 °C).For more details and discussions on such proposed systems, please refer to our published paper [48].

Description of Bi-Evaporator Ejection-Compression Refrigeration Cycle
Figure 1a illustrates the proposed bi-evaporator ejection-compression refrigeration cycle (BECRC) plant layout and corresponding P-h diagram.In such refrigeration systems, each loop operates at a different pressure level and is compressed to a common pressure (condenser pressure) with differing compression ratios.The VCC stage involves compressing the vapor that emerges from the low-pressure evaporator (state 12) to high pressure (state 13) using a compressor, and then mixing it with the outlet of the ejector (state 4).In the condenser, the mixed flow (state 5) is condensed into a saturated liquid by rejecting heat to the surroundings (state 6).The resulting liquid is then split into two streams, with one stream going to the low-pressure evaporator after throttling (state 11), completing the low-pressure loop.In the ERC stage, the remaining flow is divided into two streams.One of the streams passes through the expansion valve and is throttled to a lower-pressure, two-phase flow (state 8), while the other stream enters the generator after being pumped up to a higher pressure (state 7).The high-pressure vapor from the generator (state 1) entrains the low-pressure vapor from the high-temperature evaporator (state 7), and the resulting superheated flow leaves the ejector and enters the condenser (state 4), thereby completing the high-pressure loop.Such a dual-temperature system has potential for broad applications, where the radiant panel cooling system employs a high-temperature cold source (18/23 • C) and the air handling unit uses a low-temperature cold source (7/12 • C).For more details and discussions on such proposed systems, please refer to our published paper [48].
Working fluids for ERC are classified under two main categories: wet vapor and dry vapor, based on whether phase change occurs during the expansion process or not [49].For dry vapor, the ideal gas dynamics principle is used to model ejectors, as there is no phase change during the expansion process through the converging-diverging nozzle of ejectors.For wet vapor, ejector modeling should incorporate the phase change that occurs during the expansion process.The detailed thermophysical properties of these refrigerants are listed in Table 4. Figure 1b shows the enthalpy-entropy diagram of the working procedure in ERC.Expanding from state 1 to 2, the primary steam flow experiences a sharp pressure and temperature drop, entering the two-phase region where it becomes non-equilibrium and eventually supercooled as it expands further to the Wilson Line.A similar degree of supercooling is observed for the secondary flow, from state 9 to 10.The condensate in the diffuser evaporates as the pressure is restored from state 3 to 4. The non-equilibrium condensation model (Section 3) proposed in the article considers real gas equations of state to account for the non-ideal behavior of the refrigerants involved in the ejector-based refrigeration system.

Ejector Design
The rapid development of personal computers and commercial software has enabled more detailed and accurate analysis of ejector performance through CFD simulations.However, it is still not feasible to use CFD simulations directly for designing ejector geometries.Hence, a steam ejector is designed according to a one-dimensional theoretical model originally developed by Stoecker [50].The method requires iteration to calculate the flow states occurring inside the ejector and across the normal shock wave.At that time, calculating the flow state of the working fluid using only a calculator and conventional fluid properties table was quite challenging, especially at very low absolute pressure and near the freezing temperature.Although Stoecker's ejector model is a useful tool for designing and optimizing ejectors, its application to ERC is not yet fully established.
Table 5 summarizes the operating conditions used to run Stoecker's model, while Table 6 shows the main ejector dimensions.Reference to ESDU 94046 can provide guidance for the ejector design [51].The generator temperature is set as the motive flow inlet temperature, the evaporator temperature as the induced inlet temperature, and the condenser temperature as the mixed flow outlet temperature.Additionally, this study refers to the empirical formula reported by Zhu et al. [52] to design the initial structural parameters of the ejector.Reference to ESDU 94046 can provide guidance for the ejector design [51].The distance between nozzle exit and diffuser inlet (L 1 + L 2 ) is approximately 5-10 times the throat diameter (D 2 ); the throat length (L 2 ) should be from 2 to 4 D 2 ; the conical inlet of the mixing chamber should have an included angle of 2-10 • ; the diffuser should have a moderate divergence angle of 3-5 • ; the nozzle converging angle should be in the range of 20-40 • ; the nozzle diverging angle should be in the range of 10-12 • .Therefore, the convergent and divergent angles of the primary nozzle are equal to 20 • and 10 • , respectively; the conical inlet of the mixing chamber is equal to 6 • and the diffuser divergent angle is 4 • ; the nozzle outlet and the mixing chamber inlet are at the same cross-section.
The design procedure outlined in ESDU assumes that both primary and secondary flow exhibit ideal gas behavior with a specific heat ratio of 1.315 and a specific gas constant of 461.5 J•kg −1 •K −1 .However, this assumption is often extended to saturated steam as well under the assumption that condensation will not occur due to metastable behavior.This is because there are no nucleation sites for condensation to occur, and the steam remains in a superheated state until it encounters a surface or particle where condensation can occur.Therefore, it is important to verify the accuracy of such assumptions and consider nonequilibrium condensation models if necessary for a more accurate design and optimization of ejectors in ERC systems.

Governing Equations
The non-equilibrium condensation model includes the nucleation and droplet growth processes, in which the former process is carried out according to the classical nonisothermal nucleation theory while the latter involves the droplet growth model of Hill [54] and Yong [55].The assumption of no-slip velocity between the droplets and the gaseous phase is made since the droplet size and mass fraction are small [56].The flow field of the gas phase is modeled using the condensation phase transition theory, while the liquid phase is modeled using the compressible Navier-Stokes equations [57][58][59].
where u i , u j , u k are components of the velocity; ρ, P, E, and T refer to the mixture density, pressure, gross energy, and temperature, respectively; τ ij represents the stress tensor; λ e f f represents the effective thermal conductivity; h lv is the specific enthalpy of condensation; Γ is the generation rate of liquid mass, m −3 •m −1 ; and I is the nucleation rate, kg Two scalar equations are composed to describe the liquid fraction (β), and the droplet number per volume (N): The steam properties are measured from the truncated Virial equation of state in the third expanded term: where B and C are both functions of temperature, representing the second virial coefficient and the third virial coefficient, respectively; R is the specific gas constant for gaseous mixture of air and vapor [60].
To determine the number of droplets per unit volume (N), the vapor density (ρ v ) and the average droplet volume (V d ) are expressed as: ) where r is the average droplet radius.
The nucleation and droplet growth model are employed for the simulation of nonequilibrium condensation.The Appendix A illustrates the condensation of wet steam in supersonic flows.

Numerical Setup
Pianthong et al. [61] demonstrated that two-dimensional (2D) axisymmetric and threedimensional (3D) models predict static pressure at the center line similarly.Studies have demonstrated that using a 2D axisymmetric model can yield equivalent results to the 3D model while requiring less computational effort [62].In terms of the turbulent model, Zhang et al. [63] proposed that the k − ω SST model can predict the intensity and position of shock wave more accurately; Schlieren optical measurements were utilized by Zhu and Jiang [64] to compare the performance of Realizable, RNG, and standard k − ε and k − ω SST models and the results revealed that the RNG k − ε model had the highest level of agreement with experimental data in terms of predicting both mass flow rate and shock wave structures; Elbarghthi et al. [65] found that the k − ω SST model is more accurate than the RNG k − ε model for ejector calculations.Additionally, the k − ω SST model is widely recommended in the literature due to its superior performance in stream mixing, convergence, and flow field prediction [66].Definitions of the k − ω SST model can be found in Refs.[21,67].
ANSYS Fluent 2020 was used to perform the numerical simulation.To solve for droplet number and liquid fraction, the UDS method was employed, while the UDF method was used to model mass generation in the process of condensation in transonic flow.The inlet flows are usually low speed, which means that the total pressure and temperature are almost similar to the static magnitudes.Turbulence intensity at the inlets is measured by the rate of the root-mean-square of the velocity fluctuation in the average velocity.Because the insensitivity of the outcomes to their variability is between 1% and 10%, the parameter is determined as 5%.Although the density-based algorithm is recommended widely for supersonic flow, the pressure-based solver can effectively handle highly compressible flows, including both single-phase and two-phase flow ejectors.While the density-based algorithm is commonly recommended for supersonic flow, the pressure-based solver was deemed effective for highly compressible flows, including single-phase and two-phase flow ejectors, and provided smoother and faster convergence, making it more suitable for problems that require a large number of simulations [68].To make sure that the outcome is convergent to the expected tolerance, three convergence monitors are utilized.The residuals are lower than 10 −5 and the mass flow is still unchanged within certain times of iterations; the outlet net mass flow needs to coincide with the inlet net mass flow within a tolerance of 10 −6 kg•s −1 ; the maximum velocity value at the ejector throat is constant.

Grid Independence
Figure 2 presents the quadrilateral structured grids and outlines the settings utilized in Fluent, with a focus on refining the mesh in regions with wide velocity and pressure variations for accurate prediction of flow behavior.The axial distances of the nozzle outlet, constant section inlet, and diffuser inlet are 42 mm, 157 mm, and 232 mm, respectively.y + refers to the friction velocity ratio between the wall and the fluid velocity profile, which should typically be controlled within the range of 30 to 100 in order to obtain accurate results.Figure 3 illustrates the axial distribution of static pressure across four grid levels: 25,000 elements (coarse level), 60,000 elements (medium level), 90,000 elements (fine level), and 110,000 elements (meticulous level).Table 7 displays the secondary mass flow rate and entrainment ratio at different grid levels.The Grid Convergence Index (GCI) for a mesh containing 90,000 elements is less than 1%.Therefore, a quadrilateral grid consisting of 90,000 cells is considered a suitable reference model for this study.
cause the insensitivity of the outcomes to their variability is between 1% and 10%, the parameter is determined as 5%.Although the density-based algorithm is recommended widely for supersonic flow, the pressure-based solver can effectively handle highly compressible flows, including both single-phase and two-phase flow ejectors.While the density-based algorithm is commonly recommended for supersonic flow, the pressure-based solver was deemed effective for highly compressible flows, including single-phase and two-phase flow ejectors, and provided smoother and faster convergence, making it more suitable for problems that require a large number of simulations [68].To make sure that the outcome is convergent to the expected tolerance, three convergence monitors are utilized.The residuals are lower than 10 and the mass flow is still unchanged within certain times of iterations; the outlet net mass flow needs to coincide with the inlet net mass flow within a tolerance of 10 kg • s ; the maximum velocity value at the ejector throat is constant.

Grid Independence
Figure 2 presents the quadrilateral structured grids and outlines the settings utilized in Fluent, with a focus on refining the mesh in regions with wide velocity and pressure variations for accurate prediction of flow behavior.The axial distances of the nozzle outlet, constant section inlet, and diffuser inlet are 42 mm, 157 mm, and 232 mm, respectively. refers to the friction velocity ratio between the wall and the fluid velocity profile, which should typically be controlled within the range of 30 to 100 in order to obtain accurate results.Figure 3 illustrates the axial distribution of static pressure across four grid levels: 25,000 elements (coarse level), 60,000 elements (medium level), 90,000 elements (fine level), and 110,000 elements (meticulous level).Table 7 displays the secondary mass flow rate and entrainment ratio at different grid levels.The Grid Convergence Index (GCI) for a mesh containing 90,000 elements is less than 1%.Therefore, a quadrilateral grid consisting of 90,000 cells is considered a suitable reference model for this study.

Model Validation
The research data of Sriveerakul et al. [70,71] is compared with the numerical re obtained by dry gas and wet steam models under various back pressures and the s geometric parameters, as shown in Figure 4. Compared to the dry gas model, the steam model fits the experimental results better.The average error of the  red from 16.24% (dry gas results) to a negligible difference of 3.92% (wet steam results).indicates that the dry gas model overestimates the  values in critical modes, as served in previous studies by Yang et al [42].Furthermore, the critical back pressure o experimental data is around 50 mbar: 45 mbar for dry gas simulation and 51 mbar for steam simulation.The predicted critical pressure carried out by the wet steam mode increased; in the on-design region (critical mode zone), the wet steam model exhib wider range compared to the dry gas model (the compression effect is caused by the pact of diffuser).These results align with those observed in prior studies by Zhang e [72].On the premise of non-equilibrium condensation, the error can be considered a fective modification to the predicted entrainment ratio.In contrast, since the wet st

Model Validation
The research data of Sriveerakul et al. [70,71] is compared with the numerical results obtained by dry gas and wet steam models under various back pressures and the same geometric parameters, as shown in Figure 4. Compared to the dry gas model, the wet steam model fits the experimental results better.The average error of the ER reduces from 16.24% (dry gas results) to a negligible difference of 3.92% (wet steam results).This indicates that the dry gas model overestimates the ER values in critical modes, as observed in previous studies by Yang et al [42].Furthermore, the critical back pressure of the experimental data is around 50 mbar: 45 mbar for dry gas simulation and 51 mbar for wet steam simulation.The predicted critical pressure carried out by the wet steam model has increased; in the ondesign region (critical mode zone), the wet steam model exhibits a wider range compared to the dry gas model (the compression effect is caused by the impact of diffuser).These results align with those observed in prior studies by Zhang et al. [72].On the premise of non-equilibrium condensation, the error can be considered an effective modification to the predicted entrainment ratio.In contrast, since the wet steam model involving nonequilibrium condensation can predict the ER and compression ratio more accurately, it is more suitable for investigating complex flow phenomena in a supersonic steam ejector.
model involving non-equilibrium condensation can predict the  and compression tio more accurately, it is more suitable for investigating complex flow phenomena i supersonic steam ejector.This is due to the fact that the expansion characteristics of supersonic flow in the ejec are limited by steam condensation, which absorbs energy from the gas-phase flow, resu ing in a decrease in fluid velocity and weaker shock waves.The dry steam model indica that the shock wave appears earlier in the mixing chamber compared to that predicted wet steam simulations.Additionally, the sonic line fails to reach the ejector wall in constant section, resulting in a distinctive shape of the characteristic curve, as shown Figure 4.   [57].This is due to the fact that the expansion characteristics of supersonic flow in the ejector are limited by steam condensation, which absorbs energy from the gas-phase flow, resulting in a decrease in fluid velocity and weaker shock waves.The dry steam model indicates that the shock wave appears earlier in the mixing chamber compared to that predicted by wet steam simulations.Additionally, the sonic line fails to reach the ejector wall in the constant section, resulting in a distinctive shape of the characteristic curve, as shown in Figure 4.

Dry and Wet Steam Ejector
Figure 6a,b provides further insight into the differences between the wet steam and dry steam models through profiles and contours of static temperature.The wet steam model predicts higher temperatures compared to the dry steam model, in line with findings by Yang et al. [73].The presence of droplets releases latent heat of condensation, which heats the gas phase in the wet steam model.In the dry steam model, the minimum temperature is 130 K at x = 56 mm, with most temperatures in the nozzle, mixing chamber, and constant section below freezing point (273.15K), which is an unrealistic outcome.In contrast, the lowest temperature in the wet steam model is 198 K at x = 59 mm, and only a small portion of the temperature profile in the mixing chamber falls below freezing point.The low temperature in the nozzle outlet is due to the high Mach number at the "shock chain" region.The dry gas assumption assumes that all the internal energy of the fluid is converted into kinetic energy, and the vapor temperature is much lower than the triple point temperature of water (the hypothesis of ice crystal formation), while droplet temperature (the same as the saturation temperature at the vapor pressure) is generally well above that of the vapor.The average static temperatures for dry and wet steam are 235.48K and 288.88 K, respectively.When droplets condense in the wet steam model, latent heat is released, which heats up the gas phase and leads to a temperature jump in the primary nozzle.Figure 6a,b provides further insight into the differences between the wet steam and dry steam models through profiles and contours of static temperature.The wet steam model predicts higher temperatures compared to the dry steam model, in line with findings by Yang et al. [73].The presence of droplets releases latent heat of condensation, which heats the gas phase in the wet steam model.In the dry steam model, the minimum temperature is 130 K at  = 56 mm, with most temperatures in the nozzle, mixing chamber, and constant section below freezing point (273.15K), which is an unrealistic outcome.In contrast, the lowest temperature in the wet steam model is 198 K at  = 59 mm, and only a small portion of the temperature profile in the mixing chamber falls below freezing point.The low temperature in the nozzle outlet is due to the high Mach number at the "shock chain" region.The dry gas assumption assumes that all the internal energy of the fluid is converted into kinetic energy, and the vapor temperature is much lower than the triple point temperature of water (the hypothesis of ice crystal formation), while droplet temperature (the same as the saturation temperature at the vapor pressure) is generally well above that of the vapor.The average static temperatures for dry and wet steam are 235.48K and 288.88 K, respectively.When droplets condense in the wet steam model, latent heat is released, which heats up the gas phase and leads to a temperature jump in the primary nozzle.The temperature inside the ejector is modeled to be approximately 200 K, which is significantly lower than typical operating temperatures for refrigeration systems.However, the involved area below the triple point of water is limited, and the sudden rise in temperature could prevent ice from forming.This is because the rate of crystallization is highly dependent on the temperature and can be significantly impacted by even small changes in temperature or humidity.
In summary, the single-phase model overestimates the entrainment ratio of the ejec- The temperature inside the ejector is modeled to be approximately 200 K, which is significantly lower than typical operating temperatures for refrigeration systems.However, the involved area below the triple point of water is limited, and the sudden rise in temperature could prevent ice from forming.This is because the rate of crystallization is highly dependent on the temperature and can be significantly impacted by even small changes in temperature or humidity.
In summary, the single-phase model overestimates the entrainment ratio of the ejector due to its prediction of a higher Mach number (lower pressure) and lower temperature compared to the two-phase model.The findings demonstrate the importance of considering phase change effects in the steam ejector.

Internal Analysis of the Ejector
The ratio of the outlet and throat cross-sectional area (A 4 /A th ), outlet-throat diameter ratio (γ d ), and area ratio (AR) are critical parameters that significantly affect the flow characteristics of a converging-diverging nozzle and the Mach number at the nozzle exit.It can be calculated as follows [74]: To achieve changes in A 4 /A th and A 2 /A th , the diameter of various ejector sections needs to be adjusted.The optimization sequence follows the order of the primary nozzle outlet diameter (D 4 ) and then the diameter of the mixing chamber (D 2 ).Therefore, a detailed study of the different γ d and AR, with a fixed throat diameter D th , is conducted while taking the phase change process into consideration.
The ejector performance is measured in terms of the entrainment ratio, whilst the performance of the whole system is evaluated in terms of the COP.The change in ER in the ERC is closely related to the thermal COP (COP t = ER h 9 −h 6 h 1 −h 6 ) [75].As depicted in Figure 7, when AR reaches 6.0, the secondary flow inlet refluxes with an ER below −0.1, indicating the ejector loses its entrainment performance.It is indicated that the ERC achieves a maximum ER of 0.59 when the γ d is 2.4 and AR is 5.0, while the maximum deviation of ER for γ d and AR are 16.64% and 391.95%, respectively.AR has a great influence on the ER as compared to that of γ d .
The objective of this study is not only to facilitate the practical design of the ejector, but also to investigate the variations in flow behaviors under different parameters.The flow pattern cannot be clarified by the entrainment ratio.It is crucial to carry out research on several wet steam parameters, including Mach number, droplet average radius, nucleating rate, and liquid mass fraction, and so on.

Effect of the Primary Nozzle Outlet Diameter
Figure 8a demonstrates the variations in Mach number at different γ d .At γ d = 2.0, 2.4, and 2.8, their corresponding Mach values are 2.28, 2.52, and 2.72, respectively.For γ d = 2.0, the steam experiences a violent expansion in the mixing chamber, resulting in a Mach number of up to 4, which induces an under-expanded flow.As the γ d increases to 2.4 and 2.8, the flow structure transition slightly decreases the maximum Mach value in the mixing section, resulting in a decrease in the maximum Mach number to approximately 3.39 and 3.34, respectively.Nonetheless, the results still present a moderate under-expanded flow in the flow structure.
the ERC is closely related to the thermal COP ( = ) [75].As depicted in Figure 7, when  reaches 6.0, the secondary flow inlet refluxes with an  below −0.1, indicating the ejector loses its entrainment performance.It is indicated that the ERC achieves a maximum  of 0.59 when the  is 2.4 and  is 5.0, while the maximum deviation of  for  and  are 16.64% and 391.95%, respectively. has a great influence on the  as compared to that of  .Figure 8b depicts the non-equilibrium condensation process, which generates nanodroplets in the steam ejector.The droplet nucleation occurs inside the primary nozzle and at its exit region.As the γ d increases, the droplet average radius decreases, causing the maximum droplet average radius to move upstream.The maximum average radius of the droplet is approximately 39 nm in the mixing section for γ d = 2.0, whereas it reduces to 26 nm in the primary nozzle for γ d = 2.8.Additionally, the homogeneous nucleation process occurs for various γ d , indicating slightly different nucleating rates, as shown in Figure 8c.
Considerable heat is released into the supersonic flow during droplet growth, leading to compressive effects that bring about condensation shock.A moderate under-expanded flow with outlet-throat diameter ratios of γ d = 2.4 and 2.8 triggers a homogeneous nucleation process in both the mixing section and constant section.However, non-equilibrium condensations are absent in the constant section for γ d = 2.0.As γ d increases, the Wilson points gradually converge, and the initial position of the droplet nucleation rate moves forward.It can be inferred that the formation and evaporation of nanodroplets may be influenced by varying γ d .
2.4, and 2.8 , their corresponding Mach values are 2.28, 2.52, and 2.72, respectively.For  = 2.0, the steam experiences a violent expansion in the mixing chamber, resulting in a Mach number of up to 4, which induces an under-expanded flow.As the  increases to 2.4 and 2.8, the flow structure transition slightly decreases the maximum Mach value in the mixing section, resulting in a decrease in the maximum Mach number to approximately 3.39 and 3.34, respectively.Nonetheless, the results still present a moderate underexpanded flow in the flow structure.Figure 8b depicts the non-equilibrium condensation process, which generates nanodroplets in the steam ejector.The droplet nucleation occurs inside the primary nozzle and at its exit region.As the  increases, the droplet average radius decreases, causing the maximum droplet average radius to move upstream.The maximum average radius of the droplet is approximately 39 nm in the mixing section for  = 2.0, whereas it reduces Subsequently, Figure 8d illustrates the impact of γ d on the liquid fraction profiles.Initially, liquid droplets form in the ejector due to nucleation.As soon as the nucleation process is complete, the droplets will begin to grow.As the temperature decreases below saturation temperature, the vapor condenses, and the liquid mass fraction increases, whereas as the liquid droplets evaporate, the liquid mass fraction decreases.The conclusion can be extracted that the liquid mass fraction increases with the rising of γ d in the primary nozzle.The minimum liquid mass fraction at the nozzle outlet is 9.14% (γ d = 2.0), while the liquid mass fraction increases to 9.41% when γ d = 2.8.The maximum liquid mass fraction of 8% reported by Zhang et al. [76] and 19% reported by Sharifi [77].The maximum liquid mass fraction is attained in the constant section, followed by a gradual decrease due to droplet evaporation within the diffuser.As such, the liquid content reduces progressively.
As shown in Figure 8a, the basic driving potential of droplet vaporization and condensation is the low-pressure region caused by the transformation of the Mach number.In addition, a larger γ d not only reduces the ER but also causes accelerated droplet nucleation rate and more liquid fractions, as shown in Figure 8c,d.Figure 9b illustrates the variation trend of the droplet average radius of each AR ejectors with axial coordinates.The average radius of droplets first increases, and then decreases; the maximum droplet average radius appears at approximately 36 mm in the primary nozzle.The rise in AR does not alter the average radius of droplets.In other words, droplet growth only exists in the primary nozzle, and droplet growth no longer continues when the flow enters the mixing section.Furthermore, before the diffuser outlet, the profiles of various AR ejectors show a constant trend.In fact, there are a lot of small droplets in this region.No matter how the AR changes, the droplet average radius remains unchanged, indicating the little impact of AR on the droplet average radius.

Effect of the Mixing Chamber Throat
Figure 9c shows the profiles of droplet nucleation rates with different AR along the axis of the ejector.Once the nucleation rate reaches its peak, the droplet radius undergoes a rapid increase, compared to Figure 10b.The nucleation occurs in two regions for the AR = 4 ejector.One is in the primary nozzle, and the supercooling impact of the divergent section leads to the supersonic flow reaching the equilibrium state.The other two nucleation processes occur in the mixing section.For the AR = 5 ejector, the nucleation occurs in three regions; the difference is that nucleation occurs again in the constant section.The ending position of the droplet nucleation rate moves forward with the increase of the AR. Figure 9b illustrates the variation trend of the droplet average radius of each  ejectors with axial coordinates.The average radius of droplets first increases, and then decreases; the maximum droplet average radius appears at approximately 36 mm in the primary nozzle.The rise in  does not alter the average radius of droplets.In other words, droplet growth only exists in the primary nozzle, and droplet growth no longer continues when the flow enters the mixing section.Furthermore, before the diffuser outlet, the profiles of various  ejectors show a constant trend.In fact, there are a lot of small droplets in this region.No matter how the  changes, the droplet average radius remains unchanged, indicating the little impact of  on the droplet average radius.
Figure 9c shows the profiles of droplet nucleation rates with different  along the axis of the ejector.Once the nucleation rate reaches its peak, the droplet radius undergoes a rapid increase, compared to Figure 10b.The nucleation occurs in two regions for the  = 4 ejector.One is in the primary nozzle, and the supercooling impact of the divergent section leads to the supersonic flow reaching the equilibrium state.The other two nucleation processes occur in the mixing section.For the  = 5 ejector, the nucleation occurs in three regions; the difference is that nucleation occurs again in the constant section.The ending position of the droplet nucleation rate moves forward with the increase of the .
The distribution of liquid mass fractions is shown in Figure 9d.In the primary nozzle and mixing chamber, there is no effect of  on the liquid mass fraction distribution.The condensed phase is up to 10%, and its maximum value appears in the downstream area of the nozzle exit plane.In the constant section, the liquid mass fraction under  = 6 rapidly decreases to 0. The increase in  decreases the liquid mass fraction in the constant section and diffuser.Furthermore, with the increase of , the liquid mass fraction in the diffuser chamber decreases, and the intensity of non-equilibrium condensation in the ejector is reduced.The frequent occurrence of steam condensation and evaporation behaviors is indicated by the fluctuations of the liquid mass fraction curves in the diffuser at various .
(a)  The distribution of liquid mass fractions is shown in Figure 9d.In the primary nozzle and mixing chamber, there is no effect of AR on the liquid mass fraction distribution.The condensed phase is up to 10%, and its maximum value appears in the downstream area of the nozzle exit plane.In the constant section, the liquid mass fraction under AR = 6 rapidly decreases to 0. The increase in AR decreases the liquid mass fraction in the constant section and diffuser.Furthermore, with the increase of AR, the liquid mass fraction in the diffuser chamber decreases, and the intensity of non-equilibrium condensation in the ejector is reduced.The frequent occurrence of steam condensation and evaporation behaviors is indicated by the fluctuations of the liquid mass fraction curves in the diffuser at various AR.

Two-Phase Flow Features
For a clear view of the simulation results, Figure 10 illustrates flow structures in the steam ejector, including Mach number, subcooled vapor temperature, droplet nucleation rate, and liquid mass fraction, respectively.The supersonic flow becomes supercooled and produces vast condensation nuclei and droplets after the primary flow expands strongly.The mixed flow pressure drops significantly, producing a lot of condensed droplets.

Two-Phase Flow Features
For a clear view of the simulation results, Figure 10 illustrates flow structures in the steam ejector, including Mach number, subcooled vapor temperature, droplet nucleation rate, and liquid mass fraction, respectively.The supersonic flow becomes supercooled and produces vast condensation nuclei and droplets after the primary flow expands strongly.The mixed flow pressure drops significantly, producing a lot of condensed droplets.Meanwhile, infinitesimal droplets absorb part of the flow kinetic energy, causing the flow velocity to decrease.
As shown in Figure 10a, the subcooling degree is approximately 37 K at the nozzle throat, which leads to the first nucleation and condensation of the steam ejector with a peak nucleation rate of 5.62 × 10 20 m −3 •m −1 .The supercooled vapor temperature stays at approximately 50 K at the primary nozzle exit.The peak values of supercooled vapor temperature are found to decrease gradually in the mixing and constant section, while the trough values increase, resulting in the condensation and re-evaporation of steam.The four peak values of subcooled vapor temperature are 80 K, 55 K, 31 K, and 23 K, and the three troughs are 10 K, 17 K, and 15 K.The corresponding peak values of droplet nucleation rate are 1.51 × 10 26 m −3 •m −1 , 3.02 × 10 24 m −3 •m −1 , 2.45 × 10 16 m −3 •m −1 , and 5.13 × 10 6 m −3 •m −1 .The supercooled vapor temperature is approximately equal to 0 at downstream of the constant section (x > 220 mm), and the droplet nucleation rate is around 0 in this region (x > 175 mm).However, the liquid droplets do not disappear.Figure 10b displays that the liquid mass fraction attains its highest value of 10% at the nozzle outlet and decreases to a minimum in the diffuser due to droplet evaporation.After the droplet nucleation rate reaches its peak, the liquid mass fraction undergoes a sharp increase and reaches its peak value.Figure 10c shows that the liquid mass fraction near the ejector wall is almost zero.This is because the fluid deceleration and viscous dissipation in the boundary layer lead to some extent of heat recovery, resulting in the absence of liquid droplets in the region (represented by the blue color in the liquid mass fraction contour).The distribution of liquid mass generation rate is primarily concentrated in the mixing chamber due to the heat exchange occurring during the mixing process.
According to the experiment results, the homogeneous nucleation of ice is usually observed at a subcooling degree of about 30-40 K. Simulation results indicate that the supercooled vapor temperature reaches up to 80 K.However, the formation of ice crystals is basically a similar time-dependent phenomenon as the nucleation of liquid droplets.Although the primary flow can reach a temperature level well below the triple point, ice may not occur since the involved area is limited and the temperature suddenly rises.Alternatively, shocks may cause ice to melt suddenly if they occur immediately after steam expansion.Given the uncertain characteristics of supercooled water below the triple point temperature, (surface tension, specific heat, viscosity, etc.), the ice formation in the steam ejector should be further analyzed in future experiments.

Conclusions and Future Works
This paper proposes a novel bi-loop double-evaporator ejection-compression refrigeration cycle, which can generate refrigeration at two temperature levels and be used in temperature-and humidity-independent control systems.The cycle features a combination of vapor compression refrigeration and ejector refrigeration sub-cycles, allowing the compression of each pressure level to a common pressure (the condenser pressure) at varying compression rates.To examine the steam ejector performance and optimize the ejector geometry, homogeneous, non-equilibrium, two-phase flow computational fluid dynamics model is applied.The wet steam model is found to be more suitable than the dry gas model in predicting the entrainment ratio when compared to experimental data, with an error decreasing from 16.24% to 3.92%.Some remarkable findings are summarized below:

•
Neglecting the phase change process, the dry gas model results in an overestimation of steam expansion and lower temperature predictions compared to the wet steam model.The average Mach numbers are 1.97 and 1.76 for dry and wet steam models, respectively.The average static temperatures between dry and wet steam are 235.48K and 288.88 K, respectively.

•
The optimum values for the primary nozzle outlet diameter γ d and area ratio AR are 2.4 and 5.0, respectively.The optimized ER reaches up to 0.59.In addition, the maximum deviation of ER for γ d and AR is 16.64% and 391.95%, respectively.

•
When the γ d increases to γ d = 2.4 and 2.8, the maximum Mach number decreases to about 3.39 and 3.34.The minimum liquid mass fraction is 9.14% (γ d = 2.0) at the nozzle outlet, while there is an approximate growth of 9.41% when γ d = 2.8.

•
The AR = 5 ejector reduces 60.46% of outlet Mach number (Ma = 0.43) as compared to the AR = 4 ejector (Ma = 0.69).Droplet growth only takes place within the primary nozzle and no longer continues as the flow enters the mixing section.The nucleation occurs in two regions for the AR = 4 ejector, and the nucleation occurs in three regions for the AR = 5 ejector.
The proposed cycle has potential applications for both freezing and air-conditioning systems and could be of interest to manufacturers and researchers looking to develop more energy-efficient and environmentally friendly refrigeration systems.This study provides reference for predicting the geometrical parameter and performance of supersonic steam ejectors.However, there are limitations to this study, including its focus on a specific type of steam ejector and its reliance on numerical simulations rather than experimental data.The following areas could be explored in future research: (1) Exploring the effect of different geometries and operating conditions on the performance of supersonic steam ejectors with non-equilibrium condensing; (2) investigating the use of alternative working fluids in supersonic steam ejectors, such as refrigerants or compressed air; (3) developing new models to account for different wet fluids, as they may have varying nucleation and droplet growth rates; and (4) performing optimization and parametric studies in terms of the proposed conceptual configuration (energy, exergy, exergoeconomic, and sustainability).Additionally, experimental validation may be necessary to confirm the advantages in real-world applications.represents the droplet density; r represents the average radius of the droplet; ∂r ∂t represents the growth rate of droplets; r * refers to the Kelvin-Helmholtz critical radius.Droplet growth occurs when r exceeds r * , while droplet evaporation takes place when r is less than r * .

Nomenclature
An expression for r is given by: An expression for r * is given by: where ρ l refers to the liquid density; σ refers to the liquid surface tension; Sr refers to the super saturation ratio, that is, the ratio of vapor pressure to the equilibrium saturation pressure: where P sat (T) is the saturation pressure at temperature T.
According to the modified classical nucleation model using Kantrowitz's non-isothermal correction, the nucleation rate I can be obtained as shown below: where q c = 1 refers to condensation coefficient, M refers to the molecule mass of vapor, K b refers to the Boltzmann constant (1.3807 × 10 −23 J•K −1 ); and θ is the non-isothermal correction coefficient (Kantrowitz correction).
where γ is the ratio of specific heat capacity.Two mechanisms are involved during condensation: mass transfer from vapor to droplets and heat transfer from droplets to vapor as latent heat.Hill (1966) [54] proposed this energy transfer relation with the equation below: Gyarmathy's droplet growth model is: Young (1982) [55] modified the droplet growth model proposed by Gyarmathy (1962) [78]; this formula is tunable with two modeling parameters, α = 9 and φ = 1: (1 − r * /r) An expression for T d is given by: where ∆T is the subcooled vapor temperature, and T sat (P) is the saturation temperature at pressure P.

Figure 1 .
Figure 1.Schematic and P-h diagram of the proposed bi-evaporator ejection-compression refrigeration cycle (a); Enthalpy-entropy graphics of the working procedure for wet steam (b).

Figure 2 .
Figure 2. CFD model settings and structured meshing of the baseline geometry [69].

Figure 2 .
Figure 2. CFD model settings and structured meshing of the baseline geometry [69].

Figure 3 .
Figure 3. Axial static pressure distribution at various grid levels.

Figure 3 .
Figure 3. Axial static pressure distribution at various grid levels.

Figure 4 .
Figure 4. Validation of the wet steam ejector model.

Figure
Figure 5a,b compares the flow Mach numbers (i.e.,  =   ⁄ ) of dry and w steam.When the flow temperature is below the saturation temperature of the flow pr sure, some vapor undergoes rapid condensation in a transient process.The local ma mum Mach number is 3.28 at  = 59 mm in the wet steam model, while in the dry model, it is 3.54 at  = 56 mm.The average Mach value between wet and dry steam m els is 1.76 and 1.97, respectively, showing that the dry gas model overpredicts the M numbers by 11.93%.These results are consistent with previous studies by Sharifi et al. [5This is due to the fact that the expansion characteristics of supersonic flow in the ejec are limited by steam condensation, which absorbs energy from the gas-phase flow, resu ing in a decrease in fluid velocity and weaker shock waves.The dry steam model indica that the shock wave appears earlier in the mixing chamber compared to that predicted wet steam simulations.Additionally, the sonic line fails to reach the ejector wall in constant section, resulting in a distinctive shape of the characteristic curve, as shown Figure4.

Figure 4 .
Figure 4. Validation of the wet steam ejector model.

Figure
Figure 5a,b compares the flow Mach numbers (i.e., Ma = u/√ γRT) of dry and wet steam.When the flow temperature is below the saturation temperature of the flow pressure, some vapor undergoes rapid condensation in a transient process.The local maximum Mach number is 3.28 at x = 59 mm in the wet steam model, while in the dry gas model, it is 3.54 at x = 56 mm.The average Mach value between wet and dry steam models is 1.76 and 1.97, respectively, showing that the dry gas model overpredicts the Mach numbers by 11.93%.These results are consistent with previous studies by Sharifi et al.[57].This is due to the fact that the expansion characteristics of supersonic flow in the ejector are limited by steam condensation, which absorbs energy from the gas-phase flow, resulting in a decrease in fluid velocity and weaker shock waves.The dry steam model indicates that the shock wave appears earlier in the mixing chamber compared to that predicted by wet steam simulations.Additionally, the sonic line fails to reach the ejector wall in the constant section, resulting in a distinctive shape of the characteristic curve, as shown in Figure4.Figure6a,b provides further insight into the differences between the wet steam and dry steam models through profiles and contours of static temperature.The wet steam model predicts higher temperatures compared to the dry steam model, in line with findings by Yang et al.[73].The presence of droplets releases latent heat of condensation, which heats the gas phase in the wet steam model.In the dry steam model, the minimum temperature is 130 K at x = 56 mm, with most temperatures in the nozzle, mixing chamber, and constant section below freezing point (273.15K), which is an unrealistic outcome.In contrast, the lowest temperature in the wet steam model is 198 K at x = 59 mm, and only a small portion of the temperature profile in the mixing chamber falls below freezing point.The low temperature in the nozzle outlet is due to the high Mach number at the "shock chain" region.The dry gas assumption assumes that all the internal energy of the fluid is converted into kinetic energy, and the vapor temperature is much lower than the triple point temperature of water (the hypothesis of ice crystal formation), while droplet temperature (the same as the saturation temperature at the vapor pressure) is generally well above that of the vapor.The average static temperatures for dry and wet steam are 235.48K and 288.88 K, respectively.When droplets condense in the wet steam model, latent heat is released, which heats up the gas phase and leads to a temperature jump in the primary nozzle.

Figure 5 .
Figure 5. Axial distribution of Mach number: (a) dry steam and wet steam profiles; and (b) dry steam and wet steam contours.

Figure 5 .Figure 6 .
Figure 5. Axial distribution of Mach number: (a) dry steam and wet steam profiles; and (b) dry steam and wet steam contours.Buildings 2023, 13, x FOR PEER REVIEW 14 of 28

Figure 6 .
Figure 6.Axial distribution of static temperature: (a) dry steam and wet steam profiles; and (b) dry steam and wet steam contours.

Figure 7 .
Figure 7.The effect of key geometry parameters consists of (a) primary nozzle diameter and (b) area ratio on the entrainment ratio and the global COP.

Figure 7 .
Figure 7.The effect of key geometry parameters consists of (a) primary nozzle diameter and (b) area ratio on the entrainment ratio and the global COP.

Figure 8 .
Figure 8. Profiles of Mach number (a), droplet average radius (b), droplet nucleation rate (c), and liquid fraction (d) along the ejector center line with different primary nozzle diameter ratios.

Figure 8 .
Figure 8. Profiles of Mach number (a), droplet average radius (b), droplet nucleation rate (c), and liquid fraction (d) along the ejector center line with different primary nozzle diameter ratios.

Figure
Figure 9a depicts the Mach number along the ejector central line at AR = 4, 5, and 6.With the increase of the AR from 4 to 6, the Mach numbers maintain almost the same trend until the primary and secondary flows are mixed.The maximum value of the Mach number is unchanged, which means AR will not affect the primary flow.With the increase of AR, the Mach number after condensation shock decreases.The AR = 5 ejector reduces 60.46% of outlet Mach number (Ma = 0.43) as compared to the AR = 4 ejector (Ma = 0.69).Figure9billustrates the variation trend of the droplet average radius of each AR ejectors with axial coordinates.The average radius of droplets first increases, and then decreases; the maximum droplet average radius appears at approximately 36 mm in the primary nozzle.The rise in AR does not alter the average radius of droplets.In other words, droplet growth only exists in the primary nozzle, and droplet growth no longer continues when the flow enters the mixing section.Furthermore, before the diffuser outlet, the profiles of various AR ejectors show a constant trend.In fact, there are a lot of small droplets in this region.No matter how the AR changes, the droplet average radius remains unchanged, indicating the little impact of AR on the droplet average radius.Figure9cshows the profiles of droplet nucleation rates with different AR along the axis of the ejector.Once the nucleation rate reaches its peak, the droplet radius undergoes a rapid increase, compared to Figure10b.The nucleation occurs in two regions for the AR = 4 ejector.One is in the primary nozzle, and the supercooling impact of the divergent section leads to the supersonic flow reaching the equilibrium state.The other two nucleation processes occur in the mixing section.For the AR = 5 ejector, the nucleation occurs in three regions; the difference is that nucleation occurs again in the constant section.The ending position of the droplet nucleation rate moves forward with the increase of the AR.

Buildings 2023 ,
13, x FOR PEER REVIEW 18 of 28 number is unchanged, which means  will not affect the primary flow.With the increase of , the Mach number after condensation shock decreases.The  = 5 ejector reduces 60.46% of outlet Mach number ( = 0.43) as compared to the  = 4 ejector ( = 0.69).

Figure 9 .Figure 9 .
Figure 9. Profiles of Mach number (a), droplet average radius (b), droplet nucleation rate (c), and liquid fraction (d) along the ejector center line with different area ratios.

Figure 9 .Figure 10 .
Figure 9. Profiles of Mach number (a), droplet average radius (b), droplet nucleation rate (c), and liquid fraction (d) along the ejector center line with different area ratios.

Figure 10 .
Figure 10.Profiles and contours of Mach number, subcooled vapor temperature, droplet nucleation rate, and liquid mass fraction in the ejector (γ d = 2.4 and AR = 5).(a) Profiles of Mach number and subcooled vapor temperature.(b) Profiles of droplet nucleation rate and liquid mass fraction.(c) Contours of Mach number, subcooled vapor temperature, droplet nucleation rate, and liquid mass fraction.

Table 1 .
Effect of area ratio and nozzle exit position.

Table 2 .
Variable geometry ejector literature review.

Table 3 .
Classification of viewpoints in steam ejectors considering non-equilibrium condensation.

Table 4 .
Physical properties of the refrigerants (Critical temperature T c , Critical pressure P c , Boiling point T * [5][5].
T d is the droplet temperature; C p is the vapor isobaric specific heat capacity; λ v represents the vapor thermal conductivity; K n represents the Knudsen number; P r represents the Prandtl number; and ψ is the modeling correction coefficient: where