ORIGINAL RESEARCH article

Front. Energy Res., 20 January 2025

Sec. Solar Energy

Volume 13 - 2025 | https://doi.org/10.3389/fenrg.2025.1501335

Comparative study of parameter extractions of photovoltaic modules using analytical and numerical approaches

  • Abdelmalek Essaadi University, Tetouan, Morocco

Abstract

Developing an accurate mathematical model for parameter extraction in photovoltaic modules is a crucial endeavor in optimizing photovoltaic energy systems. This study seeks to assess and compare various analytical and numerical methods for extracting the primary five parameters of photovoltaic modules. Specifically, six established approaches based on a single diode model (SDM) are employed, including the methods introduced by Khan et al., Blas et al., Phang et al., Vika, Cubas et al., and Almonacid et al. The performance of these approaches is evaluated and compared under standard test conditions (STC) with a focus on maximum power point current and voltage. The analytical and numerical methods demonstrate their precision in predicting photocurrent-voltage (I-V) and power-voltage (U-V) curves, with the exception of the Almonacid et al. method, which tends to underestimate the I-V curve at the module’s maximum power. Among these methods, the Phang et al. approach stands out, displaying a strong agreement between experimental data and the predicted curve. This is evidenced by the lower values of root mean square error (RMSE), mean bias error (MBE), normalized RMSE (NRMSE), mean absolute percentage error (MAPE), and absolute error (AE). These findings underscore the high quality of results obtained through the Phang et al. method.

1 Introduction

The escalating costs of fossil fuels and their detrimental impact on the environment, including air pollution resulting from greenhouse gas emissions like carbon dioxide and methane, have spurred global efforts to explore and develop clean and renewable energy technologies (). Solar energy, with its abundant intensity, offers a swift and straightforward means of conversion into photovoltaic (PV) electricity. This conversion leverages the inherent properties of semiconductors, making it a promising and effective renewable energy resource (). In various energy conversion systems, it is possible to represent certain electrical photovoltaic (PV) characteristics using an equivalent electrical circuit. This approach is suitable for accurately emulating the actual behavior of solar cells, ensuring that the simulated data closely aligns with the measured current-voltage (I-V) data across all operational conditions (). Among these representation methods, some well-known equivalent electrical circuits include the single diode model (SDM), double diode model (DDM), and third diode model (TDM). These models have gained popularity and are employed by numerous researchers due to their high precision and performance. Moreover, the utilization of an increased number of parameters for extracting PV module characteristics within the single diode model enhances the accuracy of predicting the I-V behavior under varying solar radiation levels. This presents an intriguing opportunity in the research literature.

Parameter extraction methods can be categorized based on several factors, including the quantity of data samples involved in the extraction process, the specific method utilized, and the mathematical approach applied. In terms of the mathematical approach, the extraction of parameters from photovoltaic modules is typically classified into three main categories: numerical, analytical, and evolutionary methods. The primary objective of these methods is to establish a strong correlation and achieve an optimized fit between the theoretical and experimental I-V characteristics of solar cells (; ). Given that a significant number of I-V data samples are utilized during the parameter extraction process, these optimization methods are considered to collectively result in substantial improvements in obtaining the final parameter values. The relevant literature offers a wealth of references on parameter extraction. An uncomplicated method for extracting the five parameters was introduced by , Cubas et al. is one of the most cited contributors in this field. Their technique, aimed at reducing the gap between estimated and experimental peak power values, involves incorporating the series resistance as a parameter. However, a limitation of this method is its use of a fixed ideality factor value of 1.3, which makes it more accurate in the vicinity of the maximum power point but less so in other regions (). Achouby et al. have introduced a comprehensive approach that utilizes four analytical formulas, along with certain approximations, to derive the values of four parameters. Additionally, it includes an ideality factor, which is typically defined to be around 1.3 (). Zaimi et al. have developed a highly accurate numerical method for determining the five physical parameters of photovoltaic modules. This method involves fine-tuning the ideality parameter and enhancing a system of nonlinear equations. It's worth noting that the initial estimates for the four physical criteria in their method need to be precise for it to work effectively (). Stornelli et al. have introduced an innovative approach that combines both numerical and well-established methodologies. In this method, it's essential to have consistent initial values for the shunt resistance and ideality factor, which are vital for its effectiveness. This approach presents a new and promising way to determine the five parameters of photovoltaic modules (). This condensed technique allows for the identification of the ideal values for both the ideality parameter and the shunt resistance. Finding the ideal value of the ideality parameter and the shunt resistance is made possible by this condensed technique. The shunt resistance (Rsh) was advised to be disregarded by Sera et al. for the purpose of condensing the five-parameter model to a four-parameter model, proposed ignoring either Rs or Rsh. To make the computation simpler. Other analytical techniques call for more inputs, such as those in the , , and models, which call for extra inputs of estimating or computing Rsho and Rso, or the slopes at both short-circuit (SC) and open-circuit (OC) locations. For the purpose of determining the five parameters, and , also added new coefficients. Another approach was set forth by , and and is based on obtaining a closed-form equation for Rs and then numerically solving this expression for various ideality factor values. Iteratively calculating Rs, Rsh, and Iph is the foundation of Vika’s approach ().

