Determination of the Electrical Parameters of a Solar Cell in Steady State

Abstract

Photovoltaic solar energy can be obtained by using several types of technologies, including silicon solar cells. The characterization of its solar cells makes it possible to know them better. This article presents, on the one hand, the work that has been carried out on these cells. On the other hand, a theoretical study of the cell under illumination using Lamberts W function. On the basis of the electrical parameters provided by the manufacturer, the parameters such as the series and shunt resistances and the electrical quantities such as the photocurrent and the photovoltage, are determined and studied according to the ideality factor of the diode. From the results obtained the shunt resistance increases when the ideality factor increases, the series resistance decreases very weakly.

Share and Cite:

Sadio, O. , Kouyaté, M. , Traoré, P. and Barro, F. (2023) Determination of the Electrical Parameters of a Solar Cell in Steady State. Open Journal of Applied Sciences, 13, 1834-1843. doi: 10.4236/ojapps.2023.1310144.

1. Introduction

Solar photovoltaic has seen many developments and several technologies have emerged to date. The mastery of its technologies therefore becomes a necessity to consider their improvement. Research on different materials is being carried out in this direction. In this research, we distinguish the characterization of these materials in their morphological, optical and electrical aspect. Electrical characterization makes it possible to determine the electrical parameters of solar cells and to better understand their evolution. This characterization can be done in steady state (quasi-steady state), transient state or electrical or optical frequencymodulation. This present study was carried out in steady state. Theoretical studies are characterized by the presentation of a mathematical model which translates the physical behavior of the solar cell. These models are mathematical equations that are most often implicit. These equations are relations between the electrical physical quantities and their related parameters; The resolution of these equations can be done analytically or numerically.

In the literature, many authors have made the electrical characterization of solar cells by application of the analytical method. Among these methods we can cite the work of: Chan et al. [1] who used the single diode method to calculate the five electrical parameters of a solar cell. They first plotted the I-V curve of the cell. Then, from the I-V curve, they determined the short-circuit current, the open-circuit voltage, the current and voltage at the maximum power point, the slope at the open-circuit point Rso and that at the short-circuit Rsho [1] .

Chegaar et al. [2] proposed a simple conductance technique. They have, in two steps, determined the five parameters of a solar cell. First, they determined the shunt conductance Gsh. Then, they calculate the conductance G and determined the ideality factor n of the diode [2] .

Jia and Anderson [3] , considering the one-diode model of the solar cell under illumination, proposed a method for determining the series resistance and ideality factor of the diode. They considered the ideality factor as a variable. From the graph of the I-V characteristic, they fixed two values of n: the first n = 1 corresponding to the operation in open circuit and the second n = 2 corresponding to the operation in short circuit. They calculated the series resistance by making somme approximations. They then calculated the ideality factor of the diode.

Agarwal et al. [4] developed a method to determine the series resistance Rs by tuning the one-diode model of the solar cell, they considered the series resistance to be equal to infinity and the resulting short-circuit current equal to the photogenerated current under an illumination of 700 W/m2. For a low illuminance level, they found a linear dependence on the illuminance level. For an illumination level higher than 700 W/m2, they found a sublinear dependence.

Cowley and Sze [5] proposed a method for determining series resistance from semi-logarithmic I-V characteristics. They considered the one-diode model under darkness. To calculate the series resistance, they took the difference between the semi-logarithmic value and the diffusion line, on the V axis.

The numerical method is also the subject of much work. Among these works, we can mention: Mohammad Rasheed et al. [6] proposed an algorithm allowing the resolution of the Equation (1) of the solar cells. They used the one-diode model by making several iterations with an initial value of 0. The electrical parameters were determined by considering the ambient temperature. Volker Quaschning and Rolf Hanitsh [7] , proposed a general model for the description of solar generators, giving the values of voltages and currents. They used the Newton-Raphson method to characterize a partially shaded solar generator and determine the electrical parameters. LianLian Jiang et al. [8] proposed an optimization technique based on improved adaptive differential evolution (IADE) to determine the electrical parameters of solar cells. For this, they proposed new formulas for the “scaling factor” and the “crossover rate”. N Belhaouas et al. [9] proposed a numerical method based on Matlab-Simulink for the characterization of solar cells/modules under the influence of environmental parameters such as irradiance level, temperature and surface conditions. They used the electric diode model.

The equations that govern the evolution of the electrical quantities of solar cells are most often implicit. Their resolution most often requires the application of a numerical method. Numerical methods give an approximation of the calculated quantity. They require the knowledge of several variables, and they use algorithms to determine the electrical parameters. With analytical methods, the unique knowledge of the parameters given by the manufacturer makes it possible to determine the other parameters and electrical quantities of the solar cells.

