Exergoeconomic Analysis of a Variable Area Solar Ejector Refrigeration System under Hot Climatic Conditions

: The present study investigates low-grade heat utilization in ejector refrigeration systems under hot climatic conditions. A variable area ejector is used to maximize the harvested heat from the generator of the solar system at peak times. Exergy, economic, and exergoeconomic analyses are conducted to evaluate the performance of the system. A thermodynamic model of the system has been developed using Ebsilon Professional software. Available experimental and theoretical data validate the results. The effects of properties of the working ﬂuids, ejector geometry, and operation conditions are also evaluated. It was found that the coefﬁcient of performance of the system reached 0.45 at a generator pressure of 3 bars. Furthermore, it was noticed that the overall exergy efﬁciency could be increased for a ﬁxed generator temperature while increasing the ejector area ratio. A value of 21% exergetic efﬁciency was calculated for the system. The exergoeconomic analysis of the system demonstrated that heat exchangers are required to be improved thermodynamically at the expense of the capital investment cost.


Introduction
Recent years have been plagued by environmental disasters that can be directly credited to global warming effects.Many nations worldwide are making tremendous efforts toward greener and more efficient energy sources.One of the promising avenues is solar energy.Solar energy is abundant and reliable, hence its popularity [1,2].In countries located in arid climates, it has been found that one of the highest energy consumers is the residential air conditioning demand.In such areas, there is a high potential to utilize solar energy, leading to high energy savings.Research studies have been conducted to analyze the use of solar energy for refrigeration purposes [3,4].Ejector refrigeration machines have a long history of application [5][6][7].
Africa, Asia, and several parts of other continents are characterized by hot and dry climatic conditions.The average daytime temperature is usually above 30 • C, which increases the need for air conditioning [8].The ejector refrigeration systems (ERS) could be implemented in the air conditioning application.This would be advantageous because not only would it reduce reliance on the national grid, but it would also help with environmental concerns [9].An ejector is also called a "jet compressor".Ejector refrigeration systems beyond the thermally-driven refrigeration systems, i.e., that utilize low-grade heat and are attractive when driven by solar energy.[10].They can be incorporated into various system configurations and combined with power cycles [11].Ejectors are reliable because they have no moving parts.They also have relatively low capital costs and operating and maintenance expenses.Ejector refrigeration systems have much lower performance when compared to vapor-compressor systems because of their thermodynamic cycles [12].
Ejectors can be grouped based on their nozzle position, design, and number of phases.The constant pressure mixing (CPM) ejectors are more applicable than constant area mixing Energies 2022, 15, 9540 2 of 19 ejectors due to their ability to work against higher back pressures and their superior performance [13].Constant area mixing (CAM) corresponds to higher mass flow rates.A recent paper discusses combining CPM and CAM ejectors to create a constant rate of momentum change.The constant rate of momentum-change (CRMC) configuration type uses a variable area section to optimize the flow passage area, reducing the chances of thermodynamic shocks and improving ejector performance.The geometry of the nozzle can affect how the ejector operates [14,15].The nozzle geometry can be convergent, which means that the ejector can work at subsonic conditions and achieve the maximum sonic velocity at the exit of the suction.The flow through the ejector can reach supersonic speeds [16].
Ejectors are influenced by back pressure during sub-critical mode.As the pressure in the mixing chamber rises, a shockwave enters the process, interacting with the mixing and raising the pressure even higher.When the vacuum is turned off, the primary flow in the suction chamber reverses.The flow through an ejector can be either one or two phases [17].Two-phase ejectors can be classified by flow type: condensing ejectors, where the primary flow condenses in the ejector, and two-phase ejectors, where the flow is two-phase at the outlet [18].
The performance of the variable geometry ejector can be predicted by two 1-D mathematical models.The so-called Huang's model is a model that is used because it considers the ejector's critical point of operation.This model is applicable for variable geometry ejectors [19,20].Another model, Chen's model, is used as a comparison tool to study ejector operation in either the subcritical or critical.In both cases, there is some hypothetical effective throat at which the speed of the secondary flow reaches its maximum value.Just after this section of the throat, the mixing process of the primary and secondary streams begins at a uniform pressure [21].
A study of a solar ejector refrigeration system (SERS) was carried out for the Mediterranean Sea region, and the authors found that the overall efficiency had a lower performance at the solar radiation peak [22].As the flow increases, the ejector will be choked, and the harvested heat in the generator will not be fully utilized.The thermal efficiency of the SERS is decreased as the solar radiation approaches its maximum value.The ejector flow caused a stoppage, which prevented the generator from using all the solar heat it could [23].Recent research innovations have helped improve system performance, including the use of special refrigerants, the use of renewable energy, and reduced mechanical pump work [24].
Exergy and exergoeconomic analyses were carried out to study the performance of the ejector refrigeration system.An exergoeconomic analysis was conducted for series and combined double effect ejector systems to evaluate the effect of key operating conditions on the investment and product costs.It was found that the implementation of ejector to the other thermodynamic cycles has an economic advantage over the double effect system [25][26][27].An exergy study of the Kalina power-cooling cycle with an ejector was conducted to optimize the system using exergy and pinch analyses [28].The frictional losses on the ejector's wall were investigated by developing a two-dimensional model for the ejector.MATLAB software [29] was used to solve the exergy, energy, and exergoeconomic equations of the model, and a multi-objective optimization was carried out [30,31].
So-called 4E analysis was conducted for an integrated process with capturing carbon dioxide and storage, in combination with organic Rankine and absorption refrigeration cycles [32].Furthermore, an experimental analysis of a solar-assisted heat pump with energy, exergy, economic and exergoeconomic analysis was carried out [33].
In addition, it should be mentioned that hot climatic conditions significantly influence the exergetic and economic characteristics of any heat-driven refrigeration system [34].
Considering the above literature, although ejector refrigeration systems have been around for some time, many improvements to the original design have been researched.The operating principle of the ejector refrigeration systems states that the ejector's ability to expel refrigerant is limited, and a maximum flow rate is permitted, especially when the heat source is dynamic (as in a solar-driven system).The present work addresses a Energies 2022, 15, 9540 3 of 19 solar-driven ejector refrigeration system with an adjustable area ratio to cope with the variation in the generator operating conditions.The ejector as a component has been modeled in MATLAB, and simulations of the system in EbsilonProfessional.The energy, exergy, economic, exergoeconomic analyses, and a sensitivity analysis of the system were carried out to find the optimum operating conditions.The aim is to determine the optimal combination of the variable ejector dimensions.Attention has been given to applying solar ejector refrigeration systems under hot climatic conditions.