In this study, the single diode model serves as a baseline for extracting parameters related to photo-generated current (Iph), diode saturation current (I0), diode ideality factor (α), series resistance (Rs), and shunt (or parallel) resistance (Rsh) under standard conditions (STC conditions, G = 1000 W/m2, T = 25°C). Furthermore, the extracted parameters are employed to assess and compare six analytical/numerical methods, namely, the Khan et al. method, Blas et al. method, Phang et al. method, Vika method, Cubas et al. method, and Almonacid et al. method. The performance of these methods is analyzed using data sheets provided by the manufacturers of three main photovoltaic modules: the S70 and SM-210W polycrystalline, SP 75 and SPR-230 WHT-I monocrystalline, and Shell ST36 and U-EA110W thin film. All computations are conducted using Matlab software. The two sections of our study are as follows: The first section outlines the methods used by each author to extract the main five parameters of the single diode model, followed by an assessment and comparison of the effectiveness of each approach.

2 Theoretical base of analytical/numerical approaches

PV module performance accuracy is evaluated by examining the I-V characteristic under standard conditions. Various methods exist for modeling PV cells, offering different degrees of approximation to the actual device behavior. According to (), the SDM is the most commonly used approximation for PV module modeling. The SDM enables the extraction of five main parameters: the photo-generated current (Iph), diode saturation current (I0), diode ideality factor (α), series resistance (Rs), and shunt (or parallel) resistance (Rsh). These parameters are interrelated through a specific equation and are illustrated in Figure 1.

FIGURE 1

Where is the thermal voltage of the PV cell, presented by the mathematical Equation 2:

Where is the Boltzmann constant and the elementary or electron charge respectively.

Concerning performance evaluation, the I-V characterization offers valuable insights into the crucial aspects of PV cells. In this context, the precise estimation of the parameters outlined in Equation 1 is detailed through the application of six methods namely, the Khan et al. method (), Blas et al. method (), Phang et al. method (), Vika method (), Cubas et al. method (), and Femia et al. () method. The performance and effectiveness of these approaches are systematically assessed and compared under standard test conditions (STC conditions).

2.1 Khan et al. method

The Khan et al. method is an analytical technique specifically designed for parameter extraction from silicon solar cells. This method focuses on obtaining parameters of PV cells under high radiation conditions. The extraction process involves utilizing initial values from key parameters such as . The extraction is described by the following Equations 37:

2.2 Blas et al. method

The Blas et al. method is an analytical approach employed for the extraction of parameters associated with PV cell behavior. This method relies on a set of experimentally measured voltage-intensity curves conducted under elevated temperature and high radiation levels. The values obtained from these experiments are then used to extract model parameters. Similar to the previously described method, the extraction of parameters from Equation 1 is performed for the specified temperature and solar radiation levels, utilizing initial values of . The mathematical expression for this extraction process is provided as follows (Equations 811):

The calculation of

is determined as per the definition in

Equation 7

. The computational procedure for this method involves several steps:

  • 1. Initial value assignment for .

  • 2. Approximate estimation of both and values.

  • 3. Recalculation of based on the obtained estimates.

  • 4. Iterative repetition of step 3 until converges to a stable value.

The remaining measured data has been discussed in previous sections.

2.3 Phang et al. method

The Phang et al. method is an analytical approach designed for the extraction of the primary five parameters from different single-junction PV cells under standard levels of solar radiation and cell temperature. The estimation of involves implementing a linear fit of the I–V curve around both the short-circuit current and open-circuit voltage, enabling the following (Equations 12, 13):

The parameters of Equation 1 are estimated and given as follows (Equations 1416):

The value of is estimated as previously defined in Equation 6.

2.4 Vika method

The Vika method, also known as the Vika algorithm, is employed for the extraction of parameters from PV cells under standard solar radiation and cell temperature conditions. This method involves fitting the dataset adequately for a specified number of iterations, incorporating variables such as . The initial values of are computed according to the following mathematical expression Equations 1720:

2.5 Cubas et al. method

The Cubas et al. method is an analytical technique designed to extract parameters for various types of PV technology under standard levels of radiation and cell temperature. Similar to the previously discussed methods, the extraction of parameters from Equation 1 is estimated using initial values of . This estimation is accomplished through the following Equations 2123:

2.6 Almonacid et al. Method

