Energy and Exergy Analysis on Zeotropic Refrigerants R-455A and R-463A as Alternatives for R-744 in Automotive Air-Conditioning System (AACs)

: The popularity of vehicles and the increased time spent in cars with air conditioning systems has led to regulations in many countries that require the use of environmentally friendly refrigerants with minimal global warming and zero ozone depletion potential (GWP and ODP). Cars need high-performance, eco-friendly air conditioning systems to reduce their impact on the environment, lower fuel consumption, and decrease carbon emissions. The aim of the current work was to propose CO 2 -based blend zeotropic refrigerants, R-455A (R-744/32/1234yf) and R-463A (R-744/32/125/1234yf/134a), to improve the thermodynamic performance of pure CO 2 refrigerants. The thermodynamic energy and exergy analysis and system optimization of an AAC system for the new zeotropic refrigerant blends compared to carbon dioxide (R-744), using Aspen HYSYS software, were investigated. The inﬂuence of cooler/condenser pressure, average evaporator temperature, cooler/condenser outlet temperature, and refrigerant ﬂow rate on the cycles’ COP and exergy efﬁciency were conducted and are presented. The results showed that, at the same operating condition parameters, the cycle COP improved by 57.6 and 76.5% when using R455A and R463A instated of R744, respectively, with the advantage of reducing leakage problems due to the higher operating pressure of R744 (5–7 times higher than those of R455A and R463A), as well as requiring heavy equipment, but at optimal operating condition parameters, R744 and R-463A had a maximum COP of 14.58 and 14.19, respectively. The maximum COPs of R744, R455A, and R463A based on the optimal pressure of the cooler/condenser were 3.1, 4.25, and 5.4, respectively. Additionally, regarding the need for environmentally friendly air conditioning systems with acceptable performance in cars due to their impact on the environment and their contribution to global warming, the blend R455A is recommended for use as a refrigerant in AAC systems.


Introduction
Most automobiles now come equipped with automotive air conditioning (AAC) systems due to their prevalence and the amount of time people spend in their vehicles.These systems typically use refrigerant fluids, which can be harmful to the environment and contribute to global warming.To mitigate this, many countries have mandated the use of eco-friendly refrigerants with minimal global warming potential (GWP) and zero ozone depletion potential (ODP).The popularity of cars and the increasing time that people spend in them have led to the installation of air conditioning systems in most vehicles.However, the use of refrigerant fluids in these systems can be damaging to the environment, prompting many countries to enforce laws requiring the use of environmentally friendly refrigerants with zero ozone depletion potential and low global warming potential.
Due to scientific advancements and the detrimental effects of non-natural refrigerants on the ozone layer and greenhouse gas emissions, it is inevitable that current refrigerants will be phased out and replaced with new, ecologically friendly alternatives.As a result, R134a, with its 100-year GWP of 1430 [1], has been banned from air conditioning systems in the European Union since 2017 [2], being substituted with refrigerants that have a GWP of under 150 [3].The Kyoto Protocol's focus on the impact of certain greenhouse gases on global warming has led to the gradual replacement of traditional hydrofluorocarbons.Natural CO 2 (R744) has gained attention as a safe refrigerant due to its non-flammability and non-toxic properties [4].
Research on refrigerant replacements is now divided into two key approaches that are driven by environmental concerns.The first approach involves seeking innovative refrigerants with low GWP and zero ODP, such as United Signal Company's azeotropic combination R410A and DuPont Company's R1234yf and R1234ze.The second approach explores the use of natural refrigerants [5,6].CO 2 has long been used as a refrigerant by humans [7,8], offering several advantages, such as non-toxicity, environmental friendliness, safety, affordability, and excellent thermal properties.Studies by Robinson and Groll [9] from Purdue University have demonstrated that the CO 2 transcritical cycle is a cutting-edge and reliable technique for automobile air conditioning, with great potential for efficiency enhancement.Research by Yu Binbin, Wang Dandong, Xiang Wei et al. [10] has shown that CO 2 electric car air conditioning systems perform comparably to R134a.
Further studies have investigated the impact of R32 refrigerant on CO 2 water heat pump systems, revealing varying reductions in exhaust pressure when R32 was applied [11].Massuchetto LHP et al. conducted tests on three different refrigerant combinations (R744/R1270, R744/R717, and R744/RE170), determining that R744/RE170 yielded the optimal refrigeration outcome [12].Mancini et al. explored the use of CO 2 and dimethyl ether as an azeotropic combination refrigerant, finding that the addition of dimethyl ether benefited the CO 2 system [13].Adding propane to CO 2 refrigerant was also found to enhance the system's coefficient of performance by 7-9% according to the study by W. I. Mazyan et al. [14].
While CO 2 refrigerant is a popular research area for refrigerant replacement, it does have some disadvantages, including its high working pressure and relatively lower system efficiency.These factors somewhat hinder its rapid promotion and implementation [15,16].Sun et al. conducted an experimental investigation on CO 2 /R32 blends in a water-to-water heat pump system [17].Yu et al. studied the temperature difference in heat transfer using theoretical and experimental methods based on the nonlinear temperature enthalpy of R236fa/R32 mixtures, identifying distinct characteristics of the glide temperature depending on the composition [18].Liu et al. examined the detonation of R32/R1234ze(E) mixtures, analyzing the flame structure, explosion pressure, and lowest flammability boundaries of the refrigerant explosion [19].
The development of synthetic refrigerants such as R-32, R-125, R-1234yf, and R-134a with zero ODP and low GWP is hindered by safety concerns, processing complexity, and high cost [20,21].R744, an alternative to R134a, belongs to the Hydrofluoroolefin family [22].Abbood et al. [23] studied and presented environmentally friendly alternative refrigerants to R134a and compared the thermodynamic performance of an AAC system using R134a with blends of hydrocarbons (HCs) such as R290/R600a.Their theoretical analysis demonstrated that the blend (R600a/R290/134a) with a mass ratio of (43/35/22) exhibited a low global warming potential (GWP), while maintaining comparable refrigerant performance.This blend can be directly utilized in the system without requiring any modifications.A reviewed, compared, and comprehensively described outline of the main enhancements of recent advances and sustainable solutions in automobile air conditioning systems category was presented by Sagar and Rakshit [24].The thermodynamic and thermophysical properties for new mixed refrigerants, R13I1/R152a, R1234yf/R290, and R1234yf/R600a, as alternatives to R134a in automotive air conditioning systems, to solve the problem of the high global warming potential of R134a, were presented by [25,26].Savitha et al. [27] presented a review of the thermodynamic and flammability properties of the low GWP refrigerants and the compatibility of these refrigerants with construction material and lubricant.A novel design that integrated evaporative cooling with an automotive CO 2 air conditioning system was introduced and experimentally investigated by [28] to address the significant decrease in energy efficiency experienced in hot climates when using carbon dioxide (CO 2 ) refrigerant for vehicles.Lei et al. [29] presented thermotical and experimental studies to increase the heat transfer area of the evaporator and to improve the system performance of air condition systems using CO 2 in hot climate conditions.Vaccaro et al. [30] investigated and presented the effects of the addition of the second element in the CO 2 transcritical cycles to overcome the heavy expansion loss, requiring specific means for its mitigation.Luo et al. [31] theoretically studied the vapor-liquid equilibrium (VLE) properties of two eco-friendly zeotropics, CO 2 /R1234ze(Z) and CO 2 /R1336mzz(E), by using Peng-Robinson (PR) and Soave-Redlich-Kwong (SRK) models.Chen et al. [32] proposed new correlations, along with the application of deep learning-based modeling techniques, to analyze and predict saturated flow boiling heat transfer and two-phase pressure drops in the evaporating flow.Hussain et al. [33] predicted and optimized the two-phase pressure drop of refrigerant R1234yf across a diverse range of testing conditions.To achieve this, feature engineering techniques were employed to identify and select the most influential features that could accurately estimate the desired output.
However, pure R744 systems can experience leakage problems due to their higher operating pressures (7-10 times higher than conventional R134a systems) and require heavy equipment, which reduces their coefficient of performance.These issues can be mitigated by using new blends, although the properties of such blends can differ from their original constituents.Therefore, the aim of this study was to explore the potential of two R744-based blends, R-455A (R-744/32/1234yf) and R-463A (R-744/32/125/1234yf/134a), to enhance the thermodynamic performance of pure R744 refrigerants, leading to lower engine fuel consumption and hence lower carbon emissions.The need for environmentally friendly air conditioning systems with high performance in cars, due to their impact on the environment and contribution to global warming, is essentially required.To achieve this, the thermodynamic energy and exergy characteristics of an AAC system were evaluated and compared for the proposed zeotropic refrigerants and for carbon dioxide (R-744) using Aspen HYSYS software.Additionally, the effects of cooler/condenser pressure, average evaporator temperature, cooler/condenser outlet temperature, and refrigerant flow rate on the cycle coefficient of performance (COP) and exergy efficiency were studied, analyzed, and presented.The anticipated outcomes of this investigation aim to demonstrate the most advantageous refrigerant performance in comparison to R-744 and to provide suitable operating conditions for the proposed refrigerants.Furthermore, these findings will serve as a foundation for constructing an AAC system that can be mass-produced economically and in an environmentally conscious manner.    2 displays the physical and environmental properties of the studied zeotropic refrigerants.The composition of the two blends, R-455A (R-744/32/1234yf) and R-463A (R-744/32/125/1234yf/134a), by mass are (3.0/21.5/75.5)and (6.0/36.0/30.0/14.0/14.0),respectively.Furthermore, the blend R-463A has a high GWP of approximately 1386, while the GWP of the blend R-455A is less than one.

