Model Associated with the Study of the Degradation Based on the Accelerated Test : A Literature Review

Most manufacturers of solar modules guarantee the minimum performance of their modules for 20 to 25 years, and 30-year warranties have been introduced. The warranty typically guarantees that the modules will perform to at least 90% capacity in the first 10 years and to at least 80% in the following 10 15 years. Early degradation resulting from design flaws, materials or processing issues is often apparent from startup to the first few years in service. Importantly, many module failures and performance losses are the result of gradual accumulated damage resulting from long-term outdoor exposure in harsh environments, referred. Many of these processes occur on relatively long time scales and the various degradation processes may be chemical, electrical, thermal or mechanical in nature. These are either initiated or accelerated by the combined stresses of the service environment, in particular solar radiation, temperature and moisture, and other stresses such as salt air, wind and snow. Accelerated Life Testing (ALT) test methodology is normally predicated on first being able to reproduce a specific degradation or failure mode without altering it (correlation); and, second, to produce that result in less than real-time acceleration. Degradation and failure may result when an applied stress exceeds material or product strength. This may be a one-time catastrophic event, the result of cyclic fatigue, or a gradual decline in requisite properties due to ageing mechanisms. Engineers in the manufacturing industries have used accelerated test (AT) experiments for many decades. The purpose of AT experiments is to acquire reliability information quickly. Test units of a material, component, subsystem or entire systems are subjected to higher-than-usual levels of one or more accelerating variables such as temperature or stress. Then the AT results are used to predict life of the units at use conditions. The extrapolation is typically justified (correctly or incorrectly) on the basis of physically motivated models or a combination of empirical model fitting with a sufficient amount of previous experience in testing similar units. The need to extrapolate in both time and the accelerating variables