The method is a numerical approach employed for the extraction of parameters from PV cells. This method is based on a developed system with non-linear implicit equations derived from Equation 1. Within this system, the primary five unknown parameters are solved using the trust-region optimization algorithm (). This algorithm enhances robustness, especially when dealing with initial values that are far from the solution. Additionally, it can handle cases where the Jacobian becomes singular at a specific iteration. The initial values of are calculated using the following mathematical expression Equations 2427:

In the Vika, Cubas et al., and Almonacid et al. methods, the parameter α remains unidentified. To address this, multiple trials are conducted to solve the equation system within a defined interval [0.1, 1.5], with the mm parameter changing in increments of 0.1. The solutions obtained through these trials are then regarded as the extracted parameters of the I-V curve.

The selected methods were chosen for their established accuracy, reliability, and relevance to the study’s objectives. Each method offers unique strengths that make it suitable for analyzing PV systems. Phang et al. is renowned for accurately modeling I-V and P-V curves near the maximum power point, ensuring reliable performance evaluation across different PV technologies. Blas et al. provides detailed modeling that incorporates environmental factors, making it applicable in diverse climatic conditions. Khan et al. focuses on computationally efficient approaches, ideal for real-time or large-scale applications. Vika contributes valuable insights into system behaviors and environmental interactions, while Cubas et al. emphasizes parameter extraction and practical implementation for accurate and dependable results. Together, these methods offer a comprehensive and balanced evaluation framework.

3 Result and discussion

To examine the outcomes obtained from the analytical/numerical methods discussed in the preceding section, MATLAB software was employed to estimate the primary five parameters extracted using the SDM for three distinct PV module technologies. These technologies include polycrystalline (p-Si), monocrystalline (m-Si), and thin film. The performance of the selected methods is evaluated based on the electrical characteristics at STCs, manufacturing conditions, and the temperature coefficients for the maximum power point (MPP) of the utilized photovoltaic modules, as outlined in Table 1.

TABLE 1

PV modulesPolycrystalline (p-Si)Monocrystalline (m-Si)Thin film
S70SM-210WSP 75SPR-230 WHT-IShell ST36U-EA110W
Open-circuit voltage (V)21.436.121.648.7022.971
Short-circuit current (A)4.77.934.75.992.862.50
Maximum power voltage (V)16.528.9017.641.0015.854.0
Maximum power current (A)4.257.284.265.612.282.04
Maximum power (W)70.125210.47523044.6116
Number of cells in series 3660367242106
temperature coefficient (mA/°C)2423.50.3201.375
temperature coefficient (mV/°C)−76−124−76−160.7−100−276.9

Electrical properties of three PV technologies.

To evaluate and compare the selected methods, various statistical metrics are employed to gauge the precision of our analytical/numerical approaches for extracting parameters from PV modules operating under STCs and manufacturing conditions. The primary statistical metrics include root mean square error (RMSE), mean bias error (MBE), mean absolute percentage error (MAPE), normalized root mean square error (NRMSE), and Absolute error (AE). The mathematical expressions (Equations 2832) for these statistical metrics are provided as follows:

Where are the experimental and computed current, respectively. N is the number of experimental/computed currents. Table 2 below explaining the main differences between the selected statistical metrics (RMSE, MBE, MAPE, NRMSE, and AE) in the context of photovoltaic (PV) module analysis. It is noted that a lower value of these statistical indicator metrics indicates that the selected methods are suitable for extracting parameters from PV modules.

TABLE 2

MetricDefinitionPurpose and key insights
RMSEMeasures the square root of the average of squared differences between experimental and predicted valuesQuantifies the overall prediction error. Smaller values indicate better model performance. Sensitive to large errors, which are given more weight due to squaring
MAPECalculates the average of absolute percentage errors between experimental and predicted valuesEvaluates the accuracy of predictions as a percentage. Useful for comparing across datasets with different scales
Expresses error in relative terms but can be skewed by small experimental values
MBEComputes the mean of the differences between experimental and predicted valuesIndicates whether the model tends to overestimate or underestimate predictions
Positive values suggest overestimation; negative values indicate underestimation
NRMSENormalizes the RMSE by dividing it by the range or mean of experimental data, often expressed as a percentageAssesses prediction error relative to the magnitude of the data, allowing comparisons across datasets with different units or scales
Provides scale-independent insight but depends on the normalization method used
AERepresents the absolute difference between experimental and predicted values for individual data pointsFocuses on localized deviations in the data, offering point-specific error analysis
Useful for identifying outliers but does not provide an aggregate measure of model performance

Comparison of statistical metrics for evaluating PV module parameter extraction methods.

In this section, the analysis of the results obtained from the methods described in the previous section is presented. To facilitate this analysis, several I-V curves spanning the full operational range of the PV technologies were selected. Initially, the recorded I-V curves were analyzed as a function of cell temperature and effective irradiance.

