Ap1_02.vp 1 Introduction Design and performance evaluation of photovoltaic sys- tems usually involve an estimation of irradiation incident on the photovoltaic module plane. For design purposes only gloal horizontal irradiation, published by meteorological in- stitutes, is usually available. Using irradiation models, we can calculate the gain in irradiation on a tilted plane with respect to horizontal irradiation. Concerning a part of this task, most textbooks [4] recommend that the diffuse compo- nent be treated as if it were isotropically emanating from the sky vault. However, theoretical as well as experimental results have shown that this simplifying assumption is gener- ally far from reality. Dave, J., 1977 examined the validity of this isotropic distribution approximation for a sun-facing flat surface located at the bottom of plane parallel models of non-absorbing homogeneous atmospheres. Even for such idealized and somewhat unrealistic models, he showed that the use of an isotropic distribution approximation results in systematic underestimation of the diffuse energy contribu- tion to the sun-facing surfaces. That study demonstrated the need to test this approximation for more realistic, atmo- spheric conditions. Thus it appears that sky radiance should be treated as anisotropic, particularly because of the strong forward scattering effect of aerosols [2, 9, 14, 16]. In the first part of this study, some results concerning the accuracy of models to estimate irradiance on inclined planes are tested. The models chosen for discussion here are for arbitrary sky conditions and are supposed to be applicable anywhere in the world. These methods begin with measured hourly global and diffuse radiation received on a horizon- tal surface. These quantities are then transposed onto an inclined plane by a mathematical procedure. The accuracies of the models are then compared on the basis of statistical error tests, and the most accurate model is recommended. The analyses are also extended to include the optical pro- cesses occurring inside the encapsulation solar module. Solar concentrators are frequently used to increase the specific power and decrease the space of radiation degrada- tion of solar cell parameters. On the other hand, at high light concentration coefficients, a large amount of thermal energy evolves in solar cells. If cooling is not provided, the cell working temperature increases. This affects the carrier concentration and light absorption process, causing a reduc- tion in open circuit voltage, a slight increase in short circuit current, a reduction of the fill factor, and degradation of the cell power output. The degradation rate of cell output is approximately double for each 10 °C increase in temperature [7, 22]. Therefore, the modules and solar array must take full advantage of radiative, conductive and convective cooling and absorb the minimum of unused radiation. Hence, there is a need for an exact technique to calculate accurately and efficiently the temperature distribution of photovoltaic solar modules, from which we can adjust safe and proper operation at maximum ratings. The computer code was written in Fortran, a programming language quite flexible and suitable for organizing the screen with multiple windows and for drawing curves and histograms. The program structure con- sists of a main executable file, which has several external units to perform the different calculations; program maintenan- ce is therefore quite easy, as is the implementation of new models or extension to a new configuration. The derived equations are solved using an explicit technique and iteration processes [18, 20]. The models and simulation programs developed here allow us to predict the thermal performance, and they are promising tools for evaluating new photovoltaic power plants with the aim of increasing efficiency. Three main criteria are considered in evaluating the models. The first is versatility: What range of photovoltaic modules and systems can the model handle? Accuracy is the second issue: What is the physical basis and how closely can the model replicate manufacturers’ data? The final issue is computational speed. This criterion is of secondary importance for this work. 2 Apparatus and measurements Total radiation was monitored with a high precision pyranometer, which is sensitive in the wavelength range from 300 to 3000 nm. Sky diffuse radiation was measured by © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 7 Acta Polytechnica Vol. 41 No. 2/2001 Actual Optical and Thermal Performance of Photovoltaic Modules Hamdy K. Elminir, V. Benda, I. Kudláček Field testing is costly, time-consuming and depends heavily on prevailing weather conditions. Adequate security and weather protection must also be provided at the test site. Delays can be caused due to bad weather and system failures. To overcome these problems, a photovoltaic array simulation may be used. In any simulation scheme involving photovoltaic systems, one important choice is the selection of a mathematical model. In the literature several approaches to the problem have been made. Most procedures designed for this purpose are based on analytical descriptions of the physical mechanisms inside the solar cell that can be represented by a circuit diagram with discrete components, like a two-exponential model. Such simulators have some merits. However, their limited flexibility in readily simulating the influence of solar radiation, temperature and various array parameters is a serious drawback that has been noted. To get more accurate results in predicting the actual performance of photovoltaic modules, the parameters influencing incoming (optical parameters) and outgoing power flow (electrical and thermal parameters) were investigated by simulation and by some verifying experiments, to get a closer insight into the response behavior of this element, and to estimate the overall performance as well as optimization of the parameters. Keywords: solar radiation, reflection losses, radiation shape factor, radiative surface area, temperature distribution, emissivity, isotropic model, Hay and Klucher’s anisotropic models. a pyranometer equipped with a special shading device to exclude direct radiation from the sun. All data sets were subjected to various quality control tests. Three types of data checks were performed to identify missing data, data that clearly violates physical limits, and extreme data. Hours when the data was known to be ”bad” or “missing” were omitted. Second, any hour with an observation that violated a physical limit or conservation principle was eliminated from the data set, including: reported hours with a diffuse fraction greater than 1, or beam radiation exceeding the extraterrestrial beam radiation. To eliminate the uncertainty associated with radiation measurements at large incidence angles, hours with a zenith angle larger than 80� were eliminated. The final data set was constructed from the measured data that passed all of the quality control checks. 3 Optical performance of photovoltaic modules In order to make a precise representation of the actual op- tical conditions in the module, a model for an encapsulating cell was developed which determines the insolation reaching the cell from sun and sky irradiance. This was done by model- ing the optical processes occurring outside and inside the encapsulation. Most collecting devices associated with solar energy systems are tilted at some angle with respect to the horizontal. Due to the lack of measured tilted surface solar ra- diation data, models are employed to estimate the radiation incident on a collector surface from measured horizontal radiation. In the present work, methods for calculating radia- tion on a tilted surface from horizontal data are evaluated. 3.1 Evaluation of irradiation on a tilted south-oriented surface Consider a plane inclined at an angle � from the horizon- tal position. For the time being, we assume that the surface is inclined in such a way that it faces the sun; that is, direct radiation is not striking the back of this plane. Total radiation arriving on an inclined surface is composed of three com- ponents: beam radiation, sky diffuse radiation and radiation reflected from the ground. The models discussed here all share the same formulations for the beam and ground reflected components. Hourly irradiations with a solar eleva- tion angle � > 8° are tested only in order to avoid errors resulting from dividing small numbers. It is useful to write down the hourly formulation of these components. 3.1.1 Beam radiation incident on an inclined surface For a surface oriented in any direction with respect to the meridian, the trigonometric relation for the incidence angle, �I can be written in the following form: � �cos sin cos cos sin cos sin cos cos sin sin � � � � � � � � � � � I � � � � �� �cos cos cos cos sin sin sin � � � � � � � � � (1) where the meanings of the different symbols are as given in the nomenclature list. There are several commonly occurring cases for which Equation 1 is simplified. For a fixed surface sloped toward the south or north, that is, with a surface azimuth angle � of 0° or 180° (a very common situation for fixed flat plate collectors), the last term drops out. � � � � cos sin cos cos sin sin cos cos sin sin cos � � � � � � � � � � � I � � � � � � � � � cos sin sin cos cos cos . � � � � � � � � � � � (2) For horizontal surfaces, the angle of incidence, �I is the ze- nith angle of the sun, �z. Its value must be between 0° and 90° when the sun is above the horizon. For this situation � 0, and Equation 1 becomes: cos cos cos cos sin sin� � � � � ��Z � � (3) The geometric factor Rb, the ratio of beam radiation on the tilted surface, Ib, T to that on a horizontal surface at any time, can be calculated exactly by appropriate use of Equation 1. The ratio I Ib, T b, H is given by: � � � � Rb I Z � � � � �cos cos sin sin cos cos cos cos cos cos � � � � � � � � � � � � � �� sin sin . (4) The hourly beam radiation received on an inclined sur- face can be expressed as: � �I I R G I Rb, T b, H b H d, H b� � � � � . (5) It should be noted that at the grazing angles (just at sun- rise or at sunset) Rb can change rapidly and may approach infinity or zero because both the numerator and the denomi- nator are small numbers. This depends on slope, latitude and date. 3.1.2 Sky diffuse radiation incident on an inclined surface Diffuse irradiance is difficult to determine accurately with the simple parameterization methods that were used to calcu- late direct normal irradiance in the previous section, since its spatial distribution is generally unknown and time depend- ent. Three diffuse subcomponents are used to approximate the anisotropic behavior of diffuse radiation. The first is an isotropic part received uniformly from the entire sky dome. The second is circumsolar diffuse resulting from forward scat- tering 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. Several models have been pro- posed to estimate the diffuse radiation on a tilted surface (not all of which account for these three diffuse subcomponents). The Isotropic model [12] is the simplest of the tilted surface models. In this model the intensity of sky diffuse radiation is assumed uniform over the sky dome, i.e., it is independent of the azimuth and zenith angles. It approxi- mates the completely overcast sky condition. The formula for the hourly sky diffuse radiation incident on an inclined plane is given by the product of the hourly diffuse radiation incident on a horizontal surface and the configuration factor from the surface to the sky, � �1 2� cos � . I Id, T d, H� �1 2 cos � . (6) Under completely cloudy skies, the isotropic model be- comes a good approximation. As skies become clearer, the validity of the isotropic sky model deteriorates due to the presence of circumsolar and horizon brightening anisotropic effects. Klucher, T., 1979 model: This model is based on a study of clear sky conditions by Temps, R. and Coulson, K., 1977. Their model was modified by Klucher, who incorporated 8 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 41 No. 2/2001 conditions of cloudy skies. It takes into account the increase in the diffuse radiation in the vicinity of the sun (circumsolar radiation) and that near the horizon (horizon brightening). Klucher’s formulation of the hourly sky diffuse radiation inci- dent on an inclined surface is: I I F Fd, T d, H� �� � � � � � � � � � � 1 2 1 2 11 3 1 cos sin co � � � �s sin .2 3� �Z Z� � � � � � (7 ) Where F1 is the modulating function given by: F I G 1 2 1� � � � � d, H H . (8) When the skies are completely overcast, F1 0, Klucher’s model reverts to the isotropic model. Hay, J. and Davies, J., 1980 developed a model to pre- dict tilted surface diffuse radiation, which accounts for both circumsolar and isotropic diffuse radiation. Realizing that the anisotropic behavior of circumsolar diffuse radiation becomes more pronounced under clear sky conditions, Hay, J. and Davies, J., 1980 defined an ”anisotropic index, AI” to weight the circumsolar, IT, cir and isotropic, IT, iso radiation compo- nents. The anisotropy index defines a portion of the diffuse radiation to be treated as circumsolar with the remaining portion to be considered isotropic. The circumsolar diffuse radiation is projected onto the tilted surface in the same fashion as beam radiation. A I II b, H o� (9) I I A RT, cir d, H I b� � � . (10) The remaining diffuse radiation is treated as isotropic diffuse. � �I I AT, iso d, H I� � � 1 1 2 cos . � (11) The total diffuse radiation on a tilted surface is the sum of (10) and (11). � �I I A A Rd, T d, H I I b� � �� � � � � � �� � �� 1 1 2 cos . � (12) Under clear skies, the anisotropy index will be high and the circumsolar diffuse is weighted more heavily than the isotropic diffuse. Under cloudy skies, the anisotropy index goes to zero and all diffuse is treated as isotropic. The Hay model does not account for horizon brightening radiation incident on south facing surfaces. 3.1.3 Ground reflected radiation incident on an inclined surface In order to correctly evaluate the diffuse radiation inci- dent on an inclined surface of any orientation and of any tilt angle, we need to known separately the sky diffuse radiation as well as the ground reflected radiation. In this section we are going to focus on ground reflected radiation, which is also sig- nificant and which can sometimes reach values of the order of 100 W/m2 for a vertical plane. It has already been previously found, but based on restricted data, that using an albedo value measured on site and considered as constant leads to satisfactory results, Ineichen et al., 1987. � �I site Gg, H H� �� . (13) A common method for calculating the ground reflected radiation incident on a tilted surface is to assume that the foreground in the collector field of view is a diffuse reflector and that the horizon is unobstructed. Other authors have proposed anisotropic ground reflectance models [6, 8], but lack of experimental data has hampered their validation. Therefore, the ground reflected radiation is assumed to be diffuse. Under the isotropic condition, the ground reflected radiation incident on the inclined surface is given by the above quantity multiplied by the configuration factor from the ground to the inclined surface. Thus: � �I Gg, T H� � 1 2 1� �cos . (14) 3.1.4 Model evaluation Three irradiance models are considered in this work: the isotropic model, Hay’s and Klucher’s anisotropic models. For each model, measured diffuse and horizontal global values were used to calculate the irradiance on surfaces tilted at 30°, 60° and 90° above the horizon. The results were compared with the irradiances monitored and presented in terms of usual statistics: mean bias error (MBE) and root mean square error (RMSE). The results for south facing surfaces are given in Table 1 and Table 2. It is noted that the RMSEs for all three models increase as the slope of the collector increases, but remain in a domain of errors for which these relations can Acta Polytechnica Vol. 41 No. 2/2001 9 Month Isotropic Model RMSE Hay’s Model RMSE Klucher’s Model RMSE 30 ° 60 ° 90 ° 30 ° 60 ° 90 ° 30 ° 60 ° 90 ° January 18.0 22.0 22 8.0 11.0 10.0 8.0 11.0 13 February 16.0 21.0 22 7.5 11.0 12.0 7.0 11.0 13 March 9.0 12.0 12 4.0 6.0 8.5 3.0 5.0 8 April 6.0 14.0 8 4.0 12.0 8.0 3.0 12.0 8 May 3.5 5.5 12 3.0 5.5 12.0 3.0 5.0 8 June 4.0 9.0 25 3.0 9.0 21.0 4.5 8.0 24 July 3.5 6.0 19 2.0 5.0 13.0 4.0 6.5 18 August 3.0 7.0 12 2.0 5.0 10.0 3.0 5.0 12 September 11.0 15.0 17 6.5 9.0 12.0 4.5 7.0 10 October 13.0 16.0 17 6.5 7.0 9.0 5.0 6.5 9 November 24.0 31.0 36 13.0 8.0 20.0 10.0 14.0 19 December 18.0 23.0 26 9.0 13.0 14.0 9.0 11.0 13 Table 1: Root mean square errors [%] for global radiation received on inclined surfaces be applied with good accuracy. Inspecting the results, it is apparent that the models agree quite well with each other during the summer months. They deviate from each other in the winter months, when the effect of the difference in the diffuse radiation parameterization is at its maximum. The RMSE results indicate that the anisotropic models (Hay and Klucher) show similar performance on an overall basis (within 5 %) but the isotropic model exhibits much larger error. The MBE results show that the isotropic model substantially underpredicts the irradiance incident on an inclined surface, and the Klucher model considerably overpredicts irradiance incident on an inclined surface on an overall basis. Our findings confirm the observation that the Klucher model describes the irradiance on inclined planes more accurately than the isotropic model or Hay’s model. Due to lack of measured tilted surface solar radiation data, Klucher’s model is employed to estimate the radiation incident on a collector surface from measured horizontal ra- diation. The results of these calculations are plotted against the angle of tilt for summer, winter and all year round in- tended use. A summary of the effects of tilt angle on solar irradiation for south facing surfaces is shown in Fig. 1. 3.2 Optical parameters of the encapsulation module Reflection losses produced in photovoltaic modules oper- ating in real conditions in comparison with standard ones are commonly neglected or simplified to a value that does not take into account the module tilt angle. In some applica- tions of interest such as building integration of photovoltaic modules, architectural criteria or other limitations frequently force modules to have tilt angles far from their optimal values. As a consequence, optical losses increase to non-negligible quantities [17]. The different encapsulation layers cause multiple reflec- tions inside and among the slabs (see scheme in Fig. 2). The Acta Polytechnica Vol. 41 No. 2/2001 10 Month Isotropic Model MBE Hay’s Model MBE Klucher’s Model MBE 30 ° 60 ° 90 ° 30 ° 60 ° 90 ° 30 ° 60 ° 90 ° January 8.3 10.3 9.2 2.2 1.6 1.2 3.3 3.9 4.2 February 8.2 11.0 11.7 2.8 3.3 3.3 2.8 4.7 5.5 March 5.5 6.1 4.4 1.7 0.6 1.0 0.2 0.1 1.7 April 3.3 5.3 1.8 1.1 2.5 0.6 0.2 0.2 4.5 May 1.7 2.2 2.4 1.2 2.2 5.5 1.6 2.2 4.8 June 0.6 1.1 1.2 0.5 1.9 4.4 2.2 3.0 10.0 July 0.8 0.6 1.1 0.4 1.1 3.9 2.2 4.7 8.9 August 2.5 2.8 2.2 1.1 0.8 2.2 1.9 2.2 5.5 September 4.4 4.4 3.6 2.2 1.1 0.6 0.1 1.1 3.2 October 5.0 8.9 9.4 3.3 2.8 2.8 2.2 2.8 1.9 November 8.3 9.4 10.0 4.4 2.6 2.2 3.3 2.9 0.8 December 6.7 7.8 7.2 1.1 0.2 0.0 2.1 0.3 0.5 Table 2: Mean bias errors [%] for global radiation received on inclined surfaces Fig. 1: Annual variation in estimating average daily radiation for a flat plate collector facing south for a latitude of 50°, and ground reflectance of 0.2 Fig. 2: Ray tracing through the module layers abbreviation EVA stands for ethylene-vinyl-acetate, which serves as a mechanical, thermal and optical interface between the solar cells and the glass cover. ARC stands for the anti- reflective coating of the solar cell. It consists of TiO2 or another optical material to bridge over the refractive indices of silicon and the front cover. Refractive indices and trans- missivity coefficients are given in Table 3. Table 3: Refractive indices and transmissivity coefficient for the wavelength of 800 nm [5] Materials Refractive index Thickness Transmissivity coefficient Normal incidence Incidence at 60° Glass EVA ARC Silicon 1.526 1.450 1.700 3.690 3.2 mm 0.6 mm 50–60 nm 0.5 mm 0.951 0.992 1.000 – 0.941 0.990 1.000 – Light rays are modeled as objects having the properties of electromagnetic waves [1, 15]. Each ray can exist either as incident or as one of its reflected, refracted, absorbed and transmitted components, produced by the respective optical process (i.e., reflection, refraction, etc.). At the interfaces, a matrix method derived from electromagnetic theory is used to determine the probability of a photon being reflected, transmitted or absorbed in the interface layers. The calcula- tion technique is based on the model of [19, 23], which calculates the radiation transfer properties of multi-layer structures. Their model was modified to generate an equiva- lent application in the domain of photovoltaic solar cell simulation. 3.2.1 Radiation transmission through a photovoltaic module When light is incident upon the interface between two op- tically different materials, in general part of the light is reflected and part transmitted. If the medium on what is nominally chosen to be the outside (o) of the interface has a complex refractive index no and the medium on the inside (i) has a refractive index, ni then, for a given wavelength, Snell’s law states that: � �� �� �o i o i� arcsin sinn n . (15) For an incident wave of amplitude Eo � traveling across an interface from outside to inside (the superscript denotes the direction of propagation), the amplitude reflection coefficient (r) and amplitude transmission coefficient (t) are defined as: r E E t E Eoi o o + oi i + o +� �� ; . (16) Similar expressions hold for waves traveling in the opposite direction, for which it can be shown that: r r t t r rio oi oi io oi io� � � �; 1 . (17) An alternative approach is to calculate the transmission through the multi-layer using transfer matrices. Each matrix, in an ordered sequence of matrix multiplications, represents one of two possible transformations of the radiation fields. The transformations are due to: • Transmission across an interface (interface matrix), • Transmission through a layer (layer matrix). The interface matrix: The definition of the interface (M) follows from consideration of the boundary conditions at an interface. Let Eo + and Ei + represent the amplitudes of the net electric field propagation towards, and away from, an inter- face in the direction from the outside to the inside medium (Fig. 3). Eo � and Ei � are the corresponding net electric fields traveling in the reverse direction. Ei + is then composed of the transmitted component of Eo + and the reflected component of Ei �, Eo � consists of the transmitted component of Ei � and the reflected component of Eo +. Equations (16) and (17) enable this to be expressed mathematically as: E t E r E E t E r E i + oi o + io i o io i o i o + � � � � � � � (18) or E E t r t r t t E E o + o oi oi oi oi oi oi i + i � � � � � � � � � � � � � � � � � � � � 1 1 � �� � � � � � � � � � � � � ��M E Eoi E i + i . (19) When transmission across an interface is considered in terms of intensities, the following matrix equation is ob- tained: I I I I o + o oi oi oi oi oi oi i + i � � � � � � � � � � � � � � � � � � � � � � 1 1 � � � � � �� � � � � � � � � � � � � ��M I Ioi I i + i . (20) The layer matrix: The layer matrix (N) describes the attenuation of a wave as it traverses the medium within a layer. When a wave propagates in the positive direction between the interfaces of a thin layer, its initial amplitude � �Ef � 1 within the layer is attenuated to, say, � �Ef � 2 on reaching the second in- terface (Fig. 4). A wave of initial amplitude � �Ef � 2 traveling in the opposite direction through the layer is reduced to � �Ef � 1 . If the layer has thickness d then from the definition of an atten- uation factor, � � � � � � � � � � � � E d E E d E f f f f f f � � � � � � � 2 1 2 1 � � (21) or � � � � � � � � � � � � � �E E d d E E Nf f f f f f f E � � � � � � � � � � � � � � � � � � � � � � � � � � 1 1 1 0 0 2 2 � � � � � � E E f f � � � � � � � � 2 2 . (22) Matrix multiplication: The following equation repre- sents the transformation of the electric fields by a thin layer Acta Polytechnica Vol. 41 No. 2/2001 11 Fig. 3: The electric fields at an interface ( f ) sandwiched between a pair of outer (o) and inner (i) media: � � � � � �E E M N M E E o + o of E f E fi E i + i � � � � � � � � � � � � � � � � � � . (23) This involves the multiplication of three matrices, the first and third being interface materials, the second the layer ma- trix. A system is composed of a greater number of thin layers, say layers 1, 2, 3, …, q sandwiched in the transformation. � � � � � � � � � � � � E E M N M N N M N o + o o1 E 1 E 12 E 2 E q 1 E q 1, q E q � � � � � � � � � � � � � � � � � � � � � � � E qi E i + i � � � � � � ��M E E . (24) Equation (24) can be use to obtain the total amplitude reflection and transmission coefficients, �oi and �oi of the multiple thin layer structure. This follows by considering the surrounding media to be of infinite thickness, so that there is no net reflected wave in medium (i) incident upon the inner interface (i.e., Ei � � 0). By setting Eo � � 1 Equation (24) becomes: � � � � � � � � � � � � � � 1 � � � � � � � � � � � � � �� � � �� � oi o1 E 1 E 12 E 2 E q 1 E q 1, q E q E M N M N N M N � �Mqi E oi� 0 � � � � � � (25) which can be solved for �oi and �oi . 3.2.2 Absorbed solar radiation The prediction of collector performance requires infor- mation on the solar energy absorbed by the collector absorber plate. The solar energy incident on a tilted collector can be found by the methods of Section 3.1. This incident radiation has three different spatial distributions: beam radiation, dif- fuse radiation, and ground reflected radiation, and each must be treated separately. Using Klucher’s formulation, Equation (7) can be modified to give the absorbed radiation, S, by multiplying each term by the appropriate transmittance-ab- sorptance product: � � � � � �S G I R I F F � � � � � � �� � � � �� �� � � b H d, H b d d, H 1 2 1 2 11 3cos sin � � � � � � 1 2 3 1 2 1 cos sin cos . � � �� � � Z g H � �� � �� � � �G (26) The subscripts b, d, and g represent beam, diffuse, and ground, respectively. Equation (27) can be used to find the proper absorptance: � � � � � � � � � �� � � 1 a a transmitted incident o I I K L exp cos (27) where K is the proportionality constant, the extinction coeffi- cient, which is assumed to be a constant in the solar spectrum. For glass, the value of K varies from approximately 4 m 1for ”water white” glass (which appears white when viewed on the edge) to approximately 32 m 1 for poor (greenish cast of edge) glass. L is the path length in the medium. 3.2.3 Efficiency limits of photovoltaic modules The analytical results above were applied to analyze the optical behavior of an encapsulated photovoltaic solar cell. Photovoltaic modules have an angle of incidence AOI de- pendent optical behavior that can be measured and used to improve the analysis of array performance. Like absolute air mass, solar angle of incidence is time of day dependent. Its effect on the short circuit current of a photovoltaic module results from two causes. The first is familiar to solar enthusiasts as the ”cosine effect”. The ”cosine effect” is independent of the module design, and is only geometry related. Due to the cosine effect, the short circuit current from a module varies directly with the cosine of the AOI. For example, at AOI 60° the cosine effect reduces the short circuit current by one half compared to the normal incidence condition. The second way in which the short circuit current is af- fected by the AOI is dependent on the module design. The optical characteristics of the module materials located be- tween the sun and the solar cells cause the effect. Above 45° incidence angle the reflectivity starts to increase, and reaches 20% losses at 75° (see Fig. 5). The influence of a module’s op- tical characteristics on its performance is shown in Fig. 6. Acta Polytechnica Vol. 41 No. 2/2001 12 Fig. 4: Transmission through a thin absorbing layer Fig. 5: Calculated relative transmission, reflection, and absorp- tion of radiation for glass cover 4 Thermal performance of a photovoltaic module Solar radiation irradiating at the surface of the solar mo- dule is absorbed partially by various materials of the module. Almost all radiation is absorbed in the silicon. Only a very small part is absorbed in other materials of the solar module. This will produce heat inside the solar cell. There are heat loss mechanisms from the front and rear surface of the module: free convection, wind convection and radiative heat loss. Inside cell heat is transmitted by conduction. The equations describing the heat flows within the cell are given below: Radiation heat transfer equation: The total radiative heat flow from the module, Qrad [W] to the environment comprises the components of the radiation exchange of the front and rear sides of the module with the sky and the ground. These components are also functions of the module surface temper- ature, T(i, j, k), ambient air temperature, T� and the relevant emissivity � also; these components are determined by the so-called radiation shape factor, Fs the radiative surface area of the module, Am and the Stefan-Boltzmann constant �. In general this can be expressed as [11]: � �� �Q F T Trad s m i, j, k� � �� � � 4 4 . (28) Convection heat transfer equation: Convective heat flows, Qconv. cannot be treated in a closed mathematical model and, therefore have to be computed by iteration or approximation. The convective heat transfer coefficient h� [W/m 2K] is a bulky function of the actual air temperature and module surface temperature, air viscosity, thermal conductivity and heat ca- pacity, and finally wind speed. The models used are based on the following equation: � �� �Q h A T Tconv. c i, j, k� �� � . (29) Conduction heat transfer equation: can be given as: � � � �� �Q K A h T Tcond. th m i, j, k i+1, j, k� � . (30) 4.1 Problem formulation The solution starts by assuming that the photovoltaic solar module is composed of three layers. For the purposes of the 3D-thermal model the solar module is divided into a number of meshes, each mesh having six faces, (XY = Top, YX = Bot- tom, YZ = Front, ZY = Back, XZ = Right and ZX = Left), i, j, k are the coordinates of the meshes. The computer code was written in Fortran. Simulation outputs allow a complete evaluation of the temperature distribution inside the photo- voltaic module and can be of great interest both for system designers and for researchers. The technique used in this code involves calculating the temperature by an explicit method using iteration processes. The study is confined to a single silicon photovoltaic module of the following dimen- sions: thickness of silicon 0.5 mm, sandwiched between a glass layer which is taken as 3 mm and Mylar 0.5 mm, giving a total photovoltaic module thickness of 4 mm. Both length and width are taken as 40 mm. Also, it is assumed that all the heat is absorbed inside the silicon wafer homogeneously. The results are shown in Figs. 7 and 8, which show the temperature distribution along the XY-Plane at the mid Z-axis of the silicon wafer, and the temperature distribution along the XZ-Plane at mid Y-axis. © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 13 Acta Polytechnica Vol. 41 No. 2/2001 Fig. 6: Influence of angle of incidence on cell parameters Fig. 7: Temperature distribution along the XY plane at mid Z-axis Fig. 8: Temperature distribution along the XZ plane at mid Y-axis 4.2 Temperature dependence of solar module parameters Fig. 9 shows the dependence of Isc, Voc, FF, and Pm on tem- perature. As the cell temperature increases, the short circuit current, Isc somewhat increases, and the maximum power, Pm the open circuit voltage, Voc and the fill factor, FF de- crease. The short circuit current increase with temperature is attributed to light absorption variation. Under the conditions considered, the solar absorption edges shift due to a decrease in the crystal forbidden gap width, and the number of pairs generated in the bulk increases. Alternatively, as the tempera- ture increases, the intrinsic light absorption factor changes and the short circuit current increases. To explain the varia- tion of Voc with temperature, the effects of all voltage control- ling parameters (mobility, lifetime, energy band gap and absorption coefficient) have to be included in the analysis. In the general case, open circuit voltage is determined by the relation � �V I Ioc sc d� ln , by which Voc must grow as Isc increases. Depending on the current flowing mechanism, the Acta Polytechnica Vol. 41 No. 2/2001 14 For the meshes at the corners of the cell For i = 1, j = 1, k = 1 � � � �Q F A T T F A Ti, j, k s xz xz i, j, k xz 4 s xy xy i, j, � �� � � ��� � � � 4 � � � �k xy 4 s yz yz i, j, k yz 4 xz x 4 4�� � � � � � � � �� � �T F A T T h A� � � �� � � �� � � � z i, j, k xz xy xy i, j, k xy yz yz i, j, k y T T h A T T h A T T � � � � � � � � � � �� � � � � �� � � � z th zx y i, j, k i, j+1, k th yx z i, j, k i, j, k m m � � � � � K A h T T K A h T T� �� � � � � �� �+1 th zy x i, j, k i+1, j, km� �K A h T T (31) For the meshes at the edges of the cell For i = 1, k = 1, j = 2, …, (ny 1) � � � � � �Q F A T T h A Ti, j, k s yz yz i, j, k yz 4 yz yz i, j, k� � � � � � �� �� � 4 � � � � � �T F A T T h A T T� � � �� �� � � � �yz s xy xy i, j, k xy 4 xy xy i, j, k� � 4 � � � � � �� � � � xy th yx z i, j, k i, j, k+1 th zy x i, j, k i+1,m m � � � � �K A h T T K A h T T� �� � � � � �� � � � j, k th xz y i, j, k i, j+1, k th zx y i, j, k i, j m m � � � � � K A h T T K A h T T� �� ��1, k (32) For the meshes inside the cell body For i = 2, …, (nx 1), j = 2, …, (ny 1), k = 2, …, (nz 1) � � � � � �� � � �Q K A h T T K A h T Ti, j, k th xy z i, j, k i, j, k 1 th yx z i, j, k im m� � � �� � �� � � � � �� � � � , j, k+1 th zy x i, j, k i 1, j, k th yz x i, j, k m m � � � � � �K A h T T K A h T T� �� � � � � �� � � �i+1, j, k th xz y i, j, k i, j 1, k th zx y i, j, km m� � � ��K A h T T K A h T T� �� �i, j+1, k (33) For the meshes at the surface For i = 2, …, (nx 1), j = 2, …, (ny 1), k= 1 � � � � � �Q F A T T h A Ti, j, k s xy xy i, j, k xy 4 xy xy i, j, k� � � � � � �� �� � 4 � � � � � �� � � � T K A h T T K A h T T � � � � � � � xy th yx z i, j, k i, j, k+1 th yz x i, j, k i m m � �� � � � � �� � � �1, j, k th zy x i, j, k i+1, j, k th zx y i, j, k i,m m� � � �K A h T T K A h T T� �� � � � � �� � j 1, k th xz y i, j, k i, j+1, km � � � �K A h T T (34) For the meshes at the gap � � � � � �� � � �Q K A h T T K A h T Ti, j, k th yz x i, j, k i 1, j, k th zy x i, j, k im m� � � �� � �� � � � � �� � � � +1, j, k th xz y i, j, k i, j+1, k th zx y i, j, km m � � � � �K A h T T K A h T T� �� � � � � �� � � � i, j 1, k th xy z i, j, k i, j, k 1 gap yx i, j, k im � � � � � � �K A h T T h A T T� �� �, j, k+1 (35) Where: hx width of mesh in x, hy width of mesh in y, hz width of mesh in z, Axz hxhz, Axy hxhy, Ayz hyhz, m material number, nx total No. of meshes in x-direction, ny total No. of meshes in y-direction and nz total No. of meshes in z-direction. value of Id is proportional to ni 2. As the temperature increases ni grows exponentially and causes exponential growth of Id. Moreover, the crystal forbidden gap narrowing with tempera- ture increase leads to an increase in the dark current, Id and minimizes the positive effect of the absorption factor. As a re- sult, Voc decreases. The maximum power across the optimal load and the efficiency depend strongly on temperature; the shape of this dependence is close to linear. Conclusions A simple and accurate close-form solution for the optical and thermal behaviors of a solar cell array were investigated by simulation and by some verifying experiments, in order to obtain more accurate results in predicting the actual per- formance of photovoltaic modules, without resorting to a lengthy, time-consuming iterative solution which has to be repeated for any change in the parameters of the array. Nomenclature GH Hourly global radiation incident on a horizontal sur- face, kWh/m2/day GT Hourly global radiation incident on an inclined sur- face, kWh/m2/day Ib, H Hourly beam radiation incident on an inclined sur- face, kWh/m2/day Ib, T Hourly beam radiation incident on a horizontal sur- face, kWh/m2/day Id, H Hourly sky diffuse radiation incident on a horizontal surface, kWh/m2/day Id, T Hourly sky diffuse radiation incident on an inclined surface, kWh/m2/day Io Extraterrestrial hourly radiation incident on a hori- zontal surface, kWh/m2/day Ig, H Hourly ground reflected radiation incident on a hori- zontal surface Ig, T Hourly ground radiation incident on an inclined sur- face, kWh/m2/day IT, cir Circumsolar radiation component IT, iso Isotropic radiation component n Refractive index AI Anisotropic index Rb The ratio of beam radiation on the tilted surface to that on a horizontal surface r Reflection of incident � Surface slope, degrees � Declination, degrees � Latitude, degrees � Ground albedo �1 Incidence angle, the angle between beam radiation and the vertical, degrees �o Angles of refraction � Sunrise hour angle, degrees � Surface azimuth angle, degrees E Eo i � �, Amplitudes of the reflected and transmitted waves � Reflectance Qrad Total radiative heat flow from the module, W Qconv Convective heat flow, W Qcond Quantity of heat transfer by conduction from one mesh to another, W Kth Thermal conductivity of different materials of solar cell, W/mK T(i, j, k) Module surface temperature, °C T� Ambient temperature, °C � Emissivity Fs Radiation shape factor Am Radiative surface area of the module Ac The surface area of the cell, m 2 References [1] Born, M., Wolf, E.: Principles of optics. 5th Edition, Pergamon, New York, 1975 [2] Bugler, J.: The determination of hourly insolation on a tilted plane using a diffuse irradiance model based on hourly mea- sured global horizontal insolation. Solar Energy, 1977, Vol. 19, No. 5, pp. 477 [3] Dave, J.: Validity of the isotropic distribution approximation in solar energy estimations. Solar Energy, 1977, Vol. 19, pp. 331 [4] Duffie, J., Beckman, W.: Solar engineering of thermal pro- cesses. Wiley, New York, 1980 [5] Fraidenraich, N., Vilela, O.: Exact solutions for multiplayer optical structures: Application to PV modules. Solar Energy, 2000, Vol. 69, No. 5, pp. 357 [6] Gardner, C., Nadeau, C.: Estimating south slope irradiance in the Arctic – A comparison of experimental and modeled values. Solar Energy, 1988, Vol. 41, pp. 227 [7] Green, M.: High efficiency silicon solar cells. Trans Tech Publications, Switzerland, 1987 [8] Gueymard, C.: An anisotropic solar irradiance model for tilted surfaces and its comparison with selected engineering algorithms. Solar energy, 1987, Vol. 38, p. 367 [9] Hay, J.: Calculation of monthly mean solar radiation for hori- zontal and tilted surfaces. Solar Energy, 1979, Vol. 23, p. 301 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 15 Acta Polytechnica Vol. 41 No. 2/2001 Fig. 9: Variation of key parameters of a silicon solar cell at high temperatures 16 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 41 No. 2/2001 [10] Hay, J., Davies, J.: Calculation of the solar radiation incident on an inclined surface. Proceedings First Canadian Solar Radiation Data Workshop, 1980, p. 59 [11] Holman, J.: Heat Transfer. Mc Graw Hill, 1982 [12] Hottel, H. C., Woertz, W.: Performance of flat plate solar heat collectors. Transactions of American Society of Mechanical Engineers, 1942, Vol. 64, pp. 91 [13] Ineichen, P., Perez, R., Seals, R.: The importance of correct albedo determination for adequately modeling energy received by tilted surfaces. Solar Energy, 1987, Vol. 39, pp. 221 [14] Iqbal, M.: An Introduction to solar radiation. Academic Press, Toronto, 1983 [15] Klein, M., Furtak, T.: Optics. Wiley, New York, 1993 [16] Klucher, T.: Evaluation of models to predict insolation on tilted surfaces. Solar Energy, 1979, Vol. 23, p. 114 [17] Martin, N.: Comparative study of the angular influence on different PV module technologies. Proceeding 14th E.C. PVSEC, 1997 [18] Nahumura, S.: Applied Numerical Methods with Software. Prentice Hall, New York, 1991 [19] Perommer, P., Lomas, K., Seale, C., Kupke, C.: The radiation transfer through coated and tinted glazing. Solar Energy, 1995, Vol. 54, No. 5, p. 287 [20] Smith, G.: Numerical solution of partial differential equations: Finite Difference Method. Oxford, 1979 [21] Temps, R. C., Coulson, K. L.: Solar radiation incident upon slopes of different orientations, Solar Energy, 1977, Vol. 19, p. 179 [22] Wenham, S., Green, M., Watt, M.: Applied photovoltaics. University of New South Wales, Center of Photovoltaic Systems and Devices, Australia, 1994 [23] Harbeke, B.: Coherent and Incoherent Reflection and Trans- mission of Multi-layer Structures. Applied Physics, 1989, Vol. 39, pp. 165–170 Eng. Hamdy Kamal Elminir e-mail: ehamdy@hotmail.com phone: +420 2 2435 2212 fax: +420 2 2435 3949 Prof. Ing. Vítězslav Benda, CSc. Ing. Ivan Kudláček, CSc. Czech Technical University in Prague Faculty of Electrical Engineering Technická 2, 166 27 Praha 6, Czech Republic << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends false /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) /CZE >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice