Pressure Drop and Heat Transfer Characteristics of TiO 2 /R1234yf Nanorefrigerant: A Numerical Approach

: Global warming is one of the most dangerous ecological issues facing the globe. Re-frigerants are a major contributor to global warming. This investigation mainly focuses on the analysis of a greener nanorefrigerant. Nanorefrigerant can improve the efﬁciency of refrigeration and air conditioning systems that use vapor compression. In the present investigation, mathematical and computational methods are used to assess the heat transfer and pressure drop properties of TiO 2 /R1234yf. In order to analyze the heat transfer characteristics and the transport features of the innovative nanorefrigerant, appropriate mathematical predictive models were adapted from earlier investigations. The models are validated by the experiments using TiO 2 /POE nanolubricant as a test ﬂuid. The investigation was conducted with a temperature range of 10 ◦ C to 40 ◦ C and a volume percentage of nano-sized TiO 2 particles in R1234yf refrigerant ranging from 0.2 to 1%. According to the research, the introduction of nanoparticles increases viscosity, thermal conductivity, and density. However, as the amount of nanoparticles rises, the speciﬁc heat capacity of the nano-enhanced refrigerant decreases. The nanorefrigerant’s heat transfer coefﬁcient and pressure drop are improved by 134.03% and 80.77%, respectively. The outcomes observed from the predictive technique and the simulation approach had an average absolute variation of 9.91%.


Introduction
One of the most significant ecological problems influencing our planet at the present time is global warming.The prior production of refrigerant gases and carbon dioxide emissions is one of the reasons for them [1,2].These gases are released into the atmosphere, where they remain for a longer time and contribute to global warming [3].The life cycle analysis (LCA) of these gases [4] and the conversion effectiveness of refrigeration systems [5] help to quantify how seriously these gases have an influence on the environment.Over the duration of the next 100 years, Earth's temperature is predicted to increase by 2.5 to 10 • F as a result of global warming [6].Due to the impact on agriculture, there will be more heat waves, strong precipitation patterns, and rainfalls, and by the end of the next century, the sea level may increase by 1-4 feet.Therefore, in order to safeguard the ozone layer and limit net planet warming by 0.5 • C by the year 2100, the member countries of the Montreal Protocol decided in October 2016 to phase down HFCs and HCFCs.Chlorofluorocarbons (CFCs) and hydrochlorofluorocarbons (HCFCs), which were already thinning the ozone layer [7], were eliminated following the Montreal and Kyoto Protocols, respectively [8].Currently, HFCs, including R134a, R23, R404A, R407A, R410A, and R507A, make up the majority of the various refrigerants used in residential, automobile, commercial, and industrial refrigeration and air-conditioning systems [9,10].HFCs have a high global warming potential (GWP) despite having no ozone depletion potential (ODP).
Refrigeration systems are a significant energy consumption industry for the majority of advanced and some emerging countries worldwide.As a result, researchers continue to hunt for effective ways to improve the energy efficiency of refrigeration systems.An inventive technique for cooling systems employs nanorefrigerants to increase efficiency.One type of emerging nanofluid utilized in refrigeration systems to enhance heat transfer efficiency and system performance is nanorefrigerants [11].This improvement is made achievable by the nanoscale particles' greater thermal conductivity.The performance of the system's heat transmission is being improved [12], the solubility of the refrigerant and lubricant is being improved [13], and the friction coefficient and wear heat are being optimized [14].The basic refrigerant's ability to transmit heat more effectively may be increased by mixing in nanoscale particles, which boosts the system's overall performance [15].Two-phase heat transfer fluids are used as refrigerants.Refrigerant characteristics have an impact on the heat transfer capabilities and energy efficiency of the complete refrigeration system.The globally adopted refrigerant for different applications of refrigeration systems is R134a.It comes under the 'A' safety rating, and it has no flame propagation and low toxicity; meanwhile, it has about 1430 GWP [16].Due to the high GWP of R134a, researchers are tried to replace it with low GWP refrigerants.HFOs are one of the suitable candidates for substituting the R134a with or without modifications in the refrigeration systems [8,17].
When Jankovic et al. [18] examined the operation of a small power refrigeration system, they found that the cooling efficiency and COP of R1234yf were 10% and 9% lesser than those of R134a, respectively.Similar findings were made by Sieres et al. [19], who discovered that R1234yf's compressor power was 3% higher while cooling power was 6% lower than R134a.Chen et al. [20] and Li et al. [21] investigated R1234yf's potential in a refrigeration system without oil.Both came to the conclusion that R1234yf had a lower volumetric efficiency and cooling capacity than R152a and R134a.According to Navarro-Esbr et al. [22], R1234yf had 9% and 5% less cooling power and volumetric efficiencies than R134a.Investigations by Sanchez et al. [23] also showed that R1234yf's cooling impact was 4.4% to 8.5% less than R134a's.Mota-Bobiloni et al. [24] reported comparable outcomes as well.
Nanorefrigerants perform refrigeration systems better than traditional refrigerants, according to recent studies.The R140a refrigerant's efficiency was increased by 40% by the addition of Al 2 O 3 nanoparticles [25].The COP of the system containing Al 2 O 3 /R134a was determined by research by Ambhore et al. [26] to be 2.03.Similar findings were made by Subhedar et al. [27], who discovered that using Al 2 O 3 /R134a as a refrigerant increased COP by almost 85%.According to Payyala et al. [25], adding Al 2 O 3 to R140a improved the COP, pressure ratio, and energy efficiency ratio.The coefficient of the performance was increased by 14.55%, and the friction coefficient was reduced by 9.9% when the CuO nanomaterial was combined with R134a [28].CuO/R134a refrigerant's heat transfer coefficient is enhanced by 42-82% at the concentration of 1% by volume, according to Bartelt et al. [29].When CuO/R113a nano-enhanced refrigerant was examined by Katoch et al. [30] in the refrigeration system, they discovered that the energy consumption decreased by 19.82% at 0.5% nanoparticle addition.CuO/LPG-based systems improve efficiency and heat transfer rate by 36% and 46%, respectively [31].According to Dhamneya et al. [32], the evaporatively cooled condenser with TiO 2 /R134a performed substantially better and had a 51% increase in COP.Rahman et al. [33] discovered that the compressor energy utilization was reduced by 34%; however, the COP was enhanced by 4.59%.The results show that the nano-enhanced refrigerants contribute to increasing the whole system performance and heat transfer effectiveness and that the COP rises as the concentration of nanoparticles does as well [34].
The key drawback of the R1234yf refrigerant is its lower performance than R134a.Therefore, the emphasis of this study is on increasing system efficiency using R1234yf by employing nanoparticles.Few investigations have been done on the heat transfer capabilities, pressure drop, and thermo-transport characteristics of R1234yf-based nano-enhanced refrigerants.Analysis of the pressure drop properties and HTC of the R1234yf-based nanorefrigerants requires more investigation.Through the use of simulation techniques and mathematical models, this study examines the heat transfer properties and pressure drop of the TiO 2 /R1234yf refrigerant.The findings of the simulation approach are com-pared with numerical results, and the predictive models are validated using various types of prior literature.This analysis leads to the sustainable development of the refrigeration system with greener refrigerants.

