Aoues 2011
Aoues 2011
To cite this article: Kamel Aoues , Noureddine Moummi , Miloud Zellouf & Adel Benchabane (2011) Thermal performance
improvement of solar air flat plate collector: a theoretical analysis and an experimental study in Biskra, Algeria,
International Journal of Ambient Energy, 32:2, 95-102, DOI: 10.1080/01430750.2011.584469
Taylor & Francis makes every effort to ensure the accuracy of all the information (the “Content”) contained
in the publications on our platform. However, Taylor & Francis, our agents, and our licensors make no
representations or warranties whatsoever as to the accuracy, completeness, or suitability for any purpose of the
Content. Any opinions and views expressed in this publication are the opinions and views of the authors, and
are not the views of or endorsed by Taylor & Francis. The accuracy of the Content should not be relied upon and
should be independently verified with primary sources of information. Taylor and Francis shall not be liable for
any losses, actions, claims, proceedings, demands, costs, expenses, damages, and other liabilities whatsoever
or howsoever caused arising directly or indirectly in connection with, in relation to or arising out of the use of
the Content.
This article may be used for research, teaching, and private study purposes. Any substantial or systematic
reproduction, redistribution, reselling, loan, sub-licensing, systematic supply, or distribution in any
form to anyone is expressly forbidden. Terms & Conditions of access and use can be found at http://
www.tandfonline.com/page/terms-and-conditions
International Journal of Ambient Energy
Vol. 32, No. 2, June 2011, 95–102
This article presents the results of the first experimental investigation of the performance of solar air flat
plate collector at Biskra (latitude 34 480 N and longitude 5 440 E), Algeria. The thermal efficiency between
absorber plate and air in flat plate solar collector has been enhanced by introducing obstacle rows in the dynamic
air vein of the collector. For this objective, a flat plate solar collector, of 1.73 m2 area and 25 mm air gap, has
been designed and constructed. These obstacles formed with two parts: first part is perpendicular to fluid flow
and the second part is inclined, they are mounted in a staggered pattern, oriented perpendicular to the fluid flow
Downloaded by [Anadolu University] at 10:41 26 December 2014
and soldered to the back plate. The solar air heater was mounted on a stand facing south at inclination angle, and
it was tested under the environmental conditions. Moreover, a theoretical approach is employed for
determination of the thermal performances of this collector where the temperatures of all components of the
collector and outlet air are predicted. Comparisons among the experimental and theoretical results considered are
reported.
Keywords: solar energy; flat plate collector; obstacles; heat transfer; oriented flow
duced rows of fin obstacles in dynamic air vein, which obstacles. The latter is shown in Figure 1. The
are mounted in staggered pattern, oriented perpendic- obstacles are ach ¼ 70 mm in length and are
ular to fluid flow and soldered to the back plate. The mounted perpendicular to the flow.
influence of obstacles on the thermal performance has
been investigated. The experimental setup is schematised by a solar
In the theoretical study, a mathematical method, collector, as shown in Figure 1.
using the approach of Hottel, Whiller and Bliss, is
developed to determine the thermal performances of
flat plate collector. It determines the temperature
profile of the absorber and the outlet air. These results 3. Theoretical analysis
are compared with those obtained by the experiments. The air solar flat plate collector used is with a simple
pass between the absorber and the back in a galvanised
metal plate placed on the insulator; the obstacles
increase the thermal performances of the collector,
2. Solar collector and experimental setup consequently increasing its temperature at the exit.
In this study, a flat plate solar collector is constructed The method selected for modelling and studying
with a length of Lc ¼ 1:95 m and a width of lc ¼ 0:89 m the performances of this collector is the total method
(Figure 1). These collectors are made up of a number which supposes that all the components of this section
are at a constant average temperature, the average
temperature between the inlet and the outlet of this
collector.
The collectors operate under quasi steady-state
conditions (Figure 2). In these conditions, the perfor-
mance of a solar collector is described by an energy
balance that indicates the distribution of incident solar
energy into useful energy gain, energy stored and An empirical equation for the loss coefficient
energy losses. through the top of the collector Ut was developed by
Klein following the basic procedure of Hottel and
g ¼ u þ l þ st ð1Þ Woertz (1942) (Duffie and Beckman 1991):
By neglecting the thermal energy stored (thermal 2 31
inertia) in the collector, we obtain 6 N 17
Ut ¼ 4 h ie þ 5
g ¼ u þ l ð2Þ C ðTabs Ta Þ hw
Tabs ðNþfÞ
2 3
The useful heat gain by a collector can be
expressed as 6 7
6 ðTabs þ Ta Þ T2abs þ T2a 7
6
þ 6( )7 ð10Þ
1 7
_ P Tfs Tfe
u ¼ mc ð3Þ 4 ð"abs þ 0, 00591 Nhw Þ 5
þ 2Nþ f 1þ0:133"
"v
abs
N
While introducing the collector, the overall loss
coefficient between the absorber and the ambient air
UL, the useful energy gain provided by the collector is f ¼ ð1 þ 0:089hw 0:1166hw "abs Þð1 þ 0:07866NÞ
given by Equation (4) (Duffie and Beckman 1991):
Downloaded by [Anadolu University] at 10:41 26 December 2014
C ¼ 520 1 0:0000512
u ¼ Sabs FR IG ðv abs Þ UL Tfe Ta ð4Þ 1 100
e ¼ 0:430
The heat removal factor, FR is defined by Tabs
Equation (5). In the above equation, Ta is the ambient temper-
ature (K) and Tabs the mean absorber plate tempera-
_ p
mc F0 UL Sabs
FR ¼ 1 exp ð5Þ ture (K).
Sabs UL _ p
mc For 70 hh90 , use ¼ 70 , then the loss coefficient
The collector efficiency factor F 0 and collector through the bottom of the collector is
overall loss coefficient UL for this studied configura- 1
tion, shown in Figure 1, are obtained from the energy Uar ¼ ð11Þ
eis
is þ ebb þ h1vv
balances on the absorber plate, the fluid and the back
plate like in Duffie and Beckman (1991), as shown in The outlet air temperature of the collector can be
(Figure 2), where ‘Ssup ’ represents the additional area obtained from an energy balance as (Duffie and
of artificial roughness Beckman 1991):
habsf Uar þ Ssup hrabspl þ hrabspl þ Ssup habsf IG ðv abs Þeff UL Tfs Ta S F 0 UL
F0 ¼ ¼ exp abs
Uav IG ðv abs Þeff UL Tfe Ta _ f
mcp
þ hrabspl þ habsf
Uar þ Ssup habsf þ hrabspl h2rabspl ð12Þ
ð6Þ
The mean temperatures of the absorber plate are
8 9 obtained by solving the energy balance equations on
> ðUav þ Uar Þ > these plates, where the quantity ‘Ssup ’ takes into
< S h2
> >
=
sup absf þ habsf hrabspl account the fin area, which is supposed with the
>
> þS sup h
absf hrabspl
>
> galvanised metal plate temperature.
: ;
þUar Uav Ssup þ 1 habsf
UL ¼ ð7Þ 1 þ Ssup Tf þ Sabs hðabsf
u
Tpl
habsf hrabspl þ Ssup habsf Uav Tabs ¼ Þ
ð13Þ
þSsup habsf hrabspl þ Ssup h2absf Ssup
this assumption does not reflect the reality recorded for 5. Conclusion
the plate absorbing because of its capacity of significant In this study, we made theoretical and experimental
storage and which has a considerable impact on the studies on a flat plate solar collector. For those, a flat
results obtained. This quantity of stored heat which plate air solar collector was constructed of a surface of
represents the thermal inertia of the system compensates collecting equal to 1.73 m2 with the dynamic air vein
the fall of the solar radiation at the end of the day what furnished with thin metal obstacles and tested at the
is quite visible in the curves of Figure 5(b) as well during laboratory under the climatic conditions of the area of
the disturbances in Figures 4(a), (b), 5(a) and (b). Biskra.
These curves provide some information on the The comparison between the results obtained from
quality of the thermal transfer in the solar collector the theoretical and those from the experimental one
and allow us to evaluate the coefficients of convective makes it possible to make the following conclusions:
transfer what makes it possible to calculate the
(1) The theoretical approach presented in this
quantities of useful heat and lost by the system. This
article translates in a satisfactory way the
is very significant because it is starting from these data
thermal performances of the flat plate air
that one will be able to act on such parameter to
solar collector used obstacles in the dynamic
optimise the performances of such a converter of
vein.
Downloaded by [Anadolu University] at 10:41 26 December 2014
energy.
(2) The theoretical approach does not hold in
In Figure 6, we present the variations of the
account of the energy stored which is also true
efficiency of the collector during the test day. The
for the air with its physical properties poor but
average points of the experimental data are shown in
with regard to the components of the collector
the figure. Scatter of the data around the line of
(mainly the absorber), this assumption has a
average interpolation is mainly attributed to the significant and considerable impact on the
disturbance of the radiation incident and the speed of results obtained.
the wind. Figure 6 shows that there are two transient
regimes relating to the sunrise and the sunset. This is
due to the transient behaviour of the solar flat plate
Acknowledgements
collector as was already observed by other authors
The authors acknowledge the suggestions and the technical
(Esen 2008). After the system stabilisation, the effi- assistance of Belhi Guerira who was responsible for the
ciency will be located around a mean value depending technological hall at the Department of Mechanical
on environmental conditions. The efficiency calculated Engineering in Mohamed Khider university of Biskra,
for the configuration of the tested collector is accept- Algeria.
able and reached 70%. Otherwise, the efficiency
evolution is in agreement with other authors’ result,
some of them were mentioned above (Karsli 2007). Nomenclature
ach length of obstacles (m)
bch absolute height of the obstacles (m)
cp specific heat of air (J kg1 K1)
Dh hydraulic diameter (m)
e height of air tunnel in solar collector
(m)
eb thickness of the insulating (wood) (m)
eis thickness of the insulating (polysty-
rene) (m)
FR heat removal factor of solar collector
F0 efficiency factor of solar collector
G air mass flow rate (kg h1)
habsf convection heat transfer coefficient
between the absorber plate and air
(W m2 K1)
hplf convection heat transfer coefficient
between the channel back and air
(W m2 K1)
hrabspl radiation heat transfer coefficient
Figure 6. Efficiency of the collector during the time of day between the absorber plate and the
tested. channel back (W m2 K1)
International Journal of Ambient Energy 101
hvv ¼5:67þ3:86Vv convection heat transfer coefficient f cinematic viscosity of air (kg m1 s1)
caused by wind (W m2 K1) constant of Stefan–Boltzmann
IG global irradiance incident on solar air v transparent cover transmittance (0.90)
heater collector (W m2) f dynamic viscosity of air (m2 s1)
kf thermal conductivity of fin
(W m1 K1)
Lc length of the flat plate collector (m)
lc width of the flat plate collector (m)
m_ air mass flow rate (kg s1) References
N number of cover
n1 total number of the obstacles in the Abdullah, A.H., Abou-Ziyan, H.Z., and Ghoneim, A.A.,
collector, equal to 296 2003. Thermal performance of flat plate solar collector
n2 number of the obstacles in one row, using various arrangements of compound honeycomb.
equal to 8 Energy Conversion and Management, 44, 3093–3112.
g global incidence solar energy (W) Ahmed-Zaid, A., et al., 1999. Experimental study of thermal
u useful energy gain of solar air collector performance improvement of a solar air flat plate collector
(W) through the use of obstacles: application for the drying of
l loss energy of the collector (W) ‘yellow onion’. International Journal of Energy Research,
Downloaded by [Anadolu University] at 10:41 26 December 2014
Moummi, N., et al., 2004. Energy analysis of a solar air Sacadura, J.F., 1980. Initiation aux transfert thermiques.
collector with rows of fins. Renewable Energy, 29, Cast. INSA of lyon, Paris: Technique et Documentation.
2053–2064. Sopian, K., et al., 2009. Evaluation of thermal efficiency of
Moummi, N., et al., 2010. Systematic forecasts of solar double-pass solar collector with porous–nonporous media.
collector’s performance in various sites of different Renewable Energy, 34, 640–645.
climates in Algeria. International Journal of Sustainable Webb, R.L., Eckert, E.R.G., and Goldstein, R.J., 1971.
Energy, 29 (3), 142–150. Heat transfer and friction in tubes with repeated rip
Ouard, S., 1989. Optimisation des formes et disposition roughness. International Journal Heat and Mass Transfer,
d’obstacles dans la veine mobile du fluide des capteurs 14, 601–617.
solaires plans à deux couches d’air en vue de la maximisation Yeh, H.M., Ho, C.D., and Lin, C.Y., 1998. The influence of
du couple rendement et tempe´rature. PhD thesis, University collector aspect ratio on the collector efficiency of baffled
of Valenciennes, France. solar air heaters. Energy, 23 (1), 11–16.
Ozgen, F., Esen, M., and Esen, H., 2009. Experimental Youcef-Ali, S. and Desmons, J.Y., 2006. Numerical and
investigation of thermal performance of a double-flow experimental study of a solar equipped with offset
solar air heater having aluminium cans. Renewable Energy, rectangular plate fin absorber plate. Renewable Energy,
34, 2391–2398. 31, 2063–2075.
Prasad, B.N. and Mullick, S.C., 1985. Heat transfer Zugary, M.R. and Vullierne, J.J., 1993. Amélioration des
characteristics of a solar air heater used drying purposes. performances thermiques d’un capteur solaire par l’utili-
Downloaded by [Anadolu University] at 10:41 26 December 2014
Applied Energy, 13, 83–93. sation d’une structure à lamelles. Entropie, 176, 25–30.