Thermodynamic Modelling of Supersonic Gas Ejector with Droplets

This study presents a thermodynamic model for determining the entrainment ratio and double choke limiting pressure of supersonic ejectors within the context of heat driven refrigeration cycles, with and without droplet injection, at the constant area section of the device. Input data include the inlet operating conditions and key geometry parameters (primary throat, mixing section and diffuser outlet diameter), whereas output information includes the ejector entrainment ratio, maximum double choke compression ratio, ejector efficiency, exergy efficiency and exergy destruction index. In single-phase operation, the ejector entrainment ratio and double choke limiting pressure are determined with a mean accuracy of 18% and 2.5%, respectively. In two-phase operation, the choked mass flow rate across convergent-divergent nozzles is estimated with a deviation of 10%. An analysis on the effect of droplet injection confirms the hypothesis that droplet injection reduces by 8% the pressure and Mach number jumps associated with shock waves occuring at the end of the constant area section. Nonetheless, other factors such as the mixing of the droplets with the main flow are introduced, resulting in an overall reduction by 11% of the ejector efficiency and by 15% of the exergy efficiency.


Introduction
The use of supersonic ejectors as an alternative to improve the performance and reduce the environmental impact of refrigeration systems has gained interest in recent years.Depending on the role of the device within the cycle, different configurations with varying performance and complexity have been proposed in the literature [1][2][3].An interesting choice is the Heat Driven Refrigeration Cycle (HDRC) shown in Figure 1, where the ejector completely substitutes the compressor.Under this scheme, the bottom loop (streams 1-2-6-7-8) is a refrigeration cycle with the compression work coming from the top loop (1-2-3-4-5).Although of modest performance, HDRCs can be operated by low-cost low-grade energy (e.g., industrial waste reject heat or solar thermal power), rendering a cheaper and less pollutant alternative to the standard vapor compression cycle.Moreover, such cycles can also operate with a new class of environmentally friendly refrigerants, named Hydrofluoroolefins (HFOs), without diminishing too much the performances of the system compared to R134a [4].HDRC systems with water as working fluid can operate at evaporator temperatures as low as 10 • C and generator temperatures in the range 90-110 • C, resulting in a Coefficient of Performance (COP, the ratio of cooling load to generator heat input) of 0.45 under laboratory conditions [5] and of 0.245 when powered by solar thermal power [6].A similar performance (COP of 0.3) can be attained with R134a at even lower temperatures: 3 • C at the evaporator and 90 • C at the generator [7,8].The heart of these novel refrigeration systems is the supersonic ejector, a simple device where the energy of the motive flow (stream 5 in Figure 1) is used to entrain and compress the secondary flow (stream 8).They are known for their simple design, low maintenance requirements and stable operation under single-and two-phase conditions [9].The typical geometry and main sections of a supersonic ejector assuming a constant area mixing are depicted in Figure 2.Under normal operation, the primary flow (also known as motive flow) is accelerated through the motive nozzle and enters the mixing section as a supersonic jet.At this point, shear interactions develop between both inlet flows until a uniform mixture is formed (position L5).Afterwards, a series of shock waves takes place such that the flow is subsonic at the beginning of the diffuser, where it compresses to outlet conditions.At fixed inlet conditions, two operating regimes can be recognized depending on the value of the outlet pressure P out relative to a certain threshold P lim [10].For P out < P lim , the secondary mass flow rate ṁsec is choked between the primary jet and the outer wall, such that total mass flow rate is at a maximum and independent of the exit pressure.This is known as a double-choke operation.Beyond P lim the ejector operates under the single-choke regime and ṁsec decays with increasing P out .At the limit of ṁsec = 0, the ejector is said to be in the malfunction regime.The ejector operating range is usually expressed in terms of the entrainment ratio (ω r , the ratio of the secondary to primary mass flow rates, Equation ( 1)) and the pressure ratio (P ratio , the ratio of the outlet to secondary pressures, Equation ( 2)).These parameters reflect, respectively, the mass drawing and compression capacity of the device: P ratio = P out P sec . ( It has been argued that the use of ω r and P ratio is limited when the goal is to compare the effectiveness of certain ejector alternatives in the context of a complete system.Therefore, three other parameters have been considered in this study to better reflect these differences: the ejector efficiency η Elbel as proposed by Elbel and Hrnjak [11], which gives insight on the efficiency of the ejector as a compression device; the exergy efficiency η χ , which compares the exergy output relative to the sum of the primary and secondary inlets and the ejector efficiency and the exergy destruction coefficient ξ i , which reflects the distribution of exergy destruction through the device [12].These parameters are defined in Section 4.7.The reader is referred to the review of Lawrence and Elbel [13] for a thorough analysis of different ejector performance definitions. The performance of ejector based refrigeration systems is a direct function of the entrainment ratio [4].A simple approach to relate cycle operating conditions to the ejector performance is the use of thermodynamic models, which offer a good compromise between accuracy and simplicity.For single-phase ejectors, the entrainment ratio is assumed to be a function of the Effective Area, i.e., the annular passage between the motive jet and the Constant Area Section (CAS) walls where the choking of the secondary flow occurs [14].The pioneering model applying this concept is the one proposed by Huang et al. [10], which calculates the double-choke entrainment ratio and CAS diameter of gas ejectors, given the inlet operating conditions and motive nozzle dimensions.The model assumes isentropic perfect gas behavior and divides the ejector in key regions: motive nozzle, secondary inlet, mixture before and after the shock and diffuser.Despite its simplicity, it presents an experimental deviation of up to 22.99% for a vast amount of data concerning a R141b test bench with 11 different motive nozzle designs at varying primary inlet temperatures (78 • C to 95 • C).It was later extended by Chen et al. [15] to predict the entrainment ratio decay rate under single-choke operation given P lim .Galanis and Sorin [16] proposed to determine the limiting pressure and all the ejector dimensions (given the entrainment ratio) by imposing the positive entropy generation constraint.Their model also introduces the use of polytropic efficiencies instead of the constant isentropic coefficients, which better reflect the pressure ratio variations happening in single-choke operation.The perfect-gas assumption was dropped in the model of Garcia del Valle et al. [17], which used a potential flow solution to determine the secondary flow passage and entrainment ratio of single-phase ejectors.Gas properties were calculated using the Helmholtz equation of state adjusted with experimental data.An average deviation of 7% relative to the experimental data of Huang et al. [10] was reported, with an important increase in complexity.
For two-phase systems, the applicability of the above mentioned models is limited given the complexity of determining the mixture speed of sound and the possibility of phase-change and metastable states.The former can be avoided by assuming that the choked mass flow per unit area, G = ρV, is a function of the local pressure and inlet conditions [18] (see Section 2.2), whereas the latter brings the need for some assumptions regarding the liquid-vapor interactions.The simplest choice is the Homogeneous Equilibrium Model (HEM), which assumes complete thermodynamic equilibrium between the phases.This model gives accurate results in the region close to and above the critical point [19].Using this approach, Ameur et al. [20] developed a model for determining the ejector flow properties in the context of Ejector Expansion Refrigeration Cycles, where the primary and secondary flows are, respectively, liquid and vapor and the entrainment ratio is indirectly fixed by the system through the quality at the outlet of the diffuser.Results showed a deviation of 0.21% to 7.14% in the choked flow rate prediction for CO 2 supersonic nozzles and of 0.63% to 6.14% in the compression ratio of a complete ejector when comparing to published experimental data.The reader can refer to the work of He et al. [21] for a more detailed review on the thermodynamic modeling of supersonic ejectors up to 2009.
To date, the general and internal flow characteristics of ejectors under common operating conditions have been extensively studied using both experimental and numerical approaches.
The long-term objective being to increase the efficiency of ejector systems, a first attempt is made here to quantify the effect of droplet injection before the onset of the shock train on the performance of ejectors for HDRC applications.The shock train and the mixing process are indeed the two main sources of exery losses within the ejector, being responsible for 40% of the total exergy losses as shown recently by Croquer et al. [22] for a single-phase supersonic ejector working with R134a.Although it is a common practice in the operation of these cycles to avoid the presence of liquid through the device [8], its effects have not been objectively discussed so far.The presence of droplets might act on two important sources of losses inside the device: the velocity mismatch at the mixing section [23] and the intensity of the shock train in the CAS [22].
In the mixing area, the idea is based on the premise that the high velocity difference between the motive jet and the entrained fluid at the mixing section hinders the entrainment of the secondary fluid [23].An alternative to reduce the motive flow velocity while maintaining its momentum input is the use of a droplet laden gas.This hypothesis has been tested on an air-driven ejector by augmenting the humidity of the primary inlet flow [23].Although there were no noticeable effects under double-choke operation, an increase of 13% to 98% in the single choke secondary mass flow rate was observed for inlet primary pressures within the range 107 to 446 kPa and inlet humidities between 4.3% and 11.2%.A similar behavior has been reported along with an increase of 5% in P lim for motive inlet humid fractions up to 1% by Hemidi et al. [24].Nonetheless, these results should be taken with caution since the ejector was studied as an independent unit in both cases and using an uncommon fluid pair (air-water) in refrigeration systems.For R134a, it has been observed that a reduction in the motive inlet superheat (T prim − T Sat prim ) from 11 • C to 0 • C leads to a reduction in ω r of about 8%, related to greater area occupied by the motive jet in the present of condensation [25].Although the ejector was a part of a HDRC, no results concerning the system performance or P lim were provided.
The idea of injecting droplets in the CAS of ejectors sparks from the observation that mists mitigate the propagation and intensity of explosions [26].Droplets deform and breakup at the encounter with a shock front [27], extracting energy from the gas and reducing its velocity and pressure increase rate [28].In the context of ejectors, the objective is to reduce the shock intensity in the CAS, which accounts for about 40% of the exergy destroyed through the device [22].
To this end, this investigation assesses for the first time the potential effects of injecting droplets in the CAS of a supersonic ejector.A thermodynamic model has been developed to determine the entrainment and pressure ratio given the inlet conditions, geometrical parameters (motive throat, CAS and diffuser outlet diameters) and injected droplet fractions.To the best of the authors' knowledge, it has never been considered in the literature.Moreover, the ejector efficiency η Elbel [11], the ejector exergy efficiency η χ and the exergy destruction index ξ i are presented in order to reflect the potential effects over the cycle performance of a HDRC system.The model assumptions, input/output information and fundamental equations are given in Section 2. The droplet breakup energy calculation procedure is described in Section 3. Section 4 depicts the calculation procedure for each section of the ejector as well as the performance parameters.The model is then extensively validated for single-and two-phase operations in Section 5.The effect of droplet injection on the performance of a typical ejector for refrigeration purposes is analyzed in Section 6, followed by concluding remarks in Section 7.

Numerical Modeling
The thermodynamic model determines the performance of the ejector with and without droplet injection by dividing the device into sections depicted in Figure 2 and solving the mass, momentum and energy conservation equations for each region.The required and resulting information of the model as well as the basic assumptions are described in the following subsections.

•
Thermodynamic state at both inlets: P prim , T prim , T sec and P sec .
• Diameters at the primary throat D t , constant area section D mix and the diffuser exit D L8 .
• Loss coefficients for the primary nozzle η prim , secondary inlet η sec , mixing section η mix and diffuser η di f f .
-Temperature of droplets at the injector, T inj .

Output Data
• Primary ṁprim and secondary ṁsec mass flow rates and entrainment ratio.
• Thermodynamic properties at every ejector section L i indicated in Figure 2.

Main Assumptions
• Flow is 1D and steady-state and its properties are uniform at each cross section L i .

•
The secondary flow throat (Effective Area) occurs at position L3.

•
The pressure at position L3 maximizes the secondary flow.

•
A normal shock occurs before the diffuser inlet, between positions L5 and L7.

•
Before the normal shock, both inlet flows and injected droplets are fully mixed.

•
Losses are represented using expansion and compression isentropic efficiencies.
Concerning the droplet injection, the following assumptions have been made:

•
The chemical component of the droplets is the same as the gas phase used in the primary and secondary inlets.

•
The droplets are in complete mechanical and thermal equilibrium with the gas phase (V slip = 0).
• Droplets are spherical and form a monodisperse phase.

•
Droplets are injected between positions L4 and L5, normal to the main flow.

•
Droplet breakup occurs right after injection and through the shock waves.

•
Coalescence and droplet deformation effects are neglected.

General Governing Equations
For each region, the conservation principles of mass, momentum and energy are applied.In the general case of a control volume with two inlet faces (in, 1 and in, 2) and one outlet (out), the conservation equations take the following form: Conservation of mass: Conservation of momentum: Conservation of energy: where ρ is the fluid density, V is the flow velocity, A is the cross sectional area, P is the pressure, ṁ is the mass flow rate and h is the specific enthalpy.Fluid properties are determined using the CoolProp equation library, which relies on the Helmholtz free energy formulation to provide thermodynamic and transport properties within a 1% accuracy for a wide range of fluids [29].The library covers conditions up to 30 MPa and 523 K for CO 2 and between 170 K and 455 K and up to 70 MPa for R134a.
In two-phase conditions, thermodynamic equilibrium is assumed at each cross section.The model defines any two-phase state by combining the section pressure with either the specific enthalpy h or the specific entropy s and the mixture quality x l : where b is any thermodynamical property other than P, T or ρ.The subscripts m, l and v refer to mixture and saturated liquid and vapor, respectively.The superscript sat refers to saturation conditions (evaluated at the section pressure).
Losses along the expansion (inlets) and compression (diffuser) stages of the ejector are represented using isentropic coefficients η expansion and η compression respectively for each section, whereas a mixing efficiency η mix is introduced to account for friction losses in the mixing chamber [16]: where the superscript is denotes properties evaluated for an isentropic process between the same start and end points.

Calculation of the Entrainment Ratio
The entrainment ratio is calculated by assuming that the secondary mass flow rate chokes at position L3 through the area available between the primary jet and the CAS walls (A L3,sec ) [14].The choked flow for a given set of inlet conditions is calculated by maximizing the mass flow per unit area, G [30]: Such an approach has been also successfully employed by Ameur et al. [20].It enables to avoid the use of the speed of sound, which remains an open question in the literature when dealing with two-phase flows.
Throughout an isentropic expansion with fixed inlet specific total enthalpy (h 0 ) and specific total entropy (s 0 ), the local velocity and density can be expressed in terms of the inlet conditions and the varying local pressure: Such that: The function G in Equation ( 11) exhibits a maximum G * , which occurs at sonic conditions.In two-phase homogeneous systems, most of the choked conditions locate over the saturation line [18,20].
The double choke primary and secondary mass flow rates are respectively: ṁprim = G * prim A L2 and ṁsec = G * sec A L3,sec .

Droplets Effects
Droplets injected into the CAS are expected to affect the main flow by extracting energy for breakup and by altering the main flow thermodynamic properties.It is assumed that breakup may occur at two locations: at the first encounter with the main flow (between locations L4-L5) and through the shock waves before the diffuser (L5-L7).The injected droplets are assumed to become part of the main flow mixture fast enough so there is thermodynamic equilibrium at sections L5 through L8.It alters the thermodynamic properties of the flow relative to the baseline case (without injection).The terms related to the breakup energy of a droplet are described in the following subsection.

Breakup Energy
The surface energy of a droplet can be defined as the product of its surface area S d by its surface tension σ [31].For a spherical droplet of diameter φ d , it is defined as: Figure 3 presents a schematics of the surface energy transition during breakup.The droplet breakup energy ∆E br is the difference in surface energy between the mother droplet E s,1 and a daughter E s,3 .Hence, for one droplet breaking into n 3 daughters, the droplet breakup energy writes: Note that n 3 has been eliminated from Equation ( 13) by applying mass conservation.Moreover, since φ 3 < φ 1 , then ∆E br > 0, meaning breakup always extracts energy from the main flow.The breakup mechanism and resulting droplet size depend on the counter balance between the surface tension and shear forces exerted by the surrounding gas, represented by the Weber (Equation ( 14)) number We: where ρ * g is the surrounding gas density and V S is the relative velocity, i.e., the difference between the gas phase velocity and the initial velocity of the droplets.This last one being relatively small at the injection point, it is assumed that V S is directly the gas phase velocity.
Three different breakup mechanisms can be recognized depending on We and the droplet Reynolds number (Re d = ρ * g V S φ d µ g ) [27].In particular, no breakup should be expected for local Weber numbers under a critical value We c , calculated using Equation ( 15) [32]: where Oh is the Ohnesorge number, relating the viscous to inertial and surface tension forces in a droplet: Using Equation ( 15), the final droplet diameter after breakup in Equation ( 13) is: The reader can refer to the review of Liao and Lucas [33] for more details about advanced theoretical models for droplet breakup.

Droplet Injection
At the point of injection, Equations ( 14)-( 17) are used to calculate the daughter droplet size diameter φ L5 .Droplet properties at the point of injection are determined using P L3 and the injection temperature.If the corresponding state is not subcooled liquid, saturated liquid at P L3 is assumed leading to an injection temperature equal to T = 263 K here.The droplet injection velocity stems from the injector diameter, the injected droplet fraction X inj and droplet density ρ inj .
For an injection flow rate ṁinj of droplets with constant diameter φ inj , the total breakup energy between sections L4 and L5 is: For the breakup terms across the shock, the quality of the mixture is determined first at position L5.If 0 < x L,L5 < 1, an analogous calculation to that at the injection point is carried out using the liquid mass in the mixture ( ṁLiquid = ṁ5 x L,L5 ) and the properties at L5 to determine the droplet sizes and total breakup energy ∆E br,L7 .

Ejector Calculation Procedure
The following thermodynamic model has been implemented within a Matlab environment.Starting from the primary nozzle, the flow properties at each cross section L i are determined in a advancing fashion following the general procedure shown in Figure 4. Details of the calculations at each section are given below.

Motive Nozzle
The motive nozzle flow rate is determined using the approach described in Section 2.2 for inlet properties P P0 and T P0 as follows: 1.
Repeat steps 2 through 5 reducing the guess At the end of this step, the primary mass flow rate and the properties at the motive throat are determined.

Secondary Inlet
The sonic conditions corresponding to G * sec are calculated using a procedure analogous to that described in Section 4.1 but using P S0 and T S0 .By assuming that the secondary choking pressure is uniform across L3 for both flows, the area occupied by the primary jet is calculated according to Equation ( 19): where h L3,prim is determined using h P0 , s P0 and η prim , and ρ L3,prim is determined using h L3,prim and P L3 .The effective area, or minimal cross section available for the passage of the secondary flow, A L3,sec , is: Therefore, one gets:

Mixing (L3-L4)
The complete mixture of both inlet flows takes place between sections L3 and L4.Equations ( 3)-( 5) are applied in this region to determine the mixture velocity V L4 , pressure P L4 and enthalpy h L4 using the following steps: 1.

5.
If and go back to step 1. Else, the calculation is finished and flow properties at position L4 are known.

Droplet Injection (L4-L5)
Droplets are injected before the onset of the normal shock wave.Given the high relative velocity between both phases, droplets breakup shortly after injection.By inserting the droplet related terms into the conservation balances between sections L4 and L5, the following system is obtained: It is assumed here that the droplet velocity V droplets is equal to the gas velocity just after breakup.The simultaneous solution of Equations ( 22)-( 24) gives the mixture properties just before the normal shock.

Normal Shock (L5-L7)
Between sections L5 and L7, it is assumed that a normal shock wave takes place, which creates a jump in the flow conditions (P L5 < P L7 , s L5 < s L7 , ρ L5 < ρ L7 and V L5 > V L7 ) and induces a second breakup step.The cross sectional area is assumed to be constant between both locations, and equal to A L3 .
Using Equation( 3), the velocity V L7 can be eliminated from Equations ( 4) and ( 5), giving: The flow characteristics after the shock are determined as follows: 1.

3.
Determine and go back to step 1. Else, the calculation is finished and flow properties at position L7 are known.

Diffuser (L7-L8)
The ejector outlet pressure P L8 is determined using the conditions at L7, the diffuser isentropic coefficient and the cross section at L8. Since the area is constant between L6 and L7, it is assumed that changes in properties between both locations are negligible.The iterative procedure described below is applied to determine the flow properties at the ejector outlet: Calculate h is L8 by applying Equation ( 5) between sections L7 and L8.

3.
Determine and go back to step 1.

6.
Correct h L8 with P L8 and η di f f , and determine flow properties at L8.
Once h L8 and P L8 are known, the rest of flow properties at L8 can be determined.The resulting outlet pressure corresponds to the maximum compression ratio achievable by the ejector at given operating conditions [16].

Ejector Performance Parameters
Apart from ω r and P ratio , the ejector efficiency, the exergy efficiency and the exergy destruction index have been determined for better comparing the ejector performance at varying operating conditions and droplet injection characteristics.

Ejector Efficiency
The ejector efficiency proposed by Elbel and Hrnjak [11] compares the isentropic work required in the secondary flow compression to the work liberated through the isentropic expansion of the motive flow: where h is,sec out ( resp.h is,prim out ) is the specific enthalpy at the ejector outlet following the isentropic secondary flow compression (resp.isentropic primary flow expansion) and h sec and h prim are, respectively, the specific enthalpies at the secondary and primary inlet.Note that the denominator in Equation ( 27) is the work difference existing between the isentropic and isenthalpic expansion of the motive flow.

Exergy Efficiency
From the perspective of the 2nd Law of Thermodynamics, the exergy efficiency allows to analyze the performance of the ejector by comparing the output exergy flow rate of the ejector χout to the total of the primary χprim and secondary χsec inlets: with the exergy at the i-th boundary χi coming from Equation ( 29): where the subscript 0 denotes the reference or dead state.Given that the main purpose of the ejector is to provide exergy to the secondary flow, conditions at the secondary inlet are chosen for the dead state in all following calculations.

Exergy Destruction Index
The exergy destruction index ξ i compares the exergy destroyed within a specific region of the ejector Ḋi to the total exergy destroyed through the device Ḋej [12]: where Ṡ = ṁs is the entropy flow rate at a boundary.ξ i allows for comparing the contribution to the total exergy destruction of each section.

Validation of the Thermodynamic Model
The model was validated against an extensive amount of published experimental data for various operating conditions and working fluids typical of single-and two-phase ejector based refrigeration cycles.Validation points for single-phase operation include: the entrainment ratio, limiting pressure and thermodynamical properties at different sections of the ejector.Concerning two-phase operation, comparisons were made regarding the choked mass flow through a convergent divergent nozzle and the effect of the primary inlet saturation on the entrainment ratio.

Single-Phase Ejector
Figure 5 compares the predicted entrainment ratio values against experimental data from different studies [8,10,24,[34][35][36] with varying dimensions, working gases and operating conditions.As a first step, friction and mixing losses are neglected here such that: The majority of predicted values lie within 10% of the corresponding experimental point when ω r is in the range 0.35-0.6,which is a typical range for most ejector applications in refrigeration systems.Deviations beyond 20% are found at the lower range of ω r values, often associated with high motive to outlet pressure ratios.No clear tendencies are found between the model deviation and the entrainment ratio or the working gas.

Influence of the Loss Coefficients
Figure 6 demonstrates the effect of the loss coefficients on the entrainment ratio ω r and limiting pressure P lim for a constant area ejector with gaseous R134a as working fluid.Temperatures remain in the range 79-89 • C and 5-10 • C for the primary and secondary inlets, respectively, and a 10 • C superheat is applied to ensure vapor conditions [8].The loss coefficients are shown in Table 1.These were determined in a CFD study using the same ejector and operating conditions by performing energy balances through the motive nozzle, secondary inlet, mixing section and diffuser [22].Including losses in the thermodynamic model has a negligible effect on the predicted entrainment ratio, since it is only affected by the primary and secondary expansion coefficients, which are close to 1 in the present case.In terms of the limiting pressure, a better agreement is observed when using the adjusted loss coefficients η mix and η di f f , showing the greater influence of the mixing and diffusion processes on the ejector compressing capacity.Overall, the model determines P lim with an accuracy of 10%.
The reader can refer to the works of Liu and Groll [37] and Liu [38] for detailed reviews about the influence of the loss coefficients on the performance of thermodynamic models.1. Experimental data after [8].Table 2 shows a comparison between the thermodynamic model and CFD averaged values of pressure, temperature, Mach number Ma, specific enthalpy and specific entropy at four locations along the ejector: L2, L4, L7 and L8.Operating conditions correspond to Operating Point OP2 of Table 1, with T sec = 20 • C, P sec = 414.6 kPa and a fixed outlet pressure of 826.57kPa.The CFD model is described in detail in Croquer et al. [39].An average deviation of 3% is achieved at the primary throat (L2) and after mixing (L4).The greatest differences are observed at the start of the diffuser (L7): 15%, 18.53% and 19.48% for P, T and Ma, respectively.These discrepancies may be partly explained by the 2D nature of the flow at this position highlighted by the CFD model [39].In terms of global quantities, the primary and secondary mass flow rates according to the thermodynamic model (resp.CFD model) are ṁprim = 0.03753 kg/s and ṁsec = 0.0143 kg/s (resp.ṁprim = 0.03747 kg/s and ṁsec = 0.0163 kg/s), resulting in a difference of 0.15% and 14.3%, respectively, for the primary and secondary flow rates between both models.

Two-Phase Ejectors
To assess the capability of the model in handling two-phase flows, the primary mass flow rate of CO 2 through supersonic nozzles is compared with experimental values [40,41] at various inlet conditions indicated on the P-h diagram of Figure 7. Results are presented in Figure 8 for different values of the primary nozzle efficiency η prim , showing that the adjustment of η prim improves the model deviation from beyond 20% down to less than 10%.The best agreement is obtained with η prim = 0.75 for the data in Figure 8a and with η prim = 0.85 for the data in Figure 8b.
The capacity of the present model to reproduce the effects of the primary superheat at the primary inlet (∆ OH = T prim − T prim,sat ) on the entrainment ratio is shown in Figure 9 in comparison with the experimental values of Little and Garimella [25] for R134a.An average deviation of 4% is achieved.The linear dependency of ω r with the superheat level is also well reproduced.For a 0 • C superheat, the calculated quality at the primary nozzle throat (resp.position L3) is 0.942 (resp.0.933).

Effect of Droplet Injection on the Ejector Performance
The effect of droplet injection is assessed by modeling the case of a supersonic ejector with R134a as working fluid in the context of a HDRC.Primary inlet conditions are shown in Table 1 with P sec = 414.6 kPa and T sec = 20 • C. Four injection fraction values are considered: X inj = 1%, 2%, 5% and 10%.For all cases, droplets of diameter φ inj = 500 µm (average value measured by Chauvin et al. [28] from a commercial atomizer) are injected at a temperature of T inj = −13 • C to ensure liquid phase at the typical CAS pressure values encountered.No variations are expected on the entrainment ratio due to the droplet injection at the CAS, since the model determines ω r based solely on the inlet conditions and the throat areas A L2 and A L3,sec (see Section 2.2).

Changes to Pressure and Mach Profiles
The effect of droplet injection on pressure, Mach number Ma and temperature inside the ejector is shown in Figure 10, between sections L4 (right before injection) and L8 (diffuser outlet), for working conditions OP2 and injected droplet fractions between 1% and 10% of the primary mass flow rate without injection.To compute the Mach number, it is noticeable that the speed of sound is provided by the Coolprop library and not by any relation for two-phase flows. of Table 1.
Regarding the effect of droplet injection on pressure (Figure 10a), no changes are observed at sections L4 and L5, given the fact that the injection occurs in the CAS after full mixing of the primary and secondary flows.At position L7, a pressure reduction is observed with increasing injection fraction, which suggests an effect on shock intensity.Moreover, the lower pressure at the start of the diffuser is carried on to the outlet conditions affecting the compression ratio.Concerning Ma values (Figure 10b), an increase in X inj leads to lower Ma values at position L5.An opposite behavior is observed at position L7, where the Ma number slightly increases with X inj .At the diffuser outlet, the Ma number is ∼0.03 for all cases.These Ma profiles reflect that the flow in general decelerates with the injection of droplets, although the changes are barely noticeable.The most drastic impact of droplet injection is observed in the temperature profiles (Figure 10c).For an injected fraction of 10%, the temperature after injection reduces by 12 • C, and the across-shock temperature jump reduces by 60% in comparison with the case X inj = 0.A similar temperature reduction is observed at the outlet.The saturation temperature at L8 remains close to 33 • C for all the assessed conditions, meaning that increasing droplet injection reduced the superheat at the outlet from 20 • C to 5 • C.This behavior was identical for working conditions OP1 and OP3.

Shock Intensity
The intensity of the shock wave towards the end of the CAS is assessed in terms of the pressure ∆P and Mach number ∆Ma jumps between sections L5 and L7. Figure 11 shows the effect of the droplet injection on these quantities for the conditions OP2.The pressure increase has been normalized by the value obtained in the case without any injection, ∆P o .Both quantities, ∆P and ∆Ma, decrease proportionally with the injected mass fraction.At the maximum injection fraction considered (X inj = 10%) a reduction of ∼8% is observed for both quantities.Results for R134a at the inlet conditions OP2 of Table 1.
These results show that the injection of droplets effectively attenuates the shock wave intensity in the CAS, thus affecting one of the main sources of irreversibilities inside the ejector [12,22].At the maximum injected fraction considered, the entropy generation between sections L5 and L7 is reduced by 10% relative to the case without injection.

Limiting Pressure
The effect of the droplet injection fraction X inj on the limiting pressure for the three operating conditions OP1, OP2 and OP3 is shown in Figure 12.The limiting compression ratio diminishes proportionally with the increasing injection fraction, from 1% at X inj = 1% to 5% at X inj = 10%.This reduction in the outlet pressure is related to the lower pressure observed at the diffuser inlet with increased injection fraction (Figure 10a).It is important to point out that, although this behavior was assessed for fixed secondary inlet conditions, there is no reason to believe the tendency would be any different when varying the conditions of the entrained flow.
Figure 12.Effect of the droplet injection fraction X inj on the maximum compression ratio achievable by the ejector.Results for R134a at the inlet conditions OP1, OP2 and OP3 of Table 1.

Ejector Efficiency and Exergy Performance
Figure 13 shows the variations in the ejector efficiency and exergy efficiency with injection fractions of 1% to 10% at the conditions OP1, OP2 and OP3 and for P sec = 414.6 kPa and T sec = 20 • C. Both the exergy efficiency and the ejector efficiency reduce with increasing X inj , independently of the operating conditions.Regarding η Elbel , an average performance loss of 11% is observed for the greatest injection fraction, reflecting the decrease in the outlet pressure shown in Figure 10.For lower outlet pressures at fixed inlet conditions, the secondary flow compression ratio diminishes while the primary flow expands even more, resulting in less energy recovered by the ejector.1.
Concerning the exergy efficiency, a reduction of 15% in average is observed when the injection mass fraction is 10%, confirming that droplet injection results in less exergy recovered relative to the total entering at both inlets.The location where this exergy is destroyed can be pinpointed using the exergy destruction index as shown in Figure 14.Although droplets effectively attenuate the normal shocks in the CAS and reduce their contribution to exergy destruction by about 35% for X inj = 10%, other processes associated with the droplet injection (e.g., mixing with the main flow) present an even greater contribution to the total exergy destroyed.For example, in the limiting case of X inj = 10%, 3.9 kJ/kg of exergy are destroyed at the point of injection, which is six times the exergy destroyed across the shock wave in the case without injection.Thus, although the expected shock attenuation effect is carried out, other factors damaging the effectiveness of the ejector are introduced, resulting in a lower performance.From the perspective of a HDRC, a lower performance would be expected given that the injection of droplets renders the ejector less effective for energy recovery and compression.

Conclusions
In this paper, a thermodynamic model has been developed to predict the entrainment ratio and compression ratio of supersonic gas ejectors with droplets for HDRC applications, given the inlet operating conditions and key geometry parameters (primary throat, mixing section and diffuser outlet diameter).It has been extensively validated against published experimental or numerical data for single-and two-phase flows.The model has then been extended to assess the effect of injecting droplets at the constant area section of the device over different performance parameters: ejector efficiency, exergy destruction and maximum compression ratio.The following conclusions can be drawn:

•
Under single-phase conditions and fixing all losses coefficients to 1, the model determines the entrainment ratio with an average deviation of 18% for various working fluids (Air, R141b, R245fa and R134a) and the double-choke limiting pressure for R134a with a mean accuracy of 4.5%.The inclusion of adjusted losses coefficients has a negligible effect on the determined entrainment ratio but reduces the deviation in P lim to 2.5%.

•
Under two-phase conditions, the model presents a deviation of about 10% in the choked mass flow of CO 2 across a convergent-divergent nozzle.In this case, the use of loss coefficients greatly improves the results.The effect of the primary inlet superheat over the entrainment ratio is also accurately reproduced.

•
An analysis of a R134a ejector under typical operating conditions of a HDRC shows that the injection of droplets results in a lower ejector performance.Although at an injection fraction of 10%, the shock intensity reduces by 8%, the ejector efficiency and exergy efficiency reduce by 11% and 15%, respectively.Exergy destruction profiles show that the gains achieved by the shock attenuation are overcome by the entropy generated by the droplet injection and its mixing with the main flow.

Figure 1 .
Figure 1.(a) Flow process diagram and (b) pressure-enthalpy diagram of a single-phase heat driven refrigeration cycle.The dotted red lines represent the flow inside the ejector.

Figure 2 .
Figure 2. Typical geometry of a constant area ejector.

Figure 3 .
Figure 3. Variation of the droplet surface energy during breakup.

Figure 4 .
Figure 4. General calculation procedure of the thermodynamic model.Detailed steps for each region are given in Section 4.

Figure 6 .
Figure 6.Comparison between the experimental and predicted values of (a) the entrainment ratio ω r and (b) limiting pressure P lim for a single-phase R134a ejector.The operating conditions and loss coefficient values are shown in Table1.Experimental data after[8].

Figure 7 .Figure 8 .
Figure 7. Location on a P-h diagram of the experimental inlet conditions used in the CO 2 choked mass flow rate model validation.Data sources: [40,41].

Figure 9 .
Figure 9.Comparison between the experimental and predicted values of the entrainment ratio ω r for different superheat values.Experimental data after [25].

Figure 10 .
Figure 10.Effect of droplet injection fraction X inj on the (a) pressure; (b) Mach number Ma and (c) temperature at different sections of the ejector.Results for R134a at the inlet conditions OP2 of Table1.

Figure 11 .
Figure 11.Effect of the droplet injection fraction X inj on the (a) pressure and (b) Mach number Ma change across the shock wave, relative to the case without droplets.Results for R134a at the inlet conditions OP2 of Table1.

Figure 13 .
Figure13.Effect of the droplet injection fraction X inj on the ejector efficiency η Elbel and exergy efficiency η χ .Results for R134a at the inlet conditions OP1, OP2 and OP3 of Table1.

Figure 14 .
Figure14.Effect of the droplet injection fraction X inj on the contribution of each section to the exergy destruction across the ejector.Results for R134a at the inlet conditions OP2 of Table1.

Table 1 .
[22]es of the loss coefficients deduced from the CFD analysis of[22]for different operating conditions.