Materials and Methods
An ejector refrigeration system consists of a generator, evaporator, condenser, throttling valve, pump, and ejector.The schematic of the evaluated system is shown in Figure 1a.The saturated vapor from the generator (state 1) is the primary stream at a low velocity; it enters the converging-diverging nozzle, where it is accelerated.At the exit from the ejector (state 3), the supersonic accelerated refrigerant flow exits at low pressure while creating a suction effect, and the secondary stream (low-pressure refrigerant from the evaporator) is entrained (state 2).In the critical operation mode of the ejector, a shear layer is formed due to the gradient of velocity between the secondary and primary flows.The secondary stream is accelerated to sonic speed.The mixing process begins when the flow is obstructed.The interactions of the primary and secondary streams with the ejector wall add to the mixing process complexity.After the mixing process, a shock wave could form, either in the constant area chamber or at the diffuser inlet.The operational conditions determine the location of the shock wave.Normal shock causes the flow velocity to change from supersonic to subsonic, resulting in a significant increase in pressure.Finally, in the diffuser section of the ejector, the mixed flow stream pressure rises and exits at the condensation pressure.
After the condenser (state 4), the refrigerant splits into two parts.The evaporation of the refrigerant provides the target refrigeration effect (process 8-2) and leaves the evaporator as a saturated vapor.
Considering the above literature, although ejector refrigeration systems have been around for some time, many improvements to the original design have been researched.The operating principle of the ejector refrigeration systems states that the ejector's ability to expel refrigerant is limited, and a maximum flow rate is permitted, especially when the heat source is dynamic (as in a solar-driven system).The present work addresses a solardriven ejector refrigeration system with an adjustable area ratio to cope with the variation in the generator operating conditions.The ejector as a component has been modeled in MATLAB, and simulations of the system in EbsilonProfessional.The energy, exergy, economic, exergoeconomic analyses, and a sensitivity analysis of the system were carried out to find the optimum operating conditions.The aim is to determine the optimal combination of the variable ejector dimensions.Attention has been given to applying solar ejector refrigeration systems under hot climatic conditions.

