An exergetic optimization is developed to determine the optimal performance and design parameters of a solar photovoltaic (PV) array. A detailed energy and exergy analysis is carried out to evaluate the electrical performance, exergy destruction components, and exergy efficiency of a typical PV array. The exergy efficiency of a PV array obtained in this paper is a function of climatic, operating, and design parameters such as ambient temperature, solar radiation intensity, PV array temperature, overall heat loss coefficient, open-circuit voltage, short-circuit current, maximum power point voltage, maximum power point current, and PV array area. A computer simulation program is also developed to estimate the electrical and operating parameters of a PV array. The results of numerical simulation are in good agreement with the experimental measurements noted in the previous literature. Finally, exergetic optimization has been carried out under given climatic, operating, and design parameters. The optimized values of the PV array temperature, the PV array area, and the maximum exergy efficiency have been found. Parametric studies have been also carried out.
1. Introduction
Renewable energies are going to be a main substitute for fossil fuels in the coming years for their clean and renewable nature. A solar photovoltaic (PV) array is one of the most significant and rapidly developing renewable-energy technologies, and its potential future uses are notable.
PV array is a semiconductor device, which converts light energy directly into useful electricity. The energy payback time (EPBT) of a PV system lies between 10 and 15 years depending on insulation and the performance of it. If the performance of a PV array can be increased, the energy payback time can be reduced. Therefore, the optimized performance evaluation of a PV array is important.
PV array performance parametrically depends on climatic, operating, and design parameters such as ambient temperature, solar radiation intensity, PV array temperature, overall heat loss coefficient, open-circuit voltage, short-circuit current, maximum power point voltage, maximum power point current, and PV array area. It can be evaluated in terms of energy efficiency and exergy efficiency. Its evaluation based on the first and second laws of thermodynamics is known as energy efficiency and exergy efficiency, respectively.
The energy analysis has some deficiencies [1, 2]; fundamentally, the energy concept is not sensitive with respect to the assumed direction of the process; for example, energy analysis does not object if heat is considered to be transferred spontaneously in the direction of increasing temperature. It also does not distinguish the quality of energy; for example, 1 W of heat equals 1 W of work or electricity. Energy analyses on their own incorrectly interpret some processes [1, 2]; for example, environmental air, when isothermally compressed, maintains its energy (e.g., enthalpy) equal to zero, whereas the exergy of the compressed air is greater than zero. However, exergy data are more practical and realistic in comparison with the respective energy values. Thus, exergy analysis usually provides a more realistic view of process than energy analysis; some times, they are different [1, 2].
The energy conversion factor of a solar PV system sometimes is described as efficiency, but this usage sometimes leads to some difficulties such as follows [3]. The energy efficiency of a PV array can be considered as the ratio of the electricity generated to the total, or global, solar irradiation. In this definition, only the electricity generated by a PV array is considered. The other components and properties of a PV array, such as ambient temperature, PV array, temperature, chemical potential components, and heat capacity of a PV array, are not directly taken into account [3]. However, the exergy efficiency of a PV array includes most of the climatic, geometric, and operating parameters of a PV array and involves the thermal properties and chemical potential components of it directly [3].
A little work has been carried out in field of PV array exergy analysis or exergetic optimization of it.
Ross and Hsiao [4] investigated some theoretical thermodynamic limits of photochemical solar energy conversion efficiency.
Landsberg and Markvart [5] studied the Carnot factor in PV cell theory and obtained an expression for open-circuit voltage, which is equal to the band-gap multiplied by the Carnot efficiency. Markvart and Landsberg [6] also discussed the thermodynamics and reciprocity of solar energy conversion by considering PV, photochemistry, and photosynthesis phenomena.
Würfel [7] studied thermodynamic limitations of solar energy conversion, based on an entropy concept, and calculated the upper efficiency as 0.86 for maximally concentrated solar radiation.
Smestad [8] examined the concepts of hot carrier and light converter, indicating that electrons are ejected not only as heat but also as light.
Bisquert et al. [9] presented some physical and chemical principles of PV conversion and found the relation between chemical potential and open-circuit voltage of PV cell to be dependent on Carnot and statistical factors.
Gong and Kulkarni [10] considered the optimization process of a grid connected PV system. They optimized array size and the tilt angle of array surface based on energy analysis.
Černivec et al. [11] optimized the optical and electrical performance of a CIGS solar cell from a perspective based on band-gap engineering. They examined structures with different CIGS absorber thicknesses and bulk and/or hetero-interface recombination lifetimes and discussed their optimized parameters in the light of experimental achievements.
Ghoneim [12] developed a computer simulation program in order to determine the thermal performance of a photovoltaic powered water pumping system in the Kuwait climate. He found the optimum values of PV array size, PV array orientation, and pump–motor–hydraulic system characteristics.
Badescu [13] studied the effect of latitude and climate on PV module shape based on energy analysis. He obtained the optimized value of PV cells number in series and strings number for various climatic conditions.
Sahin et al. [3] carried out the thermodynamic analysis of a PV array based on chemical potential components. They also obtained exergy components and PV array exergy efficiency. Finally, they compared energy, electrical, exergy efficiencies under given climatic and operating conditions. However, they did not fulfill the computer simulation and optimization of a PV array.
Skoplaki et al. [14] suggested a simple semiempirical correlation for PV array temperature. They investigated the effect of solar irradiance, ambient temperature, and wind speed on PV array temperature. However, they ignored free convection and radiation losses from PV array to the environment.
Abdolzadeh and Ameri [15] investigated the possibility of improving the performance of a photovoltaic water pumping system by spraying water over the front of PV array experimentally. They pointed out that the efficiency of PV array can be increased due to spraying water over the front of PV array.
Joshi et al. [16] developed a thermal model for the PV module integrated with solar air collector and validated it experimentally. They indicated that PV module temperature can be controlled and reduced in consequence of changing the mass flow rate of air in solar collector and the efficiency of PV module can be increased.
Joshi et al. [17] also studied the performance characteristics of a photovoltaic (PV) and photovoltaic-thermal (PV/T) system in terms of energy efficiency and exergy efficiency, respectively. They proposed equations for the energy, electrical, and exegry efficiency of a PV system. Finally, they calculated the energy, electrical, and exegry efficiencies under given experimental data and gave useful results. However, they did not obtain the exergy destruction components of PV array. They also did not carry out the numerical simulation and optimization of PV array.
Sarhaddi, Farahat and Ajam investigated exergetic optimization of solar collector systems [18–20].
In this paper, a detailed energy and exergy analysis will be carried out to evaluate the electrical performance, exergy destruction components, and exergy efficiency of PV array. An equation for PV array exergy efficiency will be derived based on exergy destruction. A computer simulation program will be developed to predict the electrical and operating parameters of PV array. Finally, the exergetic optimization of PV array will be carried out; also, the effect of climatic, design, and operating parameters on exergy efficiency will be studied.
PV array exergy analysis is parametrically dependent on its energy analysis. Hence, firstly PV array energy analysis will be carried out. Then the exergy destruction components and exergy efficiency of PV array will be computed and optimized.
2. Energy Analysis
The proof of governing equations on PV array energy analysis is not included in order to have a brief note.
A PV array is nonlinear device and can be represented by its current-voltage (I-V) characteristic curve. There are several mathematical models, which can describe I-V characteristic curve [21]. Five-parameter photovoltaic model (Figure 1) for I-V characteristic curve is defined as [21]
I=IL-Io[exp(V+IRsa)-1]-(V+IRs)Rsh,
where I and V represent current and voltage at load, and a, IL, Io, Rs, and Rsh are ideality factor, light current, diode reverse saturation current, series resistance, and shunt resistance, respectively.
Equivalent electrical circuit in the five-parameter photovoltaic model [21].
The second terms on the right-hand side of (1) indicate diode current (ID). In order to calculate five reference parameters (aref, IL,ref, Io,ref, Rs,ref, and Rsh,ref), five pieces of information are needed at reference conditions [21]. These five pieces of information are defined as follows [21].
At short circuit current:
I=Isc,ref,V=0.
At open circuit voltage:
I=0,V=Voc,ref.
At the maximum power point:
I=Imp,ref,V=Vmp,ref.
At the maximum power point:
[d(IV)dV]mp=0.
At short circuit:
[dIdV]sc=-1Rsh,ref.
Reference conditions or standard rated conditions (SRCs) are defined as follows [22].
The solar cell temperature at reference conditions is
Tcell,ref=25C∘.
The solar radiation intensity at reference conditions is
Gref=1000W/m2.
Substituting the above five pieces of information into (1), the following equations are obtained:
Isc,ref=IL,ref-Io,ref[exp(Isc,refRs,refaref)-1]-Isc,refRs,refRsh,ref,0=IL,ref-Io,ref[exp(Voc,refaref)-1]-Voc,refRsh,ref,Imp,ref=IL,ref-Io,ref[exp(Vmp,ref+Imp,refRs,refaref)-1]-(Vmp,ref+Imp,refRs,ref)Rsh,ref,[d(IV)dV]mp=0,[dIdV]sc=-1Rsh,ref,
where Voc, Vmp, Isc, and Imp are open-circuit voltage, maximum power point voltage, short-circuit current, and maximum power point current, respectively. The subscript “ref” indicates the value of parameters at the reference conditions.
Equations (9)–(13) are a set of nonlinear equations that can be solved with numerical methods. Solving (9)–(13) gives the value of five parameters (aref, IL,ref, Io,ref, Rs,ref, and Rsh,ref), at the reference conditions (Tcell,ref=25C∘, Gref=1000W/m2). In order to calculate the model parameters at new climatic and operating conditions (Gnew,Tcell,new), a set of translation equations is used as follows [14, 21, 23]:
Tcell=(Tamb+(GGref)(UL,NOCTUL)(Tcell,NOCT-Tamb,NOCT)×[1-ηel,ref(τα)(1+γrefTcell,ref)])/(1-γrefηel,ref(τα)(GGref)(UL,NOCTUL)×(Tcell,NOCT-Tamb,NOCT)(GGref)),aaref=TcellTcell,ref,IoIo,ref=(TcellTcell,ref)3exp(εNcaref(1-Tcell,refTcell)),IL=(GGref)[IL,ref+α(Tcell-Tcell,ref)],ΔT=Tcell-Tcell,ref,ΔI=α(GGref)ΔT+(GGref-1)Isc,ref,ΔV=βΔT-RsΔI,Inew=Iref+ΔI,Vnew=Vref+ΔV,
where Tcell, Tcell,ref,Tcell,NOCT, Tamb, Tamb,NOCT, G, and Gref are solar cell temperature, solar cell temperature at the reference conditions, nominal operating cell temperature, ambient temperature, ambient temperature at NOCT conditions, solar radiation intensity, and solar radiation intensity at the reference conditions, respectively. Further, parameters UL, UL,NOCT, ηel,ref, (τα), ε, Nc, α, β, and γref are overall heat loss coefficient, overall heat loss coefficient at NOCT conditions, electrical efficiency at the reference conditions, the effective product of transmittance–absorptance, semiconductor band gap energy (1.12 eV for silicon solar cell), cells number in series, current temperature coefficient, voltage temperature coefficient, and efficiency correction coefficient for temperature (0.004C∘-1 for silicon solar cell), respectively.
PV module manufacturers usually give temperature coefficients and NOCT conditions [22]. The new values of maximum power point voltage and maximum power point current are obtained from solved I-V characteristic curve and (12) simultaneously at new climatic and operating conditions.
In the previous studies [3–17], PV array overall loss coefficient (UL) has been assumed as a constant factor or a variable with little effect, whereas it is not constant. The overall heat loss coefficient of PV array includes convection and radiation losses:
UL=hconv+hrad.
The convective heat transfer coefficient is given by [24]
hconv=2.8+3Vw,
where Vw is wind speed. The radiative heat transfer coefficient between PV array and surroundings is obtained from [25]
hrad=εgσ(Tsky+Tcell)(Tsky2+Tcell2),
where εg and σ are PV array emissivity and the Stefan-Boltzmann’s constant, respectively, and the effective temperature of the sky (Tsky) is calculated from the following empirical relation [25]:
Tsky=Tamb-6.
The energy efficiency of a PV system can be defined as the ratio of the output energy of the system (i.e., electrical energy) to the input energy (i.e., solar energy) received on photovoltaic surface.
The maximum energy efficiency of a PV system is given by [3, 17]
ηen=VocIscS.
However, this definition of energy efficiency is restricted to theoretical cases. In (20), S is solar absorbed flux and it is given by
S=GAarr=G(NsNmAmod),
where Aarr, Nm, and Ns are PV array area, number of modules in series per string, and number of strings, respectively. PV module area (Amod) is given by
Amod=L1L2,
where L1 and L2 are the length of solar module and the width of solar module, respectively.
For PV systems in practical cases, energy efficiency measures the ability of converting solar energy into electrical energy [3, 17]. The electrical power output is the product of voltage and current of photovoltaic device. This conversion efficiency is not constant, even under constant solar irradiation. However, there is a point of maximum power, where voltage is Vmp, which is less than open-circuit voltage (Voc) but close to it, and current is Imp, which is less than short-circuit current (Isc) but close to it (Figure 2). In Figure 2, E GH stands for the highest energy level of electron at maximum solar irradiation conditions.
Representation of a general current-voltage characteristic curve and its parameters [17].
E GH is equivalent to area under the I-V characteristic curve (∫V=0VocI(V)dV). In addition, EL stands for the low-energy content of electron, which is the practical case, as shown by the rectangular area in the same figure. EL is thus equivalent to ImpVmp. The maximum power point is restricted by a term called ‘‘fill factor” as
FF=VmpImpVocIsc.
The energy efficiency of a PV system at maximum power is defined as the ratio of actual electrical output to input solar energy incident on PV surface area and it is given by [3, 16, 17]
ηel=VmpImpS.
This efficiency is also called actual electrical efficiency. The electrical efficiency of a PV array can also be defined in terms of fill factor (FF) as follows:
ηel=FF×VocIscS.
3. Exergy Analysis
Exergy analysis is a technique that uses the conservation of mass and conservation of energy principles together with the second law of thermodynamics for the analysis, design, and improvement of energy and other systems. Exergy is defined as the maximum amount of work that can be produced by a system or a flow of mass or energy as it comes to equilibrium with a reference environment [26]. Exergy efficiency is the ratio of total output exergy to total input exergy [26–29]:
ηex=ExoutExin=1-Ic·vExin,
where Exin, Exout, and Ic.v are inlet exergy, outlet exergy, and irreversibility in control volume, respectively.
The inlet exergy includes solar radiation intensity exergy. According to the Petela theorem, it is given by [1, 2]
Exin=S(1-43TambTsun+13(TambTsun)4),
where Tsun is the sun’s temperature in Kelvin.
The irreversibility in control volume includes external exergy losses from control volume and internal exergy losses (exergy destructions) in control volume [28–30]:
Ic·v=∑(Exlos+Exdes).
The external exergy losses caused by heat leakage are given by [18–20]
Exlos=ULAarr(Tcell-Tamb)(1-TambTcell).
The internal exergy losses (exergy destructions) include four terms; one is caused by optical losses in PV array surface [31, 32]:
Exdes,opt=S(1-TambTsun)(1-(τα));
the second term is caused by the temperature difference between PV array surface and the sun temperature [18–20, 31, 32]:
Exdes,ΔTsun=(τα)STamb(1Tcell-1Tsun);
the third term is caused by the temperature variation of PV array with respect to the reference environmental state [18–20]:
Exdes,ΔTarr=mcellCpTambΔt(ln(TcellTamb)-(Tcell-Tamb)Tcell),
where mcell and Δt are PV array mass and time interval, respectively. Time interval is chosen according to the time step of experiment course [3].
The specific heat capacity of silicon solar cell (Cp) is calculated from [33]
Cp=0.844+1.18×10-4Tcell-1.55×104Tcell-2.
The fourth term is electrical exergy destruction [17]:
Exdes,el=(IscVoc-ImpVmp).
The exergy efficiency of PV array: substituting (27)–(34) into (26), the exergy efficiency equation of PV array is obtained:
ηex=1-{(1-(τα))+ULAarr(Tcell-Tamb)(1-TambTcell)S(1-43TambTsun+13(TambTsun)4)+(τα)STamb(1Tcell-1Tsun)S(1-43TambTsun+13(TambTsun)4)+mcellCpTambΔt(ln(TcellTamb)-(Tcell-Tamb)Tcell)S(1-43TambTsun+13(TambTsun)4)+(IscVoc-ImpVmp)S(1-43TambTsun+13(TambTsun)4)}.
4. Formulation of Optimization Problem
The electrical and exergetic models presented in the previous sections have been inserted into a MATLAB computational program. In this program, most of the geometric parameters and operating conditions can be variables. The formulation of optimization problem, considering the quantities Tamb, Tamb,NOCT, Tcell,ref, Tcell,NOCT,Tsun, G, Gref, Δt, α, β, γref, ηel,ref, ε, (τα), L1, L2, Amod, Voc,ref, Vmp,ref, Isc,ref, Imp,ref, Vw, and so forth, as constant parameters is given by
Maximizeηex=(35),subjectto(1,9-19)and(21),300≤Tcell≤350K,Amod≤Aarr≤16Amodm2,mcell,Cp,Io,IL,Rs,Rsh,a,Isc,Voc,Imp,Vmp,S,UL,ηen,ηel,FF≥0,Ns,Nm,Nc≥1andareinteger,
where Tcell and Aarr are independent parameters and mcell, Cp, Io, IL, Rs, Rsh, a, Isc, Voc, Imp, Vmp, S, UL, ηen, ηel, FF, Ns, Nm, and Nc are dependent parameters in optimization procedure. Whereas the objective function and its constraints are nonlinear, it has been maximized numerically with MATLAB optimization toolbox. MATLAB uses the variations of sequential quadratic programming (SQP) method for constrained optimization [34].
5. Results and Discussion5.1. Experimental Validation
The experimental results of Barker and Norton [35] for a rack-mounted PV array allow the checking of results obtained by computer simulation. The measured data in [35] include the solar radiation intensity, ambient temperature, PV array temperature, open-circuit voltage, maximum power point voltage, short-circuit current, and maximum power point current. Additional information about the experiment method and its conditions is found in [35]. The experimental values of the above-mentioned parameters have been obtained from Figure 8 and Figure 9 of [35].
The simulated values of PV array temperature, open-circuit voltage, maximum power point voltage, short-circuit current, and maximum power point current in present work have been validated by their corresponding experimental values in [35].
Climatic, operating, and design parameters of PV array during validation process are described in Table 1. They correspond to the experimental system described by Barker and Norton [35], except that they did not report the wind speed observed over the course of their tests. This affects the convective heat transfer coefficient between PV array surface and the ambient air. In Table 1, a wind speed of 0.5 m/s is assumed to have a comparison with the experimental data.
Climatic, operating and design parameters of PV array [22, 35].
Solar PV module parameters
Value in validation
Value in optimization
PV module type
Siemens SM55,
Siemens SM55,
monocrystalline silicon
monocrystalline silicon
Number of modules in series per string, Nm
2
From optimization
Number of strings, Ns
6
From optimization
The solar radiation intensity at the reference conditions, Gref
1000W/m2
1000W/m2
The solar radiation intensity, G
From experimental data
700W/m2
The ambient temperature, Tamb
From experimental data
300 K
The ambient temperature at NOCT conditions, Tamb,NOCT
293.15 K
293.15 K
The cell temperature at the reference conditions, Tcell,ref
298.15 K
298.15 K
The nominal operating cell temperature, Tcell,NOCT
318.15 K
318.15 K
The PV array temperature, Tcell
(14)
From optimization
The sun temperature, Tsun
5760 K
5760 K
Wind speed, Vw
0.5 m/s
0.5 m/s
The short-circuit current at the reference conditions, Isc,ref
Ns×3.45A
Ns×3.45A
(for total array)
(for total array)
The open-circuit voltage at the reference conditions, Voc,ref
Nm×21.7V
Nm×21.7V
(for total array)
(for total array)
Maximum power point current at the reference conditions, Imp,ref
Ns×3.15A
Ns×3.15A
(for total array)
(for total array)
Maximum power point voltage at the reference conditions, Vmp,ref
Nm×17.4V
Nm×17.4V
(for total array)
(for total array)
The electrical efficiency at the reference conditions, ηel,ref
0.12
0.12
The current temperature coefficient, α
1.2mA/C∘
1.2mA/C∘
The voltage temperature coefficient, β
-0.077V/C∘
-0.077V/C∘
The efficiency correction coefficient for temperature, γref
0.004C∘-1
0.004C∘-1
The semiconductor band gap energy, ε
1.12 eV
1.12 eV
The effective product of transmittance–absorptance, (τα)
0.9
0.9
The PV array emissivity, εg
0.88
0.88
The length of solar module, L1
1.293m
1.293m
The width of solar module, L2
0.329m
0.329m
Time interval, Δt
900 seconds
3600 seconds
In order to compare the simulated results with the experimental measurements, a root mean square percentage deviation (RMS) has been evaluated by following equation [16, 21]:
RMS=∑[100×(Xexp,i-Xsim,i)/Xexp,i]2n,
where n is the number of record of experimental data over the course of experiment.
The variations of solar radiation intensity, ambient temperature, and PV array temperature during the test day are shown in Figure 3. The simulated values of PV array temperature are also shown in the same figure for comparison. It can be seen from this figure that there is a good agreement between the simulated and experimental values of PV array temperature with standard percentage deviation (RMS) = 4.82%.
The variations of solar radiation intensity, ambient temperature, experimental PV array temperature, and simulated PV array temperature during the test day.
The simulated values of open-circuit voltage, maximum power point voltage, short-circuit current, maximum power point current, and the corresponding experimentally measured data during the test day are shown in Figure 4. It can be observed from this figure that there is a good agreement between the experimental and simulated values of these parameters. Further, the root mean square percentage deviations of these parameters are 2.4%, 2.28%, 4.6%, and 2.98%, respectively.
The simulated values of open-circuit voltage, maximum power point voltage, short-circuit current, maximum power point current, and the corresponding experimentally measured data during the test day.
The simulated and experimental values of energy efficiency, exergy efficiency, and electrical efficiency during the test day are shown in Figure 5. The respective values of root mean square percent deviation are 5.88%, 3.49%, and 3.58%, respectively. It can be seen from Figure 5 that there is a good agreement between the simulated and experimental values of the above-mentioned efficiencies.
The simulated and experimental values of energy efficiency, exergy efficiency, and electrical efficiency during the test day.
The good agreement between experiment and simulation values shown in the previous figures demonstrates that the choice of a wind speed of 0.5 m/s during our calculations is reasonable (Figures 3–5).
The simulated parameters errors compared with those obtained by the experimental measurement are explained as follows.
The temperature coefficients of current and voltage are assumed constant. In practical cases, they have slight fluctuation due to the solar radiation intensity and PV array temperature variations.
The experimental measurements have been obtained from the figures of [35] by interpolation and curve fitting methods. This subject decreases the precision of measured data.
Wind speed is not constant and has a direct effect on the overall heat loss coefficient that can decrease the precision of calculated overall heat loss coefficient in the computer simulation.
The effective product of transmittance–absorptance is assumed constant while it is changing during the day with the change of solar incidence angle on PV array surface.
5.2. Optimization Results
The type of PV array, its selected environmental, design conditions, and constant parameters during the optimization procedure are described in Table 1. For the fixed values in Table 1, the optimization procedure gives the following optimized values of the objective function, the independent and dependent optimization parameters:
5.3. Parametric Studies. The values of climatic, operating, and design parameters in the parametric studies of PV array are described in Table 1, except that the parameters are changed during parametric studies that are indicated by the above of respective figures.
Figure 6 shows the behavior of the exergy efficiency as a function of the PV array temperature (Tcell) and the PV array area (Aarr); it is presented to a range of operational and design conditions where the exergy efficiency takes a global maximum value. The incremented quantities are the PV array area from 0.4254 to 6.8064 m2 and the PV array temperature from 300 to 350 K. It is observed from this figure that there is a global maximum point and the coordinate of this point shows the values of optimized parameters. The calculated values of global maximum point are Tcell,opt=301.12K, Aarr,opt=2.5524m2, and ηex,max=13.33%. According to this figure, the exergy efficiency has a slight variations (~13.3%) with respect to the PV array area changes (0.4254≤Aarr≤6.8064m2) in optimum value of the PV array temperature (Tcell,opt=301.12K). However, there is an optimized point for the PV array area that is equal with Aarr,opt=2.5524m2. On the other hand, by increasing the PV array temperature, the exergy efficiency is constant initially and then decreases quickly from ~13.3% to ~2% after the PV array temperature reaches about 302 K. This indicated the optimum range of the PV array temperature (Tamb≤Tcell≤302K) for given climatic and design parameters (Table 1).
The variations of the exergy efficiency as a function of the PV array temperature and the PV array area.
In order to have maximum exergy efficiency, PV array temperature should be kept near the ambient temperature, or in other words, PV array temperature should be controlled. In order to control PV array temperature, there are some practical methods such as spraying water on the top surface of photovoltaic modules [15] or combining of PV modules in photovoltaic/thermal (PV/T) collectors [16].
Figure 7 shows the variations of exergy efficiency with respect to solar radiation intensity. It is observed that the exergy efficiency increases from ~12.61% to ~14.68% initially, and then after attaining the solar radiation intensity of about G ≈80 W⁄m2, it decreases to ~12.95%. This indicates the optimum value of solar radiation intensity for given climatic and design parameters (Table 1). Since this parameter varies during the day, the design of PV array should be based on the daily or monthly average of this parameter.
The variations of the exergy efficiency with respect to the solar radiation intensity.
The effect of ambient temperature on the exergy efficiency.
Figure 8 shows the effect of ambient temperature on the exergy efficiency. Increasing the ambient temperature from 300 to 320 K determines the sensible decrease of exergy efficiency from ~13.33% to ~6.5%. Since the ambient temperature changes during the day, PV array design should be based on the daily average of this parameter.
The average value of a function f(t) on a time interval [0, M] is defined as follows [36]:
f̅=1M∫0Mf(t)dt=∑0Mf(t)M,
where M is time in hour.
6. Conclusion
On the basis of theoretical results obtained in the present paper, the following conclusions have been drawn.
The electrical and exergetic models of PV array presented in this study are in good agreement with the experimental results of Barker and Norton [35] (Figures 3–5).
The PV array temperature has a great effect on the exergy efficiency. The exergy efficiency of a PV array can be improved if the heat can be removed from the PV array surface. In order to remove heat from the PV array surface, there are some practical methods such as spraying water on the top surface of photovoltaic modules [15] or combining of PV modules in photovoltaic/thermal (PV/T) collectors [16] (Figure 6).
Increasing the solar radiation intensity, the exergy efficiency of PV array increases initially and then decreases after attaining the solar radiation intensity of about a maximum point (Figure 7).
While the ambient temperature is increasing, the exergy efficiency of PV array decreases (Figure 8).
The design parameters such as PV array area have a little effect on the exergy efficiency (Figure 6).
Nomenclaturea:
Ideality factor (eV)
A:
Area (m2)
Cp:
Heat capacity of the silicon material (J/g·K)
CIGS:
Copper Indium Gallium Selenide
E:
Electrical power (W)
EPBT:
Energy Payback Time
Ex:
Exergy (W)
FF:
Fill Factor
G:
Solar radiation intensity (W/m2)
h:
Heat transfer coefficient (W/m2·K)
I:
Circuit current (A)
Ic·v:
Irreversibility in control volume (J)
I-V:
Current-voltage
L:
Dimensions of solar module (m)
m:
PV array mass (g)
M:
Time (hr)
n:
Number of parameters
N c:
Number of cells in PV module
N m:
Number of modules in series per string
N s:
Number of string
NOCT:
Nominal operating Cell Temperature conditions
PV:
Photovoltaic
PV/T:
Photovoltaic/thermal collector
R:
Resistance (Ω)
RMS:
Root Mean Square percentage deviation (%)
S:
Solar absorbed flux (W)
SQP:
Sequential Quadratic Programming
SRC:
Standard Rated Conditions
t:
Time interval (second)
T:
Temperature (K)
UL:
Overall heat loss coefficient from the PV array to the environment (W/m2·K)
V:
Circuit voltage (V), wind speed (m/s)
X:
Experimental or simulated value of parameter.
Greek Symbolsα:
Current temperature coefficient (mA/°C)
β:
Voltage temperature coefficient (V/°C)
γ:
Efficiency correction coefficient for temperature (°-1)
Δ:
Difference in current, temperature, time, voltage
ε:
Global emissivity, semiconductor band gap energy (eV)
η:
Efficiency (%)
σ:
Stefan-Boltzmann’s constant (W/m2K4)
(τα):
The effective product of transmittance-absorptance.
Subscripts1:
Length
2:
Width
arr:
Array
amb:
Ambient
cell:
Cell, array
conv:
Convection
des:
Destroyed
D:
Diode
el:
Electrical
en:
Energy
ex:
Exergy
exp:
Experimental
g:
Glass
GH:
The highest energy content of the electron
i:
Ith parameter
in:
Inlet
L:
The available energy content of the electron, light current
los:
Loss
max:
Maximum
mod:
Module
mp:
Maximum power point
new:
New
NOCT:
At NOCT conditions
o:
Reverse saturation
oc:
Open-circuit
opt:
Optical, optimum
out:
Outlet
rad:
Radiative
ref:
Reference
s:
Series
sc:
Short-circuit
sh:
Shunt
sim:
Simulated
sky:
Sky
sun:
Sun
w:
Wind.
Acknowledgments
The authors acknowledge Professor. A. D. Sahin [3] for responding then inquiries and acknowledge the constructive comments of Dr. M. Abdolzadeh [15].
PetelaR.An approach to the exergy analysis of photosynthesis200882431132810.1016/j.solener.2007.09.0022-s2.0-41249086339PetelaR.Exergy of undiluted thermal radiation200374646948810.1016/S0038-092X(03)00226-32-s2.0-0141794049SahinA. D.DincerI.RosenM. A.Thermodynamic analysis of solar photovoltaic cell systems2007912-315315910.1016/j.solmat.2006.07.0152-s2.0-33750704977RossR. T.HsiaoT.-L.Limits on the yield of photochemical solar energy conversion197748114783478510.1063/1.3234942-s2.0-0017552270LandsbergP. T.MarkvartT.The carnot factor in solar-cell theory19984246576592-s2.0-0032042722MarkvartT.LandsbergP. T.Thermodynamics and reciprocity of solar energy conversion2002141-2717710.1016/S1386-9477(02)00352-12-s2.0-0036529301WürfelP.Thermodynamic limitations to solar energy conversion2002141-2182610.1016/S1386-9477(02)00355-72-s2.0-0036529311SmestadG. P.Conversion of heat and light simultaneously using a vacuum photodiode and the thermionic and photoelectric effects2004821-222724010.1016/j.solmat.2004.01.0202-s2.0-1942505294BisquertJ.CahenD.HodesG.RuhleS.ZabanA.Physical chemical principles of photovoltaic conversion with nanoparticulate, mesoporous dye-sensitized solar cells2004108248106811810.1021/jp03592832-s2.0-3042598977GongX.KulkarniM.Design optimization of a large scale rooftop photovoltaic system200578336237410.1016/j.solener.2004.08.0082-s2.0-13844264281ČernivecG.KrčJ.SmoleF.TopičM.Band-gap engineering in CIGS solar cells using Nelder-Mead simplex optimization algorithm2006511-512606510.1016/j.tsf.2005.11.0952-s2.0-33646590274GhoneimA. A.Design optimization of photovoltaic powered water pumping systems20064711-121449146310.1016/j.enconman.2005.08.0152-s2.0-33644853623BadescuV.Simple optimization procedure for silicon-based solar cell interconnection in a series-parallel PV module2006479-101146115810.1016/j.enconman.2005.06.0182-s2.0-32344448257SkoplakiE.BoudouvisA. G.PalyvosJ. A.A simple correlation for the operating temperature of photovoltaic modules of arbitrary mounting200892111393140210.1016/j.solmat.2008.05.0162-s2.0-50249096178AbdolzadehM.AmeriM.Improving the effectiveness of a photovoltaic water pumping system by spraying water over the front of photovoltaic cells2009341919610.1016/j.renene.2008.03.0242-s2.0-51649118012JoshiA. S.TiwariA.TiwariG. N.DincerI.ReddyB. V.Performance evaluation of a hybrid photovoltaic thermal (PV/T) (glass-to-glass) system200948115416410.1016/j.ijthermalsci.2008.05.0012-s2.0-54849411033JoshiA. S.DincerI.ReddyB. V.Thermodynamic assessment of photovoltaic systems20098381139114910.1016/j.solener.2009.01.0112-s2.0-67649421312SarhaddiF.FarahatS.AjamH.SobhnamayanF.Thermodynamic optimization of solar parabolic cookers and comparison with energy analysisProceedings of the 5th International Chemical Engineering Congress (IChEC '08)2008Kish Island, IranFarahatS.SarhaddiF.AjamH.Exergetic optimization of flat plate solar collectors20093441169117410.1016/j.renene.2008.06.0142-s2.0-56049100927AjamH.FarahatS.SarhaddiF.Exergetic optimization of solar air heaters and comparison with energy analysis2005841831902-s2.0-32944460097De SotoW.2004chapter 2-3Madison, Wis, USASolar Energy Laboratory, University of Wisconsin—MadisonSiemens CompanySiemens solar module SM552009http://www.elk.itu.edu.tr/~gunes/sm55.pdfLuqueA.HegedusS.2003Chichester, UKJohn Wiley & SonsWatmuffJ. H.ChartersW. W. S.ProctorD.Solar and wind induced external coefficients for solar collectors1977256SukhatmeS. P.1993New York, NY, USAMcGraw-HillKotasT. J.1995Malabar, Fla, USAKriegerWongK. F. V.2000Boca Raton, Fla, USAUniversity of Miami, CRC Press LLCBejanA.1998Chichester, UKJohn Wiley & SonsFeidtM.1996Paris, FranceTech & DocBejanA.1982Chichester, UKJohn Wiley & SonsAltfeldK.LeinerW.FiebigM.Second law optimization of flat-plate solar air heaters—part I: the concept of net exergy flow and the modeling of solar air heaters198841212713210.1016/0038-092X(88)90128-42-s2.0-0023834847AltfeldK.LeinerW.FiebigM.Second law optimization of flat-plate solar air heaters—part 2: results of optimization and analysis of sensibility to variations of operating conditions198841430931710.1016/0038-092X(88)90026-62-s2.0-0023833676RegelA. R.GlazovV. M.1980Moscow, UKNaukaThe MathWorks Inc.Optimization toolbox user's guide version 22009http://www.mathworks.com/BarkerG.NortonP.Building America system performance test practices: part 1—photovoltaic systems2003NREL/TP-550-30301Golden, Colo, USANational Renewable Energy Laboratory (NREL)http://www.nrel.gov/docs/fy03osti/30301.pdfHardyG. H.LittlewoodJ. E.PólyaG.19882ndCambridge, UKCambridge University Press