In the present study, an analytical method is used for the determination of the electrical parameters of solar cells.

Unlike numerical methods, the approach consists, from the unique knowledge of the parameters given by the manufacturer, in determining the other parameters and electrical quantities of the solar cells such as the photocurrent, the short-circuit current, the series resistance and the shunt resistor. The mathematical formulation that has been proposed has been obtained by considering the one-diode model of the solar cell under illumination. Unlike the analytical methods presented above, this mathematical formulation was solved by using Lambert's W function. The expressions of the different parameters and quantities were expressed as a function of the ideality factor α of the diode and simulations were made. These simulations made it possible to determine the typical value of the diode factor for the studied solar cells. The results obtained are then compared to those available in the literature.

2. Mathematical Formulation

The structure of a solar cell is formed by the superposition of four layers which are in contact [10] . Contacts between the four zones result in an ohmic resistance Rs. The series resistance Rs characterizes the resistive effects of the solar cell [11] .The current generatedinside a solar cell is not completely collected, there are leakage currents in the solar cell. The effects of these leakage paths are characterized by the shunt resistance Rsh [12] . In steady state, the electrical equivalent circuit of solar cell is then given by the well-known one-diode model presented in Figure 1.

The current I is given by:

I = I p h I 0 [ exp ( V + I R s α V T ) 1 ] V + I R s R s h (1)

In the Equation (1), Iph is the photocurrent, I0 is the reverse saturation current, VT the thermal andαthe ideality factor or quality factor of the diode.

Figure 1. One diode model of the solar cell.

Since the objective of this work is to determine the electrical parameters of solar cells, thereafter only the following parameters will be examined: the series resistance Rs, the shunt resistance Rsh, the reverse saturation current I0 and the photocurrent Iph.

When operating in short-circuit condition, the voltage V is removed from Equation (1) and the current I is replaced by the short-circuit current Isc:

I s c = I p h I 0 [ exp ( I s c R s α V T ) 1 ] I s c R s R s h (2)

The second term of the right-hand side (diode current) is negligible. Equation (2) becomes:

I s c = I p h I s c R s R s h (3)

From Equation (3), we get the expression for the photocurrent:

I p h = I s c ( R s h + R s R s h ) (4)

At open circuit the current becomes zero and the voltage V becomes equal to the open circuit voltage Voc:

0 = I p h I 0 [ exp ( V o c α V T ) 1 ] V o c R s h (5)

Equation (5) can be reduced if we consider V o c α V T :

exp ( V o c α V T ) 1 exp ( V o c α V T ) (6)

We then have:

I 0 = ( R s h + R s ) I s c V o c R s h exp ( V o c α V T ) (7)

In addition to the short-circuit current Isc and the open-circuit voltage Voc, the I-V curve is also characterized by the maximum power point. This point corresponds to the maximum power that the cell can deliver. At this point, in Equation (3), we replace I by Im and V by Vm:

I m = I p h I 0 [ exp ( V m + I m R s α V T ) 1 ] V m + I m R s R s h (8)

The power delivered by the solar cell is: P = I V

Its derivative with respect to the voltage is then:

P V = V I V + I (9)

At this point, the derivative of the power with respect to the voltage is zero

From Equation (9) setting the maximum power point we have:

( I V ) I m , V m = I m V m (10)

Considering Equation (1) and taking the derivative with respect to V, we get

d I d V = I 0 ( 1 α V T + R s α V T d I d V ) exp ( V + I R s α V T ) ( 1 R s h + R s R s h d I d V ) (11)

At the maximum power point, considering Equation (11) we then have:

I m V m = I 0 α V T ( 1 I m V m R s ) exp ( V m + I m R s α V T ) 1 R s h ( 1 I m V m R s ) (12)

Considering Equations (10) and (11), Equation (12) can be rewritten as:

α V T V m ( 2 I m I s c ) [ V m I s c + V o c ( I m I s c ) ] [ V m R s ] α V T [ V m I s c V o c I m ] = exp [ V m + I m R s V o c α V T ] (13)

Rearranging Equation leads to:

V m ( 2 I m I s c ) V m I s c + V o c ( I m I s c ) exp ( V o c 2 V m α V T + V m I s c V o c I s c + V o c ( I m I s c ) I m V m ) = ( I m R s V m α V T + V m I s c V o c I m V m I s c + V o c ( I m I s c ) ) exp ( I m R s V m α V T + V m I s c V o c I m V m I s c + V o c ( I m I s c ) ) (14)

Let x = I m R s V m α V T + V m I s c V o c I m V m I s c + V o c ( I m I s c ) and y = V m ( 2 I m I s c ) V m I s c + V o c ( I m I s c ) exp ( V o c 2 V m α V T + V m I s c V o c I s c + V o c ( I m I s c ) I m V m )

