Theoretical and Experimental Study of a Double Air-Pass Solar Thermal Collector with an Insulating Rod of Millet

In this article, we present the study of a double-pass air insulated by crushed millet stem mixed with gum arabic. The study is carried out based on mathematical models obtained by writing energy conservation laws in the various components of the system, which made it possible to determine the evolution of the air temperature as a function of the length of the absorber and to make a comparison with the experimental results. After comparing the results obtained with those found in the literature, the influence of some physical and geometrical parameters on the performance of the solar thermal collector is presented.


Introduction
Plan solar air collectors convert solar energy into thermal energy extracted from the air into flow in the solar collector [1].This energy is used in a variety of applications, such as drying foodstuffs (fruit, wood, etc.) [2]- [18], heating, and so on.
They have a number of practical advantages, such as the direct use of air for drying products in industry, which translates into remarkable energy savings, and this is why recent developments in the field of high environmental quality have played an important role in the design offices and particularly in the industrial sector.
To produce solar thermal energy, several types of flat solar collectors are designed for different performances [3]- [21].Various studies on their performance lector according to their components, namely: glass, heat transfer fluid, absorber and insulation system.The solar radiation captured has an impact on the efficiency of the solar collector [4], the increase in the contact surface of the absorber is proportional to the rate of solar radiation received by the absorber, but also increases the heat transfer coefficient between the cover and the environment.
The solar collector glazing shall have characteristics enabling it to optimize the energy transmitted to the absorber by the greenhouse effect [5].The triple cover reduces heat loss forward more than the double cover [6].According to Wijeysundra et al., the efficiency of the double-glazed solar collector is increased by 10% -15% compared to single-glazed glass [7]- [15].
The orientation of a collector has an influence on the solar radiation received by its surface [8]- [23].The best position of the solar collector is to point it directly south.
The location of obstacles in the various channels of the solar collector influences its efficiency [9]- [17], which can reach 75% with a mass flow rate of 0.072 kg/s and a solar radiation of 788 W/m 2 [10]- [22].The finned absorber is more efficient than the finless absorber.The thermal efficiency of a double-pass solar collector with a porous medium is 60% to 70%, 20% to 30% higher than the collector without a porous medium [11]- [20].Its efficiency is high compared to conventional solar collectors, it can exceed 75%.
The porosity of the fluid influences the temperature increase according to the mass flow rate [12].Elradi A. Musa and Al's work has shown that Reynolds and Nusselt numbers also have a significant effect on low porosities [13]- [19].
To do this, a mathematical model based on the solar collector thermal balance method is studied.The results of the numerical simulations were used to evaluate the temperature profiles of each solar collector layer (glass, absorber, fluid and insulating plate).In order to validate the theoretical results, numerical results obtained with Comsol code are used to compare the evolution of the thermo-physical parameters of the fluid in the solar collector.The comparative analysis of the theoretical and experimental results of the double-pass solar thermal collector has enabled the model to be validated.

Experimental Device for Solar Collector
The solar thermal collector studied is shown schematically in Figure 1.It consists of a single glass pane with a surface area of 4 m 2 and a thickness of 0.005 m [14] [15] [16].The greenhouse effect produced at its level causes the absorber to heat up, which is made of black painted steel sheet, 0.003 m thick.The insulation layer, made of crushed millet shank, 0.030 m thick behind the sensor and on the side, minimizes losses and therefore increases the thermal performance of the solar collector.A channel (e f1 = 0.045 m on the outward flow and e f2 = 0.070 m on the return flow) connecting the inlet to the outl et al. lows air to be guided by forced convection through a fan placed at the outlet from which a drawn flow is drawn as shown in Figure 2. A total of 36 Thermocouples are placed along the Energy and Power Engineering length of the sensor to measure the temperature evolution of the different solar collector components, as illustrated in Figure 3.The air at the solar collector inlet is at room temperature with fixed air flow rates measured with the anemometer.Daytime system performance is subject to outdoor conditions (solar radiation temperature).The measurements are made with a data acquisition (Agilent 34970A) equipped with two multiplexers, comprising 19 thermocouples for the first 34901A and 17 connections for the second 34908A.The Agilent data acquisition is placed underneath the solar collector with the computer.A potentiometer is used to adjust the fan speed which is proportional to the flow rate.

