Polynomial Expressions for the Thermal Efficiency of the Parabolic Trough Solar Collector

The parabolic trough solar collector (PTC) is the most mature solar concentrating technology, and this technology is applied in numerous thermal applications. Usually, the thermal efficiency of the PTC is expressed with the aid of polynomial expressions. However, there is not a universal expression that is applied in all cases with high accuracy. Many studies use expressions with the first-degree polynomial, second-degree, or fourth-degree polynomial expressions. In this direction, this work is a study that investigates different expressions about the thermal efficiency of a PTC with a systematic approach. The LS-2 PTC module is examined with a developed numerical model in the Engineering Equation Solver for different operating temperatures and solar beam irradiation levels. This model is validated using experimental literature data. The found data are approximated with various polynomial expressions with up to six unknown parameters in every case. In every case, the mean absolute percentage error and the R2 are calculated. According to the final results, the use of the third power term leads to the best fitting results, as well as the use of the temperature difference term (ΔΤ), something that is new according to the existing literature. More specifically, the final suggested formula has the following format: “ηcol = a0 + a3∙ΔT3/Gb + b∙ΔΤ”. The results of this work can be used by the scientists for the optimum fitting of the PTC efficiency curves and for applying the best formulas in performance determination studies.


Introduction
The parabolic trough solar collector (PTC) is one of the most mature solar concentrating technologies [1] which can be used in numerous applications such as electricity production [2], industrial heat [3], refrigeration [4], space-heating [5], chemical processes [6] and polygeneration [7]. This collector is a linear concentrating technology that operates with a single-axis tracking system in the majority of the cases. There are only a few investigations with two-axis tracking systems but these designs are complex and they have a high cost. The most usual receiver is an evacuated tube absorber in order to operate in high temperatures with a satisfying efficiency [8]. The operation with thermal oil (e.g., Therminol VP-1) is able to reach up to 400 °C [9], while with the proper molten salts (e.g., nitrate salt), there is the possibility for operation up to 600 °C [10]. The use of gas working fluids (e.g., air or CO2) gives the possibility to operate in higher temperatures up to 1000 °C [11].
The literature includes numerous studies that incorporate PTC as the main or the secondary thermal energy source. Usually, the simulation of the PTC performance is conducted by using a thermal efficiency expression formula. The thermal efficiency of the solar collector is the ratio of the useful heat production to the available solar irradiation on the collector aperture (ηcol = Qu / Qs). Numerous formulas in the literature can be found and so it is interesting to present them as below.
A great percentage of the literature uses a simple linear expression for the determination of the efficiency. For example, the study of Kalogirou [12] approximates the PTC performance with a linear expression and it compares it with other linear formulas from the literature. This idea has also been used in newer studies such as in Coccia et al. [13]. Below, the linear expression is given: The previous formula is usually used for the flat plate systems according to Duffie and Beckman [14]. Thus, this linear expression can be an accurate choice for PTC systems that operate to low-temperature levels up to 100 °C. Moreover, the previous formula is also known as the Hottel-Whillier-Bliss equation.
The next part of the literature studies uses a second-order expression which is in accordance with the ISO 9806-1 [15]. The goal of this formula is to better approximate the operating points of higher temperature levels.
This formula has been used by a great part of the recent literature studies. For example, Bellos et al. [16] selected an expression as the previous one for approximating the Eurotrough module performance. In another older work, Kalogirou and Panayiotou [17] conducted experimental work with PTC and they compared the approximation of the efficiency with a linear and a second-term formula. They concluded that both expressions are accurate and the second-order formula can be applied for the cases with higher temperatures.
An alternative expression of the previous one has been introduced by Dudley et al. [18] in an experimental work about LS-2. They tried to add an extra term with (Tin − Tam) without the incorporation of the solar beam irradiation: Kutscher et al. [19] used a third-degree polynomial in order to approximate the PTC thermal efficiency and they found high values for the R 2 . More specifically, they tested the SkyTrough module with a Schott PTR-80 receiver.
They also tried to approximate the results with an alternative formula, as below: A complex polynomial expression has been suggested by Nems and Kasperski [20] which incorporates the wind speed (Vw) as an extra parameter which includes a third-order term (Gb • Tin • Vw), as below: The next part of the literature studies includes fourth-order terms. These formulas try to take into account the radiation thermal losses in a proper way because they are dependent on the fourth power of the absorber temperature level. Sallaberry et al. [21] tried to approximate the thermal losses with an expression with a fourth power and then the following expression for the thermal efficiency has been found: They stated that this expression is in accordance with the standard IEC 62727 [22]. The use of the "cos (θ)" in the denominator is justified due to the approximation of the thermal losses for specific heat input. Practically, the thermal losses have to be normalized in every case in the absorbed irradiation, and the use of the cosine is a good choice for the approximation of the incident angle modifier. The thermal efficiency was calculated by reducing the thermal losses from the absorbed heat in every case.
Odeh et al. [23] have used a similar formula but they incorporated the wind speed inside the efficiency expression: One alternative formula has been suggested also by Sallaberry et al. [24] which includes a fourth-order term and an order which takes into account the transient phenomena. The authors have expressed this formula in terms of "Qu/Aa", but below, it is presented in terms of efficiency for having a uniform presentation style in the paper: In another work, Sallaberry et al. [25] used a formula with an extra term, as below: The previous formula is associated with the IEC 62862-3-2 standard which also suggests the use of a second-order term instead of the fourth-order term.
More specifically, the previous formula has been found by eliminating terms from the general formula according to DIS/ISO 9806 (2017) [26], which indicates approximately the following: where the "EL" is the longwave irradiance with wavelength over 3 μm. Moreover, it has to be said that Sallaberry et al. [25] compared Equation (10) with Equation (12). They stated that formula (10) is slightly better than Formula (12) but the difference is negligible.
Xu et al. [27] used Formula (11) with all the terms in order to approximate the PTC performance. Furthermore, it has to be said that there are studies that tried to determine some analytical expressions for the PTC efficiency curves [28][29].
The previous literature review indicates that there is a huge interest in the literature in the approximation of the PTC thermal efficiency with effective expressions. Numerous formulas have been suggested and applied in literature studies. Moreover, there are some studies that compare different efficiency approximation equations. So, it is obvious that the researchers use various ways to approximate the PTC efficiency and there is not only one formula to expresses the thermal efficiency of a PTC. In this direction, this work comes to investigate different polynomial expressions that can be applied for the PTC thermal efficiency determination and they include the usual literature expressions. Moreover, these equations have extended in order to cover all the possible combinations of polynomial terms in a systematic way. Totally, 31 different equations are applied and evaluated according to two criteria: the mean absolute percentage error (MAPE) and the R 2 index. To our knowledge, there is no other study in the literature like the present study, so this work adds something new to the existing literature. The results are discussed and the most important polynomial approximating formulas are presented. The final conclusions of this work indicate the most effective polynomial expressions for approximating the PTC efficiency by using the operating temperature levels and the solar irradiation, and also a clear comparison among the various polynomial formulas is given. Lastly, it is important to state that the operating points for the approximations are found by a developed model in the Engineering Equation Solver (EES) [30], which is validated with numerical experimental results.

