Numerical Performance Investigation of Parabolic Dish Solar- Assisted Cogeneration Plant Using Different Heat Transfer Fluids

Institute of Energy and Power Engineering, Zhejiang University of Technology, Hangzhou, Zhejiang 310014, China Department of Mechanical Engineering, Mirpur University of Science and Technology (MUST), Mirpur 10250 (AJK), Pakistan Department of Energy Systems Engineering, Faculty of Integrated Technologies, Universiti Brunei Darussalam, Jalan Tungku Link BE 1410, Bandar Seri Begawan, Brunei Darussalam Mechanical Engineering Department, King Fahd University of Petroleum & Minerals (KFUPM), Dhahran 31261, Saudi Arabia Department of Mechanical Engineering, University of Engineering and Technology, Taxila 47050, Pakistan


Introduction
Renewable energy resources are considered the best alternative solutions to address the severe ecological problems created due to the consumption of fossil fuels [1]. Although there are various types of renewable energy resources (geothermal, solar, and wind etc.), solar energy is a unique choice as it is available abundantly [2]. Concentrating solar systems, also known as concentrated solar power (CSP) technologies, gain more interest due to their useful heat production at elevated temperatures. Particularly, dish solar collectors are an efficient source that converts solar energy into useful thermal energy and also attain maximum thermal efficiency due to higher concentration ratio [3]. Its operation relies on the concentration of solar radiations onto a small receiver where it exchanges heat with the working fluid. The application of the cavity receiver has been noticed to be a very useful technique as solar radiations are properly centered on parabolic solar dishes [4]. The useful heat gain from solar dish collectors can be used as an input source for dryers, hydrogen production systems, and Rankine and Brayton cycles [5,6].
The literature consists of both numerical and experimental works regarding parabolic dish cavity receivers. Daabo et al. [7] performed numerical analysis of various parabolic dish cavity receivers to investigate their optical performance. Results indicated that the conical shaped receiver received and absorbed maximum solar radiations per unit area as compared to the others with optical efficiency near to 75.3%. The optimization of the receiver employing optimal methods and CFDs was conducted to examine the most suitable position of the solar receiver on a concentrator collector by Przenzak et al. [8]. The optimized value of fluid flow was calculated to be 0.6 m/s without any overheating of the receiver body, and it also decreased the convection and radiation heat losses. Jilte et al. [9] investigated the dish collector with different cavity shapes mathematically, and research was performed for five inclinations and three different isothermal wall temperatures. The conical cavity receiver showed the lowest convective heat losses as compared to the other investigated receivers. Kaushika and Reddy [10] suggested a method for design and development of a lowcost steam generation system integrated with a parabolic dish (PD) system. Solar to steam conversion efficiency was to found be 70 to 80% at 723 K. Parabolic dish using cavity receiver was numerically modeled and thermally optimized by Loni et al. [11].
Experimental work has been performed to analyze the energetic and exergetic performances of the SK-14 cylindrical cavity solar dish system [12]. Heat loss factor was observed to be 4.6 W/K with 52% optical efficiency under higher solar insulations. Mohammad [13] designed and developed a parabolic dish solar (PDS) water heater for domestic usage. Thermal efficiencies were between 52% and 56% that showed the reliability of the designed system. Reddy et al. [14] carried out an experimental work on a 20 m 2 prototype fuzzy focal dish receiver. Madadi et al. [15] theoretically and experimentally analyzed the energetic and exergetic performances of PDSC. Two different cavity receivers, namely, cylindrical and conical, were taken into account with different parameters. It was observed that almost 35% to 60% of the exergy destroyed was due to the heat transfer from the sun to the receiver. The cylindrical cavity receiver showed better results than the counterpart did.
Furthermore, utilization of heat transfer fluid (HTF) has an immense impact on the efficiency of solar collectors. The most used HTFs in concentrated solar power (CSP) technologies are oil-based, while molten salts are also tested by researchers [16,17]. Loni et al. [18] conducted a study using different cavity receivers of PDC utilizing water and Behran oil as heat transfer fluids. Another study was conducted by Loni et al. [19] using air as a HTF in a dish collector integrated with organic Rankine cycle (ORC). Khan et al. [20] examined numerically a PD cavity receiver with three different nanofluids, and the dish system was coupled with sCO 2 Brayton cycle to generate electricity. Abid et al. [21] performed an analysis of different sCO 2 recompression Brayton cycles. The heat source was a parabolic dish solar collector using pressurized water. Abid et al. [22] performed another study using different nanofluids and molten salts in parabolic trough and parabolic dish collectors to power a steam cycle.
The other HTFs utilized in the solar collector are the gasses (air, helium, nitrogen, and CO 2 ) and argon [23]. Supercritical carbon dioxide has attracted a lot of attention due to its extraordinary thermophysical properties. Its critical point is quite low (305 K and 7.8 MPa) that assists properties to enhance in the region adjacent to the triple point. Numerous studies proposed and recommended sCO 2 as a heat transfer medium for higher temperature solar thermal applications [24][25][26][27].
Furthermore, the use of nanofluids in the solar collector enhances the rate of heat transfer that also raises the convective heat transfer coefficient, useful heat gain, and thermal efficiency of the solar collector. Mebarek-Oudina [28] numerically investigated laminar stationary natural convection heat transfer of different nanofluids contained in a cylindrical annulus using a discrete heat source. The objective was to investigate the effect of volume fraction, Rayleigh number, and base fluid type on the heat transfer rate. It was concluded that Raileigh number and volume fraction of nanoparticles has an influence on the thermal efficiency and growth rate of the heat transfer. Mebarek-Oudina [29] conducted a numerical investigation of natural convection heat transfer stability in vertical annulus using different heat source lengths. It was observed that heat transfer rates for smaller heater length are higher; however, heat transfer rate decreases by increasing the heater length. Kherbeet et al. [30] reported a numerical study of nanofluid flow and heat transfer of laminar convection flow adjacent to a horizontal and three-dimensional microscale forward facing step. It showed that the forward facing step increased the heat transfer, while Nusselt number also enhanced with rise in step height. Abbassi et al. [31] carried out a parametric investigation of a nanofluid-filled incinerator with a rectangular hot block on the bottom to evaluate the influence of different parameters (nanoparticle volume fraction, Rayleigh number, external magnetic intensity, etc.). It was concluded that increasing the width and height of the heater also enhanced the entropy generation but it has been reduced by increasing the Hartmann number. Goodarazi et al. [32] performed a study on the non-Newtonian nanofluid flow of carboxymethyl cellulose-aluminum oxide in a microtube. The aim was to explore the presence of nanoparticles and phenomenon of slip and temperature jump. Olia et al. [33] stated the modern developments on the applications of nanofluids for heat transfer fluid in parabolic trough solar collectors. It was concluded that nanofluids not only enhance the exergy efficiency, energy efficiency, and convective heat transfer coefficient but also mitigate the entropy generation of the system.
Gheynani et al. [34] examined the influence of nanoparticle concentration and diameter on the different heat transfer properties of non-Newtonian carboxymethyl cellulose/CuO fluid in a microtube using numerical simulation with the finite volume method. Taner [35] presented a bulgur drying plant based on the bulgur mass, temperature, and moisture for optimum production. Uncertainty analysis was carried out for data accuracy. The drying process energy and exergy efficiencies were found to be 24.95% and 44.05%, respectively. Taner [36] carried out an experimental optimization of a PEM fuel cell to enhance the efficiency and development of the simulations and modeling of PEM fuel cell. The first and second law efficiencies were almost 47.6% and 50.4%, respectively.
Topal et al. [37] performed a thermodynamic analysis of a trigeneration system that converted a single-fuel source into power, heating, and cooling and concentrated on the 2 International Journal of Photoenergy simulation of such systems to a direct cocombustion of poultry wastes. It was found that cocombustion of poultry waste was the best environment-friendly remedy to put away wastes. Topal et al. [38] conducted a thermodynamic analysis of a can circulating fluidized bed power plant cofired with olive pits and to check the effect on exergy destruction and carbon dioxide emissions. Total exergy destruction of the plant was found to be 295 MW with 31.26% exergy efficiency. Taner and Sivrioglu [39] developed a general model for technoeconomic and cost analysis of a turbine power plant. Simple payback period and unit cost for turbine power plant were calculated as 4.32 years and 3.142 $/kW, respectively. It is presented that there are many studies about the performance investigation of parabolic dish solar collector but there is lack of studies where both liquid and gas phase substances are evaluated simultaneously for their first and second law analyses. The aim of the present research is to examine in detail the energy and exergy efficiencies of PDSC using therminol VP1, pressurized water, and sCO 2 as heat transfer fluids. The thermal energy produced by the collector is further utilized to power a regenerative steam cycle for the purpose of electricity generation and rate of process heat. Therminol VP1 is a eutectic mixture of diphenyl oxide (DPO) and biphenyl. It can be used as a liquid heat transfer fluid or as a boiling-condensing heat transfer medium up to its maximum temperature. Various inlet temperature levels and mass flow rates are investigated to evaluate the work output, process heat rate, integrated efficiencies, and collector efficiencies of the solar integrated cogeneration plant. Different combinations of temperature and flow rate are tested for each heat transfer fluid, and at the end, a unique method of solar efficiency map is presented to depict the efficiencies of the solar receiver. The map shows the optimal conditions for the three HTFs exergetically and thermally with different operating temperature levels.
1.1. System Description. The proposed parabolic solar dishintegrated Rankine cycle producing electricity and rate of process heat is presented in Figure 1. Solar radiations from the sun and incidences on the aperture are reflected to the cavity receiver. The HTFs circulating inside the receiver takes the collected heat from the receiver. The high-temperature fluid at about 650 K moves to the boiler of steam cycle (point 9). The heated fluid transfers its heat to the water entering the boiler (point 5), circulating to the collector at relatively low temperature (point 10). Steam at higher pressure and higher temperature expands after passing through the turbine of a cogeneration plant (point 6). A small fraction (1/4 th ) of steam is extracted from the turbine (point 7) for process heating, while the remaining enters into the turbine to be expanded in the turbine (state point 8). The state of the liquid at the exit of the condenser is saturated, is pressurized using pump 1, and enters into the feed water heater (FWH) at higher pressure. A fraction of the steam after process heater (point 3) merges with the pressurized water in the feed water heater (FWH) and the mixture is pumped again towards the boiler for further process.

