Study of Solar Radiation in View of Photovoltaic Systems Optimization

The meteorological data concerning solar radiation are generally not sufficient to allow quantifying all the phenomena which occur when a photovoltaic panel receives solar light. It is therefore, necessary to supplement these data by using astronomical calculation for the sun position and modelling the atmosphere. A simple method to calculate global, diffuse and direct irradiance on vertical and tilted surfaces for all uniform sky conditions (clear sky and overcast sky), developed in Constantine (Algeria) and Louvain-la-neuve (Belgium), has been compared with experimental data obtained at Ghardaia (Algeria). In spite of its simplicity, the method furnishes reasonably good predictions, in comparison with experimental data, and can be proposed as a simplified tool for design purposes. This method relies on the fact that we can calculate the irradiance on a plane with arbitrary orientation and inclination, based on the measurement of a single irradiance value on a reference plane.


Introduction
Solar energy is a sustainable and abundant energy resource.Estimating solar irradiation incident on inclined surfaces of various orientations is necessary to calculate the electric power generated by photovoltaic (PV) [1], to design solar systems and to evaluate their long-term average performance.Despite the fact that many meteorological/radiometric stations measure global and diffuse irradiation received on horizontal surfaces, the data on inclined surfaces are not available and are thus estimated with different models from those measured on horizontal surfaces.
Global radiation incident on a tilted plane consists of two components: direct radiation and diffuse radiation [2].Perez et al. [3,4] developed two new and more accurate but considerably simpler versions of the original Perez diffuse irradiance model [5].The original version has been used worldwide to estimate with a short time step (hourly or less) the irradiance on tilted surfaces based on global and direct (or diffuse) irradiances measured on a horizontal surface.Li and Lam [6] evaluated the anisotropic models of Klucher [7]; Hay [8] and Perez et al. models were applied to two-years of measured data in Hong Kong (1996Kong ( -1997)).All three models produced large errors for north-facing surfaces.
Predictions for south-facing surfaces showed reasonably good agreement with measured data.Kamali et al. [9] evaluated eight diffuse models to estimate solar irradiation on tilted surfaces using daily measured solar irradiation data in Iran and recommended the Reindl et al. [10] model for estimating solar irradiation on tilted surfaces.Notton et al. [11] evaluated the combination of some well-known models to estimate the hourly global solar irradiation on a tilted surface from those on a horizontal surface.They recommended the Klucher model separately or in combination with other models to estimate diffuse solar irradiation on surfaces tilted towards the equator.
All above mentioned authors worked on models requiring the experimental values of several radiation components, for example both global and diffuse values on a horizontal plane.However, in many cases, only one component of the irradiance is measured.Often, meteorological stations include global irradiation measurements on horizontal plane.On the other hand, in PV experimental plants, often global irradiance on the PV plane is measured.
Our paper shows that, from global irradiance measurements, we can quantify the various components solar radiation at the place considered.Then, it is possible to compute global irradiance or the effect of each irradiance component on a photovoltaic field of arbitrary orientation and inclination.
This paper describes a model which makes possible this decomposition.

Basic Assumptions for Splitting the Irradiance in Various Components
In practice, measurements on the ground allow for quantification of global radiation on a plant and it is almost impossible to have at a given place experimental data for all incident directions simultaneously.However, using some models, it is possible to split the measured radiation at a place in a number of components with simple directional distributions.A first subdivision is to distinguish direct sunlight and diffuse radiation.The first component reaches a given place without having undergone absorption or diffusion from the atmosphere or reflection from the ground.The second is known as diffuse radiation, i.e. the one that comes to that place after having suffered one or more changes of direction; several components are involved.
The model first separates sky diffuse from ground diffuse as the direction is pointing to the sky or toward the ground.The sky diffuse will be seen as the sum of three components: an isotropic term, a circumsolar term and a horizon circle term.The ground diffuse, less important, is normally treated as a single isotropic component.

Plane Orientation and Inclination, and Sun Position, in Horizontal Coordinates
The processing of solar data requires the definition of several references.Their specific use is justified in the following.The orientation and inclination of a plane are characterized by variables  and  shown in Figure 1, where is the vector perpendicular to the plane.n  The sun position at any moment is given by the azimuth of the sun a and its elevation h as shown in Figure 2. The incidence angle on the considered plane is then given by     cos sin sin sin cos cos sin cos cos sin sin cos cos cos cos cos sin sin

Calculation of Sun Position and Angle of Incidence from Sun Time Coordinates
We can express the azimuth a, the elevation h and the angle of incidence against the hour angle H and solar declination δ, i.e. the time (equatorial) coordinates defined as shown in Figure 3. Angle d is defined by the direction of the sun with respect to the equatorial plane of the earth.

Sun Time Coordinates
The declination  varies during the year between two extrema respectively equal to the ecliptic obliquity and its opposite.We can reference the approximate declination  (N), depending on the day of the year N, by: On the other hand, the angle H (in) is related to the local true solar time TSV (in hours) by the formula owing to the fact that a rotation of 360 of the earth take place in 24 hours.

Example
The hour angle H, the declination , the azimuth a and the elevation h as a function of true solar time, can be observed in Figure 4 for Constantine latitude and six different dates in the year.the mean distance between earth and sun, the intensity of this radiation is:

Sky Model
The variation of the earth-sun distance in one year alters the intensity of radiation actually received by the earth (outside atmosphere).It is therefore necessary to assign at the value of I 0 a correction factor C(N) which could be calculated using Kepler's laws, where N is the number of days of the year.
where e n s and e n g are the direct and global irradiation outside atmosphere.
An approximate expression is: This leads to: Extraterrestrial radiation received by the earth will vary as follows:

Attenuation of Radiation by Atmosphere
In the atmosphere, two phenomena will affect the radiation.A phenomenon of absorption and a phenomenon of diffusion will influence the radiation during its travel through the atmosphere.
In first approximation, the intensity of direct radiation decays exponentially in function of the relative atmospheric mass m where m' is the absolute atmospheric mass defined as the integral of air density along a solar ray and m' 1 the integral of standard air density in the vertical direction taken from the sea level.
Transmission of monochromatic radiation through the atmosphere obeys Beer's law, so that 0 e k m I I where I 0λ is the incident radiation intensity of wavelength λ outside atmosphere, I λ is the transmitted radiation intensity and m is the air mass in the path of radiation.

Reference Atmosphere
To characterize any atmosphere, it is compared to a known reference atmosphere which is defined as containing only gas (clean and dry atmosphere), including a given quantity of ozone, and being of uniform composition.[12] For the total spectrum, one defines a mean factor k 0 by

Linke's Turbidity Factor
For irradiance on a real atmosphere, it remains to introduce a factor T called total Linke's turbidity factor as below: For reference atmosphere, it is obviously clear that T equal 1.This formula also shows that the irradiance will tend to 0 when T will tend to

Absorption and Diffusion Turbidity Factors
One admits that the Linke's turbidity factor T can be split into an Linke's absorption turbidity factor T abs and a Linke's diffusion turbidity factor T', thus As an approximation, we assume [13] that T abs depends only of the variable gas, i.e. the ozone layer thickness e (in mm) and of the water vapour contents w (in g/cm²).These quantities can be obtained from meteorological data or computed using local mean seasonal values.
One defines an absorption coefficient a by:   where I 0 is the intensity of the light beam before it has passed through an absorbent medium and its intensity I output of it.The absorption coefficient due to ozone and water vapor will be  a1,w and is approximately Comparing this expression to the definition of the turbidity factor equation ( 14), taking into account only the component of the absorption turbidity factor, the following analogy may be observed at first order: This leads to the following relationship linking  a1,w to T abs : It is thus possible to calculate the atmospheric absorption turbidity factor.
This graph raises some questions about the accuracy of the model used here.Indeed, the absorbing turbidity factor is supposed to be a characteristic parameter of the atmosphere.This quantity should be independent of the elevation h.So the dependence on h of the following expression should exist only through the concept of relative atmospheric mass.[14][15][16]

Radiation Components
Knowing the theoretical extraterrestrial radiation, it is possible to determine the different radiation components for each value of the Linke's diffusion turbidity factor T'.

Direct Radiation
The Equation (13) gives the value of direct radiation on a normal surface.For a tilted surface, apparent intensity of direct radiation will be equal to:  d can be estimated by the following empirical formula: , exp 2.48 sin 4 where:

Diffuse Isotropic Radiation
The diffuse isotropic radiation consists of rays that reach the plane after one or several successive reflections.Some of the reflections take place on the ground so that the intensity of diffuse isotropic radiation depends on the ground reflectance, called albedo.

     
0.05 1 cos 0.5 1 cos A semi-empirical expression of  i is obtained by the difference where: , exp 1 1.06 log(sin ) with: And a semi empirical expression for  i ' is: where A regsoil the regional ground albedo and g -is the global irradiance on a horizontal plane.This component is:

Diffuse Horizon Circle Radiation
The diffuse circle horizon radiation can be estimated by the expression (for a regional albedo 0.2): with:

Ground Diffuse Radiation
An expression used to estimate this component is as follows: where A locsoil the local soil albedo.Finally, the diffuse radiation is written in the form of 4 components:

. Global Radiation
The intensity of global radiation incident at a given moment, on any surface characterized by angles (α,) is

Evolution of the Radiation against T'
The graph below shows the general evolution of the components of global radiation based on the Linke's diffusion turbidity factor.The graph below shows the general evolution of the components of global radiation based on the diffusion turbidity factor for the 12th hour throughout the year.
The following graph shows the evolution of radiation and its components for all hours of the day As we can see from this graph, the evolution of global radiation depends on time of day.Indeed, it is easy to understand that if the PV panels are not oriented towards the sun, the proportion of diffuse radiation is larger than in case they are oriented towards the sun.And as the intensity of the diffuse component from zero, passes through a maximum for a value of T' different from zero.

Determination of the Linke's Diffuse Turbidity Factor T'
If a measurement of the direct irradiance is available, the diffuse turbidity factor T' can be easily obtained by ( 14) (15) using a computed value of T abs .
If only a measurement of the global irradiance is available, T' can be obtained by iteration in order to obtain by calculation the measured value.In some cases, as shown in Figures 7 to 9, the problem has two solutions.However, the ambiguity can generally be avoided considering the historical evolution of T'.

Data
The measuring station was located on the Applied Research Unit of Renewable Energy, Ghardaia (32.23°N, 3.81°E, elevation 450 m).The radiometric Ghardaia station using a three-dimensional system (Sun-Tracker) has two parts:       A fixed part which is EKO pyranometer for the measurement of global radiation received on a horizontal plane, thermohygrometer TECNOEL for measuring temperature and humidity, and a solarimeter.A movable part that is capable of following the path of the sun from sunrise to sunset.This part has a pyrheliometer which is pointed toward the solar disk, for measuring the normal direct solar radiation integrated over all wavelengths (between 0.2 μm and 0.4 μm) received and EKO pyranometer for the measurement of diffuse irradiance on a hori-zontal plane with a spherical ball to hide the radiative flux coming directly from the solar disk.
All devices listed above are linked to acquisition Campbell Scientific Cr10x.The data logger is connected to a PC via RS232 for data storage.
The solar radiation data measured by one of the authors, during the first 8 days of June 2008, were used in this study.Each consisted of 5 minutes, direct, diffuse and global solar radiation, on horizontal and tilted surface (32.23˚ south-facing).

Conclusions
This paper has provided a model that allows, only on knowledge of global radiation on one plane, determining the different directional components of sunlight.This model is not free of defects, as shown by the first numerical results.Nevertheless, it will allow us to determine the effectiveness of radiation.This efficiency is dependent on several phenomena of reflection, refraction and absorption.Even if there are errors in the calculation of these components, this model still allows us to develop a computer code that takes into account the directional properties of the irradiance.

5. 1 .
Extraterrestrial Irradiance Outside atmosphere, global radiation has only one direct component.At a distance from the sun corresponding to (a) (b)

Figure 4 .
Figure 4. (a) Time coordinates H and d; (b) Horizontal coordinates a and h versus TSV (in hours).

Figure 5 .
Figure 5. Variation of the global radiation outside atmosphere during the year.

Figure 8 .
Figure 8. Radiation components evolution depending on linke's turbidity factor t' throughout the year.

Figure 9 .
Figure 9. Variation of radiation different components depending on linke's turbidity factor T'.

Figure 12 .
Figure 12.Experimental and calculated irradiance on Horizontal and tilted surface (8th of June 2008).