02.03.2014 Views

the thermal performance of a solar air heater - Eurotherm 2008

the thermal performance of a solar air heater - Eurotherm 2008

the thermal performance of a solar air heater - Eurotherm 2008

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

5th European Thermal-Sciences Conference, The Ne<strong>the</strong>rlands, <strong>2008</strong><br />

THE THERMAL PERFORMANCE OF<br />

A SOLAR AIR HEATER<br />

J. Deans, A. Weerakoon<br />

The University <strong>of</strong> Auckland, Auckland, New Zealand<br />

Abstract<br />

The use <strong>of</strong> commercial ro<strong>of</strong>ing panels to harvest <strong>the</strong> <strong>solar</strong> radiation received by a building wall is<br />

becoming a common design feature <strong>of</strong> modern buildings. This energy is used to provide<br />

background <strong>air</strong> heating within <strong>the</strong> building and experience has shown that <strong>the</strong> ambient <strong>air</strong> which<br />

flows through an array <strong>of</strong> holes pierced in <strong>the</strong> ro<strong>of</strong>ing panel can extract more than 60% <strong>of</strong> <strong>the</strong><br />

incident <strong>solar</strong> radiation. The <strong>air</strong> passing through <strong>the</strong> panel flows into <strong>the</strong> <strong>air</strong> gap behind <strong>the</strong> panel.<br />

In this heating system <strong>the</strong> panel and <strong>the</strong> gap behind it effectively form a single stream heat<br />

exchanger with variable mass flow. The <strong>air</strong> warmed by <strong>the</strong> panel is transferred to <strong>the</strong> building’s<br />

heating system.<br />

The simplicity <strong>of</strong> <strong>the</strong> heat exchanger design masks a number <strong>of</strong> coupled heat transfer processes<br />

which can make <strong>the</strong> panel design unduly complicated. This study describes <strong>the</strong> experimental and<br />

analytical investigations which were focussed on <strong>the</strong> evaluation <strong>of</strong> <strong>the</strong>se processes. A large <strong>solar</strong><br />

heated panel was used in most <strong>of</strong> <strong>the</strong> experimental studies while both finite difference and finite<br />

volume models were used in <strong>the</strong> analysis. The investigation <strong>of</strong> <strong>the</strong> heat and momentum flows on<br />

ei<strong>the</strong>r side <strong>of</strong> <strong>the</strong> panel, using <strong>the</strong> CFX finite volume model, has enabled new insights to be gained<br />

about <strong>the</strong> heat transfer mechanisms between <strong>the</strong> panel and <strong>the</strong> <strong>air</strong>.<br />

1 Introduction<br />

The use <strong>of</strong> commercial ro<strong>of</strong>ing panels to extract some <strong>of</strong> <strong>the</strong> <strong>solar</strong> energy impinging on a building’s<br />

surface and to use this energy to heat <strong>the</strong> <strong>air</strong> circulating within <strong>the</strong> building has been actively<br />

investigated. One such system, commonly known as unglazed transpired <strong>solar</strong> collectors was<br />

originally developed by Conserval Inc.,[2002] and marketed by <strong>the</strong>m using <strong>the</strong> Solarwall name<br />

[Hollick, 1996]. The arrangement <strong>of</strong> this collector is shown in Figure 1(a) where <strong>the</strong> <strong>solar</strong> heated,<br />

metal panel has an array <strong>of</strong> small holes, pierced through <strong>the</strong> panel to form a permeable surface. The<br />

<strong>air</strong> being heated by <strong>the</strong> panel flows through <strong>the</strong> holes and into <strong>the</strong> <strong>air</strong> gap behind, this <strong>air</strong> is collected<br />

from a manifold at <strong>the</strong> top <strong>of</strong> <strong>the</strong> <strong>air</strong> gap and directed towards <strong>the</strong> suction fan. This fan powers <strong>the</strong><br />

complete system by drawing ambient <strong>air</strong> from <strong>the</strong> front <strong>of</strong> <strong>the</strong> <strong>solar</strong> collector, through <strong>the</strong> <strong>air</strong> gap,<br />

before delivering it to <strong>the</strong> building’s heating system.<br />

The simplicity <strong>of</strong> <strong>the</strong> panel’s design masks a number <strong>of</strong> complex heat transfer processes and a<br />

summary <strong>of</strong> <strong>the</strong>se processes are shown in Figure1 (b). Solar radiation heats <strong>the</strong> dark metal panel<br />

and <strong>the</strong> low <strong>the</strong>rmal resistance <strong>of</strong> this panel suggests that energy will be transferred by convection<br />

to <strong>the</strong> <strong>air</strong> flowing on both sides <strong>of</strong> <strong>the</strong> panel. The <strong>air</strong> approaching <strong>the</strong> front <strong>of</strong> <strong>the</strong> panel will flow<br />

radially towards each hole and that fraction <strong>of</strong> <strong>the</strong> flow which is adjacent to <strong>the</strong> panel’s outer<br />

