Design, Simulation and Optimization of a Novel Transpired Tubular Solar Air Heater

: In this paper, a novel tubular solar air heater is introduced. In this air heater, the hot boundary layer is drawn into the absorber tube and can provide thermal energy at moderate temperatures. Several di ﬀ erent cases were simulated and a correlation was proposed to predict the col-lector’s e ﬀ ectiveness as a function Rayleigh number and Reynolds number. An equation was derived to ﬁ nd the e ﬀ ectiveness of this collector. Finally, a real case was studied with non-uniform solar ﬂ ux distribution, as well as radiation heat loss. Good agreement was found between the results and those derived by the proposed analytical method. For di ﬀ erent suction values, the ﬁ rst-law and the second-law e ﬃ ciencies were calculated. Based on the exergy analysis, exergy destruction in absorption is the dominant factor that is unavoidable in low-temperature collectors. It was shown that there is an optimum suction value at which the second-law e ﬃ ciency is maximized. At the optimum point, temperature rise can reach 54 K, which is hardly possible with a ﬂ at plate collector. Based on the exergy analysis, the relation between tube wall temperature and air outlet temperature in their dimensionless forms at the optimum working condition was derived, and it was shown that e ﬀ ectiveness at the optimum working condition is around 0.5. This means that the air temperature rise shall be half of the temperature di ﬀ erence between collector wall and the ambient temperatures. A high outlet temperature besides the low cost of construction and maintenance are the main advantages of this air heater. With such a high temperature rise, this type of collector can increase the use of solar energy in domestic applications.


Introduction
Solar energy is a clean and sustainable source of energy that must be included in our daily life.Besides the use of solar energy in industrial areas such as power generation, the use of this energy source in domestic and semi-industrial areas such as cooking, air heaters, and crop dryers should not be neglected.Several studies can be found concerning domestic solar collectors [1][2][3][4].For example, Nemati and Javanamrdi [1] optimized domestic solar cylindrical-parabolic cookers using exergy analysis.They proposed a simple correlation to find the optimum dimensionless temperature for various geometries and solar radiations.Moghimi et al. [5,6] used the MOGA algorithm to optimize a Linear Fresnel Collector.They selected seven geometrical parameters for optimization.They found that the most sensitive parameters were the top insulation thickness and the cavity depth.Lingayat and Chandramohan [7] used corrugated absorber surfaces to increase the heat transfer rate between air and absorber.They used the proposed geometry in an indirect solar dryer for drying banana samples.They concluded that the average collector and dryer efficiencies were 64.5% and 55.3%, respectively.For this solar dryer, banana samples were dried from 3.5566 to 0.2604 kg/kg on a dry basis.A similar study was conducted by El-Sebaey et al. [8].They used an indirect solar dryer to reduce the moisture content of banana slices.In their experimental studies, they compared chimney-type and fan-type collectors and found that chimney-type collector has a higher thermal efficiency.The mean thermal efficiency of the chimney-type was 14.5% in comparison to 12.76% for the fan-type solar dryer.
Among the domestic applications of solar energy, solar air heaters have a special place.Hot air can be used for space heating and crop drying [9][10][11].Because a solar air heater is simple and cost-effective for construction, it is one of the main interests of researchers in this field [12].There are different types of solar air heaters like unglazed transpired collectors (UTCs) [13][14][15], solar air heaters [16,17], and many other novel collectors [18,19].Regardless of new developments in air heaters to increase thermal performance, they are mostly flat plate types.Flat plate collectors are very simple and easy to construct.However, the temperature rise, and consequently the outlet exergy in these collectors, are very low.Moreover, there is no special control available to adjust the outlet temperature in this type of collector.
In contrast, the temperature rise in concentrating collectors is high.So, it can be mixed with cold fresh air to achieve a large amount of air with the proper temperature.It means that there is no need to transfer a large amount of air from the collector to the destination, but just enough to mix hot air with ambient air at the destination.This also provides control over temperature.On the other hand, there are some types of reflectors, such as CPC, that do not require high-precision technology for construction, or sun tracking systems, which do not require tracking without losing too much performance.However, convection loss is another issue.A high temperature means a strong natural convection loss.To eliminate natural convection, it is common to use a glass tube around the absorber to intercept its direct contact with the ambient air.But the use of a glass tube has many problems, including the cost of construction and maintenance, the weight, and the possibility of breakage.In the following, a novel type of collector is proposed to utilize the advantages of high temperature, which also allows low-cost construction and maintenance.This type of collector can be used for indirect crop drying, and heating or pre-heating houses and large buildings.