Temperature Glide, P-T Envelope
During the evaporation process, the refrigerant undergoes a transition from a saturated liquid state, known as the bubble point, to a saturated vapor state, known as the dew point.The temperature at which this transition occurs remains constant under constant pressure.However, when multiple components are present in a refrigerant mixture, a temperature glide can occur.The temperature glide refers to the difference between the dew point and bubble point temperatures.In this study, the temperature glide properties of R455A and R463A were investigated by analyzing the vapor-liquid equilibrium (VLE).The Peng-Robinson equation of state (EOS) was employed to accurately calculate the saturated vapor pressure and saturated liquid density.However, it is important to use appropriate mixing rules when evaluating the thermodynamic properties of the mixture using the EOS. Figure 1 illustrates the pressure-temperature (P-T) envelopes of the studied zeotropic refrigerants (R455A and R463A) compared to pure R744.Pure R744 does not exhibit a temperature difference between the bubble and dew point temperatures since it is a single component refrigerant (see dashed lines in Figure 1a).Azeotropic blends behave similarly to pure refrigerants, where the boiling and condensation points coincide at the same composition for both vapor and liquid phases.However, for zeotropic blends, the more volatile component starts boiling first in the evaporator, followed by the less volatile component.This sequential boiling can lead to a change in concentration as the temperature changes, resulting in a temperature glide (t glide > 2 K).If the refrigerant begins boiling at t i and ends at t f , then the temperature glide can be expressed as t glide = t f − t i (see dashed lines in Figure 1b,c).The changing composition of the evaporating liquid, caused by the varying boiling point, contributes to the temperature glide.Additionally, leakage from the system can introduce changes in the composition and properties of the refrigerant, further contributing to t glide .

System Modeling
A simulation system was constructed using Aspen HYSYS V12.1 ® Software (As-penTech, Bedford, MA, USA) [35] to assess the thermodynamic performance of zeotropic alternative refrigerants in AACs, as depicted in Figure 2. Aspen HYSYS has gained recognition among academics and engineers for its reliability and capability to evaluate complex industrial processes.The user-friendly nature of the Aspen HYSYS platform enables the optimization of conceptual design and operational parameters.This robust process simulator offers a vast library of pre-built component models and property packages.Aspen HYSYS allows for the seamless integration of multiple modules, facilitating the connection of material and energy streams.This capability enables the static and dynamic modeling of various complex chemical and hydrocarbon fluid-based processes.The simulation model developed for AACs in Aspen HYSYS offers easy integration with a range of energy systems, including compressors, condensers, gas coolers, heat exchangers, evaporators, expansion valves, and more.
By leveraging Aspen HYSYS, the thermodynamic performance of zeotropic alternative refrigerants in AACs can be thoroughly examined.The software's extensive functionality and versatility contribute to the advancement of energy-efficient and environmentally conscious air conditioning technologies.

System Modeling
A simulation system was constructed using Aspen HYSYS V12.1 ® Software (Aspen-Tech, Bedford, MA, USA) [35] to assess the thermodynamic performance of zeotropic alternative refrigerants in AACs, as depicted in Figure 2. Aspen HYSYS has gained recognition among academics and engineers for its reliability and capability to evaluate complex industrial processes.The user-friendly nature of the Aspen HYSYS platform enables the optimization of conceptual design and operational parameters.This robust process simulator offers a vast library of pre-built component models and property packages.Aspen HYSYS allows for the seamless integration of multiple modules, facilitating the connection of material and energy streams.This capability enables the static and dynamic modeling of various complex chemical and hydrocarbon fluid-based processes.The simulation model developed for AACs in Aspen HYSYS offers easy integration with a range of energy systems, including compressors, condensers, gas coolers, heat exchangers, evaporators, expansion valves, and more.
By leveraging Aspen HYSYS, the thermodynamic performance of zeotropic alternative refrigerants in AACs can be thoroughly examined.The software's extensive functionality and versatility contribute to the advancement of energy-efficient and environmentally conscious air conditioning technologies.

System Description and Assumptions
It is crucial to investigate and evaluate the performance of AACs cycles considering the physical and environmental characteristics of the newly studied zeotropic alternative refrigerants.Figure 2a illustrates the schematic diagram of the traditional vapor compression refrigeration cycle in AACs, comprising an evaporator, compressor, cooler/condenser, expansion valve, and internal heat exchanger.The following is a description of the system cycle: (i) 1′-2: non-isentropic compression; (ii) 2-3: constant pressure heat rejection from the cooler/condenser; (iii) 3-3′: subcooling in the internal heat exchanger (IHE); (iv) 3′-4: isenthalpic throttling in the expansion valve (EXV); (v) 4-1: constant pressure heat absorption in the evaporator; (vi) 1-1′: superheating in the internal heat exchanger (IHE).By adjusting the average evaporation temperature, cooler/condenser refrigerant output temperature, cooler/condenser pressure, and refrigerant mass flow rate, the operating conditions of the cycle can be modified.Because the critical point pressure of the refrigerant can be comparatively low compared to the cooler pressure, a transcritical model for the cycle was developed based on the operating conditions.In transcritical cycles, a gas cooler is utilized instead of a condenser, as in subcritical cycles, since there is no phase change above the critical point.As the pressure and temperature of the refrigerant are independent of each other, both pressure and temperature need to be defined at pressures above the critical point.The coefficient of performance of the cycle varies at the same gas cooler outlet temperature but different operating pressures.Therefore, the simulation model's gas cooler was optimized by running it at the optimum gas cooler pressure [36].Figure 2b presents the pressure-enthalpy (p-h) diagram of subcritical/transcritical cycles using R-744 in AACs.
In order to ensure accuracy and simplicity in the current simulation, the following assumptions were considered to evaluate the performance of zeotropic alternative refrigerants in AACs:  The entire system operates under steady-state conditions. The effects of kinetic energy and gravity are neglected. A reference analysis is conducted with an ambient temperature of 25 °C and atmospheric pressure of 101.325 kPa.


Heat loss and pressure drop during the heat transfer process are disregarded.

System Description and Assumptions
It is crucial to investigate and evaluate the performance of AACs cycles considering the physical and environmental characteristics of the newly studied zeotropic alternative refrigerants.Figure 2a illustrates the schematic diagram of the traditional vapor compression refrigeration cycle in AACs, comprising an evaporator, compressor, cooler/condenser, expansion valve, and internal heat exchanger.The following is a description of the system cycle: (i) 1 -2: non-isentropic compression; (ii) 2-3: constant pressure heat rejection from the cooler/condenser; (iii) 3-3 : subcooling in the internal heat exchanger (IHE); (iv) 3 -4: isenthalpic throttling in the expansion valve (EXV); (v) 4-1: constant pressure heat absorption in the evaporator; (vi) 1-1 : superheating in the internal heat exchanger (IHE).By adjusting the average evaporation temperature, cooler/condenser refrigerant output temperature, cooler/condenser pressure, and refrigerant mass flow rate, the operating conditions of the cycle can be modified.Because the critical point pressure of the refrigerant can be comparatively low compared to the cooler pressure, a transcritical model for the cycle was developed based on the operating conditions.In transcritical cycles, a gas cooler is utilized instead of a condenser, as in subcritical cycles, since there is no phase change above the critical point.As the pressure and temperature of the refrigerant are independent of each other, both pressure and temperature need to be defined at pressures above the critical point.The coefficient of performance of the cycle varies at the same gas cooler outlet temperature but different operating pressures.Therefore, the simulation model's gas cooler was optimized by running it at the optimum gas cooler pressure [36].Figure 2b presents the pressure-enthalpy (p-h) diagram of subcritical/transcritical cycles using R-744 in AACs.
In order to ensure accuracy and simplicity in the current simulation, the following assumptions were considered to evaluate the performance of zeotropic alternative refrigerants in AACs:

•
The entire system operates under steady-state conditions.

•
The effects of kinetic energy and gravity are neglected.