We then have:

y = x e x (15)

Considering the Lambert function defined as z = W ( z ) e W ( z ) , z being any complex number, we obtain:

x = f 1 ( y ) = W ( y ) (16)

This implies the application of the secondary branch W-1(x).Substituting x and y in Equation (16), we obtain:

I m R s V m α V T + V m I s c V o c I m V m I s c + V o c ( I max I s c ) = W 1 ( V m ( 2 I m I s c ) V m I s c + V o c ( I m I s c ) exp ( V o c 2 V m α V T + V m I s c V o c I s c + V o c ( I m I s c ) I m V m ) ) (17)

Let’s set: A = α V T I max , B = V m ( 2 I m I s c ) V m I s c + V o c ( I m I s c ) , C = V m I s c V o c I m V m I s c + V o c ( I m I s c ) , D = V o c 2 V m α V T and E = V max α V T . We can then rewrite Equation (17) and deduce the series resistance in the form:

R s = A [ W ( B exp ( D + C ) ) C + E ] (18)

Shunt resistance is derived as:

R s h = ( V m I m R s ) ( V m R s ( I s c I m ) α V T ) ( V m I m R s ) ( I s c I m ) α V T I m (19)

Knowing the series resistance Rs and the shunt resistance Rsh will make it possible to calculate the reverse saturation current I0 from Equation (7) and the photocurrent Iph from Equation (4).

3. Results and Discussion

To perform the simulation, the values of the four parameters Isc, Voc, Im and Vm from Table 1 were used [13] .

The results obtained are presented on Table 2.

The results that have been obtained are based on the single-diode model of the solar cell. In this model, the electrical parameters depend on the ideality factor of

Table 1. Electrical characteristics for blue color and gray color cells.

Table 2. Values of the four electrical parameters of the blue cell and the gray cell.

the diode. In this present study we varied the ideality factor in the range of 1 to 1.5. This choice is explained on the one hand by the fact that for values of α < 1, the contact would be non-ohmic. Authors have reported that this means the presence of an additional diode at the rear contact [14] . On the other hand, this choice is explained by the fact that the I-V characteristics of a solar cell must agree with the model of current transport in the solar cell [15] . In the literature, authors have presented methods for determining the ideality factor. Several methods can be found in [16] .

Figures 2(a)-(d) and Figures 3(a)-(d) present the variation of the parameters according to the ideality factor of the diode.

(a) (b) (c) (d)

Figure 2. Variation of a) Rs, b) Rsh, c) Io d) Rs and Rsh both, for different values of α of the blue cell.

(a) (b) (c) (d)

Figure 3. Variation of (a) Rs, (b) Rsh, (c) Io and (d) Rs and both Rsh, for different values of α for the gray cell.

Figure 2(a), Figure 2(b) and Figure 2(c) of the blue cell and Figure 3(a), Figure 3(b) and Figure 3(c) show the respective variations of Rs, Rsh and I0 when the diode factor varies. For α between 1 and 1.4, the shunt resistance varies slightly. The same is true for the reverse saturation current. The shunt resistance and the reverse saturation current vary greatly when α is above 1.4. From the results obtained on Rs and Rsh, it appears that the shunt resistance is relatively of the order of 1000 of the series resistance [17] .

With Figure 2(a) and Figure 2(b) of the blue cell and Figure 3(a) and Figure 3(b) of the gray cell, it is possible to determine the value of α for which the shunt resistance Rsh would be a proportion of the order of 1000 the series resistance Rs

Table 3. Parameters values at the point α = αP.

for the blue cell and of the order of 100 for the gray cell. This point is then used to calculate the value of I0 and Iph corresponding to Rsh = 1000 × Rs for the blue cell and Rsh = 1000 × Rs for the gray cell. Figure 2(c) and Figure 3(d) give the value of α corresponding to this point. Values of 1209 and 1320 were found for the blue cell and the gray cell respectively. Table 3 gives the value of the four parameters at this point.

From the results obtained, it appears that the photocurrent varies very slightly when the ideality factor varies. These variations being of the order of 103, it can therefore be said that the photocurrent is equal to the short-circuit current according to the data in Table 2.

4. Conclusion