Methodology and Thermodynamic Analysis
The present section of the work is related to the mathematical equations, used for the thermodynamic modeling of PD solar receiver and regenerative Rankine cycle. Table 1 shows the input design parameters for the above said systems. A validated numerical mathematical model is developed with the help of the engineering equation solver (EES) [40]. The parabolic dish system presented in the present research is taken from Abid et al. [41], whereas Rankine cycle is adopted from [42]. The complete detail of both the systems can be found in the respective references.

First and Second Law Analyses of Solar Parabolic Dish
Collector. The performance of the solar collector has an important role, and it is closely associated with the useful heat gain by the collector. The useful energy can be calculated using the following two equations: Amount of the solar energy available to the collector can be determined by The heat transfer coefficient (h) between the working fluid in the solar collector and the receiver tube in equation (2) can be calculated using Nusselt number: Reynolds number is given by.
Gnielinski's correlation [43] for turbulent flow (3000 ≤ R e ≤ 5 × 10 6 and 0:5 ≤ P r ≤ 2000) is The Prandtl number [44] and friction factor are evaluated as The energy efficiency of the dish collector is given by 3 International Journal of Photoenergy Another crucial parameter that is required to be determined is the output exergy of the solar receiver that is presented in equation (9). Exergy is the maximum possible work available from any process [45]. Equation (9) considers different parameters such as useful heat available, pressure losses, mass flow rate, and temperature [46].
The Petela model is the most established model to calculate the rate of solar exergy.
Finally, exergy efficiency is found using exergy ratio between useful exergy and exergy of the sun.
The receiver's exergy destruction can be found using the below relation.  International Journal of Photoenergy The heat gain of solar system is determined by the famous Hottel-Whillier equation [47].
whereas heat removal factor and collector efficiency factor are summarized as

Energetic and Exergetic Analyses of Rankine Cycle.
Enthalpy is one of the most important parameters for steam cycle analysis. It is necessary to first calculate the enthalpy at each state of the system to calculate the rest of the output parameters. The efficiency of turbines is given as Work consumed by the pumps, to increase and circulate the pressure of working fluid, is given as Work done by the turbine is given as Heat input rate to the boiler and heat rejected by the condenser can be determined as [41] The thermal energy for the process heater can be evaluated as The net power output of the cycle can be computed as The thermal or energy efficiency of the steam cycle is     International Journal of Photoenergy determined as Exergy values at all the points of the examined system are required to be determined in order to calculate the exergy destruction rate. Equation (22) gives the relation for exergy at relevant points.
The exergy supplied to the boiler and condenser is mea-sured as Finally, exergy efficiency of the steam cycle and systemintegrated efficiencies are given as

Results and Discussion
This study concentrates on the integration of the parabolic dish solar collector utilizing three different HTFs, namely, sCO 2 , therminol VP1, and pressurized water with regenerative Rankine cycle, generating electricity and process heat. The present section discusses the findings of the investigated system in details. The change in mass flow rate at different inlet temperature levels of three heat transfer fluids in the solar collector is examined. The performance parameters considered in the study are overall energetic and exergetic efficiencies, power output and rate of process heat, receiver thermal and exergy efficiencies, and convective heat transfer coefficient. The steam Rankine cycle produces 875 kW of electricity and 570 kW heat for processing, and its thermal and exergy efficiencies are found to be 46.7% and 76.3%, respectively. Table 2 below shows the useful heat gain by the solar collector, net power produced, and process heat rate for the examined HTFs at three different inlet temperature levels. It is evident from the table that for lower inlet temperature, water is a better working fluid; however its effectiveness becomes less at higher temperatures. The performance of therminol VP1 gets better at elevated temperatures of above 550 K, while sCO 2 also performs well at upper temperature levels. Thermodynamic properties of regenerative Rankine cycle are presented in Table 3. Figure 2 illustrates the variation of system-integrated efficiencies for a range of inlet temperature from 350 K to 650 K. These graphs can be explained in many ways. It can be seen that higher flow rate and lower temperature gives maximum overall energy and exergy efficiencies for the investigated HTFs. However, for a mass flow rate of above 0.15 kg/sec, the variation in the properties becomes stagnant, meaning that there is no reason for employing higher mass flow rate from a thermal point of view. According to equation (25), overall energy efficiency of the integrated system is the ratio of work output from turbine and process heat rate of the pro-cess heater to the energy of the solar. Furthermore, higher value of the numerator depends on higher amount of useful heat by the solar collector according to equation (1). Furthermore, for all the inlet temperature levels, integrated efficiencies of pressurized water are the highest as compared to the other investigated fluids. The maximum difference in efficiency values of HTFs is seen at nearly 0.01 kg/sec, but after that, this difference reduces continuously. At a mass flow rate of 0.3 kg/sec, the efficiency of all   9 International Journal of Photoenergy the investigated fluids converges on the same point. At inlet temperature of 650 K, efficiency of therminol VP1 approaches very near to water values because the performance of water becomes less efficient after 600 K [24]. Pressurized water is witnessed to have 3.85% and 10% higher efficiency than the therminol VP1-and sCO 2 -based HTFs, respectively, at a mass flow rate of 0.01 kg/sec. However, these values are reduced to 0.11% and 0.31% at higher mass flow rate of 0.3 kg/sec. Figure 3 is related to the variation of mass flow rate and its influence on net power produced and rate of process heat at four different inlet temperature levels. Both the performance parameters increase sharply between 0.01 and 0.15 kg/sec and after that remain constant until 0.3 kg/sec. It is evident from Figure 3 that higher values are obtained at lower inlet temperature, and pressurized water is observed to have the highest rate of process heat and work output than the other two investigated HTFs. Net power output utilizing pressurized water in the collector loop is increased to 2.44% for a range of flow rate from 0.01 and 0.15 kg/sec; however, after 0.15 to 0.3 kg/sec, the increment is found to be as low as 0.17%. Therminol VP1 and sCO 2 has 6% and 11.9% increase for 0.01 to 0.15 kg/sec but very less increment in work output (0.17% and 0.51%, respectively). Further rise in mass flow rate clearly indicates that after a certain point, increase in mass flow rate has no importance for thermal enhancement. The highest variation is seen for sCO 2 values at all inlet temperatures due to its behavior near critical region [20]. Highest amount of work produced and rate of process heat is obtained for pressurized water at 350 K that is 1173 kW and 764 kW, respectively. For higher inlet temperature of 650 K, almost 19% reduction is found for both performance parameters. Figure 4 depicts the influence of the mass flow rate on the thermal and exergetic efficiencies of the parabolic dish receiver at four inlet temperature levels. Thermal efficiency of all investigated fluids has the same trend as discussed above for other performance parameters but exergy efficiency behaves in an opposite way with the rise in mass flow rate. PD collector with water as a HTF has the highest thermal efficiency (74.43%) at 350 K and at a mass flow rate of 0.3 kg/sec, while lowest value has been found for sCO 2 (74.2%). From the exergetic point of view, sCO 2 has a dominant feature as compared to the other HTFs and highest exergy efficiency (37.56%) is noticed at higher temperature (650 K) and lower flow rate (0.01 kg/sec). The highest exergy efficiency value of sCO 2 is due to its lower specific heat capacity as compared to therminol VP1 and pressurized water [23]. Exergy efficiency is the ratio of collector useful exergy to the exergy available from solar [22] and higher inlet temperature of fluid in the receiver will deliver maximum exergy but solar exergy is independent of the inlet temperature and remains the same for all the values of the temperature. The collector exergy depends on useful heat gain from collector, mass flow rate, and temperatures of fluid as followed by equation (9).
Figures 5(a)-5(d) represent the simulation results of integrated system efficiencies, work output, and process heat rate and collect performance as well as the heat transfer coefficient with respect to the inlet temperature. The operational parameterðT in − T 0 Þ/G b is a parameter expressing the collector efficiency. It reflects the inlet temperature of HTFs as direct normal irradiation, and ambient temperature remained constant during the simulation. The obtained results are considered at ambient temperature of 300 K, DNI of 1000 W/m 2 , and flow rate of 0.02 kg/sec. The inlet temperature increases between 350 K and 650 K with a step 10 International Journal of Photoenergy 50 K. Higher values of inlet temperature offers greater temperature of the receiver surface that results in maximum heat losses from the receiver surface. Higher heat losses decreases the useful heat production from the collector which eventually results in lowering the net power produced and process heat rate of the system. In Figure 5(a), overall energy and exergy efficiencies of the integrated system employing water reduce from 34.4% to 28.5% and 37% to 30.6%, respectively, and it shows reduction of almost 20.8%. The same trend is noticed for the therminol VP1-and sCO 2 -based integrated systems with decrease of nearly 17.93% and 17.87%, accordingly. Pressurized water was found to be 1.95% and 4.95% efficient than the systems operated with therminol VP1 and sCO 2 , respectively, at 350 K. However, after 500 K, the performance of the system using water as HTF reduces, and at 650 K, therminol VP1 was found to be 0.52% better than water. Similarly, Figure 5(b) has a very similar trend as described above; therminol VP1 has a dominant role at elevated temperatures. Collector thermal efficiency presented in Figure 5(c) has a very similar behavior to the work output and integrated efficiencies, i.e., higher thermal efficiency is found at lower temperature of 350 K. However, for exergy efficiency, sCO 2 -based collector has the highest value among the other compared HTFs. Maximum value is recorded for sCO 2 , therminol VP1, and water with 37.37%, 36.86%, and 36.98%, respectively, at 650 K. Figure 5(d) shows the graphical representation of the convective heat transfer coefficient of the absorber and investigated HTFs. It is clear that sCO 2 has the lowest value of h as compared to the other HTFs. The heat transfer coefficient (h) of therminol VP1 is more stable at all the temperature levels, while for water, h is constant at lower temperatures and after 500 K, and it drops dramatically due to evaporation.