Materials and Methods
An ejector refrigeration system consists of a generator, evaporator, condenser, throttling valve, pump, and ejector.The schematic of the evaluated system is shown in Figure 1a.The saturated vapor from the generator (state 1) is the primary stream at a low velocity; it enters the converging-diverging nozzle, where it is accelerated.At the exit from the ejector (state 3), the supersonic accelerated refrigerant flow exits at low pressure while creating a suction effect, and the secondary stream (low-pressure refrigerant from the evaporator) is entrained (state 2).In the critical operation mode of the ejector, a shear layer is formed due to the gradient of velocity between the secondary and primary flows.The secondary stream is accelerated to sonic speed.The mixing process begins when the flow is obstructed.The interactions of the primary and secondary streams with the ejector wall add to the mixing process complexity.After the mixing process, a shock wave could form, either in the constant area chamber or at the diffuser inlet.The operational conditions determine the location of the shock wave.Normal shock causes the flow velocity to change from supersonic to subsonic, resulting in a significant increase in pressure.Finally, in the diffuser section of the ejector, the mixed flow stream pressure rises and exits at the condensation pressure.
After the condenser (state 4), the refrigerant splits into two parts.The evaporation of the refrigerant provides the target refrigeration effect (process 8-2) and leaves the evaporator as a saturated vapor.In the present study, we perform the simulation based on an ejector with variable area ratios.In Figure 1b, a schematic diagram of the variable area ejector is presented.The variable area ejector consists of five parts: primary nozzle, suction chamber, mixing chamber, diffuser, and a spindle, which is driven by a motor.With its movement, the effective flow area of the primary nozzle changes accordingly.By this, the primary flow rate of the ejector can be adjusted, which in turn leads to a variation in the secondary flow rate.This model will accommodate the variation in the solar radiation during the day and allows the harvesting of more energy and increases the performance of the system.Figure 2 represents a screenshot of the proposed system in the Ebsilon Professional software.The selection of an appropriate refrigerant is a critical issue for ejector refrigeration systems.The ecological properties of the refrigerant (i.e., ozone depletion potential (ODP), In the present study, we perform the simulation based on an ejector with variable area ratios.In Figure 1b, a schematic diagram of the variable area ejector is presented.The variable area ejector consists of five parts: primary nozzle, suction chamber, mixing chamber, diffuser, and a spindle, which is driven by a motor.With its movement, the effective flow area of the primary nozzle changes accordingly.By this, the primary flow rate of the ejector can be adjusted, which in turn leads to a variation in the secondary flow rate.This model will accommodate the variation in the solar radiation during the day and allows the harvesting of more energy and increases the performance of the system.Figure 2 represents a screenshot of the proposed system in the Ebsilon Professional software.In the present study, we perform the simulation based on an ejector with variable area ratios.In Figure 1b, a schematic diagram of the variable area ejector is presented.The variable area ejector consists of five parts: primary nozzle, suction chamber, mixing chamber, diffuser, and a spindle, which is driven by a motor.With its movement, the effective flow area of the primary nozzle changes accordingly.By this, the primary flow rate of the ejector can be adjusted, which in turn leads to a variation in the secondary flow rate.This model will accommodate the variation in the solar radiation during the day and allows the harvesting of more energy and increases the performance of the system.Figure 2 represents a screenshot of the proposed system in the Ebsilon Professional software.The selection of an appropriate refrigerant is a critical issue for ejector refrigeration systems.The ecological properties of the refrigerant (i.e., ozone depletion potential (ODP), preferable if the fluid is non-corrosive, has low toxicity, is non-explosive, chemically stable, inexpensive, and commercially available.
In the present design, the refrigerant R141b has been chosen.R141b is a pure hydrofluorocarbon (HFC) working fluid with no chlorine atoms, i.e., ODP = 0 and GWP is negligibly small.It is considered an appropriate refrigerant for ejector refrigeration systems due to its high performance in terms of coefficient of performance (COP) and entrainment ratio compared with other refrigerants [35,36].

System Modeling and Evaluation
For the simulation, the following has been assumed: steady state condition; 0.01 bar pressure drop in pipes with negligible heat losses in the generator and evaporator; the cooling air temperature of the condenser is 5 • C higher than the ambient temperature; the refrigerant at the evaporator exit is saturated vapor.In addition, an adiabatic ejector is assumed with a 1D steady flow.The changes in kinetic energy at the inlet and exit of the ejector are neglected, and ejector efficiencies are used to calculate the losses in the ejector sections.A constant pressure mixing process is assumed in the mixing chamber.The thermodynamic properties are calculated using Ebsilon Professional.

Ejector Modelling
An ejector 1D mathematical model is used to evaluate and predict its performance.The used governing equations are mass, momentum, and energy for a compressible flow.The Mach number and pressure at the exit of the primary nozzle are then calculated assuming isentropic flow using the isentropic flow equations.Then, assuming that the secondary flow is choked, the pressure of the secondary flow at that section.The energy balance is applied to the mixing process to estimate the mixing temperature, which is then used to find the mixing Mach number.
where, v m is the mixed-flow speed, R g is the ideal gas constant, γ is the gas specific ratios, T m is the mixed flow temperature.
The pressure and the Mach number at the diffuser exit are calculated after determining the mixing Mach number at the constant pressure mixing process.The pressure at the exit can be calculated as where P c is the pressure of the constant mixing process, P 2 is the pressure at the diffuser exit.
Figure 3 depicts a flowchart of the process calculation and the equations used in the energy analysis are presented in the Supplementary Materials.
The ejector efficiency definition used is based on ASHRAE, which relates the actual recovered energy of compression to the theoretical maximum energy available in the motive stream.
m g is the mass flow rate in the generator, . m e is the mass flow rate at the evaporator, h e,in is the enthalpy at the inlet of the condenser, h e,out is the enthalpy at the exit of the evaporator, h g,out is the enthalpy at the generator exit.The ejector specifications are presented in Table 1.

System Performance
The entrainment ratio, ω, the compression ratio, R c , and the COP are parameters that describe the ejector refrigeration system where m s is the secondary mass flow rate and m p is the primary mass flow rate.
Q e is the evaporator cooling capacity, .
Q g is the generator heat load, and W p is the pump work.

Exergy Analysis
Exergy analysis provides a more accurate and robust analysis of the thermodynamic system compared to energy analysis.The total exergy of the system consists of chemical, physical, potential, and kinetic exergy.Only physical exergy will be considered in this study.
The exergy analysis has been conducted using the "fuel ( .E F )/product ( .E P )" approach [37]: The subscripts F, P, D, and L indicate fuel, product, destruction, and losses, respectively.

•
On the kth component level .
The evaluation is based on: • Exergetic efficiency of the entire system • Exergetic efficiency of the kth component • The exergy destruction ratio

Economic Analysis
The total revenue requirement (TRR) method is used to examine the economic analysis of the system.To perform the TRR, the system's total capital investment (TCI) must first be estimated based on the cost of each purchased item (PEC), after which the parameters for the detailed cost calculation's economic, financial, operating, and market inputs are determined.Finally, the succession of expenses will be levelized into a constant number with a comparable financial value (annuity).Estimating the PEC of each component with a precise method is the most important and difficult portion of the economic analysis; this is particularly tough for new and unproven technologies.The range of cost estimation errors is typically between 10% and 30%.Given that the suggested system is innovative, the cost estimation may be less accurate than expected.To estimate system costs and further support system performance evaluation, the TRR approach is still a viable alternative as the foundation for the exergoeconomic study.
The TRR consists of carrying charges (CC L ) and the expenses of the fuel (FC L ) and operation and maintenance (OMC L )costs: The CC L is the capital investment (including total capital recovery, preferred stock, return on investment, income taxes, and insurance).This means that it is the value of the total capital investment cost (TCI), which is composed of the fixed capital investment (FCI) and the interest accrued: The FCI is obtained by adding the bare module cost (BMC) to the service facilities and contingencies funds.The FCI represents the total system cost at time = 0; that is before construction.The indirect system costs include construction costs, contingencies, administrative fees, and engineering.
To calculate the TCI, we need to estimate the purchase equipment cost (PEC) using the following approach: where C PE,re f is the approximate equipment costs with the size X new , C PE,known is the known equipment costs with the corresponding size X known , and α is the size exponent.The third step is to adjust the estimated equipment cost to the reference year.The reference year is 2019: where ref is the year, the equipment is to be purchased and old is the year the cost of the equipment is known.After calculating the current equipment purchase cost, the equipment's nature and characteristics must be considered in the form of factors.Material and pressure correction factors (MPF) are defined for unique materials, high pressure, designs, and materials.
The final step is the module factor (MF) cost.This accounts for labor, piping, instrumentation, and everything necessary at the installation stage.The bare module cost (BMC) can be calculated using: The total direct costs are the sum of BMC, while the indirect costs are calculated as a percentage of the total direct costs.
The constant escalation levelization factor (CELF) and the capital recovery factor (CRF) are calculated as where i e f f is the effective interest rate, and n is the economic lifetime.FC L and OMC L are determined as With where r FC is the average inflation rate of fuel and r OMC is the operation and maintenance cost.FC 0 is the fuel cost in the year of operation.

Exergoeconomic Analysis
This combined exergy and economic analysis provides more comprehensive information that is unavailable through a conventional energy, exergy, and economic analysis.Exergoeconomic analysis is based on cost balances on the component level [37], i.e., . .
where the average cost of fuel c F = .E D,k .Exergoeconomic factor f k is used in the optimization procedure to make decisions of either investing in a more efficient components to reduce the exergy destruction or to sacrifice efficiency to decrease the costs associated with the carrying charges:

Results and Discussion
The mathematical model has been solved by EBSILON Professional Software, and the system key properties and parameters are calculated.The simulation results of the proposed system are reported in Table 2.
The mathematical model has been validated by already published data in the literature [38].Comparisons between the obtained results and those published in the literature for evaporator temperatures of 12 • C are presented in Table 3.The depicted good agreement between the obtained results and the data published in the literature gives confidence in the present model and its results.

Sensitivity Energy Analysis
Various area ratios were studied to define the optimum temperature in the generator.The temperature and pressure at the ejector inlet nozzle have high values and produce a suction effect causing more secondary refrigerant flow to be entrained as the generator temperature is increased.This is accompanied by an increase in COP, reaching a maximum value of 0.209 at the optimum generator temperature of 73 • C. When the generator temperature continues to rise, the COP begins to fall because there is no corresponding increase in primary mass flow due to choking conditions in critical operating mode at the ejector nozzle.As a result, a constant refrigerant mass flow prevails, which keeps the cooling capacity constant with increasing the heat input.Thus, the COP decreases.The driving pressure difference between the evaporator and the nozzle exit increases as T g increases, which increases both COP and entrainment ratio.As T g rises, so does the critical back pressure.This increase causes the thermodynamic shock wave to move to the mixing chamber, preventing clogging of the secondary flow.As the back pressure increases, the flow through the nozzle becomes slower, and the secondary stream stops because the ejector stops working properly, and the direction of the primary flow will be reversed back to the evaporator.At the downstream of the diffuser section, a thermodynamic shock wave occurs, and thus the secondary pressure rises, and as the evaporator pressure rises, so does the critical back pressure.Figure 4 depicts the effect of evaporator temperature on the performance in T g = 73 • C and T a = 38 • C operating conditions.The system COP and driving pressure difference increase rapidly as T e increases.

Sensitivity Energy Analysis
Various area ratios were studied to define the optimum temperature in the generator.The temperature and pressure at the ejector inlet nozzle have high values and produce a suction effect causing more secondary refrigerant flow to be entrained as the generator temperature is increased.This is accompanied by an increase in COP, reaching a maximum value of 0.209 at the optimum generator temperature of 73 °C.When the generator temperature continues to rise, the COP begins to fall because there is no corresponding increase in primary mass flow due to choking conditions in critical operating mode at the ejector nozzle.As a result, a constant refrigerant mass flow prevails, which keeps the cooling capacity constant with increasing the heat input.Thus, the COP decreases.The driving pressure difference between the evaporator and the nozzle exit increases as Tg increases, which increases both COP and entrainment ratio.As Tg rises, so does the critical back pressure.This increase causes the thermodynamic shock wave to move to the mixing chamber, preventing clogging of the secondary flow.As the back pressure increases, the flow through the nozzle becomes slower, and the secondary stream stops because the ejector stops working properly, and the direction of the primary flow will be reversed back to the evaporator.At the downstream of the diffuser section, a thermodynamic shock wave occurs, and thus the secondary pressure rises, and as the evaporator pressure rises, so does the critical back pressure.Figure 4 depicts the effect of evaporator temperature on the performance in Tg = 73 °C and Ta = 38 °C operating conditions.The system COP and driving pressure difference increase rapidly as Te increases.The ejector's area ratio rises as the generator temperature increases, as shown in Figure 5.This is because as the generator temperature rises, so does the generator pressure, resulting in enhanced ejector entrainment at constant condenser and evaporator temperatures.More secondary flow can be entrained into the ejector once the area ratio is increased.The area ratio increases can be doubled for every 50 • C increase in the generator temperature.
The optimum T g for a given area ratio can be estimated if the temperature and pressure of the condenser and evaporator are known.To ensure that the ejector operates under critical choking conditions, the condenser temperature, T c , is set to 38 • C. Since the ejector back pressure is equal to the pressure of condensation, the condenser temperature affects the ejector's mode of operation.As the generator temperature approaches its optimum value, the system COP will be maximum COP.As the spindle position changes during time, the primary pressure and mass flow rate change accordingly, as shown in Figure 6.The refrigerant mass flow rate increases with the increase in the spindle movement, allowing more refrigerant to pass at reduced pressure.This will lead to an increase in the COP.
The primary operating parameters selected for this analysis are the generator pressure, and evaporator and generator temperatures.Figure 7 shows the relationship between the generator temperature and the entrainment ratio.It is noted that as the generator temperature increases, the entrainment ratio increases as well.The reason for this increase is due to the increase in the primary flow temperature, which increases its energy and leads to an increase in the flow speed in the suction chamber, and more secondary flow is entrained.
resulting in enhanced ejector entrainment at constant condenser and evaporator temperatures.More secondary flow can be entrained into the ejector once the area ratio is increased.The area ratio increases can be doubled for every 50 °C increase in the generator temperature.The optimum Tg for a given area ratio can be estimated if the temperature and pressure of the condenser and evaporator are known.To ensure that the ejector operates under critical choking conditions, the condenser temperature, Tc, is set to 38 °C.Since the ejector back pressure is equal to the pressure of condensation, the condenser temperature affects the ejector's mode of operation.As the generator temperature approaches its optimum value, the system COP will be maximum COP.As the spindle position changes during time, the primary pressure and mass flow rate change accordingly, as shown in Figure 6.The refrigerant mass flow rate increases with the increase in the spindle movement, allowing more refrigerant to pass at reduced pressure.This will lead to an increase in the COP.The primary operating parameters selected for this analysis are the generator pressure, and evaporator and generator temperatures.Figure 7 shows the relationship between the generator temperature and the entrainment ratio.It is noted that as the  The optimum Tg for a given area ratio can be estimated if the temperature and pressure of the condenser and evaporator are known.To ensure that the ejector operates under critical choking conditions, the condenser temperature, Tc, is set to 38 °C.Since the ejector back pressure is equal to the pressure of condensation, the condenser temperature affects the ejector's mode of operation.As the generator temperature approaches its optimum value, the system COP will be maximum COP.As the spindle position changes during time, the primary pressure and mass flow rate change accordingly, as shown in Figure 6.The refrigerant mass flow rate increases with the increase in the spindle movement, allowing more refrigerant to pass at reduced pressure.This will lead to an increase in the COP.The primary operating parameters selected for this analysis are the generator pressure, and evaporator and generator temperatures.Figure 7 shows the relationship between the generator temperature and the entrainment ratio.It is noted that as the generator temperature increases, the entrainment ratio increases as well.The reason for this increase is due to the increase in the primary flow temperature, which increases its energy and leads to an increase in the flow speed in the suction chamber, and more secondary flow is entrained.Figure 8 depicts the COP values in relation to generator temperature at various generator pressures.The heat delivered from the solar field causes the generator temperature to rise and the COP to decrease.As Tg increases, the entrainment ratio will increase, and because of the increased rate of heat addition, the COP decreases.Figure 8 depicts the COP values in relation to generator temperature at various generator pressures.The heat delivered from the solar field causes the generator temperature to rise and the COP to decrease.As T g increases, the entrainment ratio will increase, and because of the increased rate of heat addition, the COP decreases.Figure 8 depicts the COP values in relation to generator temperature at various generator pressures.The heat delivered from the solar field causes the generator temperature to rise and the COP to decrease.As Tg increases, the entrainment ratio will increase, and because of the increased rate of heat addition, the COP decreases.As can be seen from Figure 9, the entrainment ratio increases as the temperature of the evaporator increases.The entrained flow's kinetic energy increases, which lead to an increase in the suction efficiency of the ejector.At a predetermined generator pressure, the primary flow rate is choked.As can be seen from Figure 9, the entrainment ratio increases as the temperature of the evaporator increases.The entrained flow's kinetic energy increases, which lead to an increase in the suction efficiency of the ejector.At a predetermined generator pressure, the primary flow rate is choked.Figure 10 shows that as the evaporator temperature rises, so does the COP of the ejector cooling cycle.It is also clear that the range of the COP increases as the generator pressure decreases.Figure 10 shows that as the evaporator temperature rises, so does the COP of the ejector cooling cycle.It is also clear that the range of the COP increases as the generator pressure decreases.Figure 10 shows that as the evaporator temperature rises, so does the COP of the ejector cooling cycle.It is also clear that the range of the COP increases as the generator pressure decreases.

Exergy Analysis
The data of the specific physical exergies (Table 2) were used to calculate the value of the exergy destruction within system components via applying the exergy balances (Equation ( 8)).These data and corresponding exergy destruction ratios (Equation ( 10)) are shown in Figure 11a as well as the exergy efficiencies (Equation ( 9)) are shown in Figure 11b.The definitions of the fuel and product for the solar collectors can be found in [39,40], and for the ejector in [40].
The results demonstrate that the highest exergy destruction and associated low exergetic efficiencies are associated with the heat exchangers.

Exergy Analysis
The data of the specific physical exergies (Table 2) were used to calculate the value of the exergy destruction within system components via applying the exergy balances (Equation ( 8)).These data and corresponding exergy destruction ratios (Equation ( 10)) are shown in Figure 11a as well as the exergy efficiencies (Equation ( 9)) are shown in Figure 11b.The definitions of the fuel and product for the solar collectors can be found in [39,40], and for the ejector in [40].The effect of the autonomous variable geometry ejector on the system performance as the generator temperature changes is reported in Figure 12.The higher the temperature within the generator, the higher the performance of the entire system is observed.
From Figure 12 it can be noticed that the overall exergy efficiency can be increased for a fixed generator temperature.Up to a certain point, increasing the ejector area ratio increases efficiency before drastically decreasing again.The significant reduction can be attributed to flow breakdown, as a mixture of primary and secondary streams would be very turbulent at larger area openings.To improve system performance, an ejector refrigeration system can be automated so that the nozzle opening can be regulated using a temperature sensor connected to the generator.The results demonstrate that the highest exergy destruction and associated low exergetic efficiencies are associated with the heat exchangers.
The effect of the autonomous variable geometry ejector on the system performance as the generator temperature changes is reported in Figure 12.The higher the temperature within the generator, the higher the performance of the entire system is observed.
From Figure 12 it can be noticed that the overall exergy efficiency can be increased for a fixed generator temperature.Up to a certain point, increasing the ejector area ratio increases efficiency before drastically decreasing again.The significant reduction can be attributed to flow breakdown, as a mixture of primary and secondary streams would be very turbulent at larger area openings.To improve system performance, an ejector refrigeration system can be automated so that the nozzle opening can be regulated using a temperature sensor connected to the generator.
Energies 2022, 15,9540 From Figure 12 it can be noticed that the overall exergy efficiency can be increased for a fixed generator temperature.Up to a certain point, increasing the ejector area ratio increases efficiency before drastically decreasing again.The significant reduction can be attributed to flow breakdown, as a mixture of primary and secondary streams would be very turbulent at larger area openings.To improve system performance, an ejector refrigeration system can be automated so that the nozzle opening can be regulated using a temperature sensor connected to the generator.Figure 13 shows the variation in the system's exergetic efficiency as a function of the ejector area ratio for different generator temperatures.It can be seen that the highest exergetic efficiency could be obtained at the optimum value of the generator temperature.Yet, the maximum efficiencies are observed at an area ratio of 0.75.Figure 13 shows the variation in the system's exergetic efficiency as a function of the ejector area ratio for different generator temperatures.It can be seen that the highest exergetic efficiency could be obtained at the optimum value of the generator temperature.Yet, the maximum efficiencies are observed at an area ratio of 0.75.

Economic Analysis
For conducting the economic analysis, the following parameters were assumed: Parameters assumed for the economic analysis are plant life ()-15 years; effective interest rate ( )-10%; average general inflation rate for fuel (rFC) is zero because of solar energy; average general inflation rate for OMC (rOMC) 2.5%, and total annual operation time −7500 hr./a.The calculated levelized total capital investment is 2251.34$/h, and the levelized cost of the product is 0.009 $/kWh.Figure 14 shows the distribution of the total capital investment.The evaporator, generator, air cooler, and solar system accounted for 90% of

Economic Analysis
For conducting the economic analysis, the following parameters were assumed: Energies 2022, 15, 9540 Parameters assumed for the economic analysis are plant life (n)-15 years; effective interest rate (i e f f )-10%; average general inflation rate for fuel (r FC ) is zero because of solar energy; average general inflation rate for OMC (r OMC ) 2.5%, and total annual operation time −7500 hr./a.The calculated levelized total capital investment is 2251.34$/h, and the levelized cost of the product is 0.009 $/kWh.Figure 14 shows the distribution of the total capital investment.The evaporator, generator, air cooler, and solar system accounted for 90% of the system cost, whereas the ejector, pumps, and motors accounted for 10%.

Economic Analysis
For conducting the economic analysis, the following parameters were assumed: Parameters assumed for the economic analysis are plant life ()-15 years; effective interest rate ( )-10%; average general inflation rate for fuel (rFC) is zero because of solar energy; average general inflation rate for OMC (rOMC) 2.5%, and total annual operation time −7500 hr./a.The calculated levelized total capital investment is 2251.34$/h, and the levelized cost of the product is 0.009 $/kWh.Figure 14 shows the distribution of the total capital investment.The evaporator, generator, air cooler, and solar system accounted for 90% of the system cost, whereas the ejector, pumps, and motors for 10%.

Exergoeconomic Analysis
The main results obtained from the exergoeconomic analysis are reported in Figures 15 and 16

Exergoeconomic Analysis
The main results obtained from the exergoeconomic analysis are reported in Figures 15 and 16

Conclusions
The ejector refrigeration cycle is subjected to thermodynamic modeling, simulation, energy, exergy, economic and exergoeconomic analyses, in addition to various sensitivity studies in order to useful information on design specifications.The entrainment ratio, spindle position, and COP were used for the system performance analysis, while the effect of the area ratio on system performance was investigated.The COP and entrainment ratio are two complex functions of inlet temperature, back pressure, and flow path geometry.It is also worth noting that different spindle positions can yield higher COP values due to differences in operating conditions.This demonstrates the advantage of using

Conclusions
The ejector refrigeration cycle is subjected to thermodynamic modeling, simulation, energy, exergy, economic and exergoeconomic analyses, in addition to various sensitivity studies in order to provide useful information on design specifications.The entrainment ratio, spindle position, and COP were used for the system performance analysis, while the effect of the area ratio on system performance was investigated.The COP and entrainment ratio are two complex functions of inlet temperature, back pressure, and flow path geometry.It is also worth noting that different spindle positions can yield higher COP values due to differences in operating conditions.This demonstrates the advantage of using variable geometry over fixed geometry.It is especially important in the study case because solar energy is a volatile and unpredictable source.As a result, the variable geometry design can maintain COP and cooling capacities even when the generator temperature varies.It was found that the system's coefficient of performance increased with the decrease in the generator pressure and increase in the evaporator temperature.It reached a value of 0.45 at a generator pressure of 3 bars and an evaporator temperature of 20 • C. Furthermore, it was noticed that the overall exergy efficiency increased as the generator temperature increased and approached a constant value of 21% for a generator temperature of 85 • C, while increasing the ejector area ratio.A value of 21% exergetic efficiency was recorded for the system.The exergoeconomic analysis of the system demonstrated that heat exchangers require improvement thermodynamically at the expense of the capital investment cost.The economic analysis revealed that the calculated levelized total capital investment is 2251.34$/h, and the levelized cost of the product is 0.009 $/kWh.

Figure 1 .
Figure 1.(a) Schematic of the evaluated system.(b) Schematic of the variable are ejector.

Figure 1 .
Figure 1.(a).Schematic of the evaluated system.(b).Schematic of the variable are ejector.

Figure 2 .
Figure 2. Ebsilon Professional screenshot.The selection of an appropriate refrigerant is a critical issue for ejector refrigeration systems.The ecological properties of the refrigerant (i.e., ozone depletion potential (ODP), and global warming potential (GWP)) should also be taken into account.In addition, it is
; and the exergy destruction cost rates within the component .C D,k = c F .E D,k , and the entire system .C D,tot = cF tot ∑ .

Figure 4 .Figure 4 .
Figure 4. COP vs. ejector secondary flow temperature at Tg = 73 °C and Tc = 38 °C operating conditions.The ejector's area ratio rises as the generator temperature increases, as shown in Figure 5.This is because as the generator temperature rises, so does the generator pressure,

Figure 6 .
Figure 6.Ejector pressure and mass flow vs. spindle position.

Figure 6 .
Figure 6.Ejector pressure and mass flow vs. spindle position.

Figure 6 .
Figure 6.Ejector pressure and mass flow vs. spindle position.

Figure 7 .
Figure 7. Entrainment ration vs. generator temperature at T e = 5 • C and T a = 38 • C operating conditions.

Figure 8 .
Figure 8. COP vs. generator temperature at T e = 5 • C and T c = 38 • C operating conditions.

Figure 9 .
Figure 9. Entrainment ratio vs. evaporator temperature at T g = 73 • C and T a = 38 • C operating conditions.

Figure 10 .
Figure 10.COP vs. evaporator temperature at T g = 73 • C and T a = 38 • C operating conditions.

Figure 12 .
Figure 12.Comparison of entrainment ratio and exergy efficiency vs. generator temperature.

Figure 12 .
Figure 12.Comparison of entrainment ratio and exergy efficiency vs. generator temperature.

Figure 13 .
Figure 13.Exergy efficiency vs. ejector area ratio at T e = 5 • C and T a = 38 • C operating conditions.

Figure 14 .Figure 14 .
Figure 14.Distribution of the total capital investment.
. The values Z k and C D,k for each system component and corresponding exergoeconomic factors (f k ) demonstrate that the heat exchangers require improvement thermodynamically with the expenses of the capital investment cost, i.e., the values of Z k should be increased and C D,k decreased.Energies 2022, 15, x FOR PEER REVIEW 17 of 19 . The values Zk and CD,k for each system component and corresponding exergoeconomic factors (fk) demonstrate that the heat exchangers require improvement thermodynamically with the expenses of the capital investment cost, i.e., the values of Zk should be increased and CD,k decreased.

Figure 15 .
Figure 15.Comparison of Z_k and C_D,k for each component.

Figure 15 .
Figure 15.Comparison of Z_k and C_D,k for each component.

Figure 15 .
Figure 15.Comparison of Z_k and C_D,k for each component.

Table 1 .
Ejector and entire system specifications.

Table 2 .
Simulation results of the proposed system.