This article presents an analytical method, based on Lambert’s W function, to determine the electrical parameters of solar cells on the basis of information provided by the manufacturer. The series resistance, the shunt resistance, the reverse saturation current of the diode and the photocurrent were determined and compared with those obtained in the literature with good agreements. In order to better understand these parameters, the influence of environmental parameters such as temperature and level of illumination on these parameters can be studied in perspective.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Chan, D.S.H. and Phang, J.C.H. (1984) A Method for the Direct Measurement of Solar Cell Resistance. IEEE Transactions on Electron Devices, 31, 381-383.
https://doi.org/10.1109/T-ED.1984.21532
[2] Chegaar, M., Ouennoughia, Z. and Guechi, F. (2004) Extracting DC Parameters of Solar Cells under Illumination. Vacuum, 75, 367-372.
https://doi.org/10.1016/j.vacuum.2004.05.001
[3] Jia, Q. and Anderson, W.A. (1988) A Novel Approach for Evaluating the Series Resistance of Solar Cells. Solar Cells, 25, 311-318.
https://doi.org/10.1016/0379-6787(88)90069-5
[4] Agarwal, S.K., Muralidharan, R., Agarwala, A., Tewary, V.K. and Jain, S.C. (1981) A New Method for the Measurement of Series Resistance of Solar Cells. Journal of Physics D: Applied Physics, 14, 1643-1646.
https://doi.org/10.1088/0022-3727/14/9/011
[5] Cowley, A.M. and Sze, S.H. (1965) Surface States and Barrier Height of Metal—Semiconductor Systems. Journal of Applied Physics, 30, 3212-3220.
https://doi.org/10.1063/1.1702952
[6] Rasheed, M., Alabdali, O. and Shihab S. (2021) A New Technique for Solar Cell Parameters Estimation of the Single-Diode Model, Journal of Physics: Conference Series, 1879, Article ID: 032120.
https://doi.org/10.1088/1742-6596/1879/3/032120
[7] Quaschning V. and Hanitsh R. (1996) Numerical Simulation of Current-Voltage Characteristics of Photovoltaic Systems with Shaded Solar Cells, Solar Energy, 56, 513-520.
https://doi.org/10.1016/0038-092X(96)00006-0
[8] Jiang, L.L., Maskell, D.L., Patra, J.C. (2013) Parameter Estimation of Solar Cells and Modules Using an Improved Adaptive Differential Evolution Algorithm. Applied Energy, 112, 185-193.
https://doi.org/10.1016/j.apenergy.2013.06.004
[9] Belhaouas, N., Ait Cheikh, M.S., Malek. A. and Larbes, C. (2013) Matlab-Simulink of Photovoltaic System Based on a Two-Diode Model Simulator with Shaded Solar Cells, Revue des Energies Renouvelables, 16, 65-73.
[10] Tawada, Y., Tsuge, K., Kondo, M., Okamoto, H. and Hamakawa, Y. (1982) Properties and Structure of a-SiC: H for High-Efficiencya-Si Solar Cell, Journal of Applied Physics, 53, 5273-5281.
https://doi.org/10.1063/1.331363
[11] Araujo, G.L. and Sanchez, E. (1982) A New Method for Experimental Determination of the Series Resistance of Solar Cell. IEEE Transactions on Electron Devices, 29, 1511-1513.
https://doi.org/10.1109/T-ED.1982.20906
[12] Chan, D.S.H. and Phang, J.C.H. (1984) A Method for the Direct Measurement of Solar Cell Shunt Resistance, IEEE Transactions on Electron Devices, 31, 381-383.
https://doi.org/10.1109/T-ED.1984.21532
[13] Cubas, J., Pindado, S. and De Manuel, C. (2014) Explicit Expressions for Solar Panel Equivalent Circuit Parameters Based on Analytical Formulation and the Lambert W-Function. Energies, 7, 4098-4115.
https://doi.org/10.3390/en7074098
[14] Bowden, S. and Rohatgi, A. (2001) Rapid and Accurate Determination of Series Resistance and Fill Factor Losses in Industrial Silicon Solar Cells. Proceeding of 17th European Photovoltaic Solar Energy Conference, Munich, 22-26 October, 2001, 1802-1806.
https://repository.gatech.edu/server/api/core/bitstreams/c6ae9f4e-f5d4-4a56-b871-7f8f1f113827/content
[15] Hamdy, M.A. and Call, R.L. (1987) The Effect of the Diode Ideality Factor on the Experimental Determination of Series resistance of Solar Cells. Solar Cells, 20, 119-126.
https://doi.org/10.1016/0379-6787(87)90036-6
[16] Bashahu, M. and Nkundabakura, P. (2007) Review and Tests of Methods for the Determination of the Solar Cell Junction Ideality Factors. Solar Energy, 81, 856-863.
https://doi.org/10.1016/j.solener.2006.11.002
[17] Aichouba, M.E.A. and Rahli, M. (2019) Solar Cell Parameters Extraction Optimization Using Lambert Function, Przegląd Elektrotechniczny, 95, 227-231.
https://doi.org/10.15199/48.2019.04.43

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.