APE Advances in Power Electronics 2090-1828 2090-181X Hindawi Publishing Corporation 10.1155/2015/654092 654092 Research Article Comprehensive Analysis and Experimental Validation of an Improved Mathematical Modeling of Photovoltaic Array Bal Satarupa 1 Anurag Anup 2 Nanda Mrutyunjaya 3 Sourav Suman 4 Pomilio Jose Antenor 1 Department of Electrical and Computer Engineering, National University of Singapore, 4 Engineering Drive 3 Singapore 117583 nus.edu.sg 2 Department of Electrical Engineering and Information Technology, ETH Zurich, Rämistrasse 101, 8092 Zurich Switzerland ethz.ch 3 Department of Electrical Engineering, Asian Institute of Technology, Pathumthani 12120 Thailand ait.ac.th 4 School of Computing, National University of Singapore, 13 Computing Drive Singapore 117417 nus.edu.sg 2015 1312015 2015 26 08 2014 30 11 2014 01 12 2014 1312015 2015 Copyright © 2015 Satarupa Bal et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper proposes a simple, accurate, and easy to model approach for the simulation of photovoltaic (PV) array and also provides a comparative analysis of the same with two other widely used models. It is highly imperative that the maximum power point (MPP) is achieved effectively and thus a simple and robust mathematical model is necessary that poses less mathematical complexity as well as low data storage requirement, in which the maximum power point tracking (MPPT) algorithm can be realized in an effective way. Further, the resemblance of the P-V and I-V curves as obtained on the basis of experimental data should also be taken into account for theoretical validation. In addition, the study incorporates the root mean square deviation (RMSD) from the experimental data, the fill factor (FF), the efficiency of the model, and the time required for simulation. Two models have been used to investigate the I-V and P-V characteristics. Perturb and Observe method has been adopted for MPPT. The MPP tracking is realized using field programmable gate array (FPGA) to prove the effectiveness of the proposed approach. All the systems are modeled and simulated in MATLAB/Simulink environment.

1. Introduction

Electrical energy from photovoltaic is currently regarded as the prerequisite sustainable resource for both stand-alone as well as grid connected applications, since it is abundant and clean, offers zero input fuel cost, and is distributed throughout the earth . In practical cases, photovoltaic modules operate over a highly intermittent nature of temperature and irradiance but the electrical parameters provided in the datasheet are only for the standard test conditions (STC). Moreover, in power generation from PV, optimal utilization of the available solar energy is imperative due to the high costs of PV modules. It is also seen that mathematical models of few individual components of PV system are represented and simulated for better understanding of their performances .

This calls for a simple, accurate, and easy to model approach for the simulation of photovoltaic (PV) module to track the maximum power point and to predict PV energy production under varying atmospheric conditions . In order to increase the accuracy, the following can be incorporated, but it leads to the increase in complexity of the modeling :

temperature dependence of the diode saturation current,

temperature dependence of the photo current,

inclusion of series resistance: for more accurate shape between the MPP and the open circuit (OC) voltage,

inclusion of shunt resistance in parallel with the diode,

variability of diode quality factor,

introduction of two or more parallel diodes.

The accuracy of the simulation of a PV model largely depends on the estimation of the characteristic I-V and P-V curves. Furthermore, factors such as efficiency, field factor, and simulation time affect the effectiveness of the model. A simplistic and easy to model approach is preferred so as to avoid unwanted complexity due to additional parameters.

So far, among the mathematical models of PV array proposed in the literature, the simplest is the ideal single diode model which involves only three parameters, namely, short circuit current, open circuit voltage, and the diode ideality factor . Improvement has been made, with the simplified single diode model (SSDM) being proposed in  which takes the effect of the series resistance ( R s ) which is the sum of several types of structural resistance of the device into consideration . The influence of R s only becomes dominant when the PV device operates in the voltage source region. Also, it lacks the accuracy when subjected to large temperature variations . Since, the value of R s is very low, some authors neglect its effect [5, 1719]. Further improvement has been done by the introduction of the single diode model (SDM) which includes the additional shunt resistance ( R p ) along with the series resistance . The shunt resistance exists mainly due to the leakage current of the p-n junction. The effect of R p is dominant when the PV device operates at current source region of operation. However, since the value of R p is very high many authors [4, 8, 9, 20, 21] neglect it in order to simplify the model. Although it is much more accurate than the previous models, it is not preferred on account of its computational complexity. It is also reported in  that the accuracy of this model deteriorates at low irradiance levels. In order to mitigate the inaccuracies offered by the previous models, the two-diode model was proposed in . However, this leads to more model complexity and thus more simulation time due to the involvement of a greater number of parameters. A new mathematical PV model has also been proposed in  that includes the advantages of previous models combining the three main considerations, namely, simplicity, ease of modeling, and accuracy. However, it doesnt take into consideration the effect of diode saturation current on temperature which results in model errors at the vicinity of open-circuit voltage and consequently at other regions.