Introduction
Today's manufacturers face strong pressure to develop new, higher-technology products in record time, while improving productivity, product field reliability and overall quality.This has motivated the development of methods like concurrent engineering and encouraged wider use of designed experiments for product and process improvement.The lack of accurate information about degradation rate increases the financial risk [1].The reliability and lifetime of a photovoltaic system depend mainly on the energy performance of the modules and their different degradation mode [2].A common misconception is since PV has a track record dating from the 1970's, that module reliability is now a given [3].According to literature, results from experimental tests are given.Different techniques and testing methods are used.Most of these results obtained from laboratory tests are not always accurate and consistent on the degradation rates of PV modules announced.In 2002, the National Renewable Energy Laboratory (NREL) estimated the annual performance degradation of monocrystalline and polycrystalline PV modules at 0.7% [4].John H. Wohlgemuth, formerly with BP Solar and now with NREL, describes the fundamental problem: Outdoor testing is a must, but it takes much too long to be of much use as decision maker.We clearly cannot wait 25 years or even a significant fraction of 25 years to introduce a new product.Therefore, we must develop and utilize accelerated tests to qualify these new products.Akira Terao, Chief Reliability Engineer for Sun power Corporation goes further to identify some of the life-testing roadblocks: 25 year warranty, Ill-defined field conditions, Harsh and varied outdoor conditions, Materials used near their limits; Limited acceleration factors ⇒ long tests, Large samples, small sample sizes; Subtle polymer chemistry; Cumulative effects, positive feedback loops.Clearly, assessing PV module lifetime performance is not a simple task.The LEEETISO (Laboratory of Energy, Ecology and Economy Solar-Ticino), test center of photovoltaic modules in Switzerland, stated that the power degradation rate of crystalline silicon PV modules could go from 0.7% to 9.8% during the first exposure year and 0.7% to 4.9% during the second one [5].Also, that early problems, such as encapsulate yellowing or TCO corrosion are now all resolved.However, since circa 2005, silicon wafer thickness has been greatly reduced, new lower cost materials have emerged, and other changes have been introduced making today's c-Si modules very different from previous designs.And non-crystalline silicon PV technologies have different failure modes where prior silicon history isn't useful.Jordan [6] states that in the 2000s, degradation rate measured on individual modules was 5% on average.There are only a few long-term studies on the degradation of PV modules published.Skoczek [7] have measured the performance of 204 field-aged crystalline Si based PV modules (53 module types).Exposing started in 1983 at the Joint Research Center in north Italy with a moderate subtropical climate (−10 to 35˚C, >90% RH).Today, this performance warranty is frequently >90% of rated power after 10 years and >80% after 20 or 25 years although 30 year guarantees have now appeared.These high performance losses (>20%) are related to losses of the fill factor, caused by an increased series resistance.The moderate performance losses (<20%) can be related to losses of the short-circuit current, caused by degradation of the optical properties.The long term losses are determined to be between 0.2% and 1.0% per annum.
However, there is little information on PV modules degradation modes in terms of frequency, speed of evolution and degree of impact on module lifetime and reliability.
Research on photovoltaic modules is rather focused on the race to develop new technologies without sufficient experience feedback on already operational technologies [8] [9].The requirements for higher reliability have increased the need for more up-front testing of materials, components and systems.This is in line with the modern quality philosophy for producing high-reliability products: achieve high reliability by improving the design and manufacturing processes; move away from reliance on inspection (or screening) to achieve high reliability, as described in Meeker and Hamada [10] and Meeker and Escobar [11] estimating the failure-time distribution or long-term performance of components of high-reliability products is particularly difficult.Most modern products are designed to operate without failure for years, decades or longer.Thus few units will fail or degrade appreciably in a test of practical length at normal use conditions.
For such applications, Accelerated Tests (ATs) are used in manufacturing industries to assess or demonstrate component and subsystem reliability, to certify components, to detect failure modes so that they can be corrected, to compare different manufacturers, and so forth.ATs have become increasingly important because of rapidly changing technologies, more complicated products with more components, higher customer expectations for better reliability and the need for rapid product development.There are difficult practical and statistical issues involved in accelerating the life of a complicated product that can fail in different ways.Generally, information from tests at high levels of one or more accelerating variables (e.g., use rate, temperature, voltage or pressure) is extrapolated, through a physically reasonable statistical model, to obtain estimates of life or longterm performance at lower, normal levels of the accelerating variable(s).
The present article proposes a bibliographic review on the degradation modes of PV modules [12]- [14] as well as the associated models for the study of the degradation of their power based on di-paid whose accelerated approaches used to determine the parameters.
The paper is organized in two parts, the main types of degradation of the PV module identified in the literature are exposed in the first part.A review of models associated with the degradation of the PV modules is presented in the second part.

Modes and Causes of Degradation of a PV Module
Manufacturer's guarantees of up to 20 years indicate the quality of bulk silicon PV modules currently being produced.Nevertheless, there are several failure modes and degradation mechanisms which may reduce the power output or cause the module to fail.Nearly all of these mechanisms are related to water ingress or temperature stress.Degradation mechanisms may involve either a gradual reduction in the output power of a PV module over time or an overall reduction in power due to failure of an individual solar cell in the module.They can deteriorate or become defective even when they operate on actual site for extended periods [15].
The degradation of PV modules are caused by static and dynamic mechanical loads, thermal cycling, exposure to radiation, moisture, hail impact (sand) dust accumulation, partial occultation causes found above are responsible for the degradation of the modules, which manifests itself in several ways: Corrosion of the PV module, the discoloration, delamination, breakage of the front panel (glass) module, Potential Induced Degradation (PID), hot spots (hot spots), Bubbles, dust deposits.

Method for Determination of the Degradation
The outdoor exhibition of photovoltaic modules includes several factors (wind, rain, light, heat ....) That causes damage as mentioned above.Thus these factors have a negative effect on the module power output.
To better understand these degradations, called accelerated tests can reproduce they are conducted in climatic chambers in extensive conditions of temperature and pressure to reproduce the degradation in a limited time.But it is also interesting to have a model to simulate numerically the evolution of degradation.Many authors have tried to make the degradation assessment models

