Heat-exchange relations for unglazed transpired solar collectors with circular holes on a square or triangular pitch

https://doi.org/10.1016/S0038-092X(01)00014-7Get rights and content

Abstract

The unglazed transpired solar collector is now an established solar air heater for heating outside air directly. Sample applications include pre-heating ventilation air and heating air for crop drying. The outside air in question is drawn straight from ambient, uniformly through the whole surface of a perforated blackened plate (the absorber plate) exposed to the sun. An important parameter in fixing the collector’s efficiency is the heat exchange effectiveness, ϵ. Once ϵ is known, finding the collector’s efficiency is straightforward. This paper presents measurements of ϵ for the case where the plate is perforated with circular holes on either a square or triangular layout, covering a range of wind speeds extending from zero to 5 m/s. These data extend the earlier measurements of Kutscher to a wider range of plate thicknesses, hole spacings (pitch), suction velocities, and to include a square layout of the holes. In the region where the two experiments overlap, agreement between the two is excellent. A new predictive model is developed that is based on breaking down the total heat transfer into contributions from each of the plate sections: the front, the hole and the back. Excellent agreement was found between modelled and measured ϵ; the new measured data were predicted with a 4.2% root mean squared error (RMSE) and Kutscher’s measured data is predicted within 6.3% RMSE.

Introduction

Unglazed, transpired solar collectors (Hollick and Peter, 1990, Kutscher et al., 1991, Kutscher et al., 1993) have been the subject of a number of recent investigations (Gunnewiek et al., 1997, Dymond and Kutscher, 1997; Arulanandam et al., 2000). They are effective devices for applications where outside air is to be heated directly, such as in heating ventilation air for buildings and crop drying. The outside air in question is drawn straight from ambient, through the whole surface of a transpired, dark-coloured plate (the absorber plate), which has been perforated with holes, typically with a porosity of about 0.5%. Tests conducted on several installations indicate that the unglazed transpired collector (UTC) works very well with annual solar collection efficiencies reaching 72% (Carpenter and Kokko, 1991). Over 70 large systems, each having collector areas between 500 and 10,000 m2, have been installed and are successfully operating for fresh-air heating in Canada, the United States, Germany, and Japan and heating process air for crop drying in countries throughout the world.

The Hottel-Whillier (H-W) equation (Duffie and Beckman, 1991), which is customarily used to model standard flat plate solar collectors, can be made to model the UTC as well. For an unglazed collector, the H-W equation takes on the form: η=FRαsFRUL(Ti−T)/G, where η is the efficiency, αs is the plate solar absorptivity, UL is the heat loss coefficient, Ti is the fluid inlet temperature, T is the ambient air temperature, G is the incident solar irradiance, and FR is the heat removal factor. For the UTC, the fluid to be heated is air taken straight from ambient, so Ti=T, and the equation reduces to η=FRαs. Current theoretical modelling is aimed at determining the heat removal factor, FR, appropriate to the UTC.

Such modelling started with the work of Kutscher et al. (1991), and the more recent modelling work has been summarized by Hollands (1998). Kutscher et al. (1991) argued that, under the range of wind and suction speeds relevant for the UTC operating with moderate temperature rises, there should be no significant convective heat loss from the UTC. Long-wave radiative heat transfer from the plate to the surroundings remains as the only important heat-loss mechanism (it should be noted, however, that convection is still important, in that it is the convective heat transfer between the plate and the air that is decisive in establishing the plate temperature and the radiant loss). Assuming an isothermal plate, Kutscher et al. (1991) developed an expression for the collector efficiency that (with an appropriate interpretation of hr) can be made to reduce toη=αs/(1+hr/(ϵρCpVs))where hr is the radiative heat loss coefficient from plate to ambient, Vs is the superficial suction velocity (rate at which air is sucked through the plate, per unit plate area), ρ and Cp are the density and specific heat of the air, respectively, and ϵ is the ‘heat exchange effectiveness’, to be defined presently. This means that the H-W equation applies to the UTC provided that one takes 1/(1+hr/(ϵρCpVs)) to be the appropriate UTC expression for FR. In practice, values of η ranging from 50 to 80% have been commonly achieved, with Vs ranging from 0.03 to 0.08 m/s.

The key item in the above equations for η and FR is the effectiveness, ϵ, defined byϵ≡T0−TTp−T,where T0 is the mean air temperature leaving the plate at the back-side, T is the ambient air temperature, and Tp is the plate temperature. (A similar equation for ϵ can be readily derived for the case where the plate temperature varies from the isothermal state assumed by Kutscher et al.; in this case, Tp is to be interpreted as the mean temperature of the outside surface of the plate. A further discussion on this point is given in Appendix A.) Effectiveness, ϵ, has to be evaluated, either from experimental data, or from analysis, but, once it has been evaluated, determining η and FR is straightforward.

The expression for ϵ will depend on how the holes are laid out on the plate, and it will also depend on the nature of the wind. Cao et al. (1993), Golneshan (1994), and Golneshan and Hollands, 1998, Golneshan and Hollands, 2000 reported numerical and experimental correlation equations for ϵ for a plate with holes in the form of an array of slits, the wind being assumed to be transverse to the slits but parallel to the plate. Kutscher (1994) presented an empirical model for relatively thin plates, with circular holes on a triangular layout, with the wind parallel to the plate (grazing wind). Using computational fluid dynamics (CFD) methods, Arulanandam et al. (2000) analysed a plate having circular holes on a square layout, obtaining a correlation equation for the corresponding ϵ, but only under no wind conditions, and only accounting for the heat transfer on the front face of the plate and in the holes (i.e., it neglected the heat transfer at the back of the plate). Gawlik (1995) investigated the effect of having a low plate conductivity and the effect of corrugating the plate.

This paper presents new experimental data for thick and thin plates with circular holes on a square or triangular layout over a range of typical suction velocities and wind speeds, including zero wind. It extends the data of Kutscher to a wider range of plate thicknesses and hole pitch, as well as to the square pitch layout. In the region where the two studies overlap, the agreement between the two sets of data is very good, generally within about 3%, which is comparable to the accuracy of the experiments. A new predictive model is proposed for the effectiveness; one that gives a breakdown of the contribution to the heat transfer on each of the parts of the plate: that is, on the outside face, the hole, and the back of the plate. The additional information contained in this breakdown should be useful for plate-design purposes.

Section snippets

EXPERIMENTAL

Nine test plates, each 0.6 by 0.6 m and with circular holes, were constructed according to the properties listed in Table 1. The pitch P shown represents the shortest distance between adjacent holes. Plates numbered 2 to 9 had the holes on a square pitch, whereas Plate 1 had the holes laid out on a triangular pitch. (In other words, if straight lines joined the centres of adjacent circles, the lines would form a set of squares for Plates 2–9 and a set of equilateral triangles for plate 1.) The

RESULTS

For the range of parameters studied, the effectiveness, ϵ, was found to range from a low of 0.32 to a high of 0.91. Plots were prepared (Van Decker, 1996) of ϵ against each of the relevant parameters: suction velocity, Vs, wind speed, Uw, hole pitch, P, plate thickness, t, hole diameter, D, and plate thermal conductivity, k. Effectiveness, ϵ, was found to decrease with increasing Vs, P, and D, and to increase with increasing Uw and with t. The strongest effects were found to be attributable to V

MODEL DEVELOPMENT

The model development presented in this paper assumes that the plate is isothermal. In reality, mainly because the local convective heat transfer coefficient is much greater in and near the holes than elsewhere, there will be some non-isothermality, the magnitude of which will depend upon the plate conductivity k (see Appendix A). Each aluminum plate was found to behave, for practical purposes, like an isothermal plate, but with the PVC plates, some modest degree of temperature non-uniformity

CONCLUSIONS

Eq. (17) provides a suitable model for the thermal effectiveness, ε, for UTCs, including the case of no-wind. With suitable interpretation of the measure of pitch, it is useful for plates laid out with square and triangular hole placements. It captures the effect of a range of variables, such as, suction velocity, Vs, wind velocity, U, hole pitch, P, hole diameter, D, and plate thickness, t, having been tested over the following ranges of these variables: 0.028 m/s≤Vs≤0.083 m/s; 0 m/s≤U≤5.0

NOMENCLATURE

    a

    constant in Eq. (10), chosen so as to give the best fit to experimental data

    Ad

    dimensionless group, known as the plate admittance, Ad=kst/kD

    Arep

    the surface area of a representative element of the plate

    c

    constant in Eq. (16), equal to 0.0047

    c

    constant in Eq. (15)

    Cp

    specific heat of air, J/kgK

    D

    diameter of hole, m

    e

    constant in Eq. (7), chosen to give the best fit to experimental data

    f

    constant in Eq. (11), chosen to give best fit to experimental data

    FR

    solar collector heat removal factor

    G

    solar

Acknowledgements

The authors would like to thank the Natural Sciences and Engineering Research Council Canada for providing the financial support for this work through a research grant to the University of Waterloo. We would also like to thank Dr. Ali Golneshan for all his excellent and painstaking work in designing, setting up, and breaking in the apparatus. Thanks also to Professor Ewart Brundrett for giving us access to the wind tunnel. The experiments were performed at the University of Waterloo.

References (20)

  • C. Dymond et al.

    Development of a flow distribution and design model for transpired solar collectors

    Solar Energy

    (1997)
  • Arulanandam S. J. (1995) A numerical investigation of unglazed transpired-plate solar collectors under zero-wind...
  • S.J. Arulanandam et al.

    A CFD heat transfer analysis of the transpired solar collector under no-wind conditions

    Solar Energy

    (2000)
  • Carpenter S. C. and Kokko J. P. (1991) Performance of solar preheated ventilation air systems, Proceedings Annual...
  • Cao S., Hollands K. G. T. and Brundrett E. (1993) Heat exchange effectiveness of unglazed transpired-plate solar...
  • Duffie J. A. and Beckman W. A. (1991) Solar Engineering of Thermal Processes, 2nd edn, pp. 278–279. John Wiley & Sons...
  • M.A. Ebadian et al.

    Forced convection, internal flow in ducts

  • Gawlik K. M., (1995) A Numerical and Experimental Investigation of Heat transfer Issues in the Practical Utilization of...
  • Golneshan A. A. (1994) Forced Convection Heat Transfer from Low Porosity Slotted Transpired Plates, Ph.D. Thesis,...
  • Golneshan A. A. and Hollands K. G. T. (1998) Experiments on forced convection heat transfer from slotted transpired...
There are more references available in the full text version of this article.

Cited by (123)

View all citing articles on Scopus
View full text