Materials
The performance of the refrigeration systems depends on the heat transfer fluids used in it.In this analysis, TiO 2 nano-sized particles and R1234yf refrigerant are used.The properties of TiO 2 nanoparticles, POE lubricant oil and R1234yf refrigerant are tabulated in Tables 1 and 2, respectively.The thermal properties of TiO 2 added to the R1234yf nanoenhanced refrigerant have been determined by adopting the predictive models from the literature.To analyze the impact of nanoparticles on the properties of refrigerant, the concentration of TiO 2 ranges between 0.2 and 1% by volume, and the mean diameter of the particles was 30 nm.Likewise, to study the effect of temperature on the properties, it changes from 10  The polyol ester (POE) is the lubricant oil used in commercial refrigeration systems.The fluid (TiO 2 /POE) utilized for the validation of the predictive models is prepared by a two-step technique.The lubricant oil and TiO 2 nano-sized particles were purchased from Kenizol and Sigma Alrich, respectively.The scanning electron microscopy (SEM) image shows the size of the nanoparticle, which is shown in Figure 1.The SEM was carried out by a Quanta device with a magnification of ×120,000.It indicated that the average size of the nanoparticle is 30 nm and spherical in shape.The average size of the nanoparticle chosen below 30 nm was to avoid the sedimentation of the particle in the base fluid.Initially, the TiO 2 nano-sized particles with different volume fractions (0.2 to 1%) were dispersed into the POE lubricant with the aid of the magnetic stirrer for 2 h.After that, the blend was sonicated using an ultra sonicator for 5 h.This helped to obtain a homogeneous nanoliquid.

Predictive Models
Predictive models were utilized to identify the thermo-physical properties of the TiO2/R1234yf nanorefrigerant, such as thermal conductivity, viscosity, specific heat and density, as well as predictions of the heat transfer coefficient and pressure drop of the evaporator tubes.The major assumptions considered for the predictive analysis were:  The nanoparticles are spherical in nature;  The nanorefrigerant is treated as homogeneous;  No surfactant is used;  Single-phase approach.

Thermal Conductivity
The temperature and concentration of the nanoparticles in the working fluid, the size of the particle and the interfacial layer thickness impact the thermal conductivity of the nano-enhanced refrigerant.Stiprasert et al. [36] developed a predictive model based on these factors.The thermal conductivity of TiO2/R1234yf nano-enhanced refrigerant was identified by:

Viscosity
The viscosity induces the shear stress in the flow of the working fluids.It affects the flow of the fluid.Due to viscosity, a pressure drop in the flow took place.The Brinkman model [37] was adopted to determine the viscosity of the nano-enhanced refrigerant. (2)

Density
One of the important parameters influencing the refrigeration system is the density of the working fluid.It directly affects the pressure in the system.The density is analyzed by the Pak and Cho model [38].

Predictive Models
Predictive models were utilized to identify the thermo-physical properties of the TiO 2 /R1234yf nanorefrigerant, such as thermal conductivity, viscosity, specific heat and density, as well as predictions of the heat transfer coefficient and pressure drop of the evaporator tubes.The major assumptions considered for the predictive analysis were: The nanoparticles are spherical in nature; The nanorefrigerant is treated as homogeneous; No surfactant is used; Single-phase approach.

Thermal Conductivity
The temperature and concentration of the nanoparticles in the working fluid, the size of the particle and the interfacial layer thickness impact the thermal conductivity of the nano-enhanced refrigerant.Stiprasert et al. [36] developed a predictive model based on these factors.The thermal conductivity of TiO 2 /R1234yf nano-enhanced refrigerant was identified by:

Viscosity
The viscosity induces the shear stress in the flow of the working fluids.It affects the flow of the fluid.Due to viscosity, a pressure drop in the flow took place.The Brinkman model [37] was adopted to determine the viscosity of the nano-enhanced refrigerant.

Density
One of the important parameters influencing the refrigeration system is the density of the working fluid.It directly affects the pressure in the system.The density is analyzed by the Pak and Cho model [38].