Problem Description
Figure 1 shows a temperature contour plot for natural convection around a horizontal tube (simulated by author).Based on this figure, a narrow thermal boundary layer starts at the bottom of the tube and it gradually gets thicker toward the top of the tube.At the top of the tube, a thick layer of air is separated from the solid wall in the form of a plume.Since in a round tube, the shape of the tube wall changes from a vertical plane at both sides (right and left sides) to a flat horizontal plane at the top of the tube, the boundary layer is thick and wide at the top of the tube.So, it is a great source of high-temperature thermal energy or exergy.The solid black line in this figure shows the isothermal line of 380 K (80% of the temperature difference between wall and ambient).Based on the rulers in the figure, for the current case, the isothermal line is around 12 mm in width (from −6 mm to +6 mm).This rising plume carries a considerable amount of high-temperature thermal energy that is lost.So, if a narrow slit is installed at the top of the tube, with proper suction the hot air accumulated on the top of the tube can be drawn in. Figure 2 clarifies this idea.In Figure 2b, the thermal plume vanishes and hot air is sucked in.Since a tubular absorber is a proper shape for collectors with a concentration ratio of more than one, this Unglazed Transpired Tubular Absorber (UTTA) can be used in concentrating collectors like trough collectors, linear Fresnel collectors, and CPCs.High-temperature air can later be mixed with fresh cold air to achieve an adequate amount of warm air at a proper temperature.Figure 3 schematically shows the schematic of a UTTA.

D
In a common concentrating collector, there is a glass tube around the absorber tube to reduce the natural convection heat loss.But, due to the different thermal expansion coefficients for the absorbent tube and the glass, and also due to the significant difference in temperature between these two tubes, breaking the glass cover is not unlikely to occur.In a UTTA, since the natural convection heat loss is reduced and controlled by proper suction, there is no reason to keep the glass tube cover around the absorber.Eliminating the glass tube significantly reduces construction and maintenance costs.On the other hand, since the temperature in the circumferential direction of the absorber tube varies considerably, thermal stress is one of the main concerns of designers [20][21][22].However, the narrow slit at the top of the absorber tube circumferentially provides a free displacement for the tube, and thermal stress is released.In the following, the method of simulation of a UTTA and effective non-dimensional parameters are presented.Nearly fifty different cases related to a wide range of non-dimensional parameters were simulated, and a correlation was proposed to predict the effectiveness of UTTA.Furthermore, the effects of radiation loss and non-uniformity in solar heat flux were studied, and it was shown that the proposed analytical method effectively predicts UTTA effectiveness.

Governing Equations
For a narrow slit where the pressure drop of the flow passing through the gap is much higher than the pressure drop of the air flow along the pipe, the suction is uniform along the slit and the problem can be simplified to a two-dimensional case.For steady and laminar natural convection over the outer surface of a tube, the governing equations, as well as the boundary conditions, are presented in the following.
Continuity equation: Momentum equations: Energy equation in the air zone: The energy equation in the solid zone: Air properties including density were assumed as a function of temperature.The dependency of properties on temperature is shown in Appendix A.
The diameter of the air domain around the wall was considered to be ten times the outer diameter of the tube.It was large enough to avoid numerical errors due to the false impression of boundary conditions [23][24][25].The boundary conditions at the outer surface of the tube were  =  = 0,  =  (6) and at the far field they were: At the solid-fluid interface, there was no heat flux and temperature jump.Moreover, the no-slip condition was assumed.So: where  is the normal direction.Because of the geometry and flow symmetry, only onehalf of the domain was considered in the simulation.Finally, the air mass flow rate per unit length ( ) was specified over the hole on the tube.

Solution Methodology
A very fine quadrilateral structured grid was used to discretize the domain.To capture the thermal boundary layer, the distance of the first gird to the tube wall was set to 3.95 × 10 −5 m.The total number of grids was 248,000 (See Figure 4).Equations were solved numerically, using Ansys-Fluent 2020.A couple algorithm was used for pressure-velocity coupling and the second-order upwind method was used for momentum and energy discretization.Residuals were set to 10 −6 and 10 −9 for continuity and energy equations, respectively, and the Pseudo Transient method was used to solve equations.For the grid study, natural convection over a horizontal tube with a 0.06 m diameter was simulated by the above settings.The results are presented in Figure 5. Based on this figure, the difference between  for the highest and lowest grid numbers is less than 0.17.Therefore, the selected number of the grid is satisfactory.Near fifty different cases were simulated using the aforementioned settings, and their results are verified and studied in the next parts.The ranges of contributing parameters in the simulation are presented in Table 1.Before starting the study, the method was validated by simulating the natural convection around a horizontal tube.The results were compared for a wide range of  , against two experimental correlations, i.e., Churchill and Chu [26] (Equation ( 9)), and Morgan [27] (Equation ( 10)), with a maximum uncertainty of ±5%, as reported by [27].The comparison in Figure 6 shows good agreement between the numerical and experimental data and shows the validity of the CFD approach.The average relative error with Equation ( 9) is 0.9% and with Equation ( 10) is 3.4%. where in which the convective heat transfer coefficient, ℎ , is: and All contributing parameters were calculated at the film temperature:

Data Analysis
To analyze the results, it is necessary to introduce the governing dimensionless parameters.The main parameter, i.e., effectiveness is [13]: which shows the value of thermal energy gained at the slit to the maximum theoretically accessible thermal energy. is the outlet air temperature at the slit (see Figure 3).The two other contributing parameters are: where  is the slit width and  is the air velocity at the slit entrance. =   is the air mass flow rate per unit length.Note: the thermo-physical condition required for Equation ( 16) is film temperature and for Equation ( 17) is the temperature of air at the entrance to the slit,  .
For a series of solutions, the variation in  with respect to  is shown in Figure 7.As is expected, the effectiveness decreases with increasing suction value.Generally, three zones can be assumed in this curve, i.e., low suction, high suction, and moderate suction, and three data points are selected for each zone, to be discussed in detail.Figure 8 shows the temperature contour plot of three selected points (column (a)) as well as the velocity vector plot around the slit (column (b)).For point P1 associated with a low suction, a narrow layer in the vicinity of the wall is sucked in.So, the suction layer is located entirely beneath the boundary layer thickness.
It is known that in natural convection, the thermal boundary layer thickness is proportional to  (over a vertical flat plate,  = 1 4 ⁄ [28]).Since for near-to-zero suction, only a thin wall-adjacent layer of gas is sucked in, the effectiveness has a weak dependency on thermal boundary layer thickness, and consequently on .By reducing the suction value to zero, almost all of the thermal energy is wasted by natural convection, and nothing is recovered.For moderate suction (point P2), part of the thermal boundary layer is drawn into the slit.So, some part of thermal energy is wasted and the rest is recovered.The momentum of the flow in natural convection and forced convection is proportional to  and , respectively.So, the amount of flow that is trapped and sucked in involves a competition between the momentum of natural convection flow and forced convection flow.Consequently, the effectiveness is expected to depend on the ratio of  and  .From one point onwards, the entire boundary layer as well as fresh air is drawn in (point P3).So, almost nothing is wasted.But because of the fresh air, the temperature of the outlet from the slit and also effectiveness decrease.Under this condition, in the limiting case, i.e., high , .ℎ ( −  ) =  ( −  ) Using Equations ( 11), ( 15) and ( 17), the non-dimensional form of the above equation (Equation ( 18)) is: =  −  is the heat transfer area per unit length of the tube.Finally, using Equation (10), the final form of effectiveness at high  will be Based on the above explanation and through trial and error, the following correlation is proposed for  ≤ 800 and 1.8 × 10 ≤  ≤ 7.8 × 10 .In the above equation, for  = 0 , the effectiveness approaches   ⁄ , which is slightly less than unity.At a very low suction, the velocity of the rising plume in the middle of the slit is dominant over the suction velocity, and consequently, some amount of the sucked air returns back to the ambient.
For the two limiting cases, i.e., very low and very high , Equation ( 21) reduces to: For very low , the form of the equation shows a very weak dependency on , and for very high , the equation approaches the analytical result (Equation ( 20)).The value of R 2 for Equation ( 21) is 0.998 and the RMSE is 1.1%, which shows the accuracy of the proposed correlation.The comparison of Equation ( 21) with numerical results is presented in the form of a scatter plot in Figure 9. Figure 10 shows the variation in effectiveness for various  and .Based on this figure, effectiveness decreases monotonically with increasing suction ().Moreover,  is higher at higher .

Energy Analysis of a UTTA
There are two main sources of energy loss, i.e., natural convection and radiation.So, in the absence of other sources, the total absorbed inlet energy (. ) shall be equal to energy losses.The energy balance for a UTTA is: This means that the absorbed solar heat input (. ) is wasted by radiation and convection [29].
In the above equation,  is the Stefan-Boltzmann constant and  is tube emissivity, which is equal to its absorptivity as per Kirchhoff's law of thermal radiation [29].
The first-law efficiency can be defined as: This can be presented as a function of dimensionless parameters.Starting from Equation (17), the air mass flow rate per unit length is: Alternatively, using Equations ( 11) and ( 15 Replacing the above equations as well as Equation (24) in Equation ( 26) yields: Knowing that  = .  ⁄ and using Equation (28), the above equation is simpli- where

Exergy Analysis of a UTTA
The exergy balance for a UTTA can be written as: Equation ( 33) shows that the incoming exergy from the sun is partly recovered via the suction, a part is reflected and not absorbed, a part is lost because of the heat loss, and finally, the rest is degraded due to heat transfer with a finite temperature difference.
The inlet exergy from the sun is obtained by multiplying the inlet energy from the sun (Equation ( 24)) by the maximum efficiency [30]: In the above equation,  ~4500 K is the effective sun temperature [31].
The exergy gain at the outlet of the collector (on the slit) is [13]: The exergy loss due to the reflection from the absorber is: The exergy loss from the absorber due to the radiation and convection is: The first term is the exergy loss due to radiation heat loss, while the second term is exergy loss due to convection heat loss.
The exergy destruction is [13]: where the first term is destruction due to convection heat transfer with a finite temperature difference and the second term is exergy destruction due to absorption from sun temperature to wall temperature.The second-law efficiency is defined as the ratio of gained exergy to inlet exergy from the sun [1].So: Replacing Equations ( 11), ( 27) and ( 32), the second-law efficiency will be: Using Equation ( 29), the above equation is simplified to: Alternatively, using Equation ( 31), where  is the ratio of temperature to ambient temperature:  =   ⁄ .

Analysis of a Real Case
In a real case, instead of temperature, the heat flux over the absorber tube is known, and is not uniform.Furthermore, besides natural convection, radiation heat loss shall be considered.In a real case, a parabolic trough collector with a focal length  = 0.1 m and a width of 0.4 m is considered.The absorber tube in this collector is made of aluminum with an outer diameter of  = 0.06 m and a thickness of 4 mm.The reflectivity of the parabolic reflector is  = 0.95 and the absorptivity of the absorber tube is  = 0.9.All materials are gray.Geometrical parameters are summarized in Table 2.For a parabolic collector with an aperture plane normal to incident sun rays, the flux distribution over the tube is [32]: In the above equation,  = 1000 W/m 2 is the solar intensity.The other parameters used are illustrated in Figure 11.For one meter of the tube, the absorbed solar radiation (per Equation ( 44)) is . = 333.34W/m.For the above collector, the absorber tube was modeled based on the previous procedure.To study the collector analytically, radiation and natural convection heat loss shall be calculated first.Radiation and natural convection heat loss ( ,  ) are calculated using Equation (45).
To calculate the tube wall temperature, the procedure shown in Figure 12 is used.Finally, one should guess  and then calclulate  at  .Effectiveness is then calculated using Equation ( 21), and  is calculated using Equation (15).This loop shall be repeated until the difference between the calculated and the guessed values approaches zero.The results are presented in Table 3.Furthermore, the intermediate results are also presented in Table 4 for more clarity. . = 0.003 (1/K) ℎ = 7.98 (W/m 2 K)  = 1.17Based on the results in Table 3, the maximum difference in the value of  between analytical and numerical methods is 2.2 K, which shows the accuracy of the analytical method.Figure 13 shows the temperature rise relative to ambient temperature, as well as energy gain.Based on this figure, for the current selected mass flow rates, the temperature rise varies between 85.8 K and 19.4 K, while energy gain varies between 8.8 W/m and 88.6 W/m.At a suction flow rate of around 0.002 kg/s.m, the useful energy reaches its saturated value.At this stage, the thermal boundary layer is completely sucked in, and after that, fresh air is drawn in in addition to it.Therefore, for a higher flow rate, the energy gain does not change significantly, and the rise in temperature is also reduced.For the above case, the variations in effectiveness, as well as first-and second-law efficiencies, are plotted in Figure 14.Based on this figure, we can see that the second-law efficiency is much smaller (one order of magnitude) than the first-law efficiency.This conclusion was observed before in other research [1].Furthermore, there is an optimum point at which the value of second-law efficiency is at its maximum.The optimum values are presented in Table 5.Based on Table 5, at the optimum working point, the outlet temperature is 354.2K.This shows that the ambient temperature has risen more than 54 K.This considerable temperature rise shows the superiority of this air heater in comparison to the common flattype air heaters.It must be emphasized here that the elimination of the glass cover over the tube absorber reduces the maintenance and the capital cost in comparison to concentrating collectors.
To clarify the existence of an optimal working point as well as the contribution of each exergy loss or destruction mechanisms, the dimensionless forms of exergy losses and destruction are presented: In this case,  , = 64.7% is the main source of exergy destruction.However, this is inevitable in low-temperature solar applications, because there is a large gap between the sun's temperature and that of the solar absorber.This form of exergy destruction is the main reason for low second-law efficiency.Apart from the  , , other forms is the part of the incoming exergy that is not absorbed and reflected.Since it depends only on absorptivity, it is constant and does not vary with .For a constant heat flux, absorber wall temperature, and consequently radiation loss, are constant.So,  , is also constant and independent of .However, the behaviors of  , and  , are basically different from those of the others.
Based on Figure 15, by increasing the suction value and consequently  (See Equation ( 17)), the part of energy that is recovered increases, and as a result the convection loss as well as  , decreases.However, the temperature of air outlet from the collector decreases, resulting in the increase in  , or  .A summation of these two factors is also shown in Figure 15.As expected, there is a minimum value for the summation of these two factors.In the next section, the optimum value will be presented.The last parameter that should be studied is the solar heat flux.For the above collector, at the optimum working point, i.e.,  = 120.4or  = 0.0013 kg/s.m, the heat flux was varied between 1000 W/m 2 and 2000 W/m 2 .Figure 16 shows the variation in temperature rise, with respect to heat flux.Based on this figure, as the received solar heat flux increases, the temperature rise at the slit increases almost linearly.Consequently, the thermal effectiveness of the collector increases also.Furthermore, as is expected, the secondlaw efficiency increases too.This is because, at a constant air mass flow rate, by increasing the air temperature at the slit, the exergy gain increases according to Equation (34).The variations in the effectiveness as well as the second-law efficiency with respect to heat flux are shown in Figure 17.The variation in the second-law efficiency can also be observed in Figure 17.Unlike the second-law efficiency, the first-law efficiency is reduced with the increase in the solar heat flux.The main reason is that at higher solar heat flux and thus higher tube wall temperature, the radiation heat loss is higher, and as a result, the efficiency decreases.

UTTA Optimization
For a known geometry, the main challenge is determining the optimum working condition, i.e., optimum effectiveness.According to Figure 15, it is known that there is a suction value at which the second-law efficiency is maximized.Unfortunately, there is no analytical solution for the optimum working condition.Therefore, the optimization was performed numerically for 100 different cases for  = 0.02 to 0.1 m and  = 50 to 2000 W/m.A subroutine was written in MATLAB to calculate the second law efficiency for a wide range of Re and the specified different working conditions, and for each Re range, the maximum value was derived using the adaptive Newton-Raphson method.The optimum working condition can be approximated by the following curve-fitted equation with R 2 = 99.94% and RMSE = 0.0052.
The interesting result of Equation ( 51) is that the effectiveness at the optimum working condition is around 0.5, with a weak dependency on  .
Equation ( 51) is an applicable equation that can be helpful in the design of a UTTA.For example, for a tube with a diameter of 0.03 m and a slit width of 0.008 m, it is known that the absorbed heat is 400 W/m and the ambient temperature is 300 K.Under this condition, by following the procedure shown in Figure 12,  = 1.64.It just remains to find the best working condition.Based on Equation (15), to reach the optimum working condition,  shall be 1.32.This means that the outlet air temperature shall be 396 K.As a result, the collector's effectiveness is 0.5.Using Equations ( 21) and ( 27), the fluid capacity of the collector is  = 0.00193.

Summary
In this paper, a novel tubular solar air heater was introduced.The thickness of the thermal boundary layer in natural convection is thick atop a horizontal tube.So, this offers a great source of energy that can be sucked into the tube and used.Furthermore, in this type of collector: - The glass tube is eliminated, which results in a considerable reduction in construction and maintenance costs; -Hot air can be mixed with fresh air to derive an adequate amount of warm air at the desired temperature.So, the air outlet temperature is under the control; - The mixing can be done at the destination, and therefore, smaller ducts are required to transfer hot air; - The temperature rise in this collector is relatively high.A higher temperature for thermal energy means a higher level of availability and exergy.
Several different cases were simulated, and based on the results, the suction of the boundary layer can be categorized into three zones: -Low suction-Only a narrow layer of gas adjacent to the wall is drawn in and the effectiveness has a weak dependency on  .The effectiveness approaches very slowly to a value less than unity when the suction value reduces toward zero.Consequently, most of the thermal energy is wasted and nothing is recovered; -Moderate suction-A part of the thermal boundary layer is sucked.So, only a part of thermal energy is recovered and the rest is wasted.The effectiveness depends on the ratio of  and ; -High suction-Fresh air in addition to the boundary layer is sucked.So, almost nothing is wasted.From this point on, the amount of total energy is fixed, and by increasing the air flow rate, temperature rise proportionally decreases.
A correlation was proposed to calculate effectiveness as a function of the Reynolds number and Rayleigh number.In the end, the first-law and second-law efficiencies for the case under consideration were calculated for different suction values.At the optimum point where the second-law efficiency is maximized, temperature rise can reach 54 K, which is much higher than the working zone of a flat plate collector.Based on the exergy analysis, exergy destruction in absorption is the dominant factor.This mode of destruction is unavoidable in low-temperature collectors, and is the main source of low efficiency.
Finally, a correlation was proposed to find the optimum working condition.It was shown that in the optimum working condition, the effectiveness is around 0.5.
This new type of solar collector can produce air at higher temperatures than flat plate collectors.Therefore, it can extend the use of solar energy in daily life.Since the performance of different solar collectors as well as the proposed solar collector in this research is affected by wind, it is suggested to investigate the effects of wind on the performance of such collectors in future research.

Conflicts of Interest:
The authors declare no conflict of interest.

Figure 1 .
Figure 1.Temperature contour plot for air natural convection around a horizontal tube.

Figure 2 .
Figure 2. Comparison of temperature contour plot for air natural convection around a horizontal tube, (a) without suction, (b) with suction (a narrow slit is installed at the top of the tube and the hot air is sucked in).

Figure 4 .
Figure 4.The structured mesh used and a zoomed-in view (red square) around the tube.

Figure 5 .
Figure 5. Variation of Nu versus elements number.

Figure 6 .
Figure 6.Comparison of Nu of CFD and experimental results (Error bars show ±5% uncertainty).

Figure 7 .
Figure 7. Variation of  with respect to  for  = 1.32 × 10 and its zoomed-in view for small  (dashed square).Three sample points P1, P2, and P3 show small, moderate, and high suction condition.

Figure 8 .
Figure 8. Contour plot of selected points: (a) temperature contour plot, (b) zoomed-in view of temperature contour plot around slit and velocity vectors.

Figure 9 .
Figure 9.The comparison between numerical results and the present correlation.Data points are presented with hollow circles.

Figure 11 .
Figure 11.Parameters used to calculate heat flux.

Figure 13 .
Figure 13.Temperature rise and useful energy gain.

Figure 14 .
Figure 14.Variations of effectiveness, as well as first-and second-law efficiencies, with respect to Re.

Figure 15 .
Figure 15.Contributions of different forms of exergy destruction and loss.

Figure 16 .
Figure 16.Variations in temperature rise with respect to heat flux.

Figure 17 .
Figure 17.Variations of effectiveness, as well as first-and second-law efficiencies, with respect to heat flux.

Funding:
This research received no external funding.Data Availability Statement: Data available on request from the author.

Table 1 .
Ranges of simulating parameters.

Table 2 .
Geometrical parameters of the trough collector.
* Area of one meter of absorber tube excluding slit.

Table 3 .
Calculation of  using analytical and numerical methods.

Table 4 .
Intermediate results in the calculation of  using the analytical method.
Absorbed solar heat flux over the tube (W/m 2 ) Q Heat transfer per unit length of the tube (W/m)