This paper proposes a new, simple, accurate, and easy to model approach for the simulation of PV array and also provides a comparative analysis of the same with the conventional single diode model and the improved ideal single diode model. As the PV systems are generally integrated with specific control algorithms in order to extract the maximum possible power, it is highly imperative that the MPP is achieved effectively and thus it is needed to design a model from which the MPPT algorithm can be realized in an effective way. Some MPPT techniques have been proposed in [1, 3, 4, 10]. However, for simplicity, this paper adopts the Perturb and Observe (P&O) method for MPPT.

The proposed theoretical model is verified and validated with experimental data of commercial PV array. RMSD from the experimental data, maximum efficiency of the design, the fill factor (FF), and the simulation time has also been calculated. In addition, the MPP tracking is realized in digital environment using FPGA kit to prove the effectiveness of the proposed approach. All the systems here are modeled and simulated in MATLAB/Simulink environment. The proposed modeling method can be useful for users who require simple, fast, and accurate models in simulation of PV systems.

2. Mathematical Models for a Photovoltaic Module

The major issue of real-time identification is basically the selection of a proper model. It is therefore necessary to have a proper mathematical model that can represent accurately the current-voltage characteristics of the PV array and which can be solved by analytical methods in a simplified manner . In addition to this to maximize the power extracted from a PV array with the help of MPPT control, the understanding and modeling of PV cell are also important .

Assuming the semiconductor diode equation and the Kirchhoff laws, the I - V characteristics for a PV module composed of series connected cells based on single exponential model are expressed as follows : (1) I l = I p v - I 0 exp q V + I l R s N s A k T - 1 - V + I L R s N s N S R p , where k is the Boltzmann constant ( 1.3806 × 1 0 - 23  J/K) and q the electron charge ( 1.60217 × 1 0 - 19  C). T gives the module temperature. The parameter I p v gives the photocurrent, I 0 represents the diode saturation current, and I l gives the output current. R s and R p give the series resistance and the shunt resistance. A and N s represent the diode ideality factor and the number of cells connected in series, respectively.

The first term I p v gives the photocurrent and the second part is the ideal dark current that models the emitter and base recombination. All the parameters are mostly calculated through sets of nonlinear equations .

2.1. Single Diode Model (SSDM)

The single diode model takes into account both the series resistance as well as the shunt resistance unlike the ideal single diode model or the simplified single diode model as shown in Figure 1. This resistance is the sum of several types of structural resistance of the device. R s depends mainly on four factors, namely ,

contact resistance of the metal base with p-layer,

resistance of p-layer and n-layer,

contact resistance of the metal grid with the n-layer,

resistance of the grid.

PV cell modeled as single diode model circuit.

The shunt resistance exists mainly due to the leakage current of the p-n junction. It depends basically on the fabrication method of the PV cell.

Mathematically, the current equation can be written in accordance with Kirchhoff's current law: (2) I p v = I d + I p + I l , where I p v is the photocurrent generated due to the incident light and I d is the diode current (Shockley diode equation). I l gives the output current and there is an additional term, I p which represents the leakage current of the p-n junction. In this model, the diode is given by (3) I d = I 0 exp q V + I l R s N s A k T - 1 , where the module saturation current ( I 0 ) is given by (4) I 0 T = I scr T ref + K i Δ T exp q V oc T ref + K v Δ T / A k T , where I scr is the PV module short circuit current (in A) and can be found in the product datasheet. K i is the short circuit current temperature coefficient (in A/K). Thus, from (3) and (4) the diode current can be calculated. Now, I p v can be calculated by (5) I p v = G × I scr + K i Δ T , where G is the PV module incident illumination (in kW/m2). An additional current I p is introduced here: (6) I p = V + N s I l R s N s R p . Therefore, there are three basic unknown parameters: A , R s , and R p . Here, the value of A is assumed ( 1.1 < A < 1.6 ) . The value of the resistance can be calculated from the information available in the datasheet.