Specific Heat
The specific heat of a refrigerant is the amount of energy needed to increase the unit temperature of a unit mass substance.The specific heat of the nanorefrigerant was calculated using [38]:

Flow Boiling Heat Transfer Coefficient
The HTC of the nanorefrigerant is predicted using the correlation of Peng et al. [39].
The nanoparticle impact factor, or NIF, is influenced by the particle concentration, the base refrigerant's and nanomaterials' transport characteristics, the mass flow, and the vapor quality.It is determined by:

Pressure Drop
The frictional pressure drop is evaluated by the relation from Ding et al. [40].It is expressed as: The NIF can be identified by:

Validation of Predictive Models
The predictive models for the thermo-physical properties are validated using the experiment results.For the validation, TiO 2 /POE nanolubricant is used as the test fluid.The thermal conductivity of the liquid is measured using KD 2 thermal property analyzer, and the differential scanning calorimeter (DSC) is utilized for identifying the specific heat of the TiO 2 /POE.Table 3 tabulates the deviation of the results from the prediction and experiment.From the table, it is seen that the variation between the predictive and experiment results is below 10%.Thereby, the results from the predictive models are reliable.

Simulation Method
The simulation study was carried out to identify the HTC and pressure drop of the TiO 2 nanoparticles added to the R1234yf refrigerant in a smooth evaporator tube.The analysis was conducted using Ansys Fluent (version 2020.R2) software.It identified the temperature of the wall and refrigerant, distribution of flow velocity, pressure and the vapor quality of the refrigerant in the evaporator tube to determine the heat transfer characteristic as well as the pressure drop in the test section.
The flow of the simulation analysis is shown in Figure 2. The initial stage of the simulation is the modeling of the geometry.The geometry of the evaporator tube was modeled as per the dimensions given in Figure 3.Then, the meshes were generated and the mesh quality was checked.If the mesh quality was good, it proceeded to the specified boundary conditions.After that, we chose the computation algorithm (SIMPLE) and discretization method.Then, the independency of the grid was checked, and the flow was visualized.
The working fluid (nanorefrigerant) flowed through the test section at a constant mass flux (200 kg/m 2 s) and the vapor quality varied from 0.2 to 0.9. Figure 2 shows the dimensions and boundary conditions of the horizontal evaporator for analysis.
Due to the high accuracy, consistent results and less computational time, the standard k-ε model was utilized to solve the turbulent flow problem. Figure 4 shows the mesh structure of the test section on which the analysis was performed.Therein, 108,750 elements and 112,530 nodes make upthe whole test evaporator section.
To solve the two-phase mathematical model used to simulate the flow boiling in the evaporator tube, some hypotheses were taken into account.The flow field was consider as two-dimensional and turbulent.In addition, the heat flux was uniform along the length of the test section.The fluid was considered to be incompressible, and force vectors as negligible.

Simulation Method
The simulation study was carried out to identify the HTC and pressure drop of the TiO2 nanoparticles added to the R1234yf refrigerant in a smooth evaporator tube.The analysis was conducted using Ansys Fluent (version 2020.R2) software.It identified the temperature of the wall and refrigerant, distribution of flow velocity, pressure and the vapor quality of the refrigerant in the evaporator tube to determine the heat transfer characteristic as well as the pressure drop in the test section.
The flow of the simulation analysis is shown in Figure 2. The initial stage of the simulation is the modeling of the geometry.The geometry of the evaporator tube was modeled as per the dimensions given in Figure 3.Then, the meshes were generated and the mesh quality was checked.If the mesh quality was good, it proceeded to the specified boundary conditions.After that, we chose the computation algorithm (SIMPLE) and discretization method.Then, the independency of the grid was checked, and the flow was visualized.The working fluid (nanorefrigerant) flowed through the test section at a constant mass flux (200 kg/,m 2 s) and the vapor quality varied from 0.2 to 0.9. Figure 2 shows the dimensions and boundary conditions of the horizontal evaporator for analysis.Due to the high accuracy, consistent results and less computational time, the standard k-ε model was utilized to solve the turbulent flow problem. Figure 4 shows the mesh structure of the test section on which the analysis was performed.Therein, 108,750 elements and 112,530 nodes make upthe whole test evaporator section.
To solve the two-phase mathematical model used to simulate the flow boiling in the evaporator tube, some hypotheses were taken into account.The flow field was consider as two-dimensional and turbulent.In addition, the heat flux was uniform along the length of the test section.The fluid was considered to be incompressible, and force vectors as negligible.

Grid Independency Analysis
In order to generate grids that perform the most effectively, it is necessary to take into account the grids' shape, quality, and number.A component that affects the overall processing cost and the precision of the outcomes of simulation analysis is the number of grids.A considerable spatial discretization error is produced by coarse grids, which lowers the accuracy of the analysis's findings [41].Due to this, the selection of the optimum grid number is crucial.Based on the evaluation of various grid conditions, the grid independence analysis is a procedure used to identify the best grid condition that has the smallest grids without generating a difference in the numerical results.The grid independence analysis is done to select the size of the mesh.This helps to increase the processing speed and decreases the consumption of time by optimizing the grid size.Figure 4 shows the result of the grid independency analysis.
The outlet pressure of the evaporator tube remains steady after the number of elements is 108,600, which is depicted in Figure 5. From these, it is clear that after the specific limit, the simulation results are not dependent on the size of the grid.Hence, the number of elements selected for the study was 108,750.Due to the high accuracy, consistent results and less computational time, the standard k-ε model was utilized to solve the turbulent flow problem. Figure 4 shows the mesh structure of the test section on which the analysis was performed.Therein, 108,750 elements and 112,530 nodes make upthe whole test evaporator section.
To solve the two-phase mathematical model used to simulate the flow boiling in the evaporator tube, some hypotheses were taken into account.The flow field was consider as two-dimensional and turbulent.In addition, the heat flux was uniform along the length of the test section.The fluid was considered to be incompressible, and force vectors as negligible.