surface will be convectively heated. The increasing <strong>air</strong> flow behind <strong>the</strong> panel will also be<br />

convectively heated and this transport process will be enhanced by <strong>the</strong> jets <strong>of</strong> warm <strong>air</strong> discharging<br />

from <strong>the</strong> back <strong>of</strong> <strong>the</strong> panel. The main energy losses from this heat exchanger develop at <strong>the</strong> outer<br />

panel surface and <strong>the</strong>se are surface re-radiation and <strong>the</strong> surface cooling by wind. The geometry <strong>of</strong><br />

<strong>the</strong> panel pr<strong>of</strong>ile and <strong>the</strong> variation in <strong>the</strong> rates <strong>of</strong> <strong>solar</strong> radiation give added complexity to <strong>the</strong><br />

analysis <strong>of</strong> this heat exchanger. The efficiency <strong>of</strong> energy transfer from <strong>the</strong> collector to <strong>the</strong> <strong>air</strong><br />

flowing through it is generally high (≈ 60%); Rhee and Edwards [1981] attributed this to <strong>the</strong> low<br />

convective and radiation loses from <strong>the</strong> front <strong>of</strong> <strong>the</strong> collector.


5th European Thermal-Sciences Conference, The Ne<strong>the</strong>rlands, <strong>2008</strong><br />

PANEL<br />

Solar Radiation<br />

Radiation Loss<br />

Wind Loss<br />

AIR GAP<br />

Increasing<br />

Mass Flow<br />

Rate<br />

Internal<br />

Convection<br />

External<br />

Convection<br />

(a) Physical layout.<br />

Figure 1 Arrangement <strong>of</strong> a transpired <strong>solar</strong> collector.<br />

(b) Mass and energy flow routes.<br />

The <strong>performance</strong> <strong>of</strong> transpired <strong>solar</strong> collectors has been studied by Kutscher [1993, 1994], Van<br />

Decker [2001], Gunnewiek [1996, 2002] and Arulanandam [1999] who have conducted a wide<br />

range <strong>of</strong> experimental and computational studies. A major concern with most <strong>of</strong> <strong>the</strong>se studies is<br />

that <strong>the</strong> test geometries modelled tend to over simplify <strong>the</strong> problem and consequently, limit <strong>the</strong><br />

solution. The test rig developed by Kutscher, for example, was only a small laboratory based<br />

collector which used a flat panel and no <strong>air</strong> gap. The exclusion <strong>of</strong> <strong>the</strong> <strong>air</strong> gap from this model<br />

inherently assumes that all <strong>of</strong> <strong>the</strong> energy transfer between <strong>the</strong> panel and <strong>the</strong> <strong>air</strong> will take place from<br />

<strong>the</strong> front <strong>of</strong> <strong>the</strong> panel. These studies generally demonstrate that <strong>the</strong> <strong>the</strong>rmal <strong>performance</strong> <strong>of</strong> <strong>the</strong><br />

transpired collector is dependent on <strong>the</strong> diameter <strong>of</strong> <strong>the</strong> hole through <strong>the</strong> panel, <strong>the</strong> distance between<br />

holes, <strong>the</strong> panel thickness, <strong>the</strong> suction flow rate <strong>of</strong> <strong>the</strong> <strong>air</strong> and wind speed. The model developed to<br />

describe <strong>the</strong> interaction <strong>of</strong> <strong>the</strong>se parameters assumed that <strong>the</strong> holes effectively “harvested” <strong>the</strong><br />

<strong>the</strong>rmal boundary layer from <strong>the</strong> outer surface <strong>of</strong> <strong>the</strong> panel. Both Gunnewiek [1996, 2002] and Van<br />

Decker [2001] did however recognise that heat transfer from <strong>the</strong> panel to <strong>the</strong> <strong>air</strong> would also develop<br />

within <strong>the</strong> <strong>air</strong> gap but nei<strong>the</strong>r had <strong>the</strong> test data required to develop this alternative method.<br />

The models developed by Kutscher [1993, 1994] and Van Decker [2001] do have wide recognition<br />

but Fleck [2002] has observed that “it is quite clear from our observations that such a model is a<br />

poor descriptor <strong>of</strong> <strong>the</strong> physical phenomena driving UTSC (Unglazed Transpired Solar Collector)<br />

<strong>performance</strong>”. The paucity <strong>of</strong> test information from operational unglazed transpired <strong>solar</strong> collectors<br />

is surprising since many operating systems have been commissioned, <strong>the</strong> only studies on large<br />

collectors are those published by Carpenter [1991] Fleck [2002] and Weerakoon 2006]. The test<br />

programme <strong>of</strong> Weerakoon [2007] examined <strong>the</strong> <strong>performance</strong> <strong>of</strong> transpired <strong>solar</strong> heated collectors<br />

constructed from standard size ro<strong>of</strong>ing panels and was directed towards <strong>the</strong> determination <strong>of</strong> <strong>the</strong><br />

effects which <strong>the</strong> main design variables had on <strong>the</strong> panel’s <strong>performance</strong>.<br />

2 Experiment set-up<br />

The main test rig was constructed to accommodate a standard 2.0m * 0.8m corrugated metal ro<strong>of</strong>ing<br />

panel that could be used as <strong>the</strong> <strong>solar</strong> collector. The complete test facility is shown in Figure 2<br />

where <strong>the</strong> panel was mounted to create a small uniform, sealed gap between <strong>the</strong> panel and <strong>the</strong> wall<br />

behind. The <strong>air</strong> which migrated through <strong>the</strong> panel flowed up this gap and its depth was normally set<br />

at 20 mm. The orientation <strong>of</strong> <strong>the</strong> vertical panel could be changed but was usually north facing and<br />

located on <strong>the</strong> exposed ro<strong>of</strong> <strong>of</strong> a 20 m high building.


5th European Thermal-Sciences Conference, The Ne<strong>the</strong>rlands, <strong>2008</strong><br />

Radiometer<br />

Space bars to<br />

create cavity<br />

Outlet <strong>air</strong> manifold<br />

Solar<br />

panel<br />

Wea<strong>the</strong>r<br />

shield<br />

Solar<br />

panel<br />

Thermocouple<br />

points<br />

Wheels<br />

Data logger<br />

Back wall<br />

Linear flow element<br />

Ballast<br />

Exhaust <strong>air</strong><br />

Blower fan<br />

(a) Front <strong>of</strong> <strong>the</strong> test rig.<br />

(b) Schematic rear view <strong>of</strong> <strong>the</strong> test rig.<br />

Figure 2 The Test Facility.<br />

The <strong>air</strong> transferred through <strong>the</strong> panel was collected in <strong>the</strong> outlet <strong>air</strong> manifold shown in Figure 2(b)<br />

before it was ducted to a laminar flow element and <strong>the</strong>n to a suction fan before return to <strong>the</strong><br />

environment. A differential pressure transducer within <strong>the</strong> laminar flow element measured <strong>the</strong> <strong>air</strong><br />

flow rate through <strong>the</strong> system and calibrated platinum resistance <strong>the</strong>rmometers measured both <strong>the</strong><br />

ambient and <strong>the</strong> manifold outlet <strong>air</strong> temperatures. The incident <strong>solar</strong> radiation (İ <strong>solar</strong> ) on <strong>the</strong> metal<br />

sheet was measured using an Eppley precision radiometer and <strong>the</strong> surface temperature <strong>of</strong> <strong>the</strong> panel<br />

was monitored using an Agema Thermovision 570 infrared camera.<br />

The transpired <strong>solar</strong> collectors were constructed using 0.8mm gauge, galvanised steel, ro<strong>of</strong>ing<br />

panels. The 1.6 mm diameter holes through <strong>the</strong> panels were created using a hand drill and <strong>the</strong> panel<br />

pr<strong>of</strong>iles used in this study are shown in Figure 3. The external colours <strong>of</strong> <strong>the</strong> panels are described<br />

as iron sand and ivy green; <strong>the</strong> emissivity value for both panels was determined to be 0.88.<br />

The <strong>the</strong>rmometers were calibrated using a reference <strong>the</strong>rmometer traceable to ITS-90. The error in<br />

<strong>the</strong> <strong>the</strong>rmometer readings is believed to be less than ±0.05 o C. The calibrated charts supplied by <strong>the</strong><br />

manufacturers were used to determine <strong>the</strong> <strong>air</strong> flow rate and <strong>the</strong> level <strong>of</strong> <strong>solar</strong> radiation, <strong>the</strong><br />

uncertainty in both <strong>of</strong> <strong>the</strong>se was less than 3%. The uncertainty associated with derived values such<br />

as heat flow rates, is also estimated to be 3%. The accuracy <strong>of</strong> <strong>the</strong> infrared camera was ± 1.0 o C, but<br />

its ability to determine temperature gradients was much lower.<br />

Apriori information about <strong>the</strong> overall heat transfer coefficient for <strong>the</strong> collector panel enabled its<br />

stabilisation period to be determined and <strong>the</strong>rmal response rate calculated. In most tests <strong>the</strong> test<br />

period lasted for 60 mins and <strong>the</strong> results were not accepted during <strong>the</strong> first 15 minutes. All <strong>of</strong> <strong>the</strong><br />

tests described below were conducted on days when <strong>the</strong>re were few, if any, clouds and between <strong>the</strong><br />

hours <strong>of</strong> 10.00 am and 2.30 pm.<br />

(a) Styleline Panel. .(b) Corrugated Panel. (c) Plane Panel.<br />

Figure 3 Types <strong>of</strong> collector panel’s.


5th European Thermal-Sciences Conference, The Ne<strong>the</strong>rlands, <strong>2008</strong><br />

3 Test Results<br />