•
A reference analysis is conducted with an ambient temperature of 25 • C and atmospheric pressure of 101.325 kPa.

•
Heat loss and pressure drop during the heat transfer process are disregarded.

•
The refrigerant at the evaporator outlet is assumed to be in a saturated state.

•
The operating parameter settings and modeling assumption values are presented in Tables 3 and 4.

 
The refrigerant at the evaporator outlet is assumed to be in a saturated state.The operating parameter settings and modeling assumption values are presented in Tables 3 and 4.
These assumptions are made to simplify the simulation and focus on the primary aspects of the performance evaluation of zeotropic alternative refrigerants in AACs.

Thermodynamic Model Analysis
When evaluating the system, it is important to consider several key parameters, including evaporator capacity, compressor power, pressure ratio, compressor discharge temperature, coefficient of performance (COP), and exergy efficiency.These parameters are used to analyze the energy and exergy performance of the AAC system when using different zeotropic alternative refrigerants.Furthermore, model validation is necessary to verify the accuracy of the AAC system simulation.The validation process involves


The refrigerant at the evaporator outlet is assumed to be in a saturated state.


The operating parameter settings and modeling assumption values are presented in Tables 3 and 4.
These assumptions are made to simplify the simulation and focus on the primary aspects of the performance evaluation of zeotropic alternative refrigerants in AACs.

Thermodynamic Model Analysis
When evaluating the system, it is important to consider several key parameters, including evaporator capacity, compressor power, pressure ratio, compressor discharge temperature, coefficient of performance (COP), and exergy efficiency.These parameters are used to analyze the energy and exergy performance of the AAC system when using different zeotropic alternative refrigerants.Furthermore, model validation is necessary to verify the accuracy of the AAC system simulation.The validation process involves


The refrigerant at the evaporator outlet is assumed to be in a saturated state.


The operating parameter settings and modeling assumption values are presented in Tables 3 and 4.
These assumptions are made to simplify the simulation and focus on the primary aspects of the performance evaluation of zeotropic alternative refrigerants in AACs.

Thermodynamic Model Analysis
When evaluating the system, it is important to consider several key parameters, including evaporator capacity, compressor power, pressure ratio, compressor discharge temperature, coefficient of performance (COP), and exergy efficiency.These parameters are used to analyze the energy and exergy performance of the AAC system when using different zeotropic alternative refrigerants.Furthermore, model validation is necessary to verify the accuracy of the AAC system simulation.The validation process involves


The refrigerant at the evaporator outlet is assumed to be in a saturated state.


The operating parameter settings and modeling assumption values are presented in Tables 3 and 4.
These assumptions are made to simplify the simulation and focus on the primary aspects of the performance evaluation of zeotropic alternative refrigerants in AACs.

Thermodynamic Model Analysis
When evaluating the system, it is important to consider several key parameters, including evaporator capacity, compressor power, pressure ratio, compressor discharge temperature, coefficient of performance (COP), and exergy efficiency.These parameters are used to analyze the energy and exergy performance of the AAC system when using different zeotropic alternative refrigerants.Furthermore, model validation is necessary to verify the accuracy of the AAC system simulation.The validation process involves


The operating parameter settings and modeling assumption values are presented in Tables 3 and 4.
These assumptions are made to simplify the simulation and focus on the primary aspects of the performance evaluation of zeotropic alternative refrigerants in AACs.

Thermodynamic Model Analysis
When evaluating the system, it is important to consider several key parameters, including evaporator capacity, compressor power, pressure ratio, compressor discharge temperature, coefficient of performance (COP), and exergy efficiency.These parameters are used to analyze the energy and exergy performance of the AAC system when using different zeotropic alternative refrigerants.Furthermore, model validation is necessary to verify the accuracy of the AAC system simulation.The validation process involves These assumptions are made to simplify the simulation and focus on the primary aspects of the performance evaluation of zeotropic alternative refrigerants in AACs.

Thermodynamic Model Analysis
When evaluating the system, it is important to consider several key parameters, including evaporator capacity, compressor power, pressure ratio, compressor discharge temperature, coefficient of performance (COP), and exergy efficiency.These parameters are used to analyze the energy and exergy performance of the AAC system when using different zeotropic alternative refrigerants.Furthermore, model validation is necessary to verify the accuracy of the AAC system simulation.The validation process involves comparing the model predictions with experimental or reference data to ensure the reliability of the results.The details of these parameters and the importance of model validation are discussed in the following sections.

Energy Analysis
The energy analysis of the modeled AAC system involves applying the principle of conservation of energy to each component of the system.This analysis assumes that the system is in a steady-state condition, with negligible changes in kinetic and potential energies.To calculate the evaporator capacity, the principle of conservation of energy is specifically applied to the evaporator unit, as follows: .
where ṁr is the refrigerant mass flow rate, and h is the refrigerant's specific enthalpy.
The same approach can also be used to determine the cooler/condenser capacity of the cooler/condenser unit, as follows: .
The compressor power absorbed by the refrigerant during the compression process can be calculated using the following equation, assuming the compressor operates adiabatically: .
where η is is the compressor isentropic efficiency, which is given by: where h 2,s represents the specific enthalpy of the refrigerant at the outlet of the compressor during isentropic compression, and h 2 represents the specific enthalpy of the refrigerant at the outlet of the compressor during the actual compression operation.The compressor pressure ratio is defined as: where γ is the pressure ratio of the compressor in the system, P 2 is the compressor's exhaust pressure in MPa, and P 1 is the compressor's suction pressure in MPa.The cycle's coefficient of performance:

Exergy Analysis
The study of exergy, unlike energy, reveals that exergy is not conserved and is subject to destruction in processes that involve both useful work and waste.Conducting a thermal exergy analysis is valuable in identifying the sources and magnitudes of thermodynamic inefficiencies.To achieve this objective, the following general form of the steady-state exergy rate balance equation for control volumes [38] can be applied to the components of the AAC system: .
where T 0 represents the environmental temperature dead state, .
Q j denotes the time rate of heat transfer at the boundary location, T j is the boundary temperature, .
W cv is the amount of work conducted in the control volume, .
E is the flow exergy rate, whereas in and out denote the inlet and outlet, respectively, and .I is the rate of exergy destruction in the control volume.
The exergy rate in the system at state point i is defined as follows: .
where s is the refrigerant's specific entropy, and subscript "0" indicates the reference (dead) condition.The exergy destruction in the compressor is influenced by the internal heat transfer, gas friction, and mechanical friction of the moving elements.Assuming adiabatic compression, the exergy rate balance equation leads to the following equation for the rate of exergy destruction in the compressor: .
The rate of exergy destruction in the condenser, resulting from heat transfer between the refrigerant and air streams, can be calculated using the following equation [22]: .
where T r is the refrigerant temperature.Assuming there is no heat transfer to or from the environment, the exergy destruction in the internal heat exchanger, which is caused by the internal heat transfer between the refrigerant streams, can be calculated using the following equation: .
The exergy destruction in the expansion valve, which is caused by internal friction and a rapid pressure drop, can be calculated by disregarding the heat transfer with the environment using the following equation: .
To calculate the rate of exergy destruction in the evaporator, which is caused by the heat transfer between the refrigerant and air streams, the following equation can be used [22]: .
The rate of overall exergy destruction in the AAC system's cycle can be calculated by summing up the individual exergy destruction rates in the system's components, as follows: .
Finally, the following expression for exergetic efficiency can be used to evaluate the overall exergetic performance of the AAC system:

Model Validation
The accuracy of the AAC system simulation model, generated with Aspen HYSYS, was validated through model validation.The transcritical CO 2 refrigeration cycle with a regenerator model was validated using the same boundary and operating parameters as reported by Rigola et al. [39] and Elattar and Nada [37].The experimental [39] and numerical [37] results based on COP were compared with the results obtained from the current model, and the maximum relative errors were found to be 11.24% and 2.79% for experimental and numerical results, respectively, as presented in Table 5.These results demonstrate that the simulation model technique implemented in this study using Aspen HYSYS software is sufficiently accurate.

Zeotropic Refrigerants Environmental Impacts
Considering their potential for advancing refrigeration cycles, it is important to note that R134a, R32, R125, R1234yf, and R744 have negative impacts on the depletion of the ozone layer.However, R134a and R32 also have positive impacts on global warming, whereas R1234yf and R744 do not contribute to global warming.To assess the environmental impact of the proposed refrigerant mixtures, Table 1 provides information on the ozone depletion potential (ODP) and global warming potential (GWP) of the individual pure components of refrigerants [14,40].To evaluate the environmental impact of the proposed zeotropic refrigerants, Equations ( 16) and ( 17) were used, and the results are presented in Table 2.
In the equations, w pure,i represents the mass fraction of component i in a zoetrope.The evaluated ozone depletion potential (ODP) and global warming potential (GWP) of the proposed zeotropic refrigerants are presented in Table 2.It is worth noting that all the studied zeotropic blends have an ODP of zero, as their individual components have an ODP of zero.Furthermore, based on the classification by UNIP, 2019 [41], for 100-year global warming potential levels, R-455A is classified as having a negligible GWP level, while R-463A is classified as having a high GWP level.