Heat Balance Equations
The model is based on the following assumptions:  The external and internal convective heat transfer coefficients are constant over the length of the solar collector. The thermal conduction is neglected  The pressure losses are neglected in the side walls and at the bottom of the solar collector. All surfaces of the different components are equal Under the above-mentioned hypothesis, the heat balance equation of each component of the solar collector are.

Glass
Fluid (f2) As the heat transfer fluid underneath the absorber exchanges with the plate (p) and the absorber (ab) by convection, then we have The modeling of the solar collector is based on a nodal discretization showing 5 knots.

Heat Transfer Coefficients
The system of equations presented in section 4 is based on a good knowledge of the heat transfer coefficients to take into account the heat exchanges, by conduction, convection and radiation, between the different components of the solar collector.
Heat transfer between celestial vault and glass ( ) Others use: ( ) In our calculations, we will adapt this last expression for the calculation of T c .
The heat transfer between the ambient environment and the glass is given by Mac Adams' formula: 5.7 3.8 Radiation heat transfer: Heat transfer between fluid 1 and glass: where: Heat transfer between fluid 1 and absorber: The heat transfer coefficient between glass and fluid (f1) is assumed to be equal to heat transfer coefficient between the absorber and fluid (f1).Both plates are covered by the same fluid with the same fluid velocity and the glass and absorber have the same length.
The convective and radiative heat transfer coefficients of fluid f1 were used for return fluid f2.

Thermo-Physical Properties
The properties of the individual solar collector components are given in Table 1.
The insulation used consists of crushed millet stem mixed with gum arabic binder.
It is a very accessible material in rural areas and in large quantities, often unusable.

Two-Dimensional Solar Collector Modeling
The system is modeled by Comsol multiphysics 3.5 code taking into account the couplings of conductive, convective and radiative heat transfers between the fluid flow and the solid components of the solar collector.Figure 4 shows the temperature evolution and flow direction of the heat transfer medium.The influence of the grid number on the temperature is presented in Table 2.
In the following conditions: Solar heat flux G max = 800 W/m 2 , α = 0.9, air ρ = 1.2 kg/m 3 , C p = 1.006J/(kg•K) and μ = 10 −5 Pa•s and physical properties presented in Table 1.Thus for Comsol the elements were defined with a number of degrees of freedom of 123,928, a number of mesh points equal to 9613, a number of meshelements corresponding to 18880 essentially composed of triangular meshes, with 848 delimiting elements for 12 vertices.The minimum quality of the elements is 0.652 for an element area ratio of 0.275.
The variation in mesh size has a slight influence on temperatures for numbers of elements ranging from 295 to 4720, although not optimal.It becomes regular from a number of elements equal to 18,880, which was chosen for its comparison with the thermal performance of the solar collector under working conditions (Table 1).

Numerical Results
The system of Equations ( 2) to ( 6) is discretized by finite differences method and then solved using Gauss-Seidel's iterative method.
 Air temperature evolution in the solar collector For a fixed flow rate of 0.023 kg/s, inlet air temperature of 298 K and solar heat flux ranging from 600 to 900 W/m 2 , Figure 5 and Figure 6 shows respectively the air temperature evolution along the upper channel T f1 and lower channel T f2 at 9 h, 11 h and 14 h.The maximum of air temperature flowing through the solar air collector, went respectively to 308 K, 325 K and 329 K (59 C maximum temperature).There is a temperature difference of 31C between the inlet temperature and the outlet temperature at 2 pm due to the temperature gradient of interest in the flow.