The experimental programme was structured to identify <strong>the</strong> influence which various design<br />

variables have on <strong>the</strong> collector’s <strong>performance</strong>. The general results presented in Figures 4 and 5<br />

show <strong>the</strong> influence <strong>of</strong> panel porosity, <strong>air</strong> gap depth and <strong>air</strong> suction flow rate on ei<strong>the</strong>r <strong>the</strong> overall<br />

<strong>performance</strong> <strong>of</strong> <strong>the</strong> collector or <strong>the</strong> energy transferred to <strong>the</strong> <strong>air</strong>. The influence which <strong>the</strong>se<br />

variables have on <strong>the</strong> on <strong>the</strong> design <strong>of</strong> <strong>the</strong> collector tend to support <strong>the</strong> argument that it is effectively<br />

a one stream heat exchanger.<br />

The temperature increase for <strong>the</strong> <strong>air</strong> stream passing through <strong>the</strong> panel is shown in Figure 4 where<br />

this increase is measured against <strong>the</strong> incident <strong>solar</strong> radiation and <strong>the</strong> panel’s porosity. These results<br />

were obtained for a range <strong>of</strong> collector panel porosities with a constant <strong>air</strong> mass flow rate <strong>of</strong> 0.029<br />

kg/m 2 s (based on a panel surface area) and with average wind speeds in <strong>the</strong> range 1.0 to 2 m/s. In<br />

<strong>the</strong>se tests <strong>the</strong> porosity <strong>of</strong> <strong>the</strong> panel was increased by decreasing <strong>the</strong> pitch between holes. If <strong>the</strong><br />

collector efficiency <strong>of</strong> <strong>the</strong> panel is defined as <strong>the</strong> ratio <strong>of</strong> energy extracted by <strong>the</strong> <strong>air</strong> stream and <strong>the</strong><br />

radiant energy received by <strong>the</strong> panel <strong>the</strong>n this value is effectively constant for <strong>the</strong> range <strong>of</strong> <strong>solar</strong><br />

radiation measured during <strong>the</strong>se tests. The relative unimportance <strong>of</strong> porosity is significant, since <strong>the</strong><br />

studies conducted by Kutscher [1993, 1994] and Van Decker [2001] correlated <strong>the</strong> size and <strong>the</strong><br />

spacing <strong>of</strong> <strong>the</strong> holes to <strong>the</strong> <strong>performance</strong> <strong>of</strong> <strong>the</strong> transpired <strong>solar</strong> collector.<br />

The effects which <strong>the</strong> <strong>air</strong> mass flow rate has on <strong>the</strong> <strong>performance</strong> <strong>of</strong> <strong>the</strong> plate are shown in Figure 5<br />

where <strong>the</strong> collector efficiency, increase in <strong>air</strong> temperature and <strong>the</strong> depth <strong>of</strong> <strong>the</strong> <strong>air</strong> gap are examined.<br />

These measurements were obtained using ei<strong>the</strong>r an ivy-green coloured corrugated panel or a plane<br />

ironsand coloured panel. The panel’s <strong>performance</strong> is measured in Figure 5(a) as ei<strong>the</strong>r <strong>the</strong> increase<br />

collector efficiency or a reduction in <strong>the</strong> <strong>air</strong> temperature increase when <strong>the</strong> <strong>air</strong> mass flow rate<br />

through <strong>the</strong> collector increases. This relationship evolves because <strong>of</strong> <strong>the</strong> reduction in <strong>the</strong>rmal losses<br />

from <strong>the</strong> collector when <strong>the</strong> average <strong>air</strong> temperature and <strong>the</strong> panel surface temperatures are reduced.<br />

The coupling <strong>of</strong> <strong>the</strong>se three variables permits <strong>the</strong> design <strong>of</strong> <strong>the</strong> panel to be optimised. Some<br />

confidence in <strong>the</strong> results in Figure 5(a) can be gained by extrapolating <strong>the</strong> <strong>air</strong> mass flow rates down<br />

to <strong>the</strong> values measured by Fleck (2002) who reported collector efficiencies in <strong>the</strong> range 28% to 32%<br />

when <strong>the</strong> <strong>air</strong> flow rate was 0.01 kg/m 2 s.<br />

A plain perforated panel was used in <strong>the</strong> tests which examined <strong>the</strong> influence <strong>of</strong> <strong>air</strong> gap depth behind<br />

<strong>the</strong> panel; this choice was made to ensure a well defined minimum value. The results given in<br />

Figure 5(b) show that <strong>the</strong>rmal efficiency <strong>of</strong> <strong>the</strong> collectors used in <strong>the</strong>se tests was not influenced by<br />

changes in <strong>the</strong> <strong>air</strong> gap when its depth is increased from 5 mm to 40 mm. Most <strong>of</strong> <strong>the</strong> commercial<br />

designs [www.<strong>solar</strong>wall.com] for transpired <strong>solar</strong> collectors have <strong>air</strong> gaps which are in <strong>the</strong> range 20<br />

