Renewable Energy: K. Sopian, M.A. Alghoul, Ebrahim M. Alfegi, M.Y. Sulaiman, E.A. Musa
Renewable Energy: K. Sopian, M.A. Alghoul, Ebrahim M. Alfegi, M.Y. Sulaiman, E.A. Musa
Renewable Energy: K. Sopian, M.A. Alghoul, Ebrahim M. Alfegi, M.Y. Sulaiman, E.A. Musa
Renewable Energy
journal homepage: www.elsevier.com/locate/renene
a r t i c l e i n f o a b s t r a c t
Article history: The double-pass solar collector with porous media in the lower channel provides a higher outlet
Received 17 September 2007 temperature compared to the conventional single-pass collector. Therefore, the thermal efficiency of the
Accepted 12 May 2008 solar collector is higher. A theoretical model has been developed for the double-pass solar collector. An
Available online 30 June 2008
experimental setup has been designed and constructed. The porous media has been arranged in different
porosities to increase heat transfer, area density and the total heat transfer rate. Comparisons of the
Keywords: theoretical and the experimental results have been conducted. Such comparisons include the outlet
Solar
temperatures and thermal efficiencies of the solar collector for various design and operating conditions.
Collector
Double-pass
The relationships include the effect of changes in upper and lower channel depth on the thermal effi-
Porous–nonporous media ciency with and without porous media. Moreover, the effects of mass flow rate, solar radiation, and
Thermal efficiency temperature rises on the thermal efficiency of the double-pass solar collector have been studied. In
addition, heat transfer and pressure drop relationships have been developed for airflow through the
porous media. Close agreement has been obtained between the theoretical and experimental results. The
study concluded that the presence of porous media in the second channel increases the outlet tem-
perature, therefore increases the thermal efficiency of the systems.
Ó 2008 Elsevier Ltd. All rights reserved.
0960-1481/$ – see front matter Ó 2008 Elsevier Ltd. All rights reserved.
doi:10.1016/j.renene.2008.05.027
K. Sopian et al. / Renewable Energy 34 (2009) 640–645 641
cost effective. The use of a double-pass resulted in increasing the from experimental studies used for generating the required ex-
pressure drop across the collector. However, the increase in the perimental data for validating the model are presented. It is shown
operating cost due to the increased pressure drop in the collector is that the model can predict the performance of the MNCSCD fairly
considered small. This is due to the fact that the pressure drop accurately and therefore can be used as a design tool for prototype
across the collector is small compared to the total pressure drop development.
across the system. Six different types of natural circulation air-heating solar col-
Mohamad [15], who referred the conventional double-pass lectors were designed, constructed and analysed by Koyuncu [18]
collector as a counter-current solar collector, showed that the for their performance. Each collector mainly consisted of a frame
thermal efficiency can be improved by 18% compared to the con- constructed from hardboard, vent holes, hardboard insulation,
ventional solar heater. The study also suggested an extra fan power absorbing surface made of black coated aluminium sheet and clear
of 2–3 W, which is not high. Also noted that the cost of construction plastic glazing. All solar air heaters were tested simultaneously
of the double glazing collector is comparable to the cost of the under the same environmental conditions. The experimental setup
double-pass collector. However, necessary trade-off between the was instrumented for the measurement of solar radiation, tem-
fan power and the efficiency of packed bed solar collector must be perature and relative humidity of the atmosphere air, outlet air
analyzed to obtain a cost-effective design configuration. temperature, surface temperature of the back and edge insulator
The thermal performance of a double-glass double-pass solar air and absorber plate, air speed and wind velocity. It is understood
heater with a packed bed (DPSAHPB) above the heater absorber from the results of the investigation that the performances of
plate was investigated experimentally and theoretically by Rama- Model-1, Model-2, Model-3, Model-4, Model-5 and Model-6 are
dan et al. [16]. Limestone and gravel were used as packed bed 42.11%, 45.88%, 44.23%, 39.76%, 39.05% and 36.94%, respectively,
materials. Numerical calculations were carried out, on typical and the performance of the most efficient collector (Model-2) is
summer days of 2003, to study the effect of different operational approximately 9% more than the least efficient one (Model-6). In
and configurational parameters on the heater performance. Effects addition, it is seen that unlike number of glazing sheet and air pass
of the mass flow rate of air and the mass and porosity of the packed method, the effect of the shape of the absorbing surface on the
bed material were also studied. It was inferred that for increasing performance is considerably less.
the outlet temperature Tfo of the flowing air after sunset, it is ad- In this paper, a theoretical model of a double-pass solar collector
visable to use the packed bed materials with higher masses and with porous media in the second channel has been commenced and
therefore with low porosities. It is recommended to operate the compared with indoor testing facility results. The testing facility
system with packed bed with values of equal 0.05 kg/s or lower to consists of a solar collector and a solar simulator.
have a lower pressure drop across the system.
A mathematical model for drying agricultural products in 2. Mathematical model
a mixed-mode natural convection solar crop dryer (MNCSCD) using
a single-pass double-duct solar air heater (SPDDSAH) is discussed The unsteady-state mathematical model has been developed.
by Forson et al. [17]. The model comprises the air-heating process This involves unsteady-state energy balance equations linking the
model, the drying model and the technical performance criteria outer glass cover heat transfer coefficient, and the heat transfer
model. The governing equations of the drying air temperature coefficients between the moving airstreams and surfaces forming
and humidity ratio; the material temperature and its moisture the upper and lower channels as shown in Fig. 1(a) and (b). To
content; and performance criteria indicators are derived. The simplify the analysis, the energy balance equations are written
model requires the solution of a number of interrelated non-linear under the following assumptions: (a) the temperatures of the cover
equations and a set of simultaneous differential equations. Results and plates vary only in the direction of fluid flow (x-direction);
642 K. Sopian et al. / Renewable Energy 34 (2009) 640–645
vTb
Mb Cb ¼ hrpb Tp Tb hcbf2 Tb Tf2 UR ðTb Ta Þ
vTt
v2 T
þ kb db 2b (5)
vx
The boundary conditions are obtained from the conditions that
there is no heat loss from the side of the metallic plates. One of the
boundary conditions state that at entry point, air temperature
equals ambient temperature such as
vTf1
Tf1 ð0; tÞ ¼ Ta ðtÞ; ¼ 0
vx x¼ðL;tÞ
vTg vTp vT
For x ¼ 0 and x ¼ L; ¼ ¼ b ¼ 0
vx vx vx
2=3 ! 0:14
2=3 1=3 Dh m Halogen lamps
Nu ¼ 0:116 Re 125 Pr 1þ
L mw
Heater
(12)
where m is evaluated at film temperature and mw is evaluated at wall Inlet
Collector
temperature.
(c) Turbulent flow regime (Re > 6000):
mf Cf
Pr ¼ (14)
kf
Fig. 2. Schematic of the experimental setup with the solar simulator.
The determination of the average heat transfer coefficient, h,
between the porous media and air as follow:
3. Experimental setup
T o T i
h ¼ Gp Cp (15) Fig. 2 shows the solar simulator and the collector undergoing
A T m T a
testing. The simulator uses 45 halogen lamps, each with a rated
Pressure drop or lost head is directly proportional to the length power of 300 W. The maximum average radiation of 642 W/m2 can
of the duct, proportional to the square of the flow rate, and pro- be reached. Dimmers are used to control the amount of radiation
portional to the fifth power of the duct size. Therefore, the duct- that the test collector received. The dimmers are divided into six
work designer can be relatively unconcerned about the length of scales for producing different amount of radiation values. These
the run, only moderately concerned with the circulation rate, but values have been previously measured using the pyranometer. The
must be extremely sensitive that the size of the duct is appropriate measurement errors are about 3.16% for radiation value of 277.8 W/
for the flow rate. Consequently, the pressure drop through the m2 and 4.05% for radiation value of 642 W/m2 [20]. A heater is
collectors is of highest interest since that is where the minimum placed at the inlet of the collected undergoing test to vary the inlet
dimensions are most likely to be found. We now discuss the char- temperature.
acteristics of friction factor in fluid in duct runs in the collector. An The collector consists of the glass cover, the insulated container
important fundamental relationship is the Fanning equation, given and the black painted aluminum absorber. The size of the collector
here in a modified form as is 120 cm in width and 240 cm in length. The first and second
channels can be adjusted for optimal operations. The inlet tem-
4fG2f L perature to the collector can be adjusted by heating the inlet air to
DP ¼ (16)
2gc A2x rDh the collector. Thermocouples are located strategically at the inlet,
end-of-the first pass, outlet, absorber plate and glass cover. The
where DP is the frictional loss or pressure drop, Gf is the fluid mass temperature measurements are recorded using data acquisition
flow rate, L is the length of the duct, Ax is cross-sectional area, gc is system. The flow rates are measured using the vane type
a constant (1 kg m/N s2 or 32.17 ft/s2), r is the fluid density, and f is anemometer.
a friction factor. The friction factor, f, is as empirical function of the
relative roughness of the duct.
For smooth duct, when the flow is laminar the friction factor is
given by 100 100
16 90 =70% 90
f ¼ (17)
THERMAL EFFICIENCY (%)
Re 80 Theoretical 80
4. Experimental procedures
0.25 450W/m
15 60
50 0.2 Experimental
10 40 0.15
30
0.1
5 T: Without porous media
20 Nuap Nubp Nucp Nut,
0.05
10
: Without porous media
0 0 0
0.037 0.04 0.053 0.059 0.079 5.2 6.2 7.2 8.2 9.2 10.2 11.2
MASS FLOW RATE (kg/sec) REYNOLDS NUMBER Rex10-3
Fig. 5. Effect of the solar radiation on the thermal efficiency on the double-pass solar Fig. 8. Comparison between the experimental and theoretical Nusselt number of the
collector with porous media (f ¼ 80%, Ta ¼ 33.5 C). double-pass solar collector with porous media.
K. Sopian et al. / Renewable Energy 34 (2009) 640–645 645
0.2 Fig. 9 shows the effect of the Reynolds number on the friction
factor. Porous media in the lower channel can be used to increase
the heat transfer coefficient since the friction factor depends on the
0.16
Without porous media
losses and velocity of the airflow rate.
FRICTION FACTOR (f)
0.12
7. Conclusion
quite close between the theoretical and experimental values for References
solar radiation less than 570 W/m2, but the deviation becomes
[1] Hollands KGT, Shewan EC. Optimization of flow passage geometry for air-
more pronounced for solar radiation greater than or equal to heating, plate-type solar collectors. Transactions of ASME, Journal of Solar
590 W/m2. Also, it can be seen that the thermal efficiency increases Energy Engineering 1981;103:323–30.
with an increase in the solar radiation. [2] Choudhury C, Garg HP. Design analysis of corrugated and flat plate solar air
heaters. Renewable Energy 1991;1(5/6):595–607.
Fig. 5 shows the effect of mass flow rate on the thermal effi- [3] Hachemi A. Thermal performance enhancement of solar air heaters, by fan-
ciency and temperature rise when the porosity is 80%. The collector blown absorber plate with rectangular fins. Energy Research 1995;19(7):567–
thermal efficiency increases with an increase in the mass flow rate, 78.
[4] Yeh HM, Ho CD, Hou JZ. The improvement of collector efficiency in solar air
while the temperature rise decreases with the corresponding heaters by simultaneously air flow over and under the absorbing plate. Energy
increase in the mass flow rate. The collector thermal efficiency 1999;24(10):857–71.
increases as a result of a decrease in the average temperature of the [5] Hegazy AA. Performance of flat plate solar air heaters with optimum channel
geometry for constant/variable flow operation. Energy Conversion and
absorber plate. The latter arises because of the increase in the flow Management 2000;41(4):401–17.
rate cause by an increase in the heat transfer coefficient between [6] Yeh HM, Ho CD, Lin CY. Effect of collector aspect ratio on the collector
the channel walls and the air. As is evidence in the figure, the efficiency of upward type baffled solar air heaters. Energy Conversion and
Management 2000;41(9):971–81.
efficiency of the solar collector with porous media is about 10%
[7] Swartman RK, Ogunlade O. An investigation on packed-bed collectors. Solar
better than without porous media. This is because of the increase in Energy 1966;10:106–10.
the area of the heat transfer analogous to the increase in heat [8] Collier RK, Arnold FH. Comparison of transpired beds for solar collector
applications. In: Proceedings of the annual ISES meeting; 1980. p. 451–5.
transfer by using extended surfaces such as fins and corrugated
[9] Sorour MM, Hassab MA. A screen type solar air heater. In: Proceedings of the
absorber plate. 8th international heat transfer conference; 1986. p. 3097–103.
Fig. 6 shows the effect of temperature rise on the thermal effi- [10] Lansing FL, Clarke V, Reynold R. A high performance porous flat-plate solar
ciency of double-pass solar collector with media of 90% porosity. collector. Energy 1979;4:685–94.
[11] Zaid AA, Messaoudi H, Abenne A, Ray ML, Desmons JY, Abed B. Experimental
The thermal efficiency increases with the increase in the temper- study of thermal performance improvement of a solar air flat plate collector
ature rise. through the use of obstacles: application for the drying of ‘‘yellow onion’’.
Fig. 7 shows the effect of solar incident on the temperature rise Energy Research 1999;23(12):1083–99.
[12] Moummi N, Ali SY, Moummi A, Desmons JY. Energy analysis of a solar air
and thermal efficiency. The collector thermal efficiency increases collector with rows of fins. Renewable Energy 2004;29(13):2053–64.
with an increase in the solar radiation, while the temperature rise [13] Mohamad AA. High efficiency solar air heater. Solar Energy 1997;60(2):71–6.
decreases with increasing solar radiation. Also, it can seen that the [14] Sopian K, Supranto WR, Daud W, Yatim B, Othman MY. Thermal performance
of the double-pass solar collector with and without porous media. Renewable
solar collector with porous media has a higher thermal efficiency Energy 1999;18:557–64.
than the collector without porous media for an increase in solar [15] Mohamad AA. Counter-current solar air heater. In: Proceedings of the first
radiation. Therefore, introducing porous media in the second Trabzon international energy and environmental symposium; 1996. p. 129–34.
[16] Ramadan MRI, El-Sebaii AA, Aboul-Enein S, El-Bialy E. Thermal performance of
channel increases the heat transfer area correspondingly increasing
a packed bed double-pass solar air heater. Energy 2007;8:1524–35.
the heat transfer coefficient. [17] Forson FK, Nazha MAA, Rajakaruna H. Modelling and experimental studies on
Fig. 8 shows the theoretical and experimental results of Rey- a mixed-mode natural convection solar crop-dryer. Solar Energy 2007;81(3):
346–57.
nolds number and Nusselt number of the double-pass solar col-
[18] Koyuncu Turhan. Performance of various designs of solar air heaters for crop
lector with porous media. Using porous media and for increasing drying applications. Renewable Energy 2006;31(7):1073–88.
solar radiation the effect on the performance of the solar collector is [19] Niles PN, Caxnigie EJ, Pohl JG, Cherne JM. Design and performance of an air
improved. This is because increasing the heat transfer area leads to collector for industrial crop dehydration. Solar Energy 1979;20:19–28.
[20] Supranto K, Sopian WR, Daud W, Othman MY, Yatim B. Design and con-
an increase in the heat transfer coefficient. Also the same effect is struction of a low cost solar simulator. In: Proceeding of the world renewable
seen for a decrease in the thermal conductivity. energy congress 1999 Malaysia; 1999. p. 153–6.