Grid Independency Analysis
In order to generate grids that perform the most effectively, it is necessary to take into account the grids' shape, quality, and number.A component that affects the overall processing cost and the precision of the outcomes of simulation analysis is the number of grids.A considerable spatial discretization error is produced by coarse grids, which lowers the accuracy of the analysis's findings [41].Due to this, the selection of the optimum grid number is crucial.Based on the evaluation of various grid conditions, the grid independence analysis is a procedure used to identify the best grid condition that has the smallest grids without generating a difference in the numerical results.The grid independence analysis is done to select the size of the mesh.This helps to increase the processing speed and decreases the consumption of time by optimizing the grid size.Figure 4 shows the result of the grid independency analysis.
The outlet pressure of the evaporator tube remains steady after the number of elements is 108,600, which is depicted in Figure 5. From these, it is clear that after the specific limit, the simulation results are not dependent on the size of the grid.Hence, the number of elements selected for the study was 108,750.

Grid Independency Analysis
In order to generate grids that perform the most effectively, it is necessary to take into account the grids' shape, quality, and number.A component that affects the overall processing cost and the precision of the outcomes of simulation analysis is the number of grids.A considerable spatial discretization error is produced by coarse grids, which lowers the accuracy of the analysis's findings [41].Due to this, the selection of the optimum grid number is crucial.Based on the evaluation of various grid conditions, the grid independence analysis is a procedure used to identify the best grid condition that has the smallest grids without generating a difference in the numerical results.The grid independence analysis is done to select the size of the mesh.This helps to increase the processing speed and decreases the consumption of time by optimizing the grid size.Figure 4 shows the result of the grid independency analysis.
The outlet pressure of the evaporator tube remains steady after the number of elements is 108,600, which is depicted in Figure 5. From these, it is clear that after the specific limit, the simulation results are not dependent on the size of the grid.Hence, the number of elements selected for the study was 108,750.

Governing Equations
The simulation analysis utilized the two-phase flow model and standard k-ε model in AnsysFluent software.The two-phase flow model is governed by the following equations: The n th phase continuity equation [42] is expressed as: The momentum equation [42] for the n th phase is: The n th phase energy equation [42] is given as: The standard Kk-ε model is based on the transport equations for turbulence kinetic energy and dissipation rate.This model takes less time for computation and is more reliable and accurate.The K (turbulence kinetic energy) is from the transport equation [42]as per the following: The 'ε' (dissipation rate) is identified from the transport equation as below:

Software Tools Validation
The results obtained from the experiment analysis conducted by Sun and Yang [43] are used for the validation of software tools.They determined the HTC of the CuO nanoparticles added to the R141b refrigerant.The average size of the nanoparticles for this

Governing Equations
The simulation analysis utilized the two-phase flow model and standard k-ε model in AnsysFluent software.The two-phase flow model is governed by the following equations: The n th phase continuity equation [42] is expressed as: The momentum equation [42] for the n th phase is: The n th phase energy equation [42] is given as: The standard Kk-ε model is based on the transport equations for turbulence kinetic energy and dissipation rate.This model takes less time for computation and is more reliable and accurate.The K (turbulence kinetic energy) is from the transport equation [42] as per the following: The 'ε' (dissipation rate) is identified from the transport equation as below:

Software Tools Validation
The results obtained from the experiment analysis conducted by Sun and Yang [43] are used for the validation of software tools.They determined the HTC of the CuO nanoparticles added to the R141b refrigerant.The average size of the nanoparticles for this study was 30 nm.The vapor quality ranged from 0.3 to 0.8, and the nanoparticle mass fraction ranged between 0.1 and 0.3% of the weight.They utilized a copper tube of 1 mm thickness, 10 mm inner diameter and 1400 mm length.Table 4 shows the variation of flow boiling heat transfer coefficient obtained by simulation and experiment analysis.The maximum deviation is observed at around 7%, which indicates the reliability of the present analysis.