Based on this analysis, the I-V curves presented in Table 3 were identified as a representative dataset for the various photovoltaic (PV) technologies under real operating conditions. Specifically, Curve 1 corresponds to I-V curves near the peak of the data distribution, representing typical working conditions for the PV modules. Further details on these representative I-V curves can be found in Table 3. The criteria above, as well as the I-V curves shown in Table 3, are the same previously considered by the authors to validate the models analyzed in ().

TABLE 3

G (W/m2)Tc (°C)CurvePolycrystalline (p-Si)Monocrystalline (m-Si)Thin film
S70SM-210WSP 75SPR-230 WHT-IShell ST36U-EA 110W
Isc (A)Voc (V)Isc (A)Voc (V)Isc (A)Voc (V)Isc (A)Voc (V)Isc (A)Voc (V)Isc (A)Voc (V)
1,0002514.6921.57.9336.14.721.65.9948.702.8622.92.5071

Effective normal irradiance, cell temperature, and the main electrical parameters derived from the I-V curves used for method evaluation.

The main findings of the analytical/numerical methods are presented in Tables 46, respectively. It is evident from Tables 46 that the predictions of all selected methods are approximately the same for the assessment of . Regarding , Phang et al., Vika, Cubas et al., Khan et al., and Blas et al. methods tend to predict the highest and lowest values. The range values of the considered methods are between 1.264410–09 A and 4.266407 0.10–11 A respectively, and are appropriate for p-Si, m-Si, and thin film technologies. Phang et al. indicates that the prediction results are approximately closer, while the previous methods still lead to showing the highest values (The range value is between 8.3654110–09 A and 3.65467 10−07 A).

TABLE 4

MethodsParameters extractionPolycrystalline (p-Si)
S70SM-210W
Phang et alCurve 1 (A)2.8587 10−103.62167 10−10
(A)4.690547.9376
(Ω)0.3680.237
(Ω)134141
0.9890.984
Blas et alCurve 2 (A)4.8654 10−106.1067 10−10
(A)4.704257.93657
(Ω)0.3480.373
(Ω)141144
1.0051.005
Khan et alCurve 3 (A)6.095687 0.10–107.68517 0.10–10
(A)4.699317.93741
(Ω)0.2930.309
(Ω)145.617139.240
1.0151.015
VikaCurve 4 (A)7.54797 10−109.816493 10−10
(A)4.711107.9315
(Ω)0.3960.187
(Ω)116119
1.0271.026
Cubas et alCurve 5 (A)9.377404 10−101.151188 10−9
(A)4.73655427.93546
(Ω)0.4160.216
(Ω)144.98147.01
1.0351.031
Almonacid et alCurve 6 (A)1.45876 0.10–071.85651.10–07
(A)4.720877.93768
(Ω)0.4710.487
(Ω)146.98149.15
1.3351.330

Parameters extraction of Polycrystalline (p-Si) module.

TABLE 5

MethodsParameters extractionMonocrystalline (m-Si)
SP 75SPR-230 WHT-I
Phang et alCurve 1 (A)2.5016 10−112.9320.25100.3214 10−11
(A)4.705.99
(Ω)0.3500.251
(Ω)335129
0.9811.001
Blas et alCurve 2 (A)2.557 10−112.9431 10−11
(A)4.705.99
(Ω)0.3700.321
(Ω)200137
0.9990.999
Khan et alCurve 3 (A)3.45687 0.10–113.71517 0.10–11
(A)4.75415.99
(Ω)0.3900.390
(Ω)290200
1.0051.005
VikaCurve 4 (A)2.6256 0.10–102.94517 0.10–10
(A)4.70965.9865
(Ω)0.3080.217
(Ω)298120
1.0151.019
Cubas et alCurve 5 (A)2.676229 0.10–102.966407 0.10–10
(A)4.751815.98916
(Ω)0.3000.219
(Ω)200.35119.43
1.0251.029
Almonacid et alCurve 6 (A)5.58467 10−091.698510–08
(A)4.87655.8951
(Ω)0.4500.499
(Ω)221.32200
1.3351.335

Parameters extraction of Monocrystalline (m-Si) module.

TABLE 6