The Model of Pan
Pan [16] offers the general degradation model of a photovoltaic module following ( ) ( ) where a and b are parameters that can be assessed by degradation testing.The parameter a is considered constant regardless of the level of sévérisation and b is variable reflecting its dependence on sévérisation levels.Thanks to this relationship, we can determine when the test must be stopped in order to measure the power of the module.
Wohlgemuth [17] achieved a damp heat test at a temperature of 85˚C and a relative humidity of 85%, which corresponds to the sévérisée condition that we advocate and which is that of 61215 IEC standards IEC 61646 and IEC 62108.The author followed the power loss of the photovoltaic module during the test and obtained the curve shown in Temperatures taken during the tests shall not exceed the technological limit temperature of PV module sequal to 120˚C according to Kern [18].
The 85C/85% relative humidity exposure is as accelerated as necessary.These conditions probably never happen in the real world as the modules tend to dry out at their highest temperatures, but absorb moisture at lower temperatures.It is difficult to judge what outdoor exposure the 1000 hour exposure at 85/85 represents.In a recent experiment 10 crystalline silicon modules qualified to IEC 61215 were exposed to 1250 hours of 85/85.Only 2 of the 10 types successfully passed the extended test [18].On the other hand some glass-glass encapsulated modules can easily endure more than 2000 hours of damp heat exposure.So rather than specify a particular length of time, it seems appropriate to test the control technology and the new technology through enough hours of damp heat that both begin to lose some power (say to 90% of the original) in order to verify that the new technology is no worse and has no additional failure modes than the old module technology it will replace.
The study of Wohlgemuth [17] to determine the parameters of the degradation model proposed by Pan [16]: a = 3.0868 and b = 5.762 × 10 −12 for a temperature of 85˚C and a relative humidity 85%.The parameter is considered constant whatever the test conditions (temperature and humidity) while the parameter b depends on the test conditions.For the model of Peck, b of the degradation model parameter to a temperature T and relative humidity RH can be determined by the relation   ( ) ( ) ( ) where n is the number of degradation modes.

The Exponential Model
The PV module power measured at a given time follows a Gaussian distribution.According to Reis, Sakamoto and Oshiro [18] [22], the probability density is given by: ( ) where P is the power module, μ is its mean value and σ its standard deviation.Thus, it can be calculated by Osterwald, Marion, Adelstein, and Raghuraman [23]- [25]: where P 0 is the average power at t = 0 (i.e. the nominal module power), A is a parameter which reflects annual decrease of module power, t is time in years.Of course, validity of Equation ( 6) is restricted by time (t) lower than P 0 /A.Another limitation comes from the assumption according to which A is constant in time.
( ) ( ) A/P 0 (year −1 ) is the annual degradation rate.From studies on electronic components degradation, we can consider the assumption according to which degradation rate is exponential [26] [27] ( ) ( ) Figure 2 shows PV module degradation as a function of time from the exponential model.

Using Temperature to Accelerate Failure Mechanisms
It is sometimes said that high temperature is the enemy of reliability.Increasing temperature is one of the most commonly used methods to accelerate a failure mechanism.

Arrhenius Relationship for Reaction Rates
The Arrhenius relationship is a widely used model to describe the effect that temperature has on the rate of a simple chemical reaction.This relationship can be written as  where R is the reaction rate, and temp˚K = temp ˚C + 273.15 is thermodynamic temperature in kelvin (K), k is either Boltzmann's constant or the universal gas constant and Ea is the activation energy.
The parameters Ea and 0 γ are product or material characteristics.In applications involving electronic com- ponent reliability, Boltzmann's constant k = 8.6171 × 10 −5 = 1/11605 in units of electron-volt per kelvin (eV/K) is commonly used and in this case, Ea has units of electron-volt (eV).
In most accelerated test applications, it would be more appropriate to refer to Ea in Equation ( 9) as a quasiactivation energy.

