The goal of this paper is to introduce a model to predict the maximum power of a high concentrator photovoltaic module. The model is based on simple mathematical expressions and atmospheric parameters. The maximum power of a HCPV module is estimated as a function of direct normal irradiance, cell temperature, and two spectral corrections based on air mass and aerosol optical depth. In order to check the quality of the model, a HCPV module was measured during one year at a wide range of operating conditions. The new proposed model shows an adequate match between actual and estimated data with a root mean square error (RMSE) of 2.67%, a mean absolute error (MAE) of 4.23 W, a mean bias error (MBE) of around 0%, and a determination coefficient (R2) of 0.99.
1. Introduction
High concentrator photovoltaic (HCPV) technology is based on the use of optical devices that focus the sunrays on a high efficiency solar cell with the aim of reducing the cell surface area and therefore to decrease the costs of electricity [1, 2]. A typical HCPV module is made up of multijunction (MJ) solar cells, interconnected in series or parallel, with one optical device per cell, and with a concentration factor higher than 500x. The optical devices usually consist of a primary optical reflexive or refractive element that collects and concentrates the light and a secondary optical element to homogenize the light on the solar cell surface and improve the overall acceptance angle of the module [3–7]. The efficiencies already reached by MJ solar cells and HCPV modules are noticeable and are continuously being increased [8–12], because this HCPV represents a potential alternative technology to conventional flat photovoltaic systems in the energy generation market with an installed capacity that could reach more than 1 GWp in 2020 [13].
The energy prediction and the modelling of the power output of HCPV devices are important issues to evaluate the potential and promote the market expansion of this emerging technology. However, the modelling and electrical characterization of these devices is complex due to the use of multijunction solar cells and optical elements. In order to address this issue, the scientific community has devoted considerable efforts in developing models that reproduce the electrical performance of HCPV modules [16–19]. These models present different levels of complexity and accuracy and usually need measurements of specific instruments, detailed information of the materials of the modules which is not always available, and advanced knowledge of semiconductor physics, optics, or different specific software.
A different approach for the modelling of HCPV modules based on atmospheric parameters and on a simple mathematical expression has been introduced by Fernández et al. [14]. This model predicts the maximum power of a HCPV module as a function of direct normal irradiance and temperature and introduces a spectral correction based on air mass. The air mass is widely accepted as the main parameter that affects the spectral distribution and the electrical output of HCPV modules. However, there are other parameters that play an important role in the spectral distribution of the incident irradiance and power output of HCPV modules. Among these, aerosol optical depth has been demonstrated to be a key parameter since it has a nonnegligible effect in the spectral characterization of these devices [20–25]. Bearing this in mind, the aim of this paper is to introduce an additional spectral correction based on aerosol optical depth in the mathematical model presented by Fernández et al. [14]. The approach is to find a model based on atmospheric parameters and simple mathematical relationships with an optimum balance between accuracy and simplicity in order to facilitate its application.
The paper is organized as follows. In Section 2, the mathematical model is introduced and justified. In Section 3, the experimental campaign to check the quality of the proposed model is described. In Section 4, the analysis of the results in the estimation of the output of a HCPV module is conducted. Finally, the main conclusions are outlined in Section 5.
2. Mathematical Modelling
HCPV modules only react to the direct component of the incident irradiance due to the use of optical devices. Taking this into account, the maximum power of a HCPV module can be estimated with the well-known relations used in conventional photovoltaic technology introduced by Osterwald [26] and Evans [27] as a function of direct normal irradiance (DNI) and cell temperature (Tc) as
(1)P=P*DNI*DNI1-δTc-Tc*,
where P*, DNI*, and Tc* are maximum power, direct normal irradiance, and cell temperature at standard test conditions; and δ is the cell temperature coefficient of maximum power.
HCPV modules are based on MJ solar cells and optical devices. In MJ solar cells, each junction responds to a particular band of the spectrum. At the same time, the optical devices modify the spectral distribution that strikes the solar cell surface. Because of this, HCPV modules are strongly affected by the incident spectrum, so that is necessary to introduce a spectral correction in (1). As commented, the air mass (AM) has been identified as the parameter with the largest impact on the spectral distribution and power output of HCPV modules. Because of this, (1) has been rewritten with a spectral correction based on AM by Fernández et al. [14, 28] as
(2)Pmodule=P*DNI*DNI1-δTc-Tc*×1-εAM-AMU,
where ε is the air mass coefficient of the maximum power of a HCPV module and AMU is defined as the umbral air mass. The analysis of performance of HCPV modules in outdoors has indicated that the power output of a HCPV module can be considered independent of the air mass for AM ≤ AMU and can be corrected with a linear coefficient for AM > AMU. This behaviour is shown in Figure 1 where the normalized power of a HCPV module is plotted versus AM (grey dots). For low AM values the top and middle junctions of MJ solar cells of a HCPV module almost generate the same current and therefore the power output can be considered independent of AM. However, as the AM increases the top junction strongly limits the current of the HCPV module reducing its performance. This behaviour is also shown in Figure 1 where the spectral factor (SF) of a HCPV module is plotted versus AM (red line). This factor quantifies the differential performance of a HCPV device between the incident and reference spectra as a function of the mismatch among the different junctions of the respective MJ solar cell: SF > 1 indicates spectral gains while SF < 1 indicates spectral losses [15, 29, 30]. It is important to note that the SF reaches a maximum at a certain air mass value, where the top and middle junctions generate the same current, and decreases aside this point due to limiting current by one of the junctions. On the other hand, the power output keeps almost constant and starts decreasing when the SF reaches a certain value. This can be explained since the SF quantifies the performance of a HCPV module through the mismatch among the current of the junctions of MJ solar cells and taking into account that the maximum power shows a lower spectral dependence than the current [31–33].
Normalized maximum power (grey dots) and spectral factor (red line) versus AM of a HCPV module made up of lattice-matched GaInP/GaInAs/Ge solar cells and poly(methylmethacrylate) (PMMA) Fresnel lenses [14, 15].
Equation (2) quantifies the spectral effects by the use of the AM. However, the performance of HCPV modules is also affected by aerosol optical depth (AOD). Therefore, it is adequate to introduce an additional spectral correction based on AOD in order to improve the quality of the model. In order to introduce the appropriate AOD correction, the SF of a HCPV module has been plotted versus aerosol optical depth at 550 nm (AOD550), as shown in Figure 2. As can be seen, the same behaviour to that in the previous case is observed. For low AOD550 values, the top and middle junctions of MJ solar cells of a HCPV module almost generate the same current and as the AOD550 increases the top junction limits the current of the HCPV module reducing its performance with a linear behaviour as a function of AOD550. Taking this into account, (2) can be rewritten as
(3)Pmodule=P*DNI*DNI1-δTc-Tc*×1-εAM-AMU×1-φAOD550-AOD550,U,
where φ is the aerosol optical depth coefficient of the maximum power of a HCPV module and AOD550,U is the umbral aerosol optical depth at 550 nm. This equation quantifies the AOD at 550 nm since it is the most common observed value in different aerosol data bases.
Spectral factor versus AOD550 of a HCPV module made up of lattice-matched GaInP/GaInAs/Ge solar cells and poly(methylmethacrylate) (PMMA) Fresnel lenses [15].
Equations (1), (2), and (3) need cell temperature in order to estimate the maximum power of a HCPV module. This temperature is calculated from atmospheric parameters using the equation introduced by Almonacid et al. [34] as
(4)Tcell=Tair+aDNI+bWs,
where Tair is the air temperature, Ws is the wind speed, and a and b are coefficients specific to each HCPV module. This equation has been widely analysed and has showed a good performance in the estimation of cell temperature of a HCPV module [35].
3. Experimental Campaign
In order to check the quality of the model described in the last section, a HCPV module was under study at the Centre of Advanced Studies in Energy and Environment (CEAEMA) at the University of Jaen in Southern Spain (N 37°27′36′′, W 03°28′12′′) from July 2013 to June 2014. The HCPV module was mounted on a high accuracy two-axis solar tracker located on the roof of the research centre (Figure 3(a)). Tables 1 and 2 show the main characteristics and maximum power at standard test conditions (STC) of the HCPV module. In order to measure the electrical parameters of the HCPV module a four-wire electronic load was used. Also, two four-wire PT100 were installed. One thermometer was placed close to the solar cell to measure its temperature (Figure 3(c)) and the other was placed on the back of the module to measure its heat-sink temperature. These thermometers were connected to a data logger to record both temperatures. It is important to note that the sensors were located in a receiver between the centre and the border of the module, so that the measured temperatures should be considered as the average temperature of a receiver due to the temperature distribution of a HCPV module [36]. This approach has been previously used and is considered as an adequate procedure for estimating the cell temperature of a HCPV module and for its electrical characterization [37, 38]. In addition, an atmospheric station (Figure 3(b)) recorded the main atmospheric parameters such as direct normal irradiance, air temperature, wind speed, or humidity. All the parameters were recorded daily every 5 minutes. Figure 3(d) shows the scheme of the experimental set-up. In addition, the daily average values of aerosol optical depth at 550 nm not provided by the atmospheric station were obtained from MODIS Daily Level-3 data source [39].
Characteristics of the high concentrator photovoltaic module used in the study. Every cell is protected with a bypass diode.
Geometric concentration
700
Primary optics
SOG squared flat Fresnel lens
Secondary optics
Reflexive truncated pyramid
Optical efficiency
0.80
Type of solar cells
Lattice-matched GaInP/GaInAs/Ge
Solar cells area
0.763 cm2
Number of solar cells
20
Cooling
Passive
Maximum power of the high concentrator photovoltaic module considered at standard test conditions (STC).
P (W)
DNI (W/m2)
Tcell (°C)
Spectrum
280
1000
25
AM1.5d
Experimental set-up used to carry out the study at the Centre of Advanced Studies in Energy and Environment of the University of Jaen.
Figure 4 shows the distribution of the parameters needed to apply the models described in Section 2, obtained during the experiment. Based on these parameters and performing different regression analysis, the coefficients of the models for the module considered have been obtained, as shown in Table 3.
Values of the coefficients obtained from outdoor monitored data for estimating the maximum power of the HCPV module considered.
Coefficient
Value
a (°C/Wm−2)
0.044
b (°C/ms−1)
−3.41
δ (1/ °C)
0.0016
ε (dimensionless)
0.041
AMU (dimensionless)
2.10
φ (dimensionless)
0.32
AOD550,U (dimensionless)
0.25
Frequency of the direct normal irradiance (a), the cell temperature (b), the air mass (c), and the aerosol optical depth at 550 nm (d) measured during the experimental campaign.
4. Analysis of Results
In order to check the quality of the proposed model different analyses have been conducted. Figure 5 shows the linear regression analysis between actual and predicted data using (3) and (4) and the coefficients empirically obtained for the module considered; see Table 3. As can be seen, the model presents an almost perfect fit with a slope equal to 1 and determination coefficient (R2) of 0.99. Also, different statistical parameters have been calculated: the root mean square error (RMSE), the mean absolute error (MAE), and the mean bias error (MBE). The model shows a RMSE = 2.67% and a MAE = 4.23 W which indicates the small variation of predicted data around measured data and taking into account that the HCPV module considered has a maximum power of 280 W at STC. The model shows MBE = 0.19% which neither overestimates nor underestimates the maximum power.
Linear regression analysis between measured and predicted data for the proposed model for the HCPV considered.
As commented, the output of a HCPV module is mainly given by direct normal irradiance, cell temperature, and incident spectrum. In order to evaluate the behaviour of the model versus these parameters, Figure 6 shows the RMSE for different DNI, Tc, and SF levels. The spectral factor is obtained using the procedures previously described by different authors [40–44]. Figure 6(a) shows the RMSE versus DNI. As can be seen, the model shows a good performance for all DNI levels with a RMSE almost constant ranging from 1% to 2%, although it trends to increase with DNI until a maximum value at around 3%. Figure 6(b) shows the RMSE versus Tc. The model also shows a good performance for all Tc values with a RMSE ranging from around 2.5% to 0.8%. As can also be seen, the model yields better results for high Tc values. Figure 6(c) shows the RMSE versus SF. As can be seen, the model yields good results at different spectra with a RMSE ranging from around 1% to 2.8%. Also, it can be concluded that the model yields poorest results as the SF increases.
Root mean square error (RMSE) for different DNI, Tc, and SF levels for the proposed model for the HCPV considered.
In addition to the analyses conducted, a comparison between the results obtained with the proposed model and the results obtained using models based on (1) and (2) has been carried out, as shown in Table 4. As can be seen, the proposed model significantly improves the results of the model based on (1). This can be explained since this model does not take into account any spectral correction. Regarding the comparison with the model based on (2), the proposed model also improves the results which demonstrates the better performance of the proposed model based on an additional spectral correction based on AOD550. It is important to note that Jaen is a nonindustrialized medium-size city with low-medium values of aerosol optical depth as shown in Figure 4(d). Because of this, although the proposed model has a better performance, the results obtained with the model based on (2) and on (3) are similar. However, the proposed model is expected to significantly improve the results of the model only based on AM correction at locations with high turbidity where the AOD plays an important role on the performance of HCPV modules [15].
Root mean square error (RMSE), mean absolute error (MAE), mean bias error (MBE), and determination coefficient (R2) obtained for the three models described in Section 2.
Model
RMSE (%)
MAE (W)
MBE (%)
R2
DNI, Tcell
4.96
7.46
0.42
0.92
DNI, Tcell, AM
3.53
5.58
−0.32
0.98
DNI, Tcell, AM, AOD
2.67
4.23
0.19
0.99
5. Conclusions
In this paper a mathematical model for estimating the maximum power of a HCPV module has been introduced. Following the approach presented by Fernández et al. [14], the proposed model is based on atmospheric parameters and simple mathematical equations in order to facilitate its application. The model introduces a new spectral correction based on aerosol optical depth, so that is valid for sites with high turbidity. The analysis of the results shows that the model has a good performance in the prediction of the maximum power of a HCPV module with a RMSE = 2.67%, a MAE of 4.23 W, a MBE of around 0%, and R2 of 0.99. Furthermore, the model improves the results obtained with the model previously introduced by the authors only based on an air mass correction. So the model can be considered as a new useful tool for the prediction of the electrical output of HCPV modules under diverse climate conditions.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
Partial funding for this study was provided through EPSRC funded BioCPV (EP/J000345/1) project. Also, this work is part of the project “Desenvolvemento de novos conceptos baseados en tecnoloxía de concentración fotovoltaica para a produción de enerxía eléctrica adaptados a distintas zonas climáticas,” through the program “Formación Posdoutoral do Plan Galego de Investigación, Innovación e Crecemento 2011–2015 (Plan I2C)” funded by the Xunta de Galicia and by the European Social Fund.
LuqueA.SalaG.Luque-HerediaI.Photovoltaic concentration at the onset of its commercial deployment200614541342810.1002/pip.7052-s2.0-33746080913KhamooshiM.SalatiH.EgeliogluF.Hooshyar FaghiriA.TarabishiJ.BabadiS.A review of solar photovoltaic concentrators201420141795852110.1155/2014/958521CotalH.FetzerC.BoisvertJ.KinseyG.KingR.HebertP.YoonH.KaramN.III-V multijunction solar cells for concentrating photovoltaics20092217419210.1039/b809257e2-s2.0-69249096239ZubiG.Bernal-AgustínJ. L.FracastoroG. V.High concentration photovoltaic systems applying III-V cells20091392645265210.1016/j.rser.2009.07.0022-s2.0-68749091091XieW. T.DaiY. J.WangR. Z.SumathyK.Concentrated solar energy applications using Fresnel lenses: a review20111562588260610.1016/j.rser.2011.03.0312-s2.0-79955490293VictoriaM.DomínguezC.AntónI.SalaG.Comparative analysis of different secondary optical elements for aspheric primary lenses20091796487649210.1364/oe.17.0064872-s2.0-65149095633BaigH.HeasmanK. C.MallickT. K.Non-uniform illumination in concentrating solar cells20121685890590910.1016/j.rser.2012.06.0202-s2.0-84864795481YamaguchiM.TakamotoT.ArakiK.Ekins-DaukesN.Multi-junction III-V solar cells: current status and future potential2005791788510.1016/j.solener.2004.09.0182-s2.0-20444490856Pérez-HiguerasP.MuñozE.AlmonacidG.VidalP. G.High concentrator PhotoVoltaics efficiencies: present status and forecast20111541810181510.1016/j.rser.2010.11.0462-s2.0-79851511844DimrothF.GraveM.BeutelP.FiedelerU.KarcherC.TibbitsT. N. D.OlivaE.SieferG.SchachtnerM.WekkeliA.BettA. W.KrauseR.PiccinM.BlancN.DrazekC.GuiotE.GhyselenB.SalvetatT.TauzinA.SignamarcheixT.DobrichA.HannappelT.SchwarzburgK.Wafer bonded four-junction GaInP/GaAs//GaInAsP/GaInAs concentrator solar cells with 44.7% efficiency201422327728210.1002/pip.24752-s2.0-84897866856GhosalK.BurroughsS.HeuserK.SetzD.Garralaga-RojasE.Performance results from micro-cell based high concentration photovoltaic research development and demonstration systems20132161370137610.1002/pip.22982-s2.0-84883051389GhosalK.LillyD.GabrielJ.WhiteheadM.SeelS.FisherB.WilsonJ.BurroughsS.Semprius field results and progress in system development20144270370810.1109/JPHOTOV.2013.22880262-s2.0-84897617320GlobaldataConcentrated Photovoltaics (CPV)—Global Market Size, Competitive Landscape and Key Country Analysis to 2020, UK2014FernándezE. F.AlmonacidF.RodrigoP.Pérez-HiguerasP.Model for the prediction of the maximum power of a high concentrator photovoltaic module201397121810.1016/j.solener.2013.07.0342-s2.0-84883236425FernándezE. F.AlmonacidF.Ruiz-AriasJ. A.Soria-MoyaA.Analysis of the spectral variations on the performance of high concentrator photovoltaic modules operating under different real climate conditions201412717918710.1016/j.solmat.2014.04.0262-s2.0-84901300707García-DomingoB.AguileraJ.de la CasaJ.FuentesM.Modelling the influence of atmospheric conditions on the outdoor real performance of a CPV (Concentrated Photovoltaic) module20147023925010.1016/j.energy.2014.03.1192-s2.0-84901676189ChanN. L. A.YoungT. B.BrindleyH. E.Ekins-DaukesN. J.ArakiK.KemmokuY.YamaguchiM.Validation of energy prediction method for a concentrator photovoltaic module in Toyohashi Japan20132181598161010.1002/pip.22412-s2.0-84888645129SteinerM.SieferG.HornungT.PeharzG.BettA. W.YieldOpt, a model to predict the power output and energy yield for concentrating photovoltaic modules201410.1002/pip.24582-s2.0-84892882062RodrigoP.FernándezE. F.AlmonacidF.Pérez-HiguerasP. J.Models for the electrical characterization of high concentration photovoltaic cells and modules: a review20132675276010.1016/j.rser.2013.06.0192-s2.0-84880381795ChanN. L. A.BrindleyH. E.Ekins-DaukesN. J.Impact of individual atmospheric parameters on CPV system power, energy yield and cost of energy201422101080109510.1002/pip.23762-s2.0-84878119857KurtzS.MullerM.JordanD.Key parameters in determining energy generated by CPV modules201410.1002/pip.2544FernándezE. F.Pérez-HiguerasP.Garcia LoureiroA. J.VidalP. G.Outdoor evaluation of concentrator photovoltaic systems modules from different manufacturers: first results and steps201321469370110.1002/pip.12622-s2.0-84878192388OsterwaldC. R.EmeryK. A.MullerM.Photovoltaic module calibration value versus optical air mass: the air mass function201422556057310.1002/pip.23032-s2.0-84901287615McMahonW. E.EmeryK. E.FriedmanD. J.OttosonL.YoungM. S.WardJ. S.KramerC. M.DudaA.KurtzS.Fill factor as a probe of current-matching for GaInP2/GaAs tandem cells in a concentrator system during outdoor operation200816321322410.1002/pip.7982-s2.0-42249097083MullerM.MarionB.KurtzS.RodriguezJ.An investigation into spectral parameters as they impact CPV module performance1277Proceedings of the 6th International Conference on Concentrating PhotovoltaicApril 2010Freiburg, GermanyAIP30731110.1063/1.35092182-s2.0-78649310508OsterwaldC. R.Translation of device performance measurements to reference conditions1986183-426927910.1016/0379-6787(86)90126-22-s2.0-0022777984EvansD. L.Simplified method for predicting photovoltaic array output198127655556010.1016/0038-092X(81)90051-72-s2.0-0019717697FernándezE. F.AlmonacidF.SarmahN.MallickT.SanchezI.CuadraJ. M.Soria-MoyaA.Pérez-HiguerasP.Performance analysis of the lineal model for estimating the maximum power of a HCPV module in different climate conditions1616Proceedings of the 10th International Conference on Concentrator Photovoltaic Systems (CPV '14)April 2014Albuquerque, NM, USA18719010.1063/1.4897057NofuentesG.García-DomingoB.MuñozJ. V.ChenloF.Analysis of the dependence of the spectral factor of some PV technologies on the solar spectrum distribution201411330230910.1016/j.apenergy.2013.07.0442-s2.0-84882585455Alonso-AbellaM.ChenloF.NofuentesG.Torres-RamírezM.Analysis of spectral effects on the energy yield of different PV (photovoltaic) technologies: the case of four specific sites20146743544310.1016/j.energy.2014.01.024FernándezE. F.SieferG.AlmonacidF.LoureiroA. J. G.Pérez-HiguerasP.A two subcell equivalent solar cell model for III-V triple junction solar cells under spectrum and temperature variations20139222122910.1016/j.solener.2013.03.0122-s2.0-84876702149DomínguezC.AntónI.SalaG.Multijunction solar cell model for translating I-V characteristics as a function of irradiance, spectrum, and cell temperature201018427228410.1002/pip.9652-s2.0-77952750017SieferG.BettA. W.Analysis of temperature coefficients for III-V multi-junction concentrator cells201422551552410.1002/pip.22852-s2.0-84867894675AlmonacidF.Pérez-HiguerasP. J.FernándezE. F.RodrigoP.Relation between the cell temperature of a HCPV module and atmospheric parameters201210532232710.1016/j.solmat.2012.06.0432-s2.0-84864490722FernándezE. F.AlmonacidF.RodrigoP.Pérez-HiguerasP.Calculation of the cell temperature of a high concentrator photovoltaic (HCPV) module: a study and comparison of different methods201412114415110.1016/j.solmat.2013.11.0092-s2.0-84888591771OtaY.NagaiH.ArakiK.NishiokaK.Temperature distribution in 820X CPV module during outdoor operation20121477364367FernándezE. F.RodrigoP.AlmonacidF.Pérez-HiguerasP.A method for estimating cell temperature at the maximum power point of a HCPV module under actual operating conditions201412415916510.1016/j.solmat.2014.01.0502-s2.0-84894541155StrobachE.FaimanD.KabaloS.BukobzaD.MelnichakV.GombertA.GerstmaierT.RoettgerM.Modeling a grid-connected concentrator photovoltaic system201410.1002/pip.2467MODIS Daily Level-3 data2014, http://gdata1.sci.gsfc.nasa.gov/daac-bin/G3/gui.cgi?instance_id=MODIS_DAILY_L3KingD. L.KratochvilJ.BoysonW. E.BowerW. I.Field expererience with a new performance characterization procedure for photovoltaic arraysProceedings of the 2nd World Conference and Exhibition on Photovoltaic Solar Energy ConversionJuly 1998Vienna, AustriaKingD. L.KratochvilJ. A.BoysonW. E.Measuring solar spectral and angle-of-incidence effects on photovoltaic modules and solar irradiance sensorsProceedings of the IEEE 26th Photovoltaic Specialists ConferenceOctober 1997Anaheim, Calif, USA111311162-s2.0-0031380189FanneyA. H.DoughertyB. P.DavisM. W.Short-term characterization of building integrated photovoltaic panels20031251132010.1115/1.15316422-s2.0-0037300790FanneyA.DavisM.DoughertyB.KingD.BoysonW.KratochvilJ.Comparison of photovoltaic module performance measurementsProceedings of the International Solar Energy Conference2005Orlando, Fla, USAFernándezE. F.AlmonacidF.Spectrally corrected direct normal irradiance based on artificial neural networks for high concentrator photovoltaic applications20147494194910.1016/j.energy.2014.07.075