mm to 40 mm.<br />

Air Temperature Increase [ o C]<br />

20.0<br />

16.0<br />

12.0<br />

8.0<br />

4.0<br />

Ironsand coloured Styleline panel<br />

Wind speed 1 - 2 m/s<br />

Air Flow rate 0.029 kg/m 2 s<br />

0.0<br />

400.0 500.0 600.0 700.0 800.0<br />

Solar Radiation [W/m 2 ]<br />

0.2% porosity<br />

0.4% porosity<br />

0.8% porosity<br />

Figure 4 Air temperature rise in an iron sand coloured Styleline panel.


5th European Thermal-Sciences Conference, The Ne<strong>the</strong>rlands, <strong>2008</strong><br />

Collector Efficiency [%]<br />

100.0<br />

80.0<br />

60.0<br />

40.0<br />

Ivy Green corrugated panel<br />

20.0<br />

0.4% porosity<br />

Wind speed < 1m/s<br />

4.0<br />

Solar radiation 750 - 780<br />

0.0<br />

2<br />

0.0<br />

0.0150 0.0180 0.0210 0.0240 0.0270 0.0300<br />

Mass Flow Rate [kg/m 2 s]<br />

20.0<br />

16.0<br />

12.0<br />

8.0<br />

Increase in Air Temperature [ o C]<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

5 mm gap<br />

20 mm gap<br />

40 mm gap<br />

Plane Panel<br />

Porosity 0.2%<br />

Wind speed < 1 m/s<br />

Solar Radiation 700 W/m 2<br />

0<br />

0.01 0.015 0.02 0.025 0.03<br />

Mass Flow Rate (kg/m 2 s)<br />

(a) Change in collector efficiency and outlet<br />

(b) Change <strong>of</strong> collector efficiency<br />

<strong>air</strong> temperature with <strong>air</strong> flow rate.<br />

with gap depth and <strong>air</strong> flow rate.<br />

Figure 5 Collector efficiency.<br />

4 Energy Balance for <strong>the</strong> Collector<br />

The surface temperatures <strong>of</strong> <strong>the</strong> panel during each test were obtained from infrared photographs and<br />

Figure 6 shows typical <strong>the</strong>rmal images obtained for an ivy green <strong>solar</strong> panel with 0.4% porosity.<br />

These images, obtained at three different test conditions, show that <strong>the</strong> panel’s external surface<br />

temperature is rarely uniform and that it generally decreases with increasing height. It should also<br />

be recalled that <strong>the</strong> panel has a low <strong>the</strong>rmal resistance and that <strong>the</strong> mass flow rate up <strong>the</strong> <strong>air</strong> gap also<br />

increases with height. A comparison <strong>of</strong> <strong>the</strong> photographs given in Figures 6(a) and 6(b) shows that<br />

<strong>the</strong> panel’s surface temperature increases when <strong>the</strong> <strong>air</strong> flow rate is decreased. When <strong>the</strong> results in<br />

Figures 6(a) and 6(c) are compared it can be concluded that increasing <strong>the</strong> wind speed over <strong>the</strong><br />

panel will also decrease its surface temperature. The results from this comparison confirm that both<br />

<strong>the</strong> exit <strong>air</strong> temperature and consequently <strong>the</strong> panel’s <strong>the</strong>rmal efficiency will decrease with<br />

increasing wind speed. It should be noted from <strong>the</strong>se images, that although <strong>the</strong> panel receives a<br />

constant heat flux, it is apparent that nei<strong>the</strong>r its external surface temperature nor <strong>the</strong> energy flux<br />

absorbed by <strong>the</strong> panel is constant. The low <strong>the</strong>rmal resistance <strong>of</strong> <strong>the</strong> panel implies that this<br />

approach can also be extended to <strong>the</strong> <strong>air</strong> flow in <strong>the</strong> gap behind <strong>the</strong> panel and to conclude that<br />

nei<strong>the</strong>r constant heat flux nor wall temperature can be assumed in <strong>the</strong> heating <strong>of</strong> this <strong>air</strong>.<br />

Temperature<br />

iso<strong>the</strong>rm<br />

Panel conditions (a) (b) (c)<br />

Air mass flux 0.029 kg/m 2 s 0.019 kg/ m 2 s 0.029 kg/m 2 s<br />

Solar radiation (İ <strong>solar</strong> ) 780 W/m 2 760 W/m 2 780 W/m 2<br />

Wind speed 1.5 m/s 1.5 m/s 2.5 m/s<br />

Figure 6 Thermal images for an ivy green <strong>solar</strong> panel.


5th European Thermal-Sciences Conference, The Ne<strong>the</strong>rlands, <strong>2008</strong><br />

It has been widely recognised (e.g. Kutscher [1993, 1994]) that <strong>the</strong> transpired panel and <strong>the</strong> <strong>air</strong> gap<br />

behind it form a single stream heat exchanger, with variable mass flow. Weerakoon (2007) has<br />

shown analytically that <strong>the</strong> <strong>air</strong> mass flow rate into <strong>the</strong> <strong>air</strong> gap behind <strong>the</strong> panel increases linearly<br />

with height for all <strong>of</strong> <strong>the</strong> panel geometries tested. The information presented in Figure 1(b) showed<br />

that <strong>the</strong> energy transferred to <strong>the</strong> panel surface is given by <strong>the</strong> summation <strong>of</strong> <strong>the</strong> absorbed <strong>solar</strong><br />

energy, <strong>the</strong> re-radiated energy and <strong>the</strong> convective losses from <strong>the</strong> front <strong>of</strong> <strong>the</strong> panel. All <strong>of</strong> <strong>the</strong><br />

components in <strong>the</strong> energy transfer are dependant on <strong>the</strong> surface temperature and <strong>the</strong> images<br />

presented above showed that this temperature varied with height. If it is assumed that <strong>the</strong> energy<br />

absorbed to <strong>the</strong> panel at any particular height is used to increase <strong>the</strong> local enthalpy <strong>of</strong> <strong>the</strong> <strong>air</strong>, in <strong>the</strong><br />

<strong>air</strong> gap, it is <strong>the</strong>n possible to generate a stepwise process to calculate <strong>the</strong> local energy balance and<br />

determine <strong>the</strong> local <strong>air</strong> temperature within <strong>the</strong> <strong>air</strong> gap. An element in this procedure is shown in<br />

Figure 7a. In developing this procedure <strong>the</strong> surface <strong>of</strong> <strong>the</strong> panel and <strong>the</strong> volume <strong>of</strong> <strong>the</strong> <strong>air</strong> gap<br />

behind it are sub-divided into a number <strong>of</strong> horizontal elements and a mass (equation 1) and energy<br />

balance (equation 2) for each element is produced in a stepwise procedure. The mass <strong>of</strong> ambient <strong>air</strong><br />

added to each element is known and <strong>the</strong> <strong>air</strong> temperature leaving <strong>the</strong> last element can be compared<br />

with measured values.<br />

⎛ •<br />

4 4 ⎞<br />

Q element = ⎜ I <strong>solar</strong> h loss (T s _ element T amb ) (T s _ element Tamb<br />

) ⎟<br />

α −<br />

− − ε σ −<br />

× Aelement<br />

(1)<br />

⎝<br />

⎠<br />

( m + m )C T = m C T + m C T + Q (2)<br />

in<br />

add<br />

p<br />

<strong>air</strong><br />

out _ element<br />

in<br />

p<br />

<strong>air</strong><br />

in _ element<br />

add<br />

p<br />

<strong>air</strong><br />

amb<br />

element<br />

Typical values for <strong>the</strong> derived temperature pr<strong>of</strong>iles for <strong>the</strong> <strong>air</strong>, in <strong>the</strong> <strong>air</strong> gap, are shown in Figure<br />

7(b) and were derived using equations 1 and 2. Within this calculation all <strong>of</strong> <strong>the</strong> initial values were<br />

measured; <strong>the</strong> surface temperature for each element (T s element ) was obtained using <strong>the</strong> images from<br />

<strong>the</strong> infra red camera and <strong>the</strong> external heat transfer coefficient (h loss ) was determined from an overall<br />

energy balance. The measured temperatures for <strong>the</strong> <strong>air</strong> stream at three locations within <strong>the</strong> <strong>air</strong> gap<br />

are also shown in this Figure 7(b). It is apparent from <strong>the</strong>se results that <strong>the</strong> <strong>air</strong> temperature within<br />

<strong>the</strong> <strong>air</strong> gap increases rapidly in <strong>the</strong> lower part <strong>of</strong> <strong>the</strong> gap when <strong>the</strong> <strong>air</strong> flow rate is relatively small.<br />

The <strong>air</strong> flowing up this gap can attain 75% <strong>of</strong> <strong>the</strong> exit <strong>air</strong> temperature within <strong>the</strong> first 0.3m <strong>of</strong> its<br />

height. The form <strong>of</strong> <strong>the</strong> temperature pr<strong>of</strong>ile shown in Figure 6(a) demonstrates that <strong>the</strong> energy<br />

transfer in <strong>the</strong> upper parts <strong>of</strong> <strong>the</strong> <strong>air</strong> gap is dominated by <strong>the</strong> mixing <strong>of</strong> <strong>the</strong> two <strong>air</strong> streams.<br />

m in + m add<br />

T out_element<br />

1<br />

m add<br />

T amb<br />

I <strong>solar</strong><br />

T s_element<br />

A element<br />

m in<br />

T in_element<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0<br />

Air Flow Rate Insolation<br />

0.013 kg/m 2 s 770 W/m 2<br />

0.019 kg/m 2 s 760 W/m 2<br />

0.023 kg/m 2 s 750 W/m 2<br />

0.029 kg/m 2 s 740 W/m 2<br />

Corrugated panel, 0.4% porosity<br />

Wind speed 2.0 -2.5m/s<br />

0 0.5 1 1.5 2<br />

Distance from bottom <strong>of</strong> panel (m)<br />

(a) Element from <strong>the</strong> model.<br />

(b) Calculated and measured <strong>air</strong>stream<br />

temperatures within <strong>the</strong> <strong>air</strong> gap.<br />

Figure 7 Finite difference model.


5th European Thermal-Sciences Conference, The Ne<strong>the</strong>rlands, <strong>2008</strong><br />

This conclusion can be fur<strong>the</strong>r substantiated if <strong>the</strong> local value for <strong>the</strong> efficiency (η) <strong>of</strong> <strong>the</strong> energy<br />

transfer process in each element is considered. This efficiency is defined as ratio <strong>of</strong> energy<br />

transferred to that received by radiation in each element and is shown by equation 3<br />

element<br />

η<br />

element<br />

=<br />

(3)<br />

•<br />

I<br />

Q<br />

<strong>solar</strong><br />

A<br />

element<br />

The analysis shows that <strong>the</strong> derived values for <strong>the</strong> efficiency increases with distance from <strong>the</strong> base<br />

<strong>of</strong> <strong>the</strong> panel for all mass flow rates. This increase when <strong>the</strong> temperature <strong>of</strong> <strong>the</strong> panel decreases and<br />

<strong>the</strong> temperature <strong>of</strong> <strong>the</strong> <strong>air</strong> increases fur<strong>the</strong>r suggests that convection is not <strong>the</strong> dominant energy<br />

transfer process. The importance <strong>of</strong> <strong>the</strong> mixing process within <strong>the</strong> <strong>air</strong> gap was also recognised by<br />

Belusko and Saman [2006], who studied a modified version <strong>of</strong> <strong>the</strong> <strong>solar</strong> <strong>air</strong> <strong>heater</strong> described above.<br />

Weerakoon [2007] has also examined <strong>the</strong> <strong>performance</strong> <strong>of</strong> <strong>the</strong>se collectors using a finite volume<br />

CFD package (ANSYS CFX, version 5.7.1). A range <strong>of</strong> numerical models which encompassed <strong>the</strong><br />

region in front <strong>of</strong> <strong>the</strong> panel, <strong>the</strong> gas space and a panel with 1, 2, 3, 5 or 10 holes were constructed.<br />

The models examined <strong>the</strong> energy transfer mechanisms which develop with increasing <strong>air</strong> mass flow<br />

rates through <strong>the</strong> panel. During <strong>the</strong> commissioning <strong>of</strong> <strong>the</strong>se models <strong>the</strong>ir prediction were compared<br />

with experimental values and Figure 8(a) shows <strong>the</strong> comparison <strong>of</strong> experimental values <strong>of</strong> collector<br />

efficiency with extrapolated values from <strong>the</strong> numerical models. The close agreement in <strong>the</strong>se<br />

results gives confidence to <strong>the</strong> predictions from <strong>the</strong> numerical models. The information presented<br />

in Figure 8(b) shows <strong>the</strong> vectorial velocity pr<strong>of</strong>iles which develop before and after <strong>the</strong> panel which<br />

has a 20mm <strong>air</strong> gap. In this particular condition, <strong>the</strong> hole is near <strong>the</strong> base <strong>of</strong> a panel which has an<br />

equivalent mass flux <strong>of</strong> 0.0133 kg/m 2 s and <strong>the</strong> momentum <strong>of</strong> both <strong>the</strong> <strong>air</strong> moving up <strong>the</strong> gap and<br />

<strong>the</strong> <strong>air</strong> leaving <strong>the</strong> hole are low. Figure 7(b) shows that in <strong>the</strong>se conditions <strong>the</strong>re will be a rapid<br />

increase in <strong>the</strong> <strong>air</strong> gap temperature and <strong>the</strong> results given in Figure 8(b) shows that this increase is<br />

primarily due to <strong>the</strong> mixing <strong>of</strong> <strong>the</strong> two <strong>air</strong> flows. The localised nature <strong>of</strong> <strong>the</strong> <strong>air</strong> flow entering <strong>the</strong><br />

hole is also shown in Figure 8(b); very little <strong>of</strong> this <strong>air</strong> has contact with <strong>the</strong> external wall <strong>of</strong> <strong>the</strong><br />

panel.<br />

70<br />

60<br />

50<br />

40<br />

30<br />

Eperimental<br />

One Hole Model<br />

Two Hole Model<br />

20<br />

Plain Panel<br />

Three Hole Model<br />

No Wind<br />

10 Porosity 0.2%<br />

Solar radiation 555 W/m 2<br />

0<br />

0 0.005 0.01 0.015 0.02 0.025 0.03<br />

Air Mass Flow Rate (kg./m 2 s)<br />

(a) Overall <strong>performance</strong>.<br />

(b) Air flow on both sides <strong>of</strong> <strong>the</strong> panel.<br />

Figure 8 Predictions from <strong>the</strong> CFD models.


5th European Thermal-Sciences Conference, The Ne<strong>the</strong>rlands, <strong>2008</strong><br />

5 Conclusions<br />

This experimental study <strong>of</strong> a transpired <strong>solar</strong> collector has shown that mixing <strong>of</strong> <strong>the</strong> warm <strong>air</strong> in <strong>the</strong><br />

<strong>air</strong> gap behind <strong>the</strong> panel is <strong>the</strong> controlling heat transfer mechanism. This mechanism was not<br />

considered in earlier studies. The important design variables that influence <strong>the</strong> <strong>performance</strong> <strong>of</strong> <strong>the</strong><br />

collectors were shown to be <strong>the</strong> suction flow rate and <strong>the</strong> incident <strong>solar</strong> radiation. The suction flow<br />

rate is <strong>the</strong> only controllable parameter that can be used to promote ei<strong>the</strong>r greater energy transfer from<br />

<strong>the</strong> transpired collector or to obtain a larger exit <strong>air</strong> temperature rise. The interaction <strong>of</strong> <strong>the</strong>se<br />

variables showed that <strong>the</strong> collector effectively performed as a one stream heat exchanger with<br />

variable mass flow. The experimental investigation confirmed that transpired <strong>solar</strong> collectors have<br />

<strong>the</strong> potential to transfer about 60% <strong>of</strong> <strong>the</strong> incident <strong>solar</strong> radiation.<br />

References<br />

Arulanandam, S.J., Hollands, K.G. and Brundrett, E., 1999, A CFD heat transfer analysis <strong>of</strong> <strong>the</strong><br />

transpired <strong>solar</strong> collector under no-wind conditions, Solar Energy, 67, 93-100.<br />

Belusko, M., Saman, W., 2006. Jet Impingement in an Unglazed Solar Air Collector. 13 th<br />

International Heat Transfer Conference, Sydney, Australia.<br />

Carpenter, S.C., Kokko, J.P., 1991. Performance <strong>of</strong> a Solar Preheated Ventilation Air System.<br />

Proceedings 17 th Annual Conference <strong>of</strong> Solar Energy Society <strong>of</strong> Canada, Toronto, Canada.<br />

Conserval Engineering Inc., 2002. www.<strong>solar</strong>wall.com.<br />

Fleck, B.A., Meier, R.M. and Matovic, M.D., 2002, A field study <strong>of</strong> <strong>the</strong> wind effects on <strong>the</strong><br />

<strong>performance</strong> <strong>of</strong> an unglazed transpired <strong>solar</strong> collector, Solar Energy, 73, (3), 209-216.<br />

Gunnewiek, L H, Hollands, K G and Brundrett, E, 1996, “Flow Distribution in Unglazed<br />

Transpired Plate Solar Air Heaters <strong>of</strong> Large Area” Solar Energy, 58 (4-6), 227-237<br />

Gunnewiek, L H, Hollands, K G and Brundrett, E, 2002, “Effect <strong>of</strong> Wind on Flow Distribution in<br />

Unglazed Transpired Plate Collectors” Solar Energy, 72 (4), 317-325<br />

Hollick, J.C., 1996, World's largest and tallest <strong>solar</strong> re-cladding, Conserval Engineering Inc.,<br />

Renewable Energy, 9 , 703-707.<br />

Kutscher, C.F., 1994, Heat exchanger effectiveness and pressure drop for <strong>air</strong> flow through<br />

perforated plated with and without crosswind, J. <strong>of</strong> Heat Transfer, 116, 391-399.<br />

Kutscher, C.F., Christensen, C.B. and Barker, G.M., 1993, Unglazed transpired <strong>solar</strong> collectors:<br />

heat loss <strong>the</strong>ory, Journal <strong>of</strong> Solar Energy, 115, 182-188.<br />

Rhee, S.J. and Edwards, D.K.,1981, Laminar Entrance Flow in a Flat Plate Duct with Asymmetric<br />

Suction and Heating, Numerical Heat Transfer, 4, 85-100.<br />

Van Decker, G.W.E., Hollands, K.G.T. and Brunger, A.P., 2001, Heat exchange relations for<br />

unglazed transpired <strong>solar</strong> collectors with circular holes on a square or triangular pitch, Solar Energy,<br />

71, 33-45.<br />

Weerakoon, A.N., Richards, P., McClew, I. and Deans, J., 2006. Use <strong>of</strong> Perforated Ro<strong>of</strong>ing Sheets<br />

as Solar Collectors. 13 th International Heat Transfer Conference, Sydney, Australia.<br />

Weerakoon, A.N., 2007. Use <strong>of</strong> Corrugated Panel as Transpired Solar Collectors. Ph. D. Thesis,<br />

University <strong>of</strong> Auckland, 2007.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!