Average Absolute Deviation Estimation
The average deviation of results indicates the level of prediction.It is the deviation of the results obtained from the simulation to the prediction analysis.The level of prediction is tabulated in Table 5.The AAD (average absolute deviation) is expressed by [44]: The performance of the refrigeration system influences the thermo-physical characteristics of the working fluid, such as thermal conductivity, density, specific heat, and viscosity.
The thermal conductivity of the nanorefrigerant is one of the key thermo-physical characteristics to enhance the thermal performance of the system.The heat transfer rate of the system greatly affects the thermal conductivity of the nanoparticles doped in the refrigerant [45].Figure 6 indicates the impact of temperature and nanoparticle volume concentration on the thermal conductivity of R1234yf-based nanorefrigerant.The figure shows that the thermal conductivity of the TiO 2 nanoparticle dispersed refrigerant rises with the concentration and reduces with the temperature.When 0.2% to 1% of the nanoparticles were added to the refrigerant, the thermal conductivity changed from 0.0717 to 0.0735 W/m 2 K at a temperature of 10 • C.However, when the temperature varied from 10 • C to 40 • C, it changed from 0.0735 to 0.0656 W/m 2 K at the 1% volume concentration.There is an 11.88% of augmentation in thermal conductivity obtained at the 1% volume fraction and 40 • C temperature, and when 0.2% of the nanoparticles was added to the R1234yf, a 0.6% enhancement was obtained at 10 • C. and 40°C temperature, and when 0.2% of the nanoparticles was added to the R1234yf, a 0.6% enhancement was obtained at 10°C.The thermal conductivity of the nanorefrigerant was enhanced with the doping of TiO2 nanoparticles.This is due to the Brownian motion effect of the nanometer-sized materials in the base working fluid.Photons are used to transmit heat by propagating lattice vibrations in crystalline nanoparticles immersed in the base fluid.This kind of phonon oscillates erratically.Some ballistic phonon events might lead to advances in thermal conductivity.If the ballistic phonons that originated in the particle travel through the fluid and reach neighboring particles, a significant increase in thermal conductivity occurs [46].Due to Brownian motion, heat can travel between particles more efficiently the closer they are to one another.Meanwhile, the thermal conductivity declines with the temperature.When the temperature is raised, the mean free path between the molecules also increases.This reduces the collision between them [47].While the separation between the nanoparticle grows, the Coulomb interaction (i.e., near-field radiation) degrades.This reduces the thermal conductivity [48].
Viscosity is another vital property that influences the pressure drop of the refrigeration system [49].Figure 7 depicts the variation in the viscosity of TiO2/R1234yf refrigerant with the temperature and the concentration of the nanoparticles in it.It is seen that the viscosity is increased with the addition of nanoparticles, and it decreases with the temperature.At the 0.2% concentration, the viscosity varied from 194.97 to 134.67µPas when the temperature ranged between 10 °C and 40 °C.Meanwhile, it changed from 194.97 to 198.94µPas with the 0.2% TiO2nanoparticles inclusion in the R1234yf base refrigerant at 10 °C.The highest increment in the viscosity is obtained at 1% at 2.54%, and the lowest at 0.5% of concentration at 0.5%.The high surface-to-volume area ratio of nanoparticles can be attributed as the main cause of the increase in viscosity with nanoparticles.Because nanoparticles have a larger surface area, there is more friction, which raises viscosity.On the other hand, it decreases when the temperature rises since it is hypothesized that the entropy increases; thereby, the thermal energy increases.High viscosity is not preferred for refrigeration systems, as is well known.It should be remembered that at higher temperatures, the enhanced viscosity with the nanoparticles will diminish and become closer to that of the base fluid.Therefore, for such systems, the slight increase in viscosity caused by the addition of nanoparticles may be disregarded [50].The thermal conductivity of the nanorefrigerant was enhanced with the doping of TiO 2 nanoparticles.This is due to the Brownian motion effect of the nanometer-sized materials in the base working fluid.Photons are used to transmit heat by propagating lattice vibrations in crystalline nanoparticles immersed in the base fluid.This kind of phonon oscillates erratically.Some ballistic phonon events might lead to advances in thermal conductivity.If the ballistic phonons that originated in the particle travel through the fluid and reach neighboring particles, a significant increase in thermal conductivity occurs [46].Due to Brownian motion, heat can travel between particles more efficiently the closer they are to one another.Meanwhile, the thermal conductivity declines with the temperature.When the temperature is raised, the mean free path between the molecules also increases.This reduces the collision between them [47].While the separation between the nanoparticle grows, the Coulomb interaction (i.e., near-field radiation) degrades.This reduces the thermal conductivity [48].
Viscosity is another vital property that influences the pressure drop of the refrigeration system [49].Figure 7 depicts the variation in the viscosity of TiO 2 /R1234yf refrigerant with the temperature and the concentration of the nanoparticles in it.It is seen that the viscosity is increased with the addition of nanoparticles, and it decreases with the temperature.At the 0.2% concentration, the viscosity varied from 194.97 to 134.67 µPas when the temperature ranged between 10 • C and 40 • C.Meanwhile, it changed from 194.97 to 198.94 µPas with the 0.2% TiO 2 nanoparticles inclusion in the R1234yf base refrigerant at 10 • C. The highest increment in the viscosity is obtained at 1% at 2.54%, and the lowest at 0.5% of concentration at 0.5%.The high surface-to-volume area ratio of nanoparticles can be attributed as the main cause of the increase in viscosity with nanoparticles.Because nanoparticles have a larger surface area, there is more friction, which raises viscosity.On the other hand, it decreases when the temperature rises since it is hypothesized that the entropy increases; thereby, the thermal energy increases.High viscosity is not preferred for refrigeration systems, as is well known.It should be remembered that at higher temperatures, the enhanced viscosity with the nanoparticles will diminish and become closer to that of the base fluid.Therefore, for such systems, the slight increase in viscosity caused by the addition of nanoparticles may be disregarded [50].Density influences the heat transfer properties of the refrigeration system.In ad tion, it affects the pressure drop of the system [50,51].Figure 8 indicates the deviatio the density of nano-enhanced refrigerant with various temperatures and concentrati of nanoparticles.It was observed that the density improved with the nanoparticle disp sion, and it declined with the temperature.At 10 °C temperature, the density TiO2/R1234yf refrigerant changed from 1150.2 to 1174.86 kg/m 3 when the concentratio nanoparticles ranged from 0.2 to 1% by volume.However, when the temperature va from 10 °C to 40 °C, it changed from 1150.2 to 1043.4 kg/m 3 at the 0.2% volume concen tion.The maximum enhancement in the density is 3.08%, which is obtained at 1% na particle addition and 40 °C temperature.The minimum enhancement is 0.54% at a volu concentration of 0.2% and temperature of 10 °C.The density of the nanorefrigerant is creased by the addition of nanoparticles to the base fluids.The increment in density the resulting nanorefrigerant will become increasingly apparent as the doping fractio the nanoparticles rises.It is caused by significant variations in the density between m rials in the solid and liquid phases.As the temperature rises, the volume of the fluid creases.Hence, the density of the nanorefrigerant decreases with the temperature [50,51].Density influences the heat transfer properties of the refrigeration system.In addition, it affects the pressure drop of the system [50,51].Figure 8   Density influences the heat transfer properties of the refrigeration system.In ad tion, it affects the pressure drop of the system [50,51].Figure 8 indicates the deviation the density of nano-enhanced refrigerant with various temperatures and concentrati of nanoparticles.It was observed that the density improved with the nanoparticle disp sion, and it declined with the temperature.At 10 °C temperature, the density TiO2/R1234yf refrigerant changed from 1150.2 to 1174.86 kg/m 3 when the concentration nanoparticles ranged from 0.2 to 1% by volume.However, when the temperature var from 10 °C to 40 °C, it changed from 1150.2 to 1043.4 kg/m 3 at the 0.2% volume concen tion.The maximum enhancement in the density is 3.08%, which is obtained at 1% na particle addition and 40 °C temperature.The minimum enhancement is 0.54% at a volu concentration of 0.2% and temperature of 10 °C.The density of the nanorefrigerant is creased by the addition of nanoparticles to the base fluids.The increment in density the resulting nanorefrigerant will become increasingly apparent as the doping fraction the nanoparticles rises.It is caused by significant variations in the density between m rials in the solid and liquid phases.As the temperature rises, the volume of the fluid creases.Hence, the density of the nanorefrigerant decreases with the temperature [50,51].One of the essential characteristics of the fluid is specific heat, which affects the r of heat transfer and the output temperature and so indirectly improves system per mance.Figure 9 shows the effect of temperature and concentration on the specific hea TiO2/R1234yf.From the figure, it is indicated that the specific heat capacity reduces w One of the essential characteristics of the fluid is specific heat, which affects the rate of heat transfer and the output temperature and so indirectly improves system perfor-mance.Figure 9 shows the effect of temperature and concentration on the specific heat of TiO 2 /R1234yf.From the figure, it is indicated that the specific heat capacity reduces with the rise in the concentration of nanoparticles in the refrigerant, and it improves with temperature increments.When the temperature ranged from 10 • C to 40 • C, the specific heat capacity varied from 1288.58 to 1430.95 J/kgK at the 0.2% concentration of TiO 2 nanoparticles in the R1234yf refrigerant.Meanwhile, at 10 • C, it changed from 1288.58 to 1271.36 J/kgK when the concentration ranged between 0.2 and 1%.The highest reduction in specific heat capacity was obtained at a 1% volume concentration of TiO 2 nanoparticles and a 40 • C temperature as 2.05%.The least reduction was at 0.34% at a temperature of 10 • C and the concentration of 0.2% by volume.the rise in the concentration of nanoparticles in the refrigerant, and it improves with temperature increments.When the temperature ranged from 10 °C to 40 °C, the specific heat capacity varied from 1288.58 to 1430.95 J/kgK at the 0.2% concentration of TiO2nanoparticles in the R1234yf refrigerant.Meanwhile, at 10 °C, it changed from 1288.58 to 1271.36 J/kgK when the concentration ranged between 0.2 and 1%.The highest reduction in specific heat capacity was obtained at a 1% volume concentration of TiO2nanoparticles and a 40 °C temperature as 2.05%.The least reduction was at 0.34% at a temperature of 10 °C and the concentration of 0.2% by volume.Generally, the specific heat of nanorefrigerant is reduced with increasing the concentration of nanoparticles.This is because the specific heat of nano-oxide materials is less, especially for TiO2 [52,53].However, the specific heat of nano-enhanced refrigerants improved with temperature.As the temperature rises, the liquid molecules are forced to fluctuate in their equilibrium value to a greater degree, resulting in a rise in the specific heat [54].The results indicated that the nano-enhanced refrigerants required less amount of heat to raise the temperature than the base refrigerant R1234yf.

Flow Patterns
The contours of velocity and pressure obtained by the simulation analysis are depicted in Figures 10 and 11, respectively.The flow pattern along the axial length of the test section is visualized at the heat flux of 4020W/m 2 , 1 vol.% of TiO2 nanoparticles concentration and 200 kg/m 2 s mass flux.Generally, the specific heat of nanorefrigerant is reduced with increasing the concentration of nanoparticles.This is because the specific heat of nano-oxide materials is less, especially for TiO 2 [52,53].However, the specific heat of nano-enhanced refrigerants improved with temperature.As the temperature rises, the liquid molecules are forced to fluctuate in their equilibrium value to a greater degree, resulting in a rise in the specific heat [54].The results indicated that the nano-enhanced refrigerants required less amount of heat to raise the temperature than the base refrigerant R1234yf.

Flow Patterns
The contours of velocity and pressure obtained by the simulation analysis are depicted in Figures 10 and 11, respectively.The flow pattern along the axial length of the test section is visualized at the heat flux of 4020 W/m 2 , 1 vol.% of TiO 2 nanoparticles concentration and 200 kg/m 2 s mass flux.
The velocity varied between 0 and 0.389 m/s.The maximum velocity was observed at the center, and the minimum velocity was at the walls of the test section.At the wall, the velocity is zero due to the no-slip condition.It was also observed that the velocity varied in the axial as well as the radial directions.
The pressure variation, along with the evaporation section, is depicted in Figure 11.It indicates the pressure decreased from 681 to 670 kPa.
of heat to raise the temperature than the base refrigerant R1234yf.

Flow Patterns
The contours of velocity and pressure obtained by the simulation analysis are depicted in Figures 10 and 11, respectively.The flow pattern along the axial length of the test section is visualized at the heat flux of 4020W/m 2 , 1 vol.% of TiO2 nanoparticles concentration and 200 kg/m 2 s mass flux.The velocity varied between 0 and 0.389 m/s.The maximum velocity was observed at the center, and the minimum velocity was at the walls of the test section.At the wall, the velocity is zero due to the no-slip condition.It was also observed that the velocity varied in the axial as well as the radial directions.The pressure variation, along with the evaporation section, is depicted in Figure 11.It indicates the pressure decreased from 681 to 670kPa.

Pressure Drop and Heat Transfer of TiO2/R1234yf Refrigerant
The variation of flow boiling heat transfer coefficient with vapor quality on different concentrations of nanoparticles in the refrigerant is illustrated in Figure 12.From the figure, it is seen that the HTC enhanced with the vapor quality and concentration of nanoparticles in the refrigerant.When the vapor quality change from 0.2 to 0.9, the HTC increases between 3024.56 and 7370.49W/m 2 K at 0.2% of volume concentration.Likewise, with the 1% TiO2 nanoparticles, the addition of the R1234yf refrigerant enhances the heat transfer coefficient from 6292.61 to 15,175.88W/m 2 K.When 0.2% of the nanoparticles dispersed in the refrigerant improved the coefficient by 13.35% on average, it also enhanced by 134.03% on average by mixing 1% of the TiO2 nanoparticles into the R1234yf refrigerant.The outcomes of the study indicate that the concentration of metal particles at the nanoscale affects how well the nanorefrigerant transfers heat up.The boundary layer's height is reduced, and the flow boiling HTC of the nanorefrigerant increased as a result of the molecular adsorption layer that forms on the surface of nanoparticles and their breakup [40].

Pressure Drop and Heat Transfer of TiO 2 /R1234yf Refrigerant
The variation of flow boiling heat transfer coefficient with vapor quality on different concentrations of nanoparticles in the refrigerant is illustrated in Figure 12.From the figure, it is seen that the HTC enhanced with the vapor quality and concentration of nanoparticles in the refrigerant.When the vapor quality change from 0.2 to 0.9, the HTC increases between 3024.56 and 7370.49W/m 2 K at 0.2% of volume concentration.Likewise, with the 1% TiO 2 nanoparticles, the addition of the R1234yf refrigerant enhances the heat transfer coefficient from 6292.61 to 15,175.88W/m 2 K.When 0.2% of the nanoparticles dispersed in the refrigerant improved the coefficient by 13.35% on average, it also enhanced by 134.03% on average by mixing 1% of the TiO 2 nanoparticles into the R1234yf refrigerant.The outcomes of the study indicate that the concentration of metal particles at the nanoscale affects how well the nanorefrigerant transfers heat up.The boundary layer's height is reduced, and the flow boiling HTC of the nanorefrigerant increased as a result of the molecular adsorption layer that forms on the surface of nanoparticles and their breakup [40].
Table 6 summarizes the deviation between the predictive and simulation result of flow boiling HTC for various vapor qualities at 0.2 and 1% of volume concentrations.From the simulation, the heat transfer coefficient was enhanced from 2443.86 to 6592.06 W/m 2 K and between 5699.2 and 13678.08W/m 2 K at 0.2 and 1% addition of TiO 2 nanoparticle on R1234yf refrigerant, respectively.The absolute average deviation is observed as 9.91%.Table 6 summarizes the deviation between the predictive and simulation result of flow boiling HTC for various vapor qualities at 0.2 and 1% of volume concentrations.From the simulation, the heat transfer coefficient was enhanced from 2443.86 to 6592.06 W/m 2 K and between 5699.2 and 13678.08W/m 2 K at 0.2 and 1% addition of TiO2 nanoparticle on R1234yf refrigerant, respectively.The absolute average deviation is observed as 9.91%.Figure 13 indicates the difference in pressure drop of TiO 2 /R1234yf nanorefrigerant with vapor quality for various nanoparticle concentrations.It is observed that the pressure drop rises with vapor quality and nanoparticles addition to the refrigerant.When a 0.2% volume fraction of the nanoparticles was added to the refrigerant, the pressure drop varied between 2.11 and 8.04 kPa at 0.2 and 0.9 vapor qualities, respectively.Likewise, the variation in the vapor qualities changes the pressure drop from 3.58 to 15.43 kPa at 1% of volume concentration.It was enhanced by 12.45% at 0.2% and by 80.77% at 1% by volume of TiO 2 concentration on average.Despite the fact that the mass flow rate of the heat transfer fluid was kept constant for the purposes of this investigation, the dispersion of nano-sized particles to the pure refrigerant improves the pressure decrease.When the nanoparticle concentration increases, there may be more particle-to-wall collisions, which results in a greater pressure drop for the nanorefrigerant than for the base refrigerant.The high-pressure drop of the nanorefrigerant has an effect on a refrigeration system's capacity to pump.As the nanoparticles loading rises, frictional pressure loss also rises, necessitating a significant amount of pumping energy to circulate the nanorefrigerant throughout the system.
transfer fluid was kept constant for the purposes of this investigation, the dispersion of nano-sized particles to the pure refrigerant improves the pressure decrease.When the nanoparticle concentration increases, there may be more particle-to-wall collisions, which results in a greater pressure drop for the nanorefrigerant than for the base refrigerant.The high-pressure drop of the nanorefrigerant has an effect on a refrigeration system's capacity to pump.As the nanoparticles loading rises, frictional pressure loss also rises, necessitating a significant amount of pumping energy to circulate the nanorefrigerant throughout the system.The variation between the simulation and predictive results of pressure drop is tabulated in Table 7. From the simulation, the pressure drop changes between 1.81 to 6.71 kPa, when the vapor quality ranges from 0.2 to 0.9 at 0.2% of nanoparticles concentration.In addition, it varies from 3.23 to 13.84 kPa at a 1% volume concentration of TiO2 nanoparticles on the R1234yf refrigerant.The AAD was obtained as 10.04%.

Conclusions
The HTC, pressure drop and the thermal characteristics of TiO2-doped R1234yf refrigerants with various concentrations (0.2% to 0.1% by volume) have been studied.The The variation between the simulation and predictive results of pressure drop is tabulated in Table 7. From the simulation, the pressure drop changes between 1.81 to 6.71 kPa, when the vapor quality ranges from 0.2 to 0.9 at 0.2% of nanoparticles concentration.In addition, it varies from 3.23 to 13.84 kPa at a 1% volume concentration of TiO 2 nanoparticles on the R1234yf refrigerant.The AAD was obtained as 10.04%.

Conclusions
The HTC, pressure drop and the thermal characteristics of TiO 2 -doped R1234yf refrigerants with various concentrations (0.2% to 0.1% by volume) have been studied.The temperature of the working fluid varied between 10 • C and 40 • C to identify the influence of the temperature on the properties.From the property analysis, the thermal conductivity, viscosity and density were enhanced by 7.49%, 2.54% and 2.87%, respectively, at the temperature of 25 • C and volume concentration of 1%.However, in the same conditions, the specific heat showed a 1.84% reduction.
A mathematical model and simulation approach were used to assess the nanorefrigerant's pressure drops and heat transfer characteristics.The investigation shows that the volume fraction of nanoparticles in the TiO 2 /R1234yf nano-enhanced refrigerant increases its HTC.At a 1% volume fraction of TiO 2 in R1234yf, the coefficient of heat transfer was enhanced by 134.03%.However, the presence of nanoparticles also resulted in an increase in the pressure drop.With a volume fraction of 1%, it is shown that there is an improvement of 80.77%.In the system, the improvement in heat transfer was more significant than the reduction in pressure.Therefore, adding nanoparticles to the refrigerant solution might improve the performance of the refrigeration system.The outcomes from the simulation technique and mathematical method differed by an average of 9.91%.
Numerical and simulation techniques were used in this study.The identification of nanorefrigerant performance and energy efficiency in actual cooling systems requires experimental verification, which is outside the focus of this work.Only lower concentrations and smaller sizes of nanoparticles are covered in the current research.Future analyses could expand to include higher concentrations, different nanoparticle shapes, and varied nanoparticle sizes added to the refrigerant.According to recent research, the TiO 2 /R1234yf nano-enhanced refrigerant has the potential to improve the performance of the system.It could be used to aid in human comfort as well as residential refrigeration systems and vehicular air conditioning systems.However, the performance of nano-enhanced refrigerants and their energy efficiencymust be experimentally demonstrated and must be compatible with real-time refrigeration systems.In conclusion, nano-enhanced refrigerants will soon be a key component of environmentally responsible and energy-efficient refrigeration systems.

Figure 6 .
Figure 6.Thermal conductivity of nano-enhanced refrigerant with temperature.

Figure 6 .
Figure 6.Thermal conductivity of nano-enhanced refrigerant with temperature.

Figure 7 .
Figure 7. Impact of temperature and volume concentration on the viscosity of TiO2/R1234yf.

Figure 8 .
Figure 8. Variation of volume fraction and temperature on the density of nanorefrigerant.

Figure 7 .
Figure 7. Impact of temperature and volume concentration on the viscosity of TiO 2 /R1234yf.
indicates the deviation in the density of nano-enhanced refrigerant with various temperatures and concentrations of nanoparticles.It was observed that the density improved with the nanoparticle dispersion, and it declined with the temperature.At 10 • C temperature, the density of TiO 2 /R1234yf refrigerant changed from 1150.2 to 1174.86 kg/m 3 when the concentration of nanoparticles ranged from 0.2 to 1% by volume.However, when the temperature varied from 10 • C to 40 • C, it changed from 1150.2 to 1043.4 kg/m 3 at the 0.2% volume concentration.The maximum enhancement in the density is 3.08%, which is obtained at 1% nanoparticle addition and 40 • C temperature.The minimum enhancement is 0.54% at a volume concentration of 0.2% and temperature of 10 • C. The density of the nanorefrigerant is increased by the addition of nanoparticles to the base fluids.The increment in density for the resulting nanorefrigerant will become increasingly apparent as the doping fraction of the nanoparticles rises.It is caused by significant variations in the density between materials in the solid and liquid phases.As the temperature rises, the volume of the fluid increases.Hence, the density of the nanorefrigerant decreases with the temperature rise[50,51].

Figure 7 .
Figure 7. Impact of temperature and volume concentration on the viscosity of TiO2/R1234yf.

Figure 8 .
Figure 8. Variation of volume fraction and temperature on the density of nanorefrigerant.

Figure 8 .
Figure 8. Variation of volume fraction and temperature on the density of nanorefrigerant.

Figure 9 .
Figure 9.Effect of specific heat capacity by temperature.

Figure 9 .
Figure 9.Effect of specific heat capacity by temperature.

Figure 13 .
Figure 13.Variation of pressure drops with vapor quality.

Figure 13 .
Figure 13.Variation of pressure drops with vapor quality.

Table 3 .
Summary of validation analysis.

Table 4 .
Validation of software tool.

Table 5 .
Level of prediction.

Table 6 .
Comparison of simulation and numerical results.

Table 6 .
Comparison of simulation and numerical results.

Table 7 .
Deviation of the simulation and predictive results.

Table 7 .
Deviation of the simulation and predictive results.