The Examined Collector
In this work, the module of the LS-2 PTC [18,[31][32] is used in order to apply the different performance efficiency approximation formulas. This PTC module is usually studied in the literature studies and it is illustrated in Figure 1a. Moreover, the cross-section of the receiver (cover and absorber) is given in Figure 1b Table 1 includes the main data of the examined PTC [31][32]. It is useful to state that the maximum optical efficiency, which is found for zero solar angle, it is about 75.3% by taking into consideration all the possible optical loss factors [32]. In this work, the solar angle is zero in all the cases and the ambient temperature is selected at 25 °C. At this point, it is important to state that this work investigates a typical PTC module which has been studied by many other studies in the literature. So, the found results can be adopted as reliable also for other PTCs. More specifically, the present PTC has a typical evacuated receiver which is also used in other modules and it has a vacuum between the absorber and cover. Moreover, the collector is studied in a great range of usual operating conditions and so the found results are applicable in a great range of applications. Lastly, it has to be said that the used input properties (optical and thermal in Table 1) are usual and they correspond to real commercial systems. Therefore, the found results can be used also for other PTCs.

Basic Mathematical Formulation
The present section includes the main equations which have been used in the developed model in EES for the thermal efficiency determination of the PTC. The collector thermal efficiency (ηcol) is the ratio of useful heat production (Qu) to the solar energy on the collector aperture (Qs): The solar irradiation in the collector aperture (Qs) can be written as: The useful heat production (Qu) is calculated by the energy balance in the working fluid: The inlet temperature (Tin), the outlet temperature (Tout), the oil mass flow rate (moil) and the oil specific heat capacity (cp) are used in the previous formula.
Moreover, the useful heat production (Qu) can be expressed as the heat that the absorber transfers to the working fluid: The mean receiver temperature (Tr), the mean fluid temperature (Tfm), the internal tube area (Ari) and the convection heat transfer coefficient between fluid and absorber (h) are used in the previous formula.
The energy balance on the absorber is an important step in the present methodology and it is given below: The absorbed solar energy (Qabs) can be found below: = ⋅ The thermal losses (Qloss) of the absorber to the cover can be written as below [33]: At this point, it has to be said that the convention thermal losses between absorber and cover are neglected due to the vacuum between them. Moreover, it is essential to state that the thermal losses from the absorber to the cover are the same as the losses from the cover to the ambient because the system is studied in steady-state conditions [33]. So, it can be said that: The sky temperature (Tsky) can be estimated as below by using the temperature levels in kelvin units [34]: = 0.0552 ⋅ 1.5 (22) The heat transfer coefficient between cover and ambient (hout) is found as below [35]: More specifically, the wind speed (Vw) is selected at 1 m/s, which is a typical value, and the (hout) is found to be around 10 W/m 2 K [36].
The fluid flow heat transfer coefficient (h) is estimated by exploiting the Nusselt number (Nu): The thermal conductive (k) of the thermal oil is dependent on the temperature and it is variable in this work [37].
For this work, the flow is turbulent (Re > 10,000) and there is a tubular channel. Thus, the formula of Dittus-Boelter about the Nusselt number can be used [38]: The Reynolds number (Re) is defined as: The Prandtl number (Pr) is defined as: The mass flow rate (m) in (kg/s) is connected with the volumetric flow (V) in [L/min] as below: = ⋅ 60 ⋅ 1000 (28) In the end, it has to be said that the mean fluid temperature (Tfm) can be found as [39]:

Validation of the Developed Model in EES
The validation of the developed model in EES is performed by using experimental results by the literature [18]. Totally, 8 validation cases have been examined and they are depicted in Table 2.
The inputs for every case are given on the left side of the table (Gb, Vw, Tam, Tin, V). Finally, Table 2 shows that the mean error about the outlet temperature is 1.22% and for the thermal efficiency of 1.68%. These are relatively low values which are in the range of the experimental errors. Thus, the developed model can be assumed as accurate for the present study.

Followed Methodology for the Approximation Methods
In this work, the developed model in the Engineering Equation Solver (EES) [30] was applied in order to create a dataset for comparing the different formulas for the PTC thermal efficiency. EES is a flexible program for solving quickly and with accuracy numerous non-linear equations with numerical techniques [30]. It also includes libraries for using the thermal properties of different working fluids. Moreover, this program gives the possibility of conducting optimization studies for the maximization or the minimization of an objective function.
In this work, the solar irradiation (Gb) was examined from 300 W/m 2 up to 1000 W/m 2 with a step of 50 W/m 2 , while the inlet fluid temperature (Tin) from 25 °C up to 375 °C with a step of 50 °C. In this work, the temperature difference (ΔΤ = Τin − Τam) is important parameters which are examined from 0 °C up to 350 °C with a step of 50 °C, which is an equivalent parameter with the (Tin). The working fluid is Syltherm 800 which can operate up to 400 °C [37], something that has been taken into consideration in this work. The thermal properties of this fluid have been taken from the literature [37] and they are dependent on the fluid temperature level. Totally, 120 set points have been found and they are given in Appendix A.
About the solving procedures, for every case, the inputs in the program are the following: inlet temperature (Tin), oil flow volumetric rate (V), solar beam irradiation (Gb), ambient temperature (Tam), and the wind speed (Vw). Moreover, in this work, there is a zero-incident angle and so the optical efficiency has its maximum value according to Table 1. The program uses these inputs and calculates the outputs which are mainly the oil outlet temperature, the useful heat production, the thermal losses and the collector efficiency. The solution procedure is performed by combing the Equations (13) to (29) all together and by solving this system as a matrix with all the proper unknown parameters. Practically, it is a system with 17 equations and the following 17 unknown parameters: Generally, the solution procedure is very quick with solving time lower than 1 s per case.
The literature review of this work indicates that the thermal efficiency of the PTC is expressed with polynomial expressions and thus, this strategy is investigated in this work in a detailed way. In this work, different polynomial formulas have been used in a systematic way. All these formulas can be expressed through the following general expression: In all the examined cases, the constant term "a0" is included and all the proper combinations of the coefficients have been studied. In other words, every combination includes some parameters (e.g., a0, a1, a2, a3) and the other parameters are set to be zero (e.g., a4, b). In all the cases, there are at least two coefficients, the "a0" and another. The maximum number of the coefficients is six when all the coefficients are used in the approximation. The selection of the Formula (27) is in accordance with the existing literature. The term with "ΔΤ" is an important one, especially for the LS-2 module, and it has been used also in Reference [18]. The other terms are used in various studies. The maximum power is the fourth because the thermal radiation losses include this power and also there is no study in the literature with a higher power to be used. So, there is no reason for using the fifth power for example in the approximation formula because it would be unrealistic due to the radiation thermal losses law. So, this work studies the conventional polynomial expressions with the terms "ΔΤ m /Gb" and for all the cases, the need for the extra term "ΔΤ" is evaluated. Finally, the goal is to determine the simplest expression which can lead to accurate results.
The approximation is conducted in every case with Microsoft Excel software and the command "linest". The R 2 is calculated as an effective index which has to be extremely close to 100%. Moreover, the mean absolute percentage error (MAPE) is an extra index that is calculated and it has to be close to zero. Its definition is given below: where (n) is the total number of the dataset points which is equal to 120 in this work.

Initial Analysis
The first part of this analysis is to investigate the relationship between thermal efficiency and the usually used parameters. Generally, the correlation of the thermal efficiency and the parameters "ΔΤ m /Gb" (m = 1, 2, 3, 4) is suggested by the literature. Figure 2 shows these correlations, as well as the correlation between thermal efficiency and "ΔΤ". The results are given for different solar beam irradiation levels and more specifically for 400, 600, 800 and 1000 W/m 2 . Figure 2a shows that there is not an important correlation between the thermal efficiency and the parameter "ΔΤ/Gb" because the curves are not all together but they have different trends for different solar irradiation levels. Similar results are found for Figure 2e with the parameter "ΔΤ" in the horizontal axis. Figure 2b shows that the curves are closer to each other, especially in low values of the parameter "ΔΤ 2 /Gb", but it is not enough to describe the thermal efficiency with only the parameter "ΔΤ 2 /Gb". However, the situation is different in subfigures 2c and 2d for the parameters "ΔΤ 3 /Gb" and "ΔΤ 4 /Gb" respectively. More specifically, the curves are very close to each other when there is a plotting with the parameters "ΔΤ 3 /Gb" and "ΔΤ 4 /Gb" and especially in the case with "ΔΤ 3 /Gb", which seems to be the best one. However, a deeper analysis is needed in order for the final conclusions to be extracted. In any case, it has been found that there is a need for using terms of high power and also there is a need for using more than one parameter in order to predict the thermal efficiency in an accurate way. The next step is the detailed analysis of the correlation of the thermal efficiency with simple expressions of the format "a0 + am•ΔT m /Gb" by using 120 operating points. Figure 3 depicts the obtained results. The valuation indexes (R 2 and MAPE) are also included in these figures. The vertical axis is the value of the predicted efficiency while the horizontal axis shows the real value of the efficiency which has been calculated with EES. The results show that the use of the third power term leads to the best fitting and it is a very interesting result, which is in accordance with the quality conclusion from Figure 2. More specifically, for the formula "ΔΤ 3 /Gb", the indexes are R 2 = 99.11% and MAPE = 0.47%, while for "ΔΤ/Gb" they are R 2 = 83.74% and MAPE = 1.85%, for "ΔΤ 2 /Gb" they are R 2 = 97.74% and MAPE = 0.60% and for "ΔΤ 4 /Gb" they are R 2 = 96.64% and MAPE = 0.93%. It is clear that the first power case leads to not acceptable results, while the powers of two and four lead to results with significant errors. The next step is to present approximation cases with the term "ΔΤ" to be included as an extra term to the previous studies' formulas of Figure 3. So, Figure 4 exhibits the results about the use of the "a0 + am•ΔT m /Gb + b•ΔT" approximation about the thermal efficiency. Figure 4 indicates that the combination of the third power term with the temperature difference term leads to a very good fitting with R 2 = 99.92% and MAPE = 0.13%. This high accuracy is obtained with only three parameters and this is a very interesting and promising result. Among the other cases of the "a0 + am•ΔT m /Gb + b•ΔT" formula, the use of the fourth power leads to R 2 = 99.02%, which is an acceptable result, while the use of the power of two leads to R 2 = 98.03%, which presents lower accuracy. Generally, when R 2 is over 99% and the MAPE is lower than 0.50%, then the fitting can be evaluated as relatively accurate.
To conclude, the use of the third power term leads to high accuracy, compared to the other cases, and especially when it is combined with the "ΔΤ" it leads to a very reliable model. In any case, the use of the extra term "ΔT" helps in all the cases and the expressions "a0 + am•ΔT m /Gb + b•ΔT" are more accurate than the respective "a0 + am•ΔT m /Gb", which is a reasonable result.

Total Analysis
The next stage of the present analysis is the investigation of more detailed approximation models with more terms. It is very interesting to see the approximation statistics with different linear combinations of the terms "am•ΔT m /Gb" with and without the term "ΔΤ". Figure 5 shows the results of the different linear combinations without the "ΔΤ", while Figure 6 shows results with the term "ΔΤ". The general formula that is used in Figure 5 results can be written as "ηcol = a0 + ∑(am•ΔT m /Gb)". Figure 5 indicates that there is a need for using at least three terms in order for the points to be close to the red line and so to have a relatively accurate approximation. Also, the use of more terms leads to higher R 2 and to lower MAPE, something that leads to better fitting. However, the difference between the cases with and without the term "ΔT 4 /Gb" is extremely low (R 2 = 99.35% with this term and R 2 = 99.32% without this term). This result indicates that there is a need for using the third-order term, while the fourth-order term has a very small impact on expression accuracy. The next step is the implementation of the extra term "ΔΤ" in the previous expressions and so the general formula is "ηcol = a0 + ∑(am•ΔT m /Gb) + b•ΔΤ". Figure 6 indicates that the use of this term leads to very accurate models for all the cases except the case with the first-order term "ΔT/Gb", where the R 2 is only 87.45%. The cases with the order with m = 2, m = 3 and m = 4 lead to R 2 at 99.44%, 99.94% and 99.94%, respectively. The MAPE is found to be found 1.69%, 0.37%, 0.12% and 0.11% for the cases with m = 1, m = 2, m = 3 and m = 4, respectively. Practically, the use of the fourth power term is useless because it leads to approximately the same results with the third term results. To conclude about the results of Figures 5 and 6, the use of the third-order term is very important and also the use of the "ΔΤ" is a very useful weapon for the proper prediction of the PTC thermal efficiency.
The next step in this analysis is to present all the results about the coefficients of approximation and the evaluation indexes about the 31 examined cases. Table 3 includes all these data and this table summarizes all the found results of this work. In this table, the 31 different approximations are included. The evaluation indexes (R 2 and MAPE), as well as the values of the coefficients (a0, a1, a2, a4, b) are given. In the cases with empty cells, there is no utilization of the respective term and so the empty cells indicate that the respective value is zero. In the left column of Table 3, the number terms of the used formulas are given and the equations are classified in categories with 2 to 6 terms. Among the equations with two terms, the best formula is the following: = 0.72547 − 1.4155 • 10 −6 ⋅ 3 , (R 2 = 99.11% -MAPE = 0.47%) (32) Equation (29) is a very simple formula that leads to relatively accurate results and it is a promising choice. So, the use of the third power term seems to be a very good idea for the proper approximation of the PTC performance.
In the equations family with three terms, the best fitting has been found for the following formula: Equation (30) includes the constant term, the third-order term and the term about the temperature difference. There is an excellent fitting with Equation (30) which leads to R 2 = 99.92% and MAPE = 0.13%. These results can be compared with the indexes about Equation (29) which does not include the term "ΔΤ", where R 2 = 99.11% and MAPE = 0.47%. There is a significant difference that can be obtained only by adding an extra term. This result indicates that the use of the term "ΔΤ" is an important one in the efficiency approximation study.
About the studies with four terms, practically all the formulas lead to good fitting. Generally, the cases with third-order terms and "ΔΤ" again lead to accurate fitting. For the case with four terms, the best fitting is given below: It is very interesting to state that the use of some extra terms maybe cannot increase the accuracy effectively. Comparing Equations (30) and (32), the difference between the indexes is marginal, while Equation (30) uses 3 terms and Equation (32) uses 6 terms. So, it can be said that there is great importance to the proper term selection in order to have an accurate approximation.
Moreover, it can be said that the use of more and more terms is not the only choice for having more accurate results because the computational cost and the complexity of the approximation increase without a significant increase in the fitting suitability. A characteristic example is a comparison of the Formula (30), which includes the terms (a0, a3, b) with the formula which includes the terms (a0, a1, a2, a3, a4). According to Table 3, the first one with 3 terms has R 2 = 99.92%, while the other one with 5 terms has R 2 = 99.35%. So, it is obvious that the increase in the number of terms and the power of the terms is not the proper choice for a good fitting. So, it is suggested to use terms with third-order power and the terms with "ΔΤ".
In the literature, the use of the term "ΔΤ" is not unknown and it has been used in some studies, with the most characteristic approximation of experimental results about LS-2 PTC by Dudley et al. [18]. They applied a model with (a0, a1, a2, b) in order to have a proper fitting. The fitting of this expression, in the present example, leads to R 2 = 99.44%, which is acceptable but is lower than the respective of Equation (30) with the terms model (a0, a3, b) which is 99.92%.
Another important result of this work is the determination of the third-order as the most proper one, especially in cases with a small number of terms. This result is not found in the existing literature. Usually, second or four power terms are used. Characteristic examples are the studies of References [24][25][26], where the only not-studied term is the one with the third power. This note proves that there is a scientific gap that is covered by the present work.
Another interesting comparison of the suggested model (a0, a3, b) is with the steady-state models that are suggested by the IEC 62862-3-2. Practically, the polynomial terms of this standard are taken into consideration in order to compare them with the existing optimum model. The use of the second-order terms, of the forth-order term and their combination is compared with the third-order term. The present model presents R 2 = 99.92% and MAPE = 0.13%, while the other models of the standard have not so good fitting. More specifically, the (a0, a1, a2) model has R 2 = 99.03% and MAPE = 0.43%, the (a0, a1, a4) has R 2 = 98.96% and MAPE = 0 .48% and the (a0, a1, a2, a4) has R 2 = 99.33% and MAPE = 0.35%. So, it can be said that the suggested model is a promising one for approximating the PTC thermal efficiency. However, it has to be said that the other models are also accurate models and their use leads to an efficient approximation of the thermal efficiency data. In order to give a clear image of the aforementioned comparison, Figure 7 has been added in order to directly compare the model (a0, a3, b) with the models (a0, a1, a2), (a0, a1, a4) and (a0, a1, a2, a4). In this figure, the results from the ESS for Gb = 1000 W/m 2 are plotted with the found curves for the three models. It can be said that the model (a0, a3, b) is close to the found points while the other models present some deviations especially in low and high values of (ΔΤ). In the future, there is a need for performing similar detailed and systematic approximating studies by taking into account the wind speed (especially for bare tube cases), the solar angle (through the incident angle modifier) and the flow rate. Moreover, studies about different PTCs with a non-evacuated tube, non-selective absorber, different concentration ratios, different working fluid and different materials can also be studied.

Conclusions
The objective of this work was to investigate different mathematical expressions for the determination of thermal efficiency in a parabolic trough solar collector (PTC The use of the third-order term "ΔΤ 3 /Gb" is the most important term which can lead to relatively accurate results, especially in cases with a low number of terms. - The use of the temperature difference term "ΔΤ" is important and it can lead to high improvement compared to the respective formulas without it. - The most promising formula by taking into account the use of a low number of parameters and to achieve high accuracy is the following: It is found and discussed that there is a need for using the proper terms and not using many terms in order to achieve high accuracy. - In the future, there is a need for extending the present work by taking into consideration the wind speed, the solar angle and the flow rate. Moreover, there is a need for studying extra PTC modules with different characteristics than the present one, for example with a non-selective absorber, with a non-evacuated receiver and with different working fluid.

Appendix A-Dataset for the Present Work
Appendix A includes the dataset for the approximation studies which have been found by the use of the EES model. Table A1 includes the inlet temperatures (Tin), the solar irradiation (Gb), the temperature difference (ΔΤ = Τin − Tam) and the collector thermal efficiency (ηcol). It has to be said that these points have been found for ambient temperature (Tam) at 25 °C, volumetric flow rate (V) of 100 L/min, wind speed (Vw) at 1 m/s and zero solar angle (θ).
Moreover, the thermal efficiency with three approximation formulas is given in Table A1. The three examined formulas, which have high interest, are given below with the numbering (A1), (A2) and (A3): So, Table A1 shows the "reference" values from the EES model and the approximation with three other models in order to make clear how close the efficiency is found with these three important models.