Arrhenius Relationship Time-Acceleration Factor
The Arrhenius acceleration factor is will be used to denote a time-acceleration factor.The Arrhenius-based model presents some limits.Indeed, the Arrhenius equation can be used to quantify the effect of varying temperature and irradiance on the rate of a property change.However it cannot provide a complete picture of the long-term degradation of PV modules, as other stress factors or combination of stresses are involved.These include moisture (inducing physical and chemical processes and generating mechanical stress in combination with temperature), temperature cycling (generating thermomechanical stress), electricity production (inducing electrical and electrochemical stresses), and other externally applied stresses (wind, hail, airborne pollutants.

Eyring Relationship Time-Acceleration Factor
The Arrhenius relationship Equation ( 9) was discovered by Svante Arrhenius through empirical observation in the late 1800s.Eyring [29] [30] gives physical theory describing the effect that temperature has on a reaction rate.Written in terms of a reaction rate, The Eyring relationship is: where A(temp) is a function of temperature depending on the specifics of the reaction dynamics and 0 and Ea are constants Weston and Schwarz, provides more detail [31].Applications in the literature have typically used A(temp) = (tempK) m with a fixed value of m ranging between m = 0 [32], m = 0.5 [33], to m = 1 [34] [35].The Eyring relationship temperature acceleration factor is , , AF temp temp Ea is the Arrhenius acceleration factor from Equation (10) For use over practical ranges of temperature acceleration, and for practical values of m not far from 0, the factor outside the exponential has relatively little effect on the acceleration factor and the additional term is often dropped in favor of the simpler Arrhenius relationship.

Using Humidity to Accelerate Reaction Rates
Humidity is another commonly used accelerating variable, particularly for failure mechanisms involving corrosion and certain kinds of chemical degradation.
Vapor density measures the amount of water vapor in a volume of air in units of mass per unit volume.Partial vapor pressure (sometimes simply referred to as "vapor pressure") is closely related and measures that part of the total air pressure exerted by the water molecules in the air.Partial vapor pressure is approximately proportional to vapor density.
The partial vapor pressure at which molecules are evaporating and condensing from the surface of water at the same rate is the saturation vapor pressure.
For a fixed amount of moisture in the air, saturation vapor pressure increases with temperature.
Relative humidity is usually defined as and is commonly expressed as a percent.For most failure mechanisms, physical/chemical theory suggests that RH is the appropriate scale in which to relate reaction rate to humidity especially if temperature is also to be used as an accelerating variable [33].
A variety of different humidity models (mostly empirical but a few with some physical basis) have been suggested for different kinds of failure mechanisms.
Vapor Pressure RH Saturation Vapor Pressure = (13) Much of this work has been motivated by concerns about the effect of environmental humidity on plasticpackaged electronic devices.Humidity is also an important factor in the service-life distribution of paints and coatings.In most test applications where humidity is used as an accelerating variable, it is used in conjunction with temperature.
Peck [36] presents data and models relating life of semiconductor electronic components to humidity and temperature.See also Peck and Zierdt [37] and Joyce et al. [38].Gillen and Mead [39] describe a kinetic approach for modeling accelerated aging data.LuValle, Welsher and Mitchell [40] describe the analysis of timeto-failure data on printed circuit boards that have been tested at higher than usual temperature, humidity and voltage.
They suggest ALT models based on the physics of failure.Nelson [34] and Boccaletti et al. [32] review and compare a number of different humidity models.
The Eyring/Arrhenius temperature-humidity acceleration relationship in the form of Equation ( 33) uses x 1 = 11605/temp K, x 2 = log (RH) and x 3 = x 1 x 2 where RH is relative humidity, expressed as a proportion.
An alternative humidity relationship suggested by Klinger [33], on the basis of a simple kinetic model for corrosion, uses the term x 2 = log [RH/(1−RH)] (a logistic transformation) instead.
In most applications where it is used as an accelerating variable, higher humidity increases degradation rates and leads to earlier failures.In applications where drying is the failure mechanism, however, an artificial environment with lower humidity can be used to accelerate a test.

Acceleration Model for Photodegradation
Many organic compounds degrade chemically when exposed to ultraviolet (UV) radiation.Such degradation is known as photodegradation.This section describes models that have been used to study photodegradation and that are useful when analyzing data from accelerated photodegradation tests.Many of the ideas in this section originated from early research into the effects of light on photographic emulsions [41] and the effect that UV exposure has on causing skin cancer [42].
Important applications include prediction of service life of products exposed to UV radiation (outdoor weathering) and fiber-optic systems.
It is clear that UV radiation is a major degradation factor for photovoltaic modules especially in their discoloration [43]- [45].The total dose of UV radiation can be considered as the number of photons absorbed by the material degrading and that will cause its chemical change.For PV module, this degradation is reflected by the change of the encapsulating module transmittance that induces a decrease in the PV module current-voltage characteristic.This decrease is quantified by the relative variation of module short-circuit current representing the degradation factor given by Zimmerman [46]: ) where I sc (E) is the short-circuit current of the PV module and E the dose of ultra-violet (UV).Zimmerman [46] proposes the degradation expression factor given by the following equation: ( ) P λ is the spectral power density of the sun.[ ] min max , λ λ represents the integration interval for the wavelength range in which the spectral response of the PV cell is not zero.For a Photodegradation ( ) ( ) in the encapsulant less than 70% in the range [k min , k max ], the transmittance can be written: where a cmx and b cmx are parameters of material used for PV cell [46].The relationship between UV dose E and exposure time t in the solar spectrum P(k) is: The integral extends up to 400 nm that represents a practical limit to UV photodegradation.
Combining Equation ( 15) and ( 17) and using the mean value theorem, Equation ( 18) obtained to estimate the UV degradation of PV module is given by: ( ) ( ) ( ) 0 T λ .The measurement of these characteristics requires a rigorous instrumentation without which the model accuracy is compromised.