For the calculation of the series resistance and the shunt resistance, an iterative method is employed. Some authors vary the resistance independent of each other leading to poor accuracy of results. Here, a concept of simultaneously varying both the values is shown. The idea is taken from the fact that the maximum calculated power should be equal to the power mentioned in the datasheet. From (2), (3), and (6), the current equation can be written and the equation at MPP gives the value of R p to be (7) R p = V mpp V mpp + N s I mpp R s N s V mpp I p v - I d , mpp - P max , where (8) I d , mpp = I 0 exp q V mpp + I mpp R s N s A k T .

The value of I 0 can be calculated as from (4) and I p v can be calculated from (5). However, some simplifications have to be made in order to find the unknown parameters. An initial guess is taken as (9) R s , initial = 0 , R p , initial = V mpp I scr ( T ref ) - I mpp - V oc T ref - V mpp I mpp .

Using an iterative procedure, by increasing the value of the series resistance and simultaneously updating the shunt resistance so as to match the maximum power, the output current equation can be found out. The major drawback of this model lies on its assumption of ideality factor. Also, the computational complexity and number of iterations required to obtain the output current equation are considerable.

2.2. Improved Ideal Single Diode Model

The improved ideal single diode model is basically based on the ideal single diode model as shown in Figure 2. However, the modeling involves a set of mathematical equations which produces sufficiently accurate results but with much reduced complexity.

PV cell modeled as improved ideal single diode model circuit.

The series and shunt resistances are neglected for mathematical simplicity. However, the method of deriving the parameters is of reduced complexity. The computation of these equations avoids the use of a nonlinear solver .

The current relation can be found by applying KCL: (10) I p v = I d + I l . Here, the current through the diode is given by (11) I d = I 0 exp q V A k T - 1 . The derivation of the saturation current I 0 begins by (12) V oc G , T - V oc G , T ref = - β Δ T , where V oc ( G , T ) and V oc ( G , T ref ) represent the open circuit voltages at a temperature T and at the reference temperature T ref . β gives the voltage temperature coefficient and that can be found from the product datasheet. The open circuit voltage can be found out by putting I l = 0 in (10) and equating the value of I d as in (11): (13) V oc = A k T q ln I p v I 0 + 1 . Now, finding the values of V oc at both the temperatures using (13) and replacing in (12) the following is obtained: (14) k A q T ln G I scr + K i Δ T I 0 + 1 - T ref ln G I scr I 0 T ref + 1 T ln G I scr + K i Δ T I 0 + 1 = - β Δ T . Rearranging (14), the following is found: (15) I 0 = exp q β Δ T / A k T × G × I scr + K i Δ T G I scr / I 0 T ref + 1 T ref / T    - exp q β Δ T / A k T . The I 0 ( T ref ) can be calculated according to (4). Now, writing (11) at MPP, (16) I mpp = I p v - I 0 exp q V mpp A k T - 1 , and here according to (5), at reference temperature, (17) I p v = I scr . The ideality factor A can be derived by substituting (15) and (17) into (16): (18) I mpp I scr = exp q V mpp A k T ref - I scr - I mpp I scr exp q V oc A k T ref . Thus, the unknown parameters can be found out from the above equations and from the product datasheet information.

However, the temperature dependence of the saturation current has not been considered. Also, there is a slight deviation of the I - V characteristics from the experimental curve, due to the approximation R s = 0 .

2.3. Improved PV Modeling Approach

In the proposed improved PV modeling approach, in order to avoid the complexity offered by the single diode model, it relies on the simplified single diode model, by neglecting the shunt resistance, as shown in Figure 3. Applying KCL to the above circuit, the current relation is found to be the same as (10).

PV cell modeled as improved simplified single diode model circuit.

The diode current equation incorporates the additional voltage drop across the series resistance: (19) I d = I 0 exp q V + I l R s A k T - 1 . I p v is calculated according to (5). The current equation for the simplified single diode model is given by (20) I l = I p v - I 0 exp q V + I l R s A k T - 1 . Now, there are basically four unknown parameters: I p v , I 0 , A , and R s . I p v can be determined from the information available in the manufacturer s datasheet by applying (5).

The value of the series resistance can be calculated from the information available in the product datasheet. Since, the current at maximum power point at reference temperature is available in the product datasheet, (20) is solved for maximum power point conditions: (21) I mpp = I p v - I 0 exp q V mpp + I mpp R s A k T - 1 . Thus, the series resistance can be found from (22) R s = ( A k T / q ) × ln ( ( I p v - I mpp ) / I 0 + 1 ) - V mpp I mpp . In order to find the absolute value of all these parameters, the ideality factor and the saturation current should be calculated effectively. In the earlier works, the value of the ideality factor has been assumed which leads to the degradation of the curve. This model has the novelty of calculating the value of the ideality factor from the datasheet parameters so that the simulated curves coincide with the experimental data to a larger extent. Also, an expression has been provided for finding the saturation current from the information given in the datasheet.

For the calculation of the ideality factor, advantage has been taken of the fact that the derivative of power with respect to voltage at MPP is zero.

Differentiating (20) with respect to voltage gives (23) d I l d V = - I 0 exp q v A k T × q R s A k T × exp q I l R s A k T × d I l d V + exp q I l R s A k T × q A k T × exp q V A k T exp q v A k T × q R s A k T × exp q I l R s A k T . Now, it is known at MPP the derivative power with respect to voltage is zero. And hence the derivative of load current with respect to voltage is given by (24) d I l d V MPP = - I mpp V mpp . So, at MPP, (25) - I mpp V mpp = - I 0 exp q V A k T × q R s A k T × exp q I l R s A k T × - I mpp V mpp exp q V A k T × q R s A k T × exp q I l R s A k T × - I m p p V m p p + exp q I l R s A k T × q A k T × exp q V A k T . Thus, the value of the saturation current equals (26) I 0 T = I mpp × V mpp × exp q V A k T × q R s A k T × exp q I l R s A k T × - I mpp V mpp exp q V A k T × q R s A k T × exp q I l R s A k T × - I m p p V m p p + exp q I l R s A k T × q A k T × exp q V A k T - 1 . The reverse saturation current is given by (4). Equating (4) and (26) at reference temperature (27) I scr T ref + K i Δ T exp q V oc T ref + K v Δ T / A k T = I mpp × V mpp × exp q V A k T × q R s A k T × exp q I l R s A k T × - I mpp V mpp + exp q I l R s A k T × q A k T × exp q V A k T q I l R s A k T - 1 , R s / A can be obtained from (22) and is replaced in (27). Now, the equation has only one unknown parameter A . Also, the saturation current can be calculated from (26).

This modification aims at smoothing the curve between the MPP and the open circuit point and also to match the open circuit voltages for a large range of temperature variations.

3. Validating the Model

The PV array model is simulated in the MATLAB/ Simulink environment in order to validate the claims. Tables 1 and 2 give a comparative view on the parameters on the datasheet and the parameters obtained from the proposed simulation model. It is seen that the three main points coincide with that of the datasheet. The characteristic curves developed from the model have been plotted with the experimental data at four different temperature conditions and three different irradiation values. This has been done so as to effectively verify the notion that the proposed model gives acceptable results even if the conditions are considerably far from STC as shown in Figures 4 and 5.

Parameters of TBP-1237 solar module at STC (Temp: 25°C; A.M.: 1.5; G = 1000 W/M2).

Parameter Label Value from Value from
datasheet model
Maximum power P max 20 W 20 W
Open circuit voltage V oc 20.91 V 20.91 V
Short circuit current I scr 1.3 A 1.3 A
Voltage at MPP V mpp 17.1 V 17.1 V
Current at MPP I mpp 1.17 A 1.17 A

Parameters of MSX60 solar module at STC (Temp: 25°C; A.M.: 1.5; G = 1000 W/M2).

Parameter Label Value from Value from
datasheet model
Maximum power P max 59.85 W 59.85 W
Open circuit voltage V oc 21.1 V 21.1 V
Short circuit current I scr 3.8 A 3.8 A
Voltage at MPP V mpp 17.1 V 17.1 V
Current at MPP I mpp 3.5 A 3.5 A

(a) I - V model curves with experimental data for TBP-1237 solar array at different temperatures and (b) I - V model curves with experimental data for TBP-1237 solar array at different irradiations.

(a) I - V model curves with experimental data for MSX60 solar array at different temperatures and (b) I - V model curves with experimental data for MSX60 solar array at different irradiations.

Figure 4 gives the I - V curves for TBP-1237 at different temperatures and different irradiations, respectively.

However, in order to test the validity efficiently, it is required to test the model with other solar panels. Figures 5(a) and 5(b) give the I - V curves for MSX60 at different temperatures and different irradiations, respectively.

3.1. Effect of Variation in Temperature and Irradiation

It can be seen that the values of the parameters found can be used for irradiation and temperature close to the STC. With a large change in the temperature or irradiation, a small modification of the algorithm can be put forward to account for the change in the model parameters.

Although there is an increase of the photo current with the increase in temperature increase owing to the slight decrease of the band gap energy, E g , the main output characteristics, such as efficiency, are negatively influenced by high temperature. The series resistance increases with increase in temperature. However, a decrease is observed with the increase in irradiation . On the other hand, the junction characteristic parameters such as ideality factor and saturation current are highly temperature dependent. The variation of ideality factor is almost linear with temperature. It increases with increase in temperature . This is explained from the fact that, at the increased temperatures, imperfections of basic material are more pronounced and thus create defects in the lattice structure.

The short-circuit current from a solar cell depends linearly on light intensity which results in the increase in PV output power as the solar radiation increases. The series resistance is basically derived from the slope of the I - V curve between the MPP and the open-circuit point. It is hence treated as fitting parameter rather than fixed panel s effective series resistance.

The slope is strongly dependent on the panel short-circuit current and maximum power point current and therefore also on the irradiation. Hence, it tends to increase with increase in light intensity.

It has to be pointed out that the variation of the parameters with respect to temperature can be accounted for the calculation of I p v directly from (5). The dependence of resistance with respect to temperature has been shown in (22). However, the change in the series resistance and the ideality factor of the modeling method can be accounted for. The series resistance can also be given by  (28) R s = 2 V oc I sc - Ar I sc 2 - A V t I sc . The values of V oc and I sc can be determined at any given temperature and irradiation according to . Ar gives the area under the curve. Now, an iterative procedure is done so as to find the exact values for the ideality factor and the series resistance. This can be described as Figure 6. “ e ” can be defined as any small value close to zero. The ideality factor shows a linear behavior and, thus, it can be calculated manually by the linear nature as in  for the practical working temperatures: (29) A A STC = T T STC , where A STC and T STC are the ideality factor and the temperature at STC. This model thus proves to be an efficient model as compared to the previous approaches, even though it neglects the shunt resistance, on account of the fact that all the parameters are found from the datasheet parameters, instead of taking any assumptions into account. The improved ideal single diode model also finds the parameters but it neglects the series resistance which leads to increased deviation from the experimental values. Since, none of the parameters are assumed; this approach provides a better model than its predecessors, which rely on intelligent assumptions so as to fit the curve.

Flowchart for determining the fianl values of the parameters at any temperature and irradiation.

4. Comparative Results and Analysis

With the purpose of comparison between the three different approaches of mathematical modeling of PV array and validating them experimentally, a small-scale module of has been considered. The accuracy of the model is experimentally validated using TBP-1237 and MSX60. The proposed PV circuit model is implemented using a current controlled source and simple computational blocks. The inclusion of series resistance and ideality factor as unknown parameters aids in the smoothness of the output characteristics between MPP and open-circuit voltage and, thus, the curve coincides closely with the experimental data obtained not only at the three main points but throughout the region. Table 3 gives the components used in the prototype. For experimental purpose, three arrays of bulbs with three bulbs of 200 W in each row are used as artificial sun. The conventional DC-DC boost converter with a control tracker is implemented for MPP tracking.

Components used in the study.

Parameter Label Value/equipment
Inductance L 1.7 mH
Output capacitor C out 300 μF
A to D converter ADC THS1030
Current sensor CS LEM LA-55P

4.1. Analysis on Basis of <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M118"> <mml:mi>P</mml:mi> <mml:mtext>-</mml:mtext> <mml:mi>V</mml:mi></mml:math> </inline-formula> and <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M119"> <mml:mi>I</mml:mi> <mml:mtext>-</mml:mtext> <mml:mi>V</mml:mi></mml:math> </inline-formula> Curves

The I - V and P - V output curves for different mathematical models are compared with the experimental results for determining the accuracy of the models and their closeness to the actual curve. For the single diode model, the graph obtained coincides with the experimental result approximately at I sc and V oc but the MPP deviates from the one obtained from experimental result. For the improved ideal single diode model, the simulation result at MPP points exhibit less deviation with improvement in V oc point as compared to the former. It takes advantage of the simplicity of ideal models and enhances the accuracy by deriving a mathematical representation, capable of extracting accurate estimates of the model parameters, directly related to manufacturer datasheets. The inclusion of the calculation of ideality factor affects the curvature of the I - V curve and thus expresses the ideality of the diode. The mathematical derivation of the same improves the model fitting. However, owing to the assumption of no series or shunt resistance, the model shows some deviation from the actual curve between the MPP and the open circuit point. Extracting the advantage of the above two models, an improved simplified single diode model has been proposed here, where the result closely matches not only at the three points: the open circuit point, the MPP, and the short circuit point, but also throughout the curve.

4.2. Analysis on Basis of Simulation Time

The simulation time taken by the models has been shown in Table 4. More simulation time indicates a higher complexity of the model. As expected, the single diode model takes the highest time on account of the computational complexity offered by it. It is then followed by the proposed model and the lowest time is taken by the improved ideal single diode model. However, the time difference offered by the latter ones is very less and hence no significant time is elapsed during simulation.

Parameters used for comparative analysis.

Model MPP points RMSD FF Max. eff. (in %) Simulation time
Single diode Model 17.092 V1.172 A 0.0561 0.7256 15.7366 2.14 s

Improved ideal single diode model 17.087 V1.169 A 0.0587 0.7244 15.6881 0.86 s

Proposed model 17.1 V1.17 A 0.0531 0.7263 15.7862 0.89 s

4.3. Analysis of Basis of Fill Factor

Another factor which determines the quality of a model is the fill factor. This factor gives the ratio of the power produced at the maximum power point to the maximum theoretical power that should be possible to extract from the module. For good cells, the value should generally be more than 0.7. It is known that the higher the fill factor, the better the model. Table 4 shows the fill factors for various models. Much difference is not observed between the fill factors. However, it is seen that the proposed model gives the highest fill factor. It is followed by the single diode model and the least fill factor is seen in the improved ideal single diode model. It is due to its deviation in the curve from the actual experimental data. Although the curve at the MPP deviates towards the lower side in case of single diode model, it shows a higher fill factor on account of its greater deviation towards the higher side at the other regions.

4.4. Analysis on Basis of Maximum Efficiency

The subsequent parameter used for the analysis is the maximum efficiency which is the ratio between the maximum power and the incident power. At an irradiation of 1000 W/m2 and a temperature of 35°C, the maximum efficiency offered by the models in % is illustrated in Table 4.

4.5. Analysis on Basis of RMSD

The accuracy of a particular model can be known from the RMSD it shows. Higher the RMSD value, less accurate the model. Table 4 shows the RMSD comparison between the various PV models. The least RMSD is shown by the proposed model, since it takes into account all the parameters. The improved ideal single diode model also shows a low RMSD value even if some of the parameters are ignored since it is not based on any assumption in calculation. The single diode model shows a better RMSD than the improved ideal single diode model. However, it is not true in all cases. In some cases of varying temperature and irradiation, the improved ideal single diode model proves to be a better model than the single diode model. The absolute errors for these models are compared in Figure 9(a) for STC. It is seen that the proposed model provides better results when compared with the results found by different approaches for the same module.

In order to verify the notion that the proposed model gives better accuracy, even far from STC, comparisons have been made based on the absolute values of error as shown in Figures 9(b) and 9(c).

4.6. Experimental Validation for MPP Tracking

The proposed method being an improvement over the conventional method and exhibiting closeness to the real cell at MPP is validated by comparing with the experimental result. Figure 7 depicts the prototype used with PV panel as the input to the converter and battery of 45 W as its output. The current and voltage sensed by the sensors are used to generate duty ratio in the FPGA environment. With battery as the output of the converter, the input is regulated by the P&O control algorithm where irrespective of the alteration in irradiation MPP point for each P - V curve is tracked efficiently. Figure 8 shows the experimental waveform of the MPP points for three different irradiation levels which perfectly matches with the one found from the modeling.

(a) Bulbs used as artificial sun for PV panel and (b) prototype of the conventional boost converter used.

(a) Experimental waveforms of power, voltage, and current for TBP-1237 module ( Y -Axis-V/div:10 V/div;A/div:1 A/div;W/div:10 W/div), (b) simulation PV waveforms, and (c) experimental power points from MPP tracking using P and O method.

Comparison of absolute errors among three different approaches for TBP-1237 module at (a) STC (Temp: 25°C; A.M.: 1.5; G = 1000 W/m2), (b) Temp: 35°C; A.M.: 1.5; G : 900 W/m2, and (c) Temp: 35°C; A.M.: 1.5; G =1000 W/m2.

5. Conclusions

This paper analyses the development of an effective approach to modeling of PV modules used for simulation purposes. The proposed modeling fits the mathematical characteristic model equations to the experimental curve. In addition to it, it avoids complexities and achieves better accuracy by the inclusion of the series resistance. The dependency of diode saturation current on the temperature as well as on the MPP helps in making the equation effective and allows the user to adjust the output curve at the open-circuit voltage, short circuit current point, and the MPP at different temperatures and irradiations. Also, the inclusion of the calculation of the ideality factor from the information available in the datasheet improves the accuracy of the curve and fits the curve with the experimental one. The results are experimentally validated to prove the effectiveness of the proposed method. The Perturb and Observe MPPT algorithm implemented along with the boost converter at different irradiation substantiate the closeness of the MPP points of the proposed model to the real time observed points. Along with it, the comparison of the proposed model with the conventional along with the improved single diode model in terms of the fundamental quantities such as fill factor, RMSD, and maximum efficiency gives a better insight into the advantage of the proposed model to the previous ones.

The proposed mathematical model is easy to implement in various simulation platforms for PV power system studies and it avoids the underlying complexities involved in PV parameter identification. This paper is intended to be a useful tool, serving the cause of both beginners as well as seasoned users.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

de Brito M. A. G. Sampaio L. P. Junior L. G. Canesin C. A. Evaluation of MPPT techniques for photovoltaic applications Proceedings of the IEEE International Symposium on Industrial Electronics (ISIE '11) June 2011 Gdansk, Poland 1039 1044 10.1109/isie.2011.5984303 2-s2.0-80052866518 Villalva M. G. Gazoli J. R. Filho E. R. Comprehensive approach to modeling and simulation of photovoltaic arrays IEEE Transactions on Power Electronics 2009 24 5 1198 1208 10.1109/TPEL.2009.2013862 2-s2.0-67349141404 Sokolov M. Green T. C. Mitcheson P. D. Shmilovitz D. Dynamic analysis of photovoltaic system with MPP locus emulation Proceedings of the IEEE 26th Convention of Electrical and Electronics Engineers in Israel (IEEEI '10) November 2010 Eliat, Israel 212 215 10.1109/eeei.2010.5661961 2-s2.0-78651231382 Walker G. Evaluating MPPT converter topologies using a matlab PV model Journal of Electrical & Electronics Engineering, Australia 2001 21 1 49 55 2-s2.0-0035173089 Chan D. S. H. Phang J. C. H. Analytical methods for the extraction of solar-cell single- and double-diode model parameters from I-V characteristics IEEE Transactions on Electron Devices 1987 34 2 286 293 2-s2.0-0023292133 Gow J. A. Manning C. D. Development of a photovoltaic array model for use in power-electronics simulation studies IEE Proceedings: Electric Power Applications 1999 146 2 193 200 10.1049/ip-epa:19990116 2-s2.0-0032629423 Ahmad G. E. Hussein H. M. S. El-Ghetany H. H. Theoretical analysis and experimental verification of PV modules Renewable Energy 2003 28 8 1159 1168 10.1016/s0960148102002288 2-s2.0-0037410823 Veerachary M. PSIM circuit-oriented simulator model for the nonlinear photovoltaic sources IEEE Transactions on Aerospace and Electronic Systems 2006 42 2 735 740 10.1109/TAES.2006.1642586 2-s2.0-33745727042 Celik A. N. Acikgoz N. Modelling and experimental verification of the operating current of mono-crystalline photovoltaic modules using four- and five-parameter models Applied Energy 2007 84 1 1 15 10.1016/j.apenergy.2006.04.007 2-s2.0-33749242233 Kuo Y.-C. Liang T.-J. Chen J.-F. Novel maximum-power-point-tracking controller for photovoltaic energy conversion system IEEE Transactions on Industrial Electronics 2001 48 3 594 601 2-s2.0-0035366448 10.1109/41.925586 Elhagry M. T. Elkousy A. A. T. Saleh M. Elshatter T. Abou-Elzahab E. M. Fuzzy modeling of photovoltaic panel equivalent circuit 1 Proceedings of the 40th Midwest Symposium on Circuits and Systems August 1997 60 63 Liu S. Dougal R. A. Dynamic multiphysics model for solar array IEEE Transactions on Energy Conversion 2002 17 2 285 294 10.1109/TEC.2002.1009482 2-s2.0-0036601542 Sera D. Teodorescu R. Rodriguez P. PV panel model based on datasheet values Proceedings of the IEEE International Symposium on Industrial Electronics (ISIE '07) June 2007 Vigo, Spain 2392 2396 10.1109/isie.2007.4374981 2-s2.0-47049122000 Vitorino M. A. Hartmann L. V. Lima A. M. N. Corrêa M. B. R. Using the model of the solar cell for determining the maximum power point of photovoltaic systems Proceedings of the European Conference on Power Electronics and Applications (EPE '07) September 2007 Aalborg, Denmark 1 10 10.1109/epe.2007.4417684 2-s2.0-51049122553 Patel H. Agarwal V. MATLAB-based modeling to study the effects of partial shading on PV array characteristics IEEE Transactions on Energy Conversion 2008 23 1 302 310 2-s2.0-40849142550 10.1109/TEC.2007.914308 Noroozian R. Abedi M. Gharehpetian G. Hosseini S. Modelling and simulation of microturbine generation system for on-grid and offgrid operation modes Bal S. Anurag A. Babu B. C. Comparative analysis of mathematical modeling of Photo-Voltaic (PV) array Proceedings of the Annual IEEE India Conference (INDICON '12) December 2012 Kochi, India 269 274 2-s2.0-84874124451 10.1109/indcon.2012.6420627 Tan Y. T. Kirschen D. S. Jenkins N. A model of PV generation suitable for stability analysis IEEE Transactions on Energy Conversion 2004 19 4 748 755 10.1109/TEC.2004.827707 2-s2.0-10244260364 Kajihara A. Harakawa T. Model of photovoltaic cell circuits under partial shading Proceedings of the IEEE International Conference on Industrial Technology (ICIT '05) December 2005 866 870 10.1109/icit.2005.1600757 2-s2.0-33847247911 Glass M. Improved solar array power point model with spice realization 1 Proceedings of the 31st Intersociety Energy Conversion Engineering Conference (IECEC '96) August 1996 286 291 Altas I. H. Sharaf A. M. A photovoltaic array simulation model for matlab-simulink GUI environment Proceedings of the International Conference on Clean Electrical Power (ICCEP '07) May 2007 Capri, Italy 341 345 10.1109/iccep.2007.384234 2-s2.0-35348881452 Salam Z. Ishaque K. Taheri H. An improved two-diode photovoltaic (PV) model for PV system Proceedings of the Joint International Conference on Power Electronics, Drives and Energy Systems & Power India (PEDES '10) December 2010 New Delhi, India IEEE 1 5 2-s2.0-79952725860 10.1109/pedes.2010.5712374 Mahmoud Y. Xiao W. Zeineldin H. H. A simple approach to modeling and simulation of photovoltaic modules IEEE Transactions on Sustainable Energy 2012 3 1 185 186 2-s2.0-84859042388 10.1109/tste.2011.2170776 Bonkoungou D. Koalaga Z. Njomo D. Modelling and simulation of photovoltaic module considering single-diode equivalent circuit model in MATLAB International Journal of Emerging Technology and Advanced Engineering 2013 3 3 493 502 Hatti M. Meharrar A. Tioursi M. Novel approach of maximum power point tracking for photovoltaic module neural network based Proceedings of the International Symposium on Environment Friendly Energies in Electrical Applications November 2010 Ghardaïa, Algeria 1 6 Nishioka K. Sakitani N. Uraoka Y. Fuyuki T. Analysis of multicrystalline silicon solar cells by modified 3-diode equivalent circuit model taking leakage current through periphery into consideration Solar Energy Materials and Solar Cells 2007 91 13 1222 1227 2-s2.0-34249808493 10.1016/j.solmat.2007.04.009 Bal S. Anurag A. Babu B. C. Comparative analysis of mathematical modeling of Photo-Voltaic (PV) array Proceedings of the Annual IEEE India Conference (INDICON '12) December 2012 Kochi, India 269 274 10.1109/indcon.2012.6420627 2-s2.0-84874124451 Chowdhury S. Chowdhury S. P. Taylor G. Song Y. H. Mathematical modelling and performance evaluation of a stand-alone polycrystalline pv plant with mppt facility Proceedings of the 2008 IEEE Power and Energy Society General Meeting—Conversion and Delivery of Electrical Energy in the 21st Century July 2008 1 7 10.1109/PES.2008.4596376 de Soto W. Klein S. A. Beckman W. A. Improvement and validation of a model for photovoltaic array performance Solar Energy 2006 80 1 78 88 10.1016/j.solener.2005.06.010 2-s2.0-29144478675 Sabry M. Ghitas A. E. Influence of temperature on methods for determining silicon solar cell series resistance Journal of Solar Energy Engineering 2007 129 3 331 335 10.1115/1.2735350 2-s2.0-34548205271 Yordanov G. H. Midtgård O.-M. Saetre T. O. Equivalent cell temperature calculation for PV modules with variable ideality factors Proceedings of the 38th IEEE Photovoltaic Specialists Conference (PVSC '12) June 2012 Austin, Tex, USA 505 508 10.1109/pvsc.2012.6317666 2-s2.0-84869452321