11
International Journal of Photoenergy Figure 6(a) is the portrayal of density and heat capacity of three investigated fluids. The C p values of pressurized water have an upward trend at low temperatures; however, a sudden decrease is observed after 550 K, while sCO 2 has shown an increment of nearly 23.3%. The synthetic oil is found to have better thermophysical properties, as C p of the therminol increases by 48.5 with a rise in temperature. A decrease in the density of all HTFs is found with the rise in temperature. Maximum value of density can be seen for oil at lower temperature, and a constant fall is found for sCO 2 as well, while water depicts a dramatic change after 550 K [41]. Thermal conductivity values in Figure 6(b) shows the similar tendency to that of viscosity; maximum k value can be seen at 400 K for water and a reduction can be observed afterwards. Thermal conductivity of oil reduces to 61.5%, while a small increase for sCO 2 is noticed.
The following section provides the optimum operating conditions for the investigated heat transfer fluids in solar dish collector at various inlet temperature levels. The combinations of mass flow rate and inlet temperature of the working fluid in the receiver that lead to higher thermal and exergy efficiency will be assumed the optimum conditions. Table 4 and Figure 7 consist of the outcomes for optimal values against all levels of inlet temperatures. For lowtemperature levels of up to 450 K, the optimum flow rate is 0.07 kg/sec and corresponding efficiency values are tabulated in front of them. At 550 K, the optimum value of the mass flow rate can be seen to be 0.20 kg/sec, while higher inlet temperatures have optimum features at 0.30 kg/sec. The highest value of the exergetic efficiency for all the three investigated fluids is obtained at 750 K with flow rate of 0.30 kg/sec; however, a higher value of thermal efficiency is witnessed at 350 K and 0.07 kg/sec. This table also proves that the highest exergy performance does not associate with the higher thermal efficiency, and hence, the optimal operating scenario can be selected by taking both energy and exergy into consideration. The efficiency map of the collector is very useful and is an innovative tool to evaluate the collector performance under different criteria. Furthermore, validation of the parabolic dish receiver has been conducted by comparing the thermal efficiency of the collector using Al 2 O 3 /thermal oil as a HTF with the model presented by Loni et al. [48] in Figure 8. The proposed model shows a discrepancy of 0.85% in efficiency in comparison to the reference model when evaluated against inlet temperature.