Time Scale and Model for Total Effective UV Dosage
As described in Martin et al. [47], the appropriate time scale for photodegradation is the total (i.e., cumulative) effective UV dosage, denoted by D Tot .Intuitively, this total effective dosage can be thought of as the cumulative number of photons absorbed into the degrading material and that cause chemical change.The total effective UV dosage at real time t can be expressed as ( ) ( ) where the instantaneous effective UV dosage at real time τ is Here ( ) 0 , E τ λ is the spectral irradiance (or intensity) of the light source at time τ (both artificial and natural light sources have potentially time-dependent mixtures of light at different wavelengths, denoted by λ ), ( ) ( )  is the spectral absorbance of the material being exposed (damage is caused only by photons that are absorbed into the material), and ( ) φ λ is a quasi-quantum efficiency of the absorbed radiation (allow- ing for the fact that photons at certain wavelengths have a higher probability of causing damage than others).
The functions E 0 and A in the integrand of Equation ( 20) can either be measured directly or estimated from data and the function ( ) φ λ can be estimated from data.A simple log-linear model is commonly used to describe quasi-quantum efficiency as a function of wavelength.That is: The integrals over wavelength, like that in Equation ( 20), are typically taken over the UV-B band (290 nm to 14,320 nm), as this is the range of wavelengths over which both ( ) φ λ and ( ) 0 , E t λ are importantly different from 0. Longer wavelengths (in the UV-A band) are not terribly harmful to organic materials  .Shorter wavelengths (in the UV-C band) have more energy, but are generally filtered out by ozone in the atmosphere ( )

Model for Photodegradation and UV Intensity
Degradation (or damage) D(t) at time t depends on environmental variables like UV, temp and RH, that may vary over time, say according to a multivariable profile ξ(t) = [UV, temp, RH, ⋅⋅⋅].Laboratory tests are conducted in well-controlled environments, usually holding these variables constant (although sometimes such variables are purposely changed during an experiment, as in step-stress accelerated tests).Interest often centers, however, on life in a variable environment.These sample paths might be modeled by a given functional form, where z is scaled time and g(z) would usually be suggested by knowledge of the kinetic model (e.g., linear for zeroth-order kinetics and exponential for first-order kinetics), although empirical curve fitting may be adequate for purposes where the amount of extrapolation in the time dimension is not large.As in SAFT models, μ can be modeled as a function of explanatory variables like temperature and humidity when these variables affect the degradation rate.

Inverse Power Relationship
The inverse power relationship is frequently used to describe the effect that stresses like voltage and pressure have on lifetime.Voltage is used in the following discussion.When the thickness of a dielectric material or insulation is constant, the voltage stress is proportional to the square of the voltage.Let volt denote voltage and let volt U be the voltage at use conditions.The lifetime at stress level volt is given by where 1 β , in general, is negative.The 1model has SAFT form with acceleration factor If T(volt U ) has a log-location-scale distribution with parameters μ U and σ, then T(volt) also has a log locationscale distribution with , where x U = log(volt U ), x = log(volt), 0 and σ does not depend on x.
The inverse power relationship is widely used to model life as a function of pressure-like accelerating variables (e.g., stress, pressure, voltage stress).
This relationship is generally considered to be an empirical model because it has no formal basis from knowledge of the physics/chemistry of the modeled failure modes.It is commonly used because engineers have found, over time, that it often provides a useful description of certain kinds of AT data.
Analytically, suppose that degrading dielectric strength at age t can be expressed as ( ) ( ) Then the acceleration factor for volt versus volt U is which is an inverse power relationship, as in Equation (24).To extend this model, suppose that higher voltage also leads to an increase in the degradation rate and that this increase is described with the degradation model: where ( ) ( )  .Then equating D(t) to volt and solving for failure time t gives the failure time ( ) ( ) Then the ratio of failure times at volt U versus volt is the acceleration factor which is again an inverse power relationship with β 1 = 1 − 2. This motivation for the inverse power relationship described here is not based on any fundamental understanding of what happens to the insulating material at the molecular level over time.

Generalized Eyring Relationship
The generalized Eyring relationship extends the Eyring relationship in Section 3.3.3,allowing for one or more nonthermal accelerating variables (such as humidity or voltage).For one additional nonthermal accelerating variable X, the model, in terms of reaction rate, can be written as ( ) ( ) where X is a function of the non-thermal stress.The parameters 1 Ea γ = (activation energy) and 0, 2, 3 are characteristics of the particular physical/chemical process.Additional factors like the one on the right-hand side of Equation ( 31) can be added for other non-thermal accelerating variables.
In the following sections, following common practice, we set (tempK) m = 1, using what is essentially the Arrhenius temperature-acceleration relationship.
These sections describe some important special-case applications of this more general model.If the underlying model relating the degradation process to failure is a SAFT model, the generalized Eyring relationship can be used to describe the relationship between times at different sets of conditions temp and X.In particular, the acceleration factor relative to use conditions temp U and X U is Suppose that T(temp U ) (time at use or some other baseline temperature) has a log location scale distribution with parameters μ U and σ.Then T(temp) has the same log-location-scale distribution with ( ) where

Temperature-Voltage Acceleration Models
Many different models have been used to describe the effect of the combination of temperature and voltage on acceleration.For instance, Meeker and Escobar [48] analyzed data from a study relating voltage and temperature to the failure of glass capacitors.They modeled the location parameter of log-lifetime as a simple linear function of temp ˚C and volt.The generalized Eyring relationship can also be used with X = log(volt), as done in Boyko and Gerlach [49].Klinger [33] modeled the Boyko and Gerlach [49] data by including second-order terms for both accelerating variables.To derive the time-acceleration factor for the extended Arrhenius relationship with temp and volt, one can follow steps analogous to those outlined in Section 3.6.Using Equation ( 31) with X = log(volt), one obtains Again, failure occurs when the dielectric strength crosses the applied voltage stress, that is, D(t) = volt.This occurs at time ( ) ( ) From this, one computes where ( ) ( ) . When 3 0 γ = , there is no interaction between tem- perature and voltage.In this case, AF(temp, volt) can be factored into two terms, one that involves temperature only and another term that involves voltage only.Thus, if there is no interaction, the contribution of temperature (voltage) to acceleration is the same at all levels of voltage (levels of temperature).