Parametric Studies
A transcritical/subcritical cycle was developed for R744, R455A, and R463A due to the critical temperatures of these refrigerants.In transcritical cycles, a gas cooler is used instead of a condenser since there is no phase change above the critical point.At pressures above the critical point, the pressure and temperature of the refrigerant become independent variables and both need to be specified.The coefficient of performance (COP) of the cycle is influenced by the cooler/condenser outlet temperature and the operating pressures.Figure 3 illustrates the optimum cooler/condenser pressure, P 2 , for the studied zeotropic refrigerants, R455A and R463A, compared to R744, based on the maximum COP of the cycle.The common operating conditions for the comparison are t evap = 7.5 • C, t 3 = 35 • C, and ṁr = 0.075 kg/s.According to the figure, the maximum COP values achieved with R744, R455A, and R463A were 3.1, 4.25, and 5.4, respectively.These values were obtained at the optimum pressures of 8.9, 1.65, and 2.3 MPa for R744, R455A, and R463A, respectively.
above the critical point, the pressure and temperature of the refrigerant become i ent variables and both need to be specified.The coefficient of performance (CO cycle is influenced by the cooler/condenser outlet temperature and the opera sures.Figure 3 illustrates the optimum cooler/condenser pressure, P2, for the stu tropic refrigerants, R455A and R463A, compared to R744, based on the maximu the cycle.The common operating conditions for the comparison are tevap = 7.5 °C, and ṁr = 0.075 kg/s.According to the figure, the maximum COP values achi R744, R455A, and R463A were 3.1, 4.25, and 5.4, respectively.These values were at the optimum pressures of 8.9, 1.65, and 2.3 MPa for R744, R455A, and R463 tively.W comp ), exergy efficiency (η ex ), and coefficient of performance (COP) with the cooler/condenser pressure (P 2 ) at t evap = 7.5 • C, t 3 = 35 • C, and ṁr = 0.075 kg/s.In Figure 4a, it is evident that as P 2 increased, the evaporator capacity ( .Q evap ) sharply increased for R744, while it remained nearly constant for R455A and R463A.This behavior is due to the difference in refrigerant enthalpy through the evaporator.With R744, the refrigerant enthalpy decreased at the inlet of the evaporator, while it remained approximately constant for R455A and R463A as P 2 increased.This is because the cycle is transcritical for R744 using a gas cooler, while it is subcritical for R455A and R463A using a condenser.Furthermore, R463A exhibited the highest evaporator capacity ( .Q evap ) among the three refrigerants, while R455A had a higher .Q evap than R744 until P 2 ≤ 8.9 MPa, after which R744 surpassed R455A.Figure 4b illustrates the compressor power ( .W comp ) of R744, R455A, and R463A.It was observed that .W comp increased with an increase in the cooler/condenser pressure (P 2 ).This behavior can be attributed to the thermodynamic relationship between pressure and temperature and the nature of the thermodynamic cycle.
exhibited the highest evaporator capacity ( evap Q  ) among the three refrigerants, while R455A had a higher evap Q  than R744 until P2 ≤ 8.9 MPa, after which R744 surpassed R455A.Figure 4b illustrates the compressor power ( comp W  ) of R744, R455A, and R463A.It was observed that comp W  increased with an increase in the cooler/condenser pressure (P2).This behavior can be attributed to the thermodynamic relationship between pressure and temperature and the nature of the thermodynamic cycle.Figure 4c,d illustrates the impact of cooler/condenser pressure (P2) on exergy efficiency (ex) and coefficient of performance (COP).For R744, it can be observed that as the pressure of the cooler/condenser increased, both ex and COP initially rose, reached a peak around P2 = 9 MPa, and then started to decline.This behavior is primarily attributed to the increase in evaporator capacity ( evap Q  ), which dominated over the increase in compressor power ( comp W  ) until P2 ≤ 9 MPa.However, after P2 > 9 MPa, the situation reversed, and the (a)  Figure 4c,d illustrates the impact of cooler/condenser pressure (P 2 ) on exergy efficiency (η ex ) and coefficient of performance (COP).For R744, it can be observed that as the pressure of the cooler/condenser increased, both η ex and COP initially rose, reached a peak around P 2 = 9 MPa, and then started to decline.This behavior is primarily attributed to the increase in evaporator capacity ( .Q evap ), which dominated over the increase in compressor power ( .W comp ) until P 2 ≤ 9 MPa.However, after P 2 > 9 MPa, the situation reversed, and the COP became inversely correlated with .Q evap and .W comp .On the other hand, for R455A and R463A, the exergy efficiency (η ex ) and coefficient of performance (COP) decreased with increasing P 2 .This decrease is mainly due to the increase in compressor power ( .W comp ) and the nearly constant values of evaporator capacity ( .Q evap ) as P 2 rose.At P 2 = 15 MPa, the maximum .Q evap values for R744, R455A, and R463A were 12.2 kW, 9 kW, and 12.5 kW, respectively.Furthermore, within the studied range of P 2 , the increase in evaporator capacity ( .Q evap ) was 614% for R463A compared to R744 at P 2 = 7 MPa and 2.5% at P 2 = 15 MPa.Additionally, within the studied range of P 2 , the COP of R455A and R463A decreased by 27% and 37%, respectively.

Impact of Average Evaporator Temperature, t evap
Figure 5 illustrates the impact of the average evaporator temperature (t evap ) on the following key performance parameters: evaporator capacity ( .Q evap ), compressor power ( .W comp ), exergy efficiency (η ex ), and coefficient of performance (COP).The analysis was performed at a condenser temperature of t 3 = 35 • C, mass flow rate of ṁr = 0.075 kg/s, and condenser pressure of P 2 = P opt .The optimum pressures for R744, R455A, and R463A were 8.9, 1.65, and 2.3 MPa, respectively, as shown in Table 2. Table 2 also provides the critical pressures and temperatures for R744, R455A, and R463A, indicating that the R744 cycle operates in the transcritical region, while the cycles for R455A and R463A are in the subcritical region.The .Q evap of R455A and R463A increased with rising t evap , while it decreased for R744.Additionally, the .Q evap of R744 was higher than that of R455A until t evap reached 9 • C. and COP became more pronounced when t3 exceeded 35 °C.This phenomenon occurs due to the reduction in refrigerant enthalpy difference across the evaporator with increasing t3, while the compression work remained constant.Figure 6b indicates that t3 had no significant effect compressor power ( comp W  ) since the cooler/condenser pressure (P2) and average evaporator temperature (tevap) were fixed at their optimum values of 7.5 °C.Exergy efficiency (ex) serves as a fundamental measure of the system's thermodynamic performance.Figure 6d illustrates that for R744, R455A, and R463A, as t3 increased, (a) In Figure 5b,c, it can be observed that for all refrigerants, .W comp decreased and COP increased with increasing average evaporator temperature (t evap ).This trend occurs due to the decrease in compression work with higher t evap while maintaining a constant cooler/condenser pressure (P 2 ).Furthermore, Figure 5d demonstrates the effect of t evap on exergy efficiency (η ex ), clearly showing a decrease as t evap increased.At t evap = 15 • C, the maximum .Q evap values were 10 kW and 12.5 kW for R455A and R463A, respectively, while R744 achieved a maximum .Q evap of 9.48 kW at t evap = 5 • C. Within the studied range of t evap , using R455A and R463A instead of R744 resulted in .Q evap enhancements of 16.3% and 49%, respectively, at t evap = 15 • C. The maximum COP values were 4.25, 6.7, and 7.5 for R744, R455A, and R463A, respectively, at t evap = 15 • C. Additionally, within the studied range of t evap , the COP for R744, R455A, and R463A increased by 52%, 76%, and 56%, respectively.The improvements in COP were 57.6% and 76.5% when using R455A and R463A instead of R744, respectively, at t evap = 15 W comp ), exergy efficiency (η ex ), and coefficient of performance (COP).The analysis was conducted at a fixed average evaporator temperature of t evap = 7.5 • C, mass flow rate of ṁr = 0.075 kg/s, and condenser pressure of P 2 = P opt .The optimum pressures for R744, R455A, and R463A were 8.9, 1.65, and 2.3 MPa, respectively.As shown in Figure 6a,c, the .Q evap and COP decreased as the cooler/condenser outlet temperature (t 3 ) increased for all refrigerants.Furthermore, the decline in .Q evap and COP became more pronounced when t 3 exceeded 35 • C.This phenomenon occurs due to the reduction in refrigerant enthalpy difference across the evaporator with increasing t 3 , while the compression work remained constant.Figure 6b indicates that t 3 had no significant effect on compressor power ( .W comp ) since the cooler/condenser pressure (P 2 ) and average evaporator temperature (t evap ) were fixed at their optimum values of 7.5 • C. R455A, and R463A decreased by 57%, 77%, and 67%, respectively.The maximum COP values were 4.8, 6.5, and 6.5 for R744, R455A, and R463A, respectively, at t3 = 20 °C.Moreover, within the studied range of t3, the COP for R744, R455A, and R463A decreased by 58.3%, 77%, and 66%, respectively.When compared to R744 at t3 = 20 °C, the improvements in COP were 35.4% and 76.5% when using R455A and R463A, respectively.Exergy efficiency (η ex ) serves as a fundamental measure of the system's thermodynamic performance.Figure 6d illustrates that for R744, R455A, and R463A, as t 3 increased, ηex initially rose, reached a peak, and then started to decline.Typically, this peak in η ex occurs around t 3 = 35 • C. The decrease in η ex after t 3 > 35 • C can be attributed to the increase in total irreversibility of the cycle while keeping the compressor power ( .W comp ) constant.
W comp increased with an increase in the refrigerant flow rate ( ṁr ).However, ṁr had no significant impact on the exergy efficiency (η ex ) and coefficient of performance (COP).The maximum .Q evap values were 18.8 kW, 17.7 kW, and 24.9 kW for R744, R455A, and R463A, respectively, at ṁr = 0.15 kg/s.Within the studied range of ṁr , the .Q evap for R744, R455A, and R463A increased by 208%, 200%, and 207%, respectively.Similarly, the maximum .W comp values were 6 kW, 4 kW, and 4.6 kW for R744, R455A, and R463A, respectively, at ṁr = 0.15 kg/s.Within the studied range of ṁr , the .W comp for R744, R455A, and R463A increased by 200%, 700%, and 206%, respectively.The maximum COP values were 3.15, 4.4, and 5.3 for R744, R455A, and R463A, respectively, at t 3 = 20 • C. Within the studied range of ṁr , the improvements in COP were 40% and 68% when using R455A and R463A instead of R744, respectively.It is worth noting that while         Q evap , COP, and ηex, with values of 12.35 kW, 5.39, and 47.8%, respectively, surpassing both R744 and R455A.To further understand the system behavior, the percentage of exergy destruction for each component in the studied refrigerant blends compared to R744 is illustrated in Figure 10.This provides insight into the influence of each component on the overall exergy destruction.It can be observed that the valve (EXV) contributed the highest irreversibility, accounting for 42% in the R455A cycle, followed by the compressor with 35% in the R463A cycle.The evaporator and condenser also contributed to exergy destruction in all configurations, albeit to a lesser extent.Overall, the findings from these analyses demonstrated the performance advantages and trade-offs associated with different refrigerant blends in the AACs system, highlighting the potential of R463A as the most promising option based on the evaluated performance indicators.

System Optimization
After conducting a parametric analysis, the performance of two refrigerant blend mixtures (R455A and R463A) and pure CO2 (R744) was optimized to maximize the coefficient of performance (COP) within the studied ranges of operating condition parameters.The optimal operating conditions, based on maximum COP, are presented in Table 6.Furthermore, Table 7 lists the exergy destruction and efficiency for each component, calculated based on the optimal COP case.The table provides the maximum system performance and refrigeration capacity achievable for R744, R455A, and R463A, with values of 14.21 kW, 15.05 kW, and 19.15 kW, respectively.The corresponding maximum COP values were 14.58, 12.86, and 14.19, and the exergy efficiencies (ηex) were 45.4%, 26.8%, and 26.7%, respectively.
The outcomes of the optimization study were further analyzed in terms of COP, cycle operating pressures, and refrigerant environmental impact.In terms of COP, R744 and R463A yielded the best results, with their cycle performance indicators being most similar.Regarding cycle operating pressures, R463A and R455A exhibited lower pressures compared to R744.This aspect may have implications for the materials used in cycle components, cycle leakages, longevity, and compressor lubrication.In terms of environmental impact, R455A showed the least negative impact on the environment with a global warming potential (GWP) less than 1.However, it had a slightly lower COP compared to R744 and R463A.Additionally, the average difference in performance (COP) between R455A and R744 and R463A was approximately 11%.
Considering the need for environmentally friendly air conditioning systems with acceptable performance in cars, given their impact on the environment and contribution to global warming, it is recommended to use the R455A blend as the refrigerant in the AAC system.

System Optimization
After conducting a parametric analysis, the performance of two refrigerant blend mixtures (R455A and R463A) and pure CO 2 (R744) was optimized to maximize the coefficient of performance (COP) within the studied ranges of operating condition parameters.The optimal operating conditions, based on maximum COP, are presented in Table 6.Furthermore, Table 7 lists the exergy destruction and efficiency for each component, calculated based on the optimal COP case.The table provides the maximum system performance and refrigeration capacity achievable for R744, R455A, and R463A, with values of 14.21 kW, 15.05 kW, and 19.15 kW, respectively.The corresponding maximum COP values were 14.58, 12.86, and 14.19, and the exergy efficiencies (η ex ) were 45.4%, 26.8%, and 26.7%, respectively.The outcomes of the optimization study were further analyzed in terms of COP, cycle operating pressures, and refrigerant environmental impact.In terms of COP, R744 and R463A yielded the best results, with their cycle performance indicators being most similar.Regarding cycle operating pressures, R463A and R455A exhibited lower pressures compared to R744.This aspect may have implications for the materials used in cycle components, cycle leakages, longevity, and compressor lubrication.In terms of environmental impact, R455A showed the least negative impact on the environment with a global warming potential (GWP) less than 1.However, it had a slightly lower COP compared to R744 and R463A.Additionally, the average difference in performance (COP) between R455A and R744 and R463A was approximately 11%.
Considering the need for environmentally friendly air conditioning systems with acceptable performance in cars, given their impact on the environment and contribution to global warming, it is recommended to use the R455A blend as the refrigerant in the AAC system.

Conclusions
The growing popularity of vehicles has led to increased time spent in cars equipped with air conditioning systems.However, the refrigerant fluids used in these systems often have significant negative environmental effects.To address this issue, many countries have implemented regulations mandating the use of refrigerants with minimal global warming potential (GWP) and zero ozone depletion potential (ODP) in cars.Therefore, there is a pressing need for environmentally friendly and high-performance air conditioning systems in vehicles, given their detrimental impact on the environment and contribution to global warming in addition to lower engine fuel consumption and hence lower carbon emissions.
The objective of this study was to propose the use of CO 2 -based blend zeotropic refrigerants, specifically R-455A (R-744/32/1234yf) and R-463A (R-744/32/125/1234yf/134a), to improve the thermodynamic performance of a pure CO 2 refrigerant in automotive air conditioning systems.The study involved analyzing the thermodynamic energy and exergy of an AAC system and optimizing it using Aspen HYSYS software.The investigation focused on comparing the performance of the new zeotropic refrigerant blends with that of carbon dioxide (R-744).Additionally, the study examined the impact of cooler/condenser pressure, average evaporator temperature, cooler/condenser outlet temperature, and refrigerant flow rate on the cycles' coefficient of performance (COP) and exergy efficiency.
The results of the analyses indicate that R-463A provided the best system performance among the investigated refrigerants, resulting in the longest possible driving range when used in automotive air conditioning systems compared to R744 and R-455A at the same operating conditions.The study identified optimal cooler/condenser pressures for R744, R455A, and R463A as 8.9, 1.65, and 2.3 MPa, respectively.The maximum COP values based on these optimal cooler/condenser pressures were 3.1, 4.25, and 5.4 for R744, R455A, and R463A, respectively.Furthermore, the study revealed the maximum refrigeration capacities ( .Q evap ) for R744, R455A, and R463A as 12.2, 9, and 12.5 kW, respectively, at a specific

.
Q I HE = .m r (h I HE,in − h I HE,out ) Hot = .m r (h I HE,out − h I HE,in ) cold .I I HE = ∑ .E I HE,in − ∑ .E I HE,out Heat exchanger model: Simple end point ∆P Hot stream = 0 kPa ∆P Cold stream = 0 kPa

Figure 4 .
Figure 4. Impact of cooler/condenser pressure on AAC output parameters: (a)

P 2 =Figure 5 .
Figure 5. Impact of average evaporator temperature on AAC output parameters: (a)

Processes 2023 ,
11, x FOR PEER REVIEW 15 of 23 ex initially rose, reached a peak, and then started to decline.Typically, this peak in ex occurs around t3 = 35 °C.The decrease in ex after t3 > 35 °C can be attributed to the increase in total irreversibility of the cycle while keeping the compressor power ( comp W  ) constant.At t3 = 20 °C, the maximum evap Q  values were 14.1 kW, 12.9 kW, and 15 kW for R744, R455A, and R463A, respectively.Within the studied range of t3, the evap Q  for R744,

.
significantly influenced by the refrigerant flow rate ( ṁr the exergy efficiency (η ex ) and coefficient of performance (COP) remained unaffected by changes in ṁr .
the refrigerant flow rate (ṁr), the exergy efficiency (ex) and coefficient of performance (COP) remained unaffected by changes in ṁr.

Figure 7 .
Figure 7. Impact of refrigerant flow rate on AAC output parameters: (a) evap Q  , (b) comp W  , (c)

Figure 7 .
Figure 7. Impact of refrigerant flow rate on AAC output parameters: (a)

Figure 8
Figure8provides a cycle analysis and the state properties at the inlet and exit of each component for the R744 refrigerant using Aspen HYSYS software.The purpose of this analysis was to compare and assess the performance indicators of the Advanced Adiabatic Compressed Air Energy Storage (AACs) system when charged with pure R744 and with R455A and R463A blend refrigerants as working fluids.The selected performance indicators include evaporator capacity ( .Q evap ), compressor power ( .W comp ), exergy efficiency (η ex ), and coefficient of performance (COP).

Figure 9
Figure 9 presents the system enhancements achieved with the different working fluids.It is evident that the blend R463A exhibited the highest values for.

Figure 9 .
Figure 9. Comparisons of studied refrigerant blends on AAC performance indicators: (a) evap Q  ,

Figure 10 .
Figure 10.Exergy destruction percentages of each component for the studied refrigerant blends compared to R-744.

Figure 10 .
Figure 10.Exergy destruction percentages of each component for the studied refrigerant blends compared to R-744.

Table 1
presents the key physical and environmental characteristics of individual refrigerants: R134a, R32, R125, R1234yf, and R744.As shown in Table1, all refrigerants have an ozone depletion potential (ODP) of 0. However, R125 and R134a have higher global warming potential (GWP) values compared to the other refrigerants.Notably, R744 and R1234yf have the lowest GWP values.Additionally, R744 has the highest critical pressure and the lowest normal boiling point, with values of 7.3 MPa and −78 • C, respectively.

Table 4 .
Energy and exergy equations and modelling assumptions.

Table 4 .
Energy and exergy equations and modelling assumptions.

Table 4 .
Energy and exergy equations and modelling assumptions.

Table 4 .
Energy and exergy equations and modelling assumptions.

Table 4 .
Energy and exergy equations and modelling assumptions.

Table 4 .
Energy and exergy equations and modelling assumptions.

Table 6 .
Optimal operating parameters based on maximum COP.

Table 6 .
Optimal operating parameters based on maximum COP.

Table 7 .
Exergy destruction and efficiency for each component for the optimal COP case.