MethodsParameters extractionThin film
Shell ST36U-EA110W
Phang et alCurve 1 (A)1.264410–091.25342 10−11
(A)2.862.4576
(Ω)0.3480.398
(Ω)237310
0.9830 0.999
Blas et alCurve 2 (A)1.9564. 10–091.4845. 10–11
(A)2.862.4976
(Ω)0.4150.572
(Ω)254288
1.0051.005
Khan et alCurve 3 (A)2.45687 0.10–092.08517 0.10–11
(A)2.86122.5165
(Ω)0.26150.469
(Ω)258.768275.43
1.0151.020
VikaCurve 4 (A)2.905128. 10–092.432145. 10–11
(A)2.85982.497765
(Ω)0.398760.266769
(Ω)217235
1.0251.027
Cubas et alCurve 5 (A)3.422709 0.10–092.930493 0.10–11
(A)2.849862.5243
(Ω)0.418660.4376
(Ω)246.88246.99
1.0331.035
Almonacid et alCurve 6 (A)3.65467 10−078.3654110–09
(A)2.86522.5768
(Ω)0.5680.5791
(Ω)286.89294.01
1.3351.335

Parameters extraction of thin film module.

In terms of the ideality factor, the estimated value of α is technology-dependent (p-Si, m-Pi, and thin film). It is clear from the three technologies of PV modules that all proposed methods predict a significant value of the ideality factor except for the last (Almonacid et al.) method. The range values of Phang et al., Vika, Cubas et al., Khan et al., and Blas et al. methods are between 0.983 and 1.335 for the three technologies. The Almonacid et al. method presents almost the same prediction values with a slight variation, which is higher compared to the other methods.

According to Tables 46, Phang et al., Vika, Cubas et al., Khan et al., and Blas et al. methods surpass the Almonacid et al. methods in terms of and . The range value of and for appropriate methods are respectively between 0.187 Ω and 0.510 Ω, 116 Ω and 275 Ω for the three technology.

To evaluate the performance accuracy of each method under standard conditions (GNI = 1000 W/m2, T = 25°C), the figures below depict both I-V and P-V curves for three technologies (polycrystalline, monocrystalline, and thin film). As evident from the figures, all the proposed methods exhibit well-fitted curves, indicating good agreement with experimental data at standard conditions. Notably, the Almonacid et al. method for the three technologies demonstrates the least favorable curve compared to the preceding methods.

Figures 2A–D depicts the I-V and P-V curves of S-70 and SM-210W p-Si PV modules generated by six methods under STC conditions. The predictive methods exhibit high performance across all methods, excluding the Almonacid et al. method. Considering both Table 7 and Figures 2A–D, Phang et al., Blas et al., Khan et al., and Vika methods are deemed appropriate, closely aligning with experimental data. The range values of RMSE (%), MBE (%), and MAPE are 0.080575, −0.95201166, and 0.0722083 for S-70, and 0.20266, 0.98897, and 0.10125 for SM-210W polycrystalline PV modules, respectively. In terms of I-V and P-V curves, the Cubas et al. method exhibits a slight variation compared to the preceding methods. The Almonacid et al. method records the highest value of NRMSE, approximately 0.25906 and 0.48904 for S-70 and SM-210W, respectively. The results obtained from the single diode model for polycrystalline modules indicate that the SM-210W module is more accurate compared to the S-70 PV module.

TABLE 7

APolycrystalline (P-Si)
S-70SM-210W
RMSEMBEMAPENRMSERMSEMBEMAPENRMSE
Phang et al0.05897−1.376390.031500.117940.04209−1.443280.030710.08418
Blas et al0.06378−1.338710.046440.127560.04186−1.342390.045240.08372
Khan et al0.07078−1.067390.056600.141560.05967−1.155080.055640.11934
Vika0.07946−0.918240.074060.158920.06806−0.939220.071570.13612
Cubas et al0.08093−0.571380.092010.161860.14563−0.578000.091010.29126
Almonacid et al0.12953−0.439960.132640.259060.24452−0.454310.131960.48904
BMonocrystalline (M-Si)
SP-75SPR-230 WHT-I
RMSEMBEMAPENRMSERMSEMBEMAPENRMSE
Phang et al0.03874−1,54,8190.027520.077480.02378−1,55,5120.020500.04756
Blas et al0.03772−1,34,4970.045020.075440.02148−1,35960.044500.04296
Khan et al0.03592−1,16,7960.053840.071840.02017−1,19510.049580.04034
Vika0.03144−0,96,2050.071170.062880.01958−0,97,4870.070490.03916
Cubas et al0.13093−0,58,9570.088970.261860.19033−0,65,8740.085020.38066
Almonacid et al0.23043−0,46,8910.129660.460860.20169−0,48,4990.116710.40338
CThin film
ST36U-EA110W
RMSEMBEMAPENRMSERMSEMBEMAPENRMSE
Phang et al0.07432−1.360330.039900.148640.06558−1.369030.038140.13116
Blas et al0.07975−1.232380.049550.15950.06986−1.244200.049030.13972
Khan et al0.08107−0.996860.069730.162140.07571−1.000570.067810.15142
Vika0.07979−0.841970.080170.159580.07885−0.892270.077770.1577
Cubas et al0.09337−0.529200.107900.186740.12930−0.549460.105510.2586
Almonacid et al0.16369−0.402410.138990.327380.28078−0.419100.136910.56156