Temperature-Humidity Acceleration Models
Relative humidity is another environmental variable that can be combined with temperature to accelerate corrosion or other chemical reactions.
Gillen and Mead [39] and Klinger [33] studied kinetic models relating aging with humidity.LuValle, Welsher and Mitchell [40] provided physical basis for studying the effect of humidity, temperature and voltage on the failure of circuit boards.See Boccaletti et al. [32], Nelson [34], Joyce et al. [38], Peck [36] and Peck and Zierdt [37] for ALT applications involving temperature and humidity.
The extended Arrhenius relationship Equation (31) applied to ALTs with temperature and humidity uses x 1 = 11605/tempK, x 2 = log(RH) and x 3 = x 1 x 2 where RH is a proportion denoting relative humidity.The case when 3 0 β = (no temperature-humidity interaction) is known as "Peck's relationship" and was used by Peck [36] to study failures of epoxy packing.
Klinger [33] suggested the term x 2 = log[RH/(1−RH)] instead of log(RH).This alternative relationship is based on a kinetic model for corrosion.
The acceleration model with the capacity to take into account temperature (T) and relative humidity (RH) is the Peck model [50].
According to Charki et al. [19], the Wiener process can be estimated as follows: In this case, x 1 and x 2 correspond respectively to relative humidity (RH) and module temperature (T); 0 1 x and 0 2 x respectively correspond to known relative humidity (HR 0 ) and temperature (T 0 ) in the reference conditions.The acceleration factor is expressed as [19]: ( ) , n Ea γ = .We consider that 0 a a = for each stress level.The relation Equation (36)  .
If the acceleration factor is defined by relation Equation (37), it is possible to estimate an equivalent temperature, T eq , which represents the degradation that would have occurred if the module had been aged for the same length of time but at a constant temperature.The equivalent temperature T eq can be calculated using the following relation [51]: where t is time, T module(t) is the time-dependent module temperature, and t 1 and t 2 are the integration start and end times.
Using the same methodology below, the equivalent of relative humidity H eq can be estimated as where t is time, H(t) is the time-dependent environmental relative humidity, and t 1 and t 2 are the integration start and end times.
The module temperature T module(t) and the relative humidity H(t) are simulated hour by hour for a period of 50 years, for at four cities: Paris (France), Oslo (Norway), Madrid (Spain) and Tamanrasset (Algeria).The equivalent temperature T eq and the equivalent relative humidity H eq obtained versus Ea and n respectively using relations Equations ( 39) and ( 40) are plotted in Figure 3.

Modeling Photodegradation Temperature and Humidity Effects
When modeling photodegradation, it is often necessary to account for the effect of temperature and humidity.The Arrhenius rate reaction model Equation ( 9) can be used to scale time (or dosage) in the usual manner.Humidity is also known to affect photodegradation rate.Sometimes the rate of degradation will be directly affected by moisture content of the degrading material.In this case one can use a model such as described in Burch, Martin and Van Landingham [52] to predict moisture content as a function of relative humidity.
Combining these model terms with the log of total effective UV dosage from Equation (19) gives where temp K is temperature in kelvin, MC(RH) is a model prediction of moisture content, as a function of relative humidity, k is Boltzmann's constant, E a is a quasi-activation energy, and 0 β and C are parameters that are characteristic of the material and the degradation process.

Figure 1 .
The values of a and b are estimated respectively: a = 3.0868 and b = 5.762 × 10 −12 .
the parameter b of degradation model for a temperature of 358˚K (85˚C) and a relative humidity of 85%.The parameter b of the acceleration model is equal to 1.948 × 10 −9 for a temperature of 105˚C and a relative humidity of 85%, and equalto1.614× 10 −11 for conditions of 85˚C/95%HR.The constants a and b depend on the degradation mode considered.Thus, Equation (2) and (3) give respectively the degradation model due to corrosion and discoloration: a (a corrosion and a discoloration) and b (b corrosion and b discoloration) are the parameters of the degradation model.They are determined from accelerated testing Charki, and Laronde [19]-[21].

λ 0 T
λ is the spectral response of PV cell, ( ) respectively represent the transmit- tances of the encapsulant and the glass slide of PV module, ( ) λ is the transmittance of not irradiated PV cell.
interval for the wavelength range where the spectral response of PV cell is not zero.This model presents a major constraint.In fact, its use requires knowledge of the intrinsic characteristics of the materials used for the production of PV cells such as a cmx , b cmx ,

Figure 3 .
Figure 3. Equivalent temperature versus Ea (a) and equivalent relative humidity versus n (b) for the four cities studied.