The air enters the second channel with the outlet temperature of T f1 .Thus the lower air temperature through the solar collector at the return flow went respectively up to 320 K (47 C), 338 K (65 C) and 341 K (68 C) output temperature.A temperature difference of 17 C is less important than in the first channel and remains relatively constant when the thermal equilibrium is reached (from 70 cm onwards). Influence of mass flow rate  For a solar heat flux of 900 W/m 2 at 14 h, Figure 7 and Figure 8 show the effect of mass flow on the air outlet temperature in the solar collector.It can be seen that for 0.023 kg/s, 0.035 kg/s and 0.051 kg/s, the output temperatures of the heat transfer medium are 341 K (68 C), 338 K (65 C) and 334 K (61 C) respectively.Thus, it can be seen that the more the flow rate tends towards its minimum value, the higher is the air outlet temperature.This is due to the length of time, more or less long, that the fluid stays in the solar collector, which causes more exchange between walls and air. Comparison of numerical results For a flow rate of 0.023 kg/s and solar heat flux of 900 W/m 2 , the analysis of the obtained results shows a small difference of 7C between the two models (Figure 9 and Figure 10).This difference can be explained by the very nature of the two models studied, the first one (global model) being more global based mainly on a good knowledge of the heat transfer coefficients, the second

Comparison with Experimental Results
Figure 11 and Figure 12 show the evolution of the air temperature throughout the solar collector with a flow rate of 0.023 kg/s using theoretical and experimental models.For the same rate fixed at 0.023 kg/s in the experiment as in the numerical resolution with irradiation of 600 W/m 2 and 900 W/m 2 at 9 h and 14 h respectively, we note that the pace of the experimental curves is relatively equal to that of the theoretical curves.At 9 o'clock, the experimental and numerical curves have the same output temperature 320 K (47 C) and 14 o'clock respectively 335 K (62 C) and 341 K (63 C) with a deviation of 6 C.This difference can be explained by the thermal pressure losses of the solar collector in the experiment.

Solarcollector Output
The collector efficiency is defined as the ratio of the effective power u φ extracted Energy and Power Engineering from the collector to the incident solar flux G.The effective power is evaluated using enthalpy balance.
T i is the inlet temperature of the fluid f1 and T o is the outlet temperature of fluid f2.
Figure 13 shows the performance evolution of the solar collector for a fixed flow rate of 0.023 kg/s and solar radiation of 900 W/m 2 .
Figure 13.Performance of solar collector.
We note, that at this flow rate and maximum irradiation, that the optimal efficiency is 80% for a double-pass single-glazed air solar collector with crushed millet rod as insulator.

Conclusions
In this work, we have proposed a numerical and experimental study of a double air-pass solar collector with the objective of producing hot air to supply and improve drying techniques.
For this, we realized a double pass air solar collector and developed a global A good agreement between calculated and measured temperatures is observed in the pm of the day while in the am.The results show a significant difference due to the lateral thermal losses and also of the thermal inertia of the insulation material of the solar collector that are not taken into account in our global model.
have been analyzed and have shown results on the parameters of the solar col-M.M. D. Ndiaye et al.DOI: 10.4236/epe.2018.103008107 Energy and Power Engineering

Figure 1 .
Figure 1.Picture of the double air pass solar thermal collector.

Figure 2 .
Figure 2. Diagram of the studied solar collector.

Figure 3 .
Figure 3. Position of thermocouples in the sensor collector.

Figure 4 .
Figure 4. Temperature field over the solar collector.

Figure 5 .Figure 6 .
Figure 5. Evolution of the upper air temperature Tf1 along the solar collector

Figure 7 .
Figure 7. Distribution of the upper air temperature T f1 for different mass flow rates at 14 h.

Figure 8 .
Figure 8. Distribution of the lower air temperature T f2 for different mass flow rates at 14 h.

Figure 9 and
Figure 9 and Figure 10 show a comparison of temperature changes obtained with the two approaches (global model and Comsol).For a flow rate of 0.023 kg/s and solar heat flux of 900 W/m 2 , the analysis of the obtained results shows a small difference of 7C between the two models (Figure9and Figure10).This difference can be explained by the very nature of the two models studied, the first one (global model) being more global based mainly on a good knowledge of the heat transfer coefficients, the second

Figure 9 .
Figure 9. Distribution of temperatures T f1 of the upper canal.

Figure 10 .
Figure 10.Distribution of temperatures T f2 of the lower canal.
model and a numerical code to simulate and follow its thermal behavior.The solar collector is insulated thermally by a local material composed of easily accessible crushed millet stem.The influences of masse flow rate on the air transient temperature response are presented along the solar collector at different instants and solar heat fluxes.

Table 1 .
Thermo-physical properties of materials.

Table 2 .
Influence of mesh size.