Statistical metric of all proposed methods for three different technologies.

FIGURE 2

Figures 3A–D provides a comparative analysis of six proposed methods for SP-75 (A-B) and SPR-230 WHT-I (C-D) monocrystalline PV modules under STCs. This figure illustrates the I-V and P-V curves for each method, highlighting the relationship between the prediction methods and experimental data. Referring to Table 7 and Figures 3A–D, the maximum and minimum values of RMSE (%) vary between 0.23043 and 0.03144 for SP-75 monocrystalline silicon and 0.20169 and 0.01958 for SPR-230 WHT-I monocrystalline silicon, respectively. The range values of RMSE (%) suggest that the Phang et al., Blas et al., Khan et al., Vika, and Cubas et al. methods exhibit smaller values than the Almonacid et al. method. In terms of I-V and P-V curves for SP-75 and SPR-230 WHT-I monocrystalline PV modules, the methods depicted in Figures 3A–D closely align with the experimental data, earning recognition as “the most accurate prediction methods” compared to the S-70 and SM-210W polycrystalline PV modules.

FIGURE 3

Figures 4A–D and Table 7 present the prediction accuracy of Phang et al., Blas et al., Khan et al., Vika, Cubas et al., and Almonacid et al. methods under STCs for (A-B) ST36 and (C-D) U-EA110W thin film PV modules. Among the six methods, Phang et al., Blas et al., Khan et al., Vika, and Cubas et al. exhibit the lowest values of MBE (%), MBE (%), MAPE, and RMSE. In terms of I-V and P-V curves, all five methods curve well compared to the Almonacid et al. method and closely align with the experimental data. It is evident from the table that the U-EA110W thin film PV module is more accurate compared to the ST36 thin film PV module. Similarly, consistent with the previous prediction methods, thin film PV modules are recognized as “the least accurate prediction methods” compared to monocrystalline and polycrystalline PV modules.

FIGURE 4

Table 7 provides numerical values for the selected methods, employing computational performance analysis to compare and evaluate their accuracy in predicting experimental data for monocrystalline, polycrystalline, and thin film PV modules. The results show that all methods, except for Almonacid et al., exhibit lower prediction errors for the three PV technologies at the maximum power point. It is important to note that there is no specific criterion for determining which method is more appropriate, as no technique consistently outperforms the others under all conditions. Nevertheless, the Phang et al. method demonstrated a strong alignment with experimental I-V and P-V curves, particularly near the maximum power point, based on lower RMSE, MAPE, and AE values across all PV technologies. Thus, it can be regarded as a suitable approach for all three technologies. It is worth mentioning that our findings align with those studied in (; ), where analytical methods are highlighted for their strong performance, closely matching fitting curves and numerical methods under standard STC conditions.

In this section, the optimal technology for a SDM is illustrated in Figures 5A,B. The figures depict the AE versus PV module voltage for both SP-75 and SPR-230 WHT-I monocrystalline PV modules using various methods. The maximum AE values are 0.004238 for Phang et al., 0.015832 for Almonacid et al., 0.005772 for Vika, 0.001702 for Khan et al., 0.003988 for Blas et al., and 0.004598 for Cubas et al. in the case of SP-75. Similarly, for SPR-230 WHT-I, the maximum AE values are 0.004958 for Phang et al., 0.015381 for Almonacid et al., 0.004286 for Vika, 0.000151 for Khan et al., 0.00204 for Blas et al., and 0.003514 for Cubas et al. method.

FIGURE 5

The parameter extraction for the SDM of three different photovoltaic modules has been presented, with the first five analytical/numerical methods identified as suitable models based on their lower values of several statistical metrics, particularly the Phang et al. and Blas et al. methods. The Phang et al. method, in particular, achieves high accuracy through its use of a linear fit near critical points, such as the short-circuit current and open-circuit voltage. This approach enhances its robustness in minimizing errors across different PV technologies and demonstrates strong agreement between experimental data and the predicted curves. These attributes underscore the high quality and reliability of the results obtained using the Phang et al. method for the proposed technologies.

4 Conclusion

This paper provides a comparative study of various analytical/numerical methods for estimating parameter extraction in different technologies of PV modules. The simulations were conducted under STCs using a single diode model implemented in the Matlab Environment software. Each method was introduced and modeled, and their appropriateness was assessed based on statistical metrics such as RMSE, NRMSE, MAPE, MBE, and AE.

The performance analysis of the six methods was tested across different manufacturers of monocrystalline, polycrystalline, and thin film modules. The I-V and P-V curves of all selected methods exhibited a strong relationship between experimental and predicted values, except for the Almonacid et al. method. Notably, the Phang et al., Blas et al., Khan et al., and Vika methods demonstrated lower error values for various PV modules, including S-70 and SM-210W polycrystalline, SP 75 SPR-230 WHT-I Monocrystalline, and Shell ST36 and U-EA110W thin film modules. Conversely, the Almonacid et al. method yielded the highest error values compared to the other methods.

In conclusion, the Phang et al. method is identified as the most effective for both technologies, providing high-quality results with a strong agreement between experimental and predicted values.

Statements

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

BB: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Conflict of interest

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Abbreviations

  • SDM

    Single diode model

  • STC

    Standard tests conditions

  • RMSE

    Root mean square error

  • MBE

    Mean bias error

  • NRMSE

    Mean absolute percentage error

  • MAPE

    Normalized root mean square error

  • AE

    Absolute error

  • DDM

    Double diode model

  • TDM

    Third diode model

  • p-Si

    Polycrystalline

  • m-Si

    Monocrystalline

  • MPP

    Maximum power point

  • GNI

    Global normal irradiation (W/m2)

Nomenclature

  • SDM

    Single diode model

  • STC

    Standard tests conditions

  • RMSE

    Root mean square error

  • MBE

    Mean bias error

  • NRMSE

    Mean absolute percentage error

  • MAPE

    Normalized root mean square error

  • AE

    Absolute error

  • DDM

    Double diode model

  • TDM

    Third diode model

  • p-Si

    Polycrystalline

  • m-Si

    Monocrystalline

  • MPP

    Maximum power point

  • GNI

    Global normal irradiation (W/m2)

Symbols

  • I-V

    Photocurrent–voltage

  • U-V

    Power-voltage

  • Rsh

    Shunt resistance (Ω)

  • Rs

    Series resistance (Ω)

  • Iph

    Photo-generated current (A)

  • I0

    Diode saturation current

  • α

    Diode ideality factor

  • Thermal voltage

  • Boltzmann constant (1.3806503 × 10−23 J/K)

  • Electron charge (1.60217646 × 10−19C)

  • Short circuit current at STC (A)

  • Open circuit voltage at STC (V)

  • Maximum power current of PV module (A)

  • Maximum power voltage (V)

  • Thermal voltage

References

  • 1

    AlmonacidF.RodrigoP.EnergyE. F. (2016). undefined: determination of the current–voltage characteristics of concentrator systems by using different adapted conventional techniques. Elsevier.

  • 2

    AppelbaumJ.PeledA. (2014). Parameters extraction of solar cells - a comparative examination of three methods. Sol. Energy Mater. Sol. Cells.122, 164173. 10.1016/j.solmat.2013.11.011

  • 3

    BaiJ.LiuS.HaoY.ZhangZ.JiangM.ZhangY. (2014). Development of a new compound method to extract the five parameters of PV modules. Energy Convers. Manag.79, 294303. 10.1016/j.enconman.2013.12.041

  • 4

    BatzelisE. I.PapathanassiouS. A. (2016). A method for the analytical extraction of the single-diode PV model parameters. IEEE Trans. Sustain. Energy.7, 504512. 10.1109/TSTE.2015.2503435

  • 5

    BelliaH.YoucefR.FatimaM. (2014). A detailed modeling of photovoltaic module using MATLAB. NRIAG J. Astron. Geophys.3, 5361. 10.1016/j.nrjag.2014.04.001

  • 6

    Breisnes VikaH. (2014). Modelling of photovoltaic modules with battery energy storage in simulink/matlab with in-situ measurement comparisons.

  • 7

    CannizzaroS.Di PiazzaM. C.LunaM.VitaleG. (2014). “Generalized classification of PV modules by simplified single-diode models,” in IEEE international symposium on industrial electronics (Institute of Electrical and Electronics Engineers Inc.), 22662273. 10.1109/ISIE.2014.6864971

  • 8

    CelikA. N.AcikgozN. (2007). Modelling and experimental verification of the operating current of mono-crystalline photovoltaic modules using four- and five-parameter models. Appl. Energy.84, 115. 10.1016/j.apenergy.2006.04.007

  • 9

    ChanD. S. H.PhillipsJ. R.PhangJ. C. H. (1986). A comparative study of extraction methods for solar cell model parameters. Solid State Electron29, 329337. 10.1016/0038-1101(86)90212-1

  • 10

    ChinV. J.SalamZ.IshaqueK. (2015). Cell modelling and model parameters estimation techniques for photovoltaic simulator application: a review. Appl. Energy154, 500519. 10.1016/j.apenergy.2015.05.035

  • 11

    CubasJ.PindadoS.VictoriaM. (2014). On the analytical approach for modeling photovoltaic systems behavior. J. Power Sources.247, 467474. 10.1016/j.jpowsour.2013.09.008

  • 12

    De BlasM. A.TorresJ. L.PrietoE.GarcíaA. (2002). Selecting a suitable model for characterizing photovoltaic devices. Renew. Energy.25, 371380. 10.1016/S0960-1481(01)00056-8

  • 13

    El AchoubyH.ZaimiM.IbralA.AssaidE. M. (2018). New analytical approach for modelling effects of temperature and irradiance on physical parameters of photovoltaic solar module. Energy Convers. Manag.177, 258271. 10.1016/j.enconman.2018.09.054

  • 14

    FemiaN.PetroneG.SpagnuoloG.VitelliM. (2012). Power electronics and control techniques for maximum energy harvesting in photovoltaic systems (Femia, N. et al; 2013). CRC Press, Taylor and Francis Group.

  • 15

    HusseinA. (2017). A simple approach to extract the unknown parameters of PV modules. Turk. J. Electr. Eng. Comput. Sci.25, 44314444. 10.3906/elk-1703-14

  • 16

    KarmalkarS.HaneefaS. (2008). A physically based explicit $J$– $V$ model of a solar cell for simple design calculations. IEEE Electron Device Lett.29, 449451. 10.1109/LED.2008.920466

  • 17

    KhanF.BaekS. H.ParkY.KimJ. H. (2013). Extraction of diode parameters of silicon solar cells under high illumination conditions. Energy Convers. Manag.76, 421429. 10.1016/j.enconman.2013.07.054

  • 18

    LouzazniM.BelmahdiB. (2022). “Analytical modelling approach of photovoltaic curves: analysis and comparison,” in Lecture notes in networks and systems (Springer Science and Business Media Deutschland GmbH), 562572.

  • 19

    PhangJ. C. H.ChanD. S. H.PhillipsJ. R. (1984). Accurate analytical method for the extraction of solar cell model parameters. Electron. Lett.20, 406408. 10.1049/el:19840281

  • 20

    PowellM. (1968). A Fortran subroutine for solving systems of non-linear algebraic equations. undefined.

  • 21

    SaleemH.KarmalkarS. (2009). An analytical method to extract the physical parameters of a solar cell from four points on the illuminated $J{-}V$ curve. IEEE Electron Device Lett.30, 349352. 10.1109/LED.2009.2013882

  • 22

    SeraD.TeodorescuR.RodriguezP. (2008). “Photovoltaic module diagnostics by series resistance monitoring and temperature and rated power estimation,” in IECON proceedings (industrial electronics conference) (IEEE Computer Society), 21952199.

  • 23

    Sheraz KhalidM.AbidoM. A. (2014). A novel and accurate photovoltaic simulator based on seven-parameter model. Electr. Power Syst. Res.116, 243251. 10.1016/j.epsr.2014.06.010

  • 24

    StornelliV.MuttilloM.de RubeisNardiI. (2019). A new simplified five-parameter estimation method for single-diode model of photovoltaic panels. Energies12, 4271. 10.3390/en12224271

  • 25

    VillalvaM. G.GazoliJ. R.FilhoE. R. (2009). Comprehensive approach to modeling and simulation of photovoltaic arrays. IEEE Trans. Power Electron.24, 11981208. 10.1109/TPEL.2009.2013862

  • 26

    WangG.ZhaoK.ShiJ.ChenW.ZhangH.YangX.et al (2017). An iterative approach for modeling photovoltaic modules without implicit equations. Appl. Energy.202, 189198. 10.1016/j.apenergy.2017.05.149

  • 27

    Yahya-KhotbehsaraA.ShahhoseiniA. (2018). A fast modeling of the double-diode model for PV modules using combined analytical and numerical approach. Sol. Energy.162, 403409. 10.1016/j.solener.2018.01.047

  • 28

    ZaimiM.El AchoubyH.IbralA.AssaidE. M. (2019). Determining combined effects of solar radiation and panel junction temperature on all model-parameters to forecast peak power and photovoltaic yield of solar panel under non-standard conditions. Sol. Energy.191, 341359. 10.1016/j.solener.2019.09.007

Summary

Keywords

parameter extraction, photovoltaic modules, analytical/numerical approaches, single diode model, standard test conditions, I–V characteristics

Citation

Belmahdi B (2025) Comparative study of parameter extractions of photovoltaic modules using analytical and numerical approaches. Front. Energy Res. 13:1501335. doi: 10.3389/fenrg.2025.1501335

Received

24 September 2024

Accepted

03 January 2025

Published

20 January 2025

Volume

13 - 2025

Edited by

Kaiwen Sun, University of New South Wales, Australia

Reviewed by

Yongli Lu, Massachusetts Institute of Technology, United States

Yinyan Liu, University of New South Wales, Australia

Updates

Copyright

*Correspondence: Brahim Belmahdi,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics