AVAILABLE SOLAR RADIATION THEORETICAL BACKGROUND

Similar documents
Chapter 2 Available Solar Radiation

Estimation of Hourly Solar Radiation on Horizontal and Inclined Surfaces in Western Himalayas

Chapter 1 Solar Radiation

4. Solar radiation on tilted surfaces

Observer-Sun Angles. ), Solar altitude angle (α s. ) and solar azimuth angle (γ s )). θ z. = 90 o α s

ME 476 Solar Energy UNIT THREE SOLAR RADIATION

XI. DIFFUSE GLOBAL CORRELATIONS: SEASONAL VARIATIONS

COMPUTER PROGRAM FOR THE ANGLES DESCRIBING THE SUN S APPARENT MOVEMENT IN THE SKY

Hottel s Clear Day Model for a typical arid city - Jeddah

Calculating equation coefficients

Modeling of the Optimum Tilt of a Solar Collector to Receive Maximum Radiation

Evaluation of Incident Solar Radiation on Inclined Plane by Empirical Models at Kuching, Sarawak, Malaysia

Modeling Solar Irradiance

SURFACE ORIENTATIONS AND ENERGY POLICY FOR SOLAR MODULE APPLICATIONS IN DHAKA, BANGLADESH

Motion of the Sun. View Comments

Sunlight and its Properties II. EE 446/646 Y. Baghzouz

Solar Radiation 1.1 THE SUN

ESTIMATION OF HOURLY GLOBAL SOLAR RADIATION INCIDENT ON INCLINED SURFACES IN IRAQ AT DIFFERENT SKY CONDITION

Radiation in the atmosphere

Basic Solar Geometry. Contents

BMM4753 RENEWABLE ENERGY RESOURCES

LECTURE 3 - SOLAR ENERGY AND SOLAR RADIATION -

Solar Radiation Measurements and Model Calculations at Inclined Surfaces

Effect of Solar Angles on Incident Energy of the Flat Collectors

CLASSICS. Handbook of Solar Radiation Data for India

1. SOLAR GEOMETRY, EXTRATERRESTRIAL IRRADIANCE & INCIDENCE ANGLES

ME 430 Fundamentals of Solar Energy Conversion for heating and Cooling Applications

Exercise 6. Solar Panel Orientation EXERCISE OBJECTIVE DISCUSSION OUTLINE. Introduction to the importance of solar panel orientation DISCUSSION

1/71 AES-LECTURE 1. Solar energy. solar radiation definitions incident solar energy

EGEE 437: HWK #2. Brownson. Yaqdaan Alkayyoomi, Mackenzie Ailes, Nicholas Minutillo, Sheel Vora. Group 17. Due on Thursday, Feb.

Orientation of Building

1/55. Solar energy. solar radiation definitions incident solar energy

Comparative study of two models to estimate solar radiation on an inclined surface

OPTIMIZATION OF GLOBAL SOLAR RADIATION OF TILT ANGLE FOR SOLAR PANELS, LOCATION: OUARGLA, ALGERIA

Sunlight and its Properties Part I. EE 446/646 Y. Baghzouz

Chapter Seven. Solar Energy

Page 1. Name:

Detlev Heinemann ENERGY METEOROLOGY LECTURE NOTES POSTGRADUATE PROGRAMME RENEWABLE ENERGY CARL VON OSSIETZKY UNIVERSITÄT OLDENBURG 2002

AVAILABILITY OF DIRECT SOLAR RADIATION IN UGANDA

Optimizing the Photovoltaic Solar Energy Capture on Sunny and Cloudy Days Using a Solar Tracking System

Basic information about developed Calculator for solar / wind hybrid power supply

Average Monthly Solar Radiations At Various Places Of North East India

Lecture 2: Global Energy Cycle

Simplified Collector Performance Model

Detlev Heinemann ENERGY METEOROLOGY LECTURE NOTES POSTGRADUATE PROGRAMME RENEWABLE ENERGY CARL VON OSSIETZKY UNIVERSITÄT OLDENBURG 2000

Hourly Solar Radiation Analysis of Buildings within 48 Facings in FuZhou, China

C) wavelength C) eastern horizon B) the angle of insolation is high B) increases, only D) thermosphere D) receive low-angle insolation

Direct Normal Radiation from Global Radiation for Indian Stations

Chapter 6. Solar Geometry. Contents

C2: NĂNG LƯỢNG ĐIỆN MẶT TRỜI

PV 2012/2013. Radiation from the Sun Atmospheric effects Insolation maps Tracking the Sun PV in urban environment

Estimation of Hourly Global Solar Radiation for Composite Climate

Solar Resource Mapping in South Africa

Which graph best shows the relationship between intensity of insolation and position on the Earth's surface? A) B) C) D)

DEVELOPMENT OF ANALYTICAL EQUATIONS FOR OPTIMUM TILT OF TWO-AXIS AND SINGLE-AXIS ROTATING SOLAR PANELS FOR CLEAR-ATMOSPHERE CONDITIONS

Study on the Optimal Tilt Angle of Solar Collector According to Different Radiation Types

Photovoltaic Systems Solar Radiation

Principles of Energy Conversion Part 11A. Solar Energy Insolation

Interaction between building services and environment.

Solar resource. Radiation from the Sun Atmospheric effects Insolation maps Tracking the Sun PV in urban environment

Mohammed Awwad Al-Dabbas. Mutah University, Assistance Prof, Mechanical Engineering department, Mutah University, Karak,Jordan

POSSIBILITIES OF SOLAR ENERGY APPLICATION IN RUSSIAN CITIES. Zhanna MINGALEVA 1, Natalia SHPAK 2 *

Chapter 5. Daylighting

Which Earth latitude receives the greatest intensity of insolation when Earth is at the position shown in the diagram? A) 0 B) 23 N C) 55 N D) 90 N

DETERMINATION OF OPTIMAL TILT ANGLE FOR MAXIMUM SOLAR INSOLATION FOR PV SYSTEMS IN ENUGU-SOUTHERN NIGERIA

CCMR Educational Programs

Chapter 2 Solar and Infrared Radiation

Solar Insolation and Earth Radiation Budget Measurements

ATMOSPHERIC ENERGY and GLOBAL TEMPERATURES. Physical Geography (Geog. 300) Prof. Hugh Howard American River College

The optimal tilt angle of a solar collector

Solar radiation in Onitsha: A correlation with average temperature

Correlation Between Sunshine Hours and Global Solar Radiation in Warri, Nigeria.

Evaluation of monthly total global and diffuse solar radiation in Ibi, Taraba state, Nigeria

Lecture 3: Global Energy Cycle

C) the seasonal changes in constellations viewed in the night sky D) The duration of insolation will increase and the temperature will increase.

EELE408 Photovoltaics Lecture 04: Apparent Motion of the Sum

LP PYRA Installation and Mounting of the Pyranometer for the Measurement of Global Radiation:

Solar radiation. Chapter Introduction

IMPROVED MODEL FOR FORECASTING GLOBAL SOLAR IRRADIANCE DURING SUNNY AND CLOUDY DAYS. Bogdan-Gabriel Burduhos, Mircea Neagoe *

HORIZONTAL AND VERTICAL ILLUMINANCE/IRRADIANCE FROM THE IDMP STATION IN GENEVA

The Atmosphere. Importance of our. 4 Layers of the Atmosphere. Introduction to atmosphere, weather, and climate. What makes up the atmosphere?

Available online at ScienceDirect. Energy Procedia 81 (2015 )

Optimum solar flat-plate collector slope: Case study for Helwan, Egypt

Experimental and Theoretical Study on the Optimal Tilt Angle of Photovoltaic Panels

Estimation of global solar radiation using clear sky radiation in Yemen

Why the Earth has seasons. Why the Earth has seasons 1/20/11

A) usually less B) dark colored and rough D) light colored with a smooth surface A) transparency of the atmosphere D) rough, black surface

The Spatial Analysis of Insolation in Iran

Chapter 2. Heating Earth's Surface & Atmosphere

ISESCO Science and Technology Vision

SEASONAL AND DAILY TEMPERATURES

Estimation of Seasonal and Annual Albedo of the Earth s Atmosphere over Kano, Nigeria

DELIVERABLE 5.3 Report on Predictive control system development and performance tests 1

Model 3024 Albedometer. User s Manual 1165 NATIONAL DRIVE SACRAMENTO, CALIFORNIA WWW. ALLWEATHERINC. COM

A FIRST INVESTIGATION OF TEMPORAL ALBEDO DEVELOPMENT OVER A MAIZE PLOT

Asian Journal on Energy and Environment

Development and Validation of Flat-Plate Collector Testing Procedures

Questions you should be able to answer after reading the material

Estimation of Diffuse Solar Radiation for Yola, Adamawa State, North- Eastern, Nigeria

ASSESSMENT OF SOLAR RADIATION DATA USED IN ANALYSES OF SOLAR ENERGY SYSTEMS

Transcription:

AVAILABLE SOLAR RADIATION THEORETICAL BACKGROUND DETERMINING THE OPTIMUM ORIENTATION OF A GREENHOUSE ON THE BASIS OF THE TOTAL SOLAR RADIATION AVAILABILITY

THE SOLAR CONSTANT At a distance of one astronomical unit, 1.495 10 11 m, the mean earth-sun distance, the sun subtends an angle of about 0.5. The radiation emitted by the sun and its spatial relationship to the earth result in a nearly fixed intensity of solar radiation outside of the earth s atmosphere. The solar constant Gsc is the energy from the sun per unit time received on a unit area of surface perpendicular to the direction of propagation of the radiation at mean earth-sun distance outside the atmosphere. Sun-earth relationships

EXTRATERRESTRIAL RADIATION The spectral distribution of the extraterrestrial radiation is the radiation that would be received in the absence of the atmosphere. Variation of the earth-sun distance does lead to variation of extraterrestrial radiation flux in the range of 3,3 %. A simple equation with accuracy adequate for most engineering calculations is given by first eqaution, while second provides a more accurate data: where G on is the extraterrestrial radiation incident on the plane normal to the radiation on the nth day of the year and B is given by: Variation of extraterrestrial solar radiation with time of year

DIRECTION OF BEAM RADIATION The geometric relationships between a plane of any particular orientation relative to the earth at any time (whether that plane is fixed or moving relative to the earth) and the incoming beam solar radiation, that is, the position of the sun relative to that plane, can be described in terms of several angles. Latitude angle of a place - the angular location north or south of the equator, north positive; 90 90 Declination angle - the angular position of the sun at solar noon (when the sun is on the local meridian) with respect to the plane of the equator, north positive The declination δ can be found from the approximate equation of Cooper (1969) = 23.45 sin(360 284+n 365 ) nth day of the year (starts from January 1) Maximum and minimum value of declination angle -23.45 o +23.45 o

The day of the year - n Variation in sun-earth distance and declination are all continuously varying functions of time of year. For many computational purposes it is customary to express the time of year in terms of n, the day of the year, and thus as an integer between 1 and 365. The use of integer values of n is adequate for most engineering calculations. Recommended Average Days for Months and Values of n by Months (Klein, 1977, do not use for > 66.5 o

Direction of beam radiation The hour angle - the angular displacement of the sun east or west of the local meridian due to rotation of the earth on its axis at 15 per hour: =0 at noon negative in the morning positive in the afternoon Surface azimuth angle - the deviation of the projection on a horizontal plane of the normal to the surface from the local meridian. In northern hemisphere: =0 o for south facing surfaces =180 o for north facing surfaces = -90 o for east facing surfaces = +90 for west facing surfaces Slope the angle between the plane of the surface and the with horizontal ( 90 o means tha the surface has a downward facing component): 0 o 180 o

Additional angles are defined that describe the position of the sun in the sky Zenith angle z - the angle between the vertical and the line to the sun the angle of incidence of beam radiation on a horizontal surface cos θz = cos φ cos δ cos ω + sin φ sin δ Solar altitude angle s - the angle between the horizontal and the line to the sun the complement of the zenith angle Solar azimuth angle s - the angular displacement from south of the projection of beam radiation on the horizontal plane Displacements east of south are negative and west of south are positive The required values of these parameters have been computed (Duffie and Beckman, 1991). Zenith angle of sun on an inclined surface of the greenhouse i is:

Extraterrestrial radiation on a horizontal surface At any point in time, the solar radiation incident on a horizontal plane outside of the atmosphere is the normal incident solar radiation given as: G o = G SC (1 + 0.033 cos 360n ) cos 365 z G sc =1367 W/m 2 is solar constant It is often necessary for calculation of daily solar radiation to have the integrated daily extraterrestrial radiation on a horizontal surface H o. Thus is obtained by integrating G o over the period from sunrise to sunset (J/m 2 ): H o G o = G SC (1 + 0.033 cos 360n ) (cos φ cos δ cos ω + sin φ sin δ) 365 24 3600 GSC 360n = + s 1 0. 033cos cos cos sin s + sin sin 365 180 It is also of interes to calculate the extraterrestrial radiation on a horizontal surface for an hour period. Thus is obtained by integrating G o for a period between hour angles 1 and 2 which define an hour (where 2 is the larger): I o = 12 x 3600 G SC (1 + 0.033 cos 360n ) [cosφ cosδ (sinω 365 2-sinω 1 ) + π( ) 2 1 sinφ sinδ) 180

Available solar radiation The primary radiation fluxes on a surface at or near the ground that are important in connection with solar thermal processes. It is convenient to consider radiation in two wavelength ranges: Solar or short-wave radiation is radiation originating from the sun, in the wavelength range of 0.3 to 3 µm. In the terminology used here, solar radiation includes both beam and diffuse components unless otherwise specified. Long-wave radiation is radiation originating from sources at temperatures near ordinary ambient temperatures and thus substantially all at wavelengths greater than 3 µm (longwave radiation is emitted by the atmosphere, by a collector, or by any other body at ordinary temperatures). The radiant energy fluxes of importance in solar thermal processes. Short-wave solar radiation is shown by Long-wave radiation is shown by

Clearness index K T Before reaching the surface of the earth, radiation from the sun is attenuated by the atmosphere and the clouds. The ratio of solar radiation at the surface of the earth to extraterrestrial radiation is called the clearness index. The clearness index gives a measure of the atmospheric effects at a place on the insolation. It is a stochastic parameter, which is a function of time of year, season, climatic condition and geographic location. A daily clearness index is the ratio of a particular day s radiation to the extraterrestrial radiation for that day: An hourly clearness index can defined as: k T = I I o K T = H H o Annual Monthly Averaged Clearness Index in Serbia

Solar radiation data Solar radiation data are available in several forms. The following information about radiation data is important in their understanding and use: whether they are instantaneous measurements (irradiance) or values integrated over some period of time (irradiation)(usually hour or day); the time or time period of the measurements; whether the measure-ments are of beam, diffuse, or total radiation; the instruments used; the receiving surface orientation (usually horizontal, sometimes inclined at a fixed slope, or normal to the beam radiation); if averaged, the period over which they are averaged (e.g., monthly averages of daily radiation). Most radiation data available are for horizontal surfaces were measured with thermopile pyranometers. Average daily solar radiation data are also available from maps that indicate general trends. In some geographical areas where climate does not change abruptly with distance (i.e., away from major influences such as mountains or large industrial cities), maps can be used as a source of average radiation if data are not available. However, large-scale maps must be used with care because they do not show local physical or climatological conditions that may greatly affect local solar energy availability.

Hourly Global Solar Radiation on Horizontal Surfaces Global solar radiation on horizontal surfaces can be measured with a pyranometer, which is an instrument that measures global solar radiation from all directions. The global solar radiation on horizontal surfaces can be categorized as follows: Direct beam solar radiation (I b ) Diffuse solar radiation (I d ) Solar radiation on a horizontal surface is the sum of the horizontal direct and diffuse radiation: I = I d + I b Instruments for measuring solar radiation are of two basic types: A pyrheliometer is an instrument using a collimated detector for measuring solar radiation from the sun and from a small portion of the sky around the sun (i.e., beam radiation) at normal incidence. A pyranometer is an instrument for measuring total hemispherical solar (beam plus diffuse) radiation, usually on a horizontal surface. If shaded from the beam radiation by a shade ring or disc, a pyranometer measures diffuse radiation. CMP3 (KIP&ZONEN) Pyranometer Automatic meteorological station in Čačak

Hourly Diffuse Radiation on Horizontal Surface Diffuse solar radiation is a part of the sunlight that passes through the atmosphere and is consumed, scattered or reflected by water vapor, dust particles or pollution. This type of radiation cannot throw sharp shadows and therefore cannot be focused. Large numbers of meteorological/radiometric stations normally measure the received global irradiation on horizontal surfaces. It is not an easy task to collect such measurements, mainly due to the high price of the measuring equipment. For this reason, a number of mathematical models have been defined to estimate diffuse radiation on horizontal surfaces. The models that can determine diffuse radiation on horizontal surfaces can be classified into two types as follows: Parametric models Decomposition models

Parametric models for diffuse radiation on horizontal surfaces Specific information of environmental conditions are necessitates of parametric models, such as: - atmospheric turbidity - fractional sunshine - cloud cover - perceptible water content I = I bn cos z + I d I bn = A exp( B cos z ) I d = CI z - the angle of incidence of beam radiation on a horizontal surface θz = cos -1 ( cos φ cos δ cos ω + sin φ sin δ) I bn is normal solar beam radiation on horizontal surfaces A, B, C are values of constants and are given for different model

Parametric models for diffuse radiation on horizontal surfaces The ASHRAE algorithm offers a simpler method, which is widely utilized by the engineering communities Machler and Iqbal s Model studied estimating hourly diffuse irradiation from hourly global solar radiation measured on horizontal surfaces. They recommended investigating solar altitude ranges above 40. They improved the invariables A, B and C (of the ASHRAE model), which considers the improvements in solar radiation research over the previous decades. Parishwad s Model assessed the constants of the ASHRAE model (A, B, and C) using statistical indicator analysis of measured solar radiation data of six cities from different locations in India. Nijegorodov s Model revised the constants of ASHRAE s model, using a computer program to predict hourly and daily global solar radiation in Botswana based on recorded solar radiation components in Botswana, Namibia, and Zimbabwe

Decompoition models for diffuse radiation on horizontal surfaces Decomposition models typically utilize only data pertaining to global radiation to estimate diffuse radiation from global solar radiation data. Decomposition models are based on a correlation between the diffuse and total radiation on a horizontal surface (k T = clearness index). This correlation is defined as a function of the hourly clearness index (ratio of hourly global horizontal I to hourly extraterrestrial radiation. Orgill and Hollands (1977) have used data of this type from Canadian stations, Erbs et al. (1982) have used data from four U.S. and one Australian station, and Reindl et al.(1990) have used an independent data set from the United States and Europe. The three correlations are essentially identical, although they were derived from three separate databases. The ratio Id /I as function of hourly clearness index k T showing the Orgill and Hollands (1977) Erbs et al. (1982) and Reindl et al. (1990) correlations

Most used models for diffuse radiation on horizontal surfaces The Erbs et al. correlation: The Orgill and Hollands correlation has been widely used, produces results that are for practical purposes the same as those of Erbs et al., and is represented by the following equations: The ratio Id/I as a function of hourly clearness index k T From Erbs et al. (1982)

RADIATION ON SLOPED SURFACES The general problem of calculation of radiation on sloped surfaces is when only the total radiation on a horizontal surface is known. For this we need the directions from which the beam and diffuse components reach the surface in question. The direction from which diffuse radiation is received is a function of conditions of cloudiness and atmospheric clarity, which are highly variable. Clear-day data suggest a diffuse radiation model as being composed of three parts: The first is an isotropic part, received uniformly from the entire sky dome. The second is circumsolar diffuse, resulting from forward scattering of solar radiation and concentrated in the part of the sky around the sun. The third, referred to as horizon brightening, is concentrated near the horizon and is most pronounced in clear skies. Schematic of the distribution of diffuse radiation over the sky dome showing the circumsolar and horizon brightening components added to the isotropic component

Solar radiation on a sloped surface ( ) The incident solar radiation is the sum of a set of radiation streams including beam radiation, the three components of diffuse radiation from the sky, and radiation reflected from the various surfaces seen by the tilted surface: I = I,b + I,d,iso + I,d,cs + I,d,hz + I,ref where the subscripts iso, cs, hz, and ref refer to the isotropic, circumsolar, horizon, and reflected radiation streams. Generally, diffuse radiation models for tilted surfaces can be classified into two groups: isotropic and anisotropic models. They differ in the division of the sky into regions with normal and elevated diffuse radiation intensities. Isotropic models assume there is uniformity in the distribution of diffuse radiation intensity over the sky. Anisotropic models include appropriate modules for representing areas of elevated diffuse radiation.

Solar radiation on a sloped surface ( ) For a surface (a collector) of area A c, the total incident radiation can be expressed in terms of the beam and diffuse radiation on the horizontal surface and the total radiation on the surfaces that reflect to the tilted surface. I A c = I b R b A c + I d,iso A s F s-c + I d,cs R b A c + + I d,hz A hz F hz-c + I i i A i F i-c 1. The first term is the beam contribution. 2. The second is the isotropic diffuse term, which includes the product of sky area A s (an undefined area) and the radiation view factor from the sky to the collector F s c. 3. The third is the circumsolar diffuse, which is treated as coming from the same direction as the beam. 4. The fourth term is the contribution of the diffuse from the horizon from a band with another undefined area A hz. 5. The fifth term is the set of reflected radiation streams from the buildings, fields, and so on, to which the tilted surface is exposed. The symbol i refers to each of the reflected streams: a. I i is the solar radiation incident on the ith surface, b. ρ i is the diffuse reflectance of that surface c. F i c is the view factor from the ith surface to the tilted surface.

Solar radiation on a sloped surface ( ) useful form In general, it is not possible to calculate the reflected energy term in detail, to account for buildings, trees, and so on, the changing solar radiation incident on them, and their changing reflectances. Standard practice is to assume that there is one surface, a horizontal, diffusely reflecting ground, large in extent, contributing to this term. In this case, I i is simply I and ρ i becomes ρ g, a composite ground reflectance. Previous equation can be rewritten in a useful form by interchanging areas and view factors (since the view factor reciprocity relation requires that, e.g., A s F s c = A c F c s ). This eliminates the undefined areas A s and A hz (The area A c appears in each term in the equation and cancels): I = I b R b + I d,iso F s-c + I d,cs R b + I d,hz F hz-c + I g F c-g The result is an equation that gives I in terms of parameters that can be determined either theoretically or empirically. Many models have been developed, of varying complexity, as the basis for calculating I. The differences are largely in the way that the diffuse terms are treated.

The geometric factor R For purposes of solar process design and performance calculations, it is often necessary to calculate the hourly radiation on a tilted surface of a collector from measurements or estimates of solar radiation on a horizontal surface. The most commonly available data are total radiation for hours or days on the horizontal surface, whereas the need is for beam and diffuse radiation on the plane of a collector. The geometric factor R b, the ratio of beam radiation on the tilted surface to that on a horizontal surface at any time, can be calculated from: For horizontal surfaces, the angle of incidence is the zenith angle of the sun z - its value must be between 0 and 90 when the sun is above the horizon: Equation relating the angle of incidence of beam radiation on a surface to the other mentioned angles is:

RADIATION ON SLOPED SURFACES: ISOTROPIC SKY It can be assumed (suggested by Hottel and Woertz) that the combination of diffuse and groundreflected radiation is isotropic. With this assumption, the sum of the diffuse from the sky and the ground-reflected radiation on the tilted surface is the same regardless of orientation. The total radiation on the tilted surface is the sum of the beam contribution calculated as I b R b and the diffuse on a horizontal surface, Id. An improvement on this model, the isotropic diffuse model, was derived by Liu and Jordan (1963). The radiation on the tilted surface was considered to include three components: beam, isotropic diffuse, and solar radiation diffusely reflected from the ground: I β = I b R b + I d A surface tilted at slope from the horizontal has a view factor to the sky F c s = 1 cos β The tilted surface has a view factor to the ground F c g = 2 Total solar radiation on the tilted surface for an hour and by the definition of R: R= I b R I b + I d I 1+cos β 2 1+cos β 2 + I g 1 cos β 2 + g 1 cos β 1+cos β The isotropic diffuse model is easy to understand, is conservative (i.e., it tends to underestimate I ), and makes calculation of radiation on tilted surfaces easy. 2 2.

RADIATION ON SLOPED SURFACES: ANISOTROPIC SKY The improved models have been developed which take into account the circumsolar diffuse and/or horizon-brightening components on a tilted surface that are shown schematically in the figure. Hay and Davies (1980) estimate the fraction of the diffuse that is circumsolar and consider it to be all from the same direction as the beam radiation; they do not treat horizon brightening. Reindl et al. (1990) add a horizon brightening term to the Hay and Davies model, as proposed by Klucher, giving a model to be referred to as the HDKR model. Perez et al. (1987,1988, 1990) treat both circumsolar diffuse and horizon brightening in some detail in a series of models. Beam, diffuse, and ground-reflected radiation on a tilted surface

RADIATION ON SLOPED SURFACES: ANISOTROPIC SKY THE HAY-AND-DAVIES model The Hay-and-Davies model is based on the assumption that all of the diffuse can be represented by two parts, the isotropic and the circumsolar. The diffuse radiation on a tilted collector is written as: I β,d = I β,d,iso + I β,d,cs = I d 1 A i 1 + cos β 2 A i is an anisotropy index which is a function of the transmittance of the atmosphere for beam radiation. A i = I b I o + A i R b The anisotropy index determines a portion of the horizontal diffuse which is to be treated as forward scattered; it is considered to be incident at the same angle as the beam radiation. Under clear conditions, the A i will be high, and most of the diffuse will be assumed to be forward scattered. When there is no beam, A i will be zero, the calculated diffuse is completely isotropic. The Hay-and-Davies method for calculating I is not much more complex than the isotropic model and leads to slightly higher estimates of radiation on the tilted surface.

RADIATION ON SLOPED SURFACES: ANISOTROPIC SKY HDKR model - Hay, Davies, Klucher, Reindl model Temps and Coulson account for horizon brightening on clear days by applying a correction factor of 1 + sin 3 ( /2) to the isotropic diffuse. Klucher modified this correction factor by a modulating factor f so that it has the form 1 + f sin 3 ( /2) to account for cloudiness. Reindle found that with increasing overcast sky, there is a decrease in the diffuse radiation intensity originating from the given region. Reindl et al. have modified the Hay-and-Davies model by the addition of a term like Klucher and give the diffuse on the tilted surface like: I β,d = I d 1 A i 1 + cos β 2 1 + f sin 3 β 2 + A i R b f = I b I When the beam and ground-reflected terms are added, the HDKR model (the Hay, Davies, Klucher, Reindl model) results. The total radiation on the tilted surface is then: I β = I b + I d A i R b + I d 1 A i 1 + cos β 2 1 + f sin 3 β 2 + Iρ g 1 cos β 2

RADIATION ON SLOPED SURFACES: ANISOTROPIC SKY PEREZ ET AL. model The Perez et al. (1990) model is based on a more detailed analysis of the three diffuse components. The diffuse on the tilted surface is given by: I β,d = I d 1 F i 1 + cos β 2 + F 1 a b + F 2 sin β where F 1 and F 2 are circumsolar and horizon brightness coefficients and a and b are terms that account for the angles of incidence of the cone of circumsolar radiation on the tilted and horizontal surfaces. The circumsolar radiation is considered to be from a point source at the sun. The brightness coefficients F 1 and F 2 are functions of three parameters that describe the sky conditions, the zenith angle z (in degrees), a clearness ε, and a brightness, where ε is a function of the hour s diffuse radiation I d and normal incidence beam radiation I b,n : ε = I d + I b,n I d + 5.535 x 10 6 θ z 3 1 + 5.535 x 10 6 θ z 3 = m I d I on where m is the air mass (m=1/cos z) and I on is the extraterrestrial normal-incidence radiation. These parameters are all calculated from data on total and diffuse radiation.

RADIATION ON SLOPED SURFACES: ISOTROPIC OR ANISOTROPIC SKY MODEL? The next question is which of these models for total radiation on the tilted surface should be used. The isotropic model is the simplest, gives the most conservative estimates of radiation on the tilted surface, and has been widely used. The HDKR model is almost as simple to use as the isotropic and produces results that are closer to measured values. For surfaces sloped toward the equator, the HDKR model is suggested. The Perez model is more complex to use and generally predicts slightly higher total radiation on the tilted surface; For surfaces with (surface azimuth angle) far from 0 in the northern hemisphere or 180 in the southern hemisphere, the Perez model is suggested.

SOLAR GREENHOUSE The greenhouses are generally tilted towards the sun in order to maximize the solar irradiance on the surfaces. The main purpose is to provide an environment conducive to plant production on a year-round basis or to extend the growing season. The environment inside a greenhouse is dependen on many factors including : - the time of year, - the amount and duration of natural sunlight, - the relative humidity, - the size and type of equipment and structure and - the type of plants growing in the house. Total solar radiation received by a greenhouse at a particular time and locations depends upon its shape as well as orientation, which ultimately determines the inside air temperature, one of the most dominant parameters affecting the plant growth. It is already established that inside air temperature of a passive greenhouse directly depends upon the ambient air temperature, the solar radiation intensity, the overall heat transfer coefficient, the cover material and the wind velocity.

Greenhouse shapes Various researchers have used different greenhouse shapes at different latitudes for raising off-season vegetables and ornamental plants around the world. The most commonly used single span shapes of greenhouses for the agriculture applications are even-span, uneven-span, vinery, modified arch and quonset type. These greenhouses are either constructed along east-west (E-W) or north-south (N-S) orientation (longer axis). A greenhouse receives: most of the beam radiation at its floor (which is responsible for the increase in inside air temperature) receives diffuse and ground reflected radiation from each wall and roof (thus shape and orientation of the greenhouse has some bearing on the greenhouse air temperature). The selection of optimum shape and orientation of a greenhouse can lower the heating and cooling loads of the installed systems thereby saving a lot of operating cost.

Analytical approach of the model This study presents an approach to determining the optimum orientation of a greenhouse for year round applications for different climatic conditions. The most commonly used uneven-span single shape of greenhouse in E-W and N-S orientation have been selected for comparison. Some previous research results shows that uneven-span shape greenhouse receives the maximum solar radiation during each month of the year at all latitudes, and air temperature remains the highest inside an uneven-span shape as compared to other shapes during different months of the year. Selected greenhouse shape in E-W orientation

Total solar radiation availability on greenhouse cover Total solar radiation input to each wall, included surfaces and roofs, is computed for both solar greenhouse orientation and compared for each month of the year at different latitudes at the northern hemisphere. Greenhouse is divided into various sections along E-W and N-S orientation - the details of each section are shown in Table. Experimental validation is carried out for measured global solar radiation data for horizontal surface for location (44 o N and 20 o E), Belgarde Serbia. L=12 m, W=4 m, H=3 m Table 1. Sectional details of selected geometry greenhouses E - W orientation ( o C) A (m 2 ) SW 90 24 SR 18.5 37.9 NW 90 24 NR 135 16.8 EW 90 10 WW 90 10 N - S orientation ( o C) A (m 2 ) EW 90 24 ER 18.5 37.9 WW 90 24 WR 45 16.8 NW 90 10 SW 90 10 Section details: SW- south wall, SR - south roof, NW - north wall, NR - north roof, EW - east wall, WW - west wall, ER - east roof, WR - west roof; Sectional details of selected geometry greenhouses

Total solar radiation availability on greenhouse cover The hourly total solar radiation incident on different inclined and vertical surface of the greenhouse is total sum of solar radiation falling on different surfaces (each wall and roof) of the greenhouse. Adopting the isotropic diffuse model, the solar radiation incident on greenhouse surfaces can be broken down into three main components: 1. direct (beam) radiation emanating from the sky 2. diffuse sky radiation 3. radiation scattered or reflected by the ground. Total solar radiation of the greenhouse has been computed for all the sections of each shape of the greenhouse in both E-W and N-S orientations: t i i i = 1 * A i and I i are the area of ith section and total solar radiation available on ith section For E-W orientation the hourly total solar radiation is given by: S = 6 A I Solar irradiance components S = A I + A I + A I + A I + A I + t SW SW SR SR NW NW NR NR EW EW A WW I WW

Total radiation [kwm - 2] Model validation Model validation is carried out for the measured data of solar radiation at the specific location (Belgrade, Serbia). Monthly average global solar radiation data measured on the horizontal surface is compared with that theoretically computed on the same surface. It is observed that both the set of values are closely matched indicating that the developed model is validated well. 8 7 Predicted Measured 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 10 11 12 Month Predicted and measured monthly average daily global solar radiation on horizontal surface

Solar radiation (W/m 2 ) Irradiation [Whm -2 ] Total solar radiation availability on greenhouse cover The hourly total solar radiation on greenhouse cover surfaces is calculated for a typical summer and winter day. The hourly variation of total solar radiation on different walls and roofs in E-W orientation greenhouse at 44 o N latitude, for a typical summer day (17-07-05) is shown in figure. 1200 NR SW SR EW WW NW 1050 900 750 600 450 300 150 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 Time of the day

Irradiation [Wh/m -2 ] Solar radiation (W/m 2 ) Total solar radiation availability on greenhouse cover Hourly variation of the total solar radiation on different walls and roofs in E-W orientation greenhouse for a typical winter day (17-01-05) at 44 o N is shown in figure. 1000 900 800 700 600 500 400 300 200 100 0 NR SW SR EW WW NW 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 Time of the day

A comparison of total solar radiation availability on the selected typical day of each month of the year in E-W and N-S orientations 700 600 E-W N-S 24 o N Total solar radiation availability, kw 500 400 300 200 100 24 o N latitude 0 J17 F16 M16 A15 M15 J11 J17 A16 S15 O15 N14 D10 Month N-S orientation receives more solar radiation during all months of the year as compared to E-W orientation. In summer eastern and western section in N-S orientation receive more radiation as compared to northern and southern section in E-W orientation. The difference is small in winter months than in summer months. It can be concluded that due to higher ambient air temperatures throughout the year an orientation that receives less solar radiation should be preferred as this would help in maintaining lower inside air temperature i.e. E-W orientation is more suited.

Total solar radiation availability, kw A comparison of total solar radiation availability on the selected typical day of each month of the year in E-W and N-S orientations 700 E-W N-S 600 34 o N 500 400 300 200 100 34 o N latitude 0 J17 F16 M16 A15 M15 J11 J17 A16 S15 O15 N14 D10 Month N-S orientation of greenhouse receives less solar radiation in winter months (7.7% less in January) but more in summer months (24% more in July) as compared to E-W orientation. In winter the eastern and western side sections in N-S orientation receive less radiation as compared to northern and southern side sections in E-W orientation. However in summer months, eastern and western side sections in N-S orientation receive more solar radiation as compared to the eastern and western side sections in N-S orientation. It can be concluded that E-W orientation should be preferred, as it would provide more radiation in winter and less in summer, which is desirable in composite climates.

Total solar radiation availability, kw A comparison of total solar radiation availability on the selected typical day of each month of the year in E-W and N-S orientations Total solar radiation availability, kw 600 500 400 E-W N-S 44 o N 600 500 400 E-W N-S 54 o N 300 300 200 200 100 100 0 J17 F16 M16 A15 M15 J11 J17 A16 S15 O15 N14 D10 Month 0 J17 F16 M16 A15 M15 J11 J17 A16 S15 O15 N14 D10 Month With increase in the latitude angle, the difference in the radiation received during winter months further increases. N-S orientation of greenhouse receives much less solar radiation in winter months (less in January 14.2% at 44 o N, and 14.6% at 54 o N) but more in summer months (more in July 11.8% at 44 o N, and 9.9% at 54 o N) as compared to E-W orientation. In can be concluded that N-S orientation receives more radiation in summer, but significantly less in winter months. As winters are severe and longer here, E-W orientation should be preferred due to greater solar radiation availability in winter months as this would help in maintaining higher inside air temperatures.

Conclusions The developed model for the solar radiation availability is validated as the predicted and measured values are in close agreement. East-west orientation is best suited for year round greenhouse applications at all latitudes as this orientation receives greater radiation in winter but less in summer except near the equator. This effect is more significant at higher latitudes.