Conclusions
Thermodynamic performance investigation of PD solarassisted cogeneration plant is presented in this study. Three working fluids/HTFs, namely, sCO 2 , pressurized water, and therminol VP1, are used to collect the heat produced by the solar collector. The collected heat is further used to drive the regenerative steam Rankine cycle producing power and process heat simultaneously. Firstly, PDSC is thermodynamically modeled by considering three HTFs (sCO 2 , pressurized water, and therminol VP1) and those modeling results are further utilized to evaluate the performance of integrated system. Different combinations of flow rates and inlet temperatures have been tested and an optimal scenario is presented that leads to maximum energy and exergy efficiencies of the solar collector. The major outcomes are presented as follows: (i) Higher flow rate and lower temperature offer maximum overall energy and exergy efficiencies for the investigated HTFs. Pressurized water is observed to have an integrated efficiency of 10% and 3.85% higher than sCO 2 -and therminol VP1-based HTFs, respectively, at a flow rate of 0.01 kg/sec. However, these values are reduced to 0.11% and 0.31% at a higher flow rate of 0.3 kg/sec (ii) Net power output of the proposed integrated system for water as a HTF in the collector loop is increased to 2.44% for flow rate from 0.01 to 0.15 kg/sec; however, after 0.15 to 0.3 kg/sec, this increment is as low as 0.17%. Therminol VP1 and sCO 2 has 6% and 11.9% increase for flow rate between 0.01 and 0.15 kg/sec but very less increment in work output values (0.17% and 0.51%, respectively) for further rise in the mass flow rate. Highest amount of net power produced (1173 kW) and the rate of process heat (764 kW) are obtained for pressurized water at 350 K, respectively. For higher inlet temperature of 650 K, almost 19% reduction is found for both performance parameters (iii) The PD collector with water as a HTF has the highest thermal efficiency (74.43%) at 350 K and at a flow rate of 0.3 kg/sec, while lowest value has been found for sCO 2

12
International Journal of Photoenergy (iv) Pressurized water was found to be 1.95% and 4.95% more efficient than the systems operating with therminol VP1 and sCO 2 , working fluids at a temperature of 350 K. However, the performance of the system using water as HTF starts to decrease with rise in temperature, and at 650 K, therminol VP1 is observed to surpass water with 0.52% better efficiency (v) An innovative solar efficiency map is proposed with an aim to find out the optimal operating region with respect to collector energy and exergy efficiencies. Greek Letters η pump : Pump efficiency (%) η tur : Turbine efficiency (%) ρ: Density (kg/m 3 ) η pet : Petela efficiency η en : Energy efficiency (%) _ ψ : Rate of exergy destruction (kW).

Data Availability
Data will be provided on request.

Conflicts of Interest
The authors declare that they have no conflicts of interest.