Effect of the Dynamic Resistance on the Maximum Output Power in Dynamic Modelling of Photovoltaic Solar Cells ()

Siaka Touré^{}

Solar Energy Laboratory, Félix Houphouët-Boigny University of Abidjan Cocody, Abidjan, Côte d’Ivoire.

**DOI: **10.4236/ojmsi.2022.101003
PDF HTML XML
97
Downloads
430
Views
Citations

Solar Energy Laboratory, Félix Houphouët-Boigny University of Abidjan Cocody, Abidjan, Côte d’Ivoire.

Several studies on PV solar cells are found in the literature which use static models. Those models are mainly one-diode, two-diode or three-diode models. In the dynamic modelling, a variable parallel capacitance is incorporated. Unlike the previous studies which do not clearly establish a relationship between the capacitance and the voltage, in the present paper, the link between the capacitance and the voltage is investigated and established. In dynamic modelling investigated in this paper, the dynamic resistance is introduced in the modelling of the solar cell. It is introduced in the current-voltage characteristic. The value of the dynamic resistance is evaluated at the maximum power point and its effect on the maximum power is investigated. The study shows for the first time, that the dynamic resistance must be introduced in the current-voltage characteristic, because it has an influence on the PV cell output.

Keywords

Photovoltaic Solar Cells, Dynamic Modelling, Dynamic Resistance, Diffusion Capacitance, Transition Capacitance

Share and Cite:

Touré, S. (2022) Effect of the Dynamic Resistance on the Maximum Output Power in Dynamic Modelling of Photovoltaic Solar Cells. *Open Journal of Modelling and Simulation*, **10**, 48-57. doi: 10.4236/ojmsi.2022.101003.

1. Introduction

The world is facing an energy crisis due to the growing energy demand and the environmental damages linked to fossil fuel energy sources. Photovoltaic conversion of solar energy is an alternative source which is clean, sustainable and very promising, because of the researches for higher efficiency solar cells. Several PV models have been studied. Most of them are static models based on a single diode and some parasitic parameters [1] [2] [3]. Other simulations are based on a two-diode model [4] [5] [6] [7]. The one-diode, two-diode and three-diode models have also been used in several works on hybrid Photovoltaic/Thermal panels [8] [9]. Comparatively to the works about the static models, few studies are made about dynamic modelling. The dynamic modelling incorporates a variable parallel capacitance. The different parameters of those models, such as the PV cell impedance, are frequency dependant. In [10], the effect of capacitance on the output characteristics of solar cells is discussed. The model takes the parallel capacitance and the parasite resistances Rsh and Rs into account. The I - V curves are drawn for different frequencies. A split of the characteristics is observed, which is due to the charge or discharge of the internal capacitance. The model is built in a MATLAB environment. In [11], the dynamic model also includes a parallel capacitance as well as Rsh and Rs. The bases of the Impedance Spectroscopy (IS) are explained and applied to detect the Potential-Induced Degradation (PID) of c-Si PV panels. In [12], a dynamic PV model is studied, which incorporates capacitance and reverse-bias characteristics. The model also includes a series inductance. Some cell parameters are obtained from measurements of the impedance at various frequencies. The nonlinear equation of the I - V characteristics is numerically solved using MATLAB Simulink Toolbox. Some work has also been made to study the effect of solar array capacitance on the performance of switching shunt voltage regulator [13]. Some other work has been performed for the characterization of solar cells by Impedance Spectroscopy [14]. The dynamic modelling is also a tool for the study of Negative Capacitance (NC) that has some effects in Organic Light-Emitting Diodes (OLEDs) [15] [16] [17]. This literature review shows that, in the case of dynamic modelling of the PV cell, the parallel capacitance is a very important parameter. But the previous studies do not clearly establish a relationship between the capacitance and the voltage; moreover, in those previous dynamic modellings, where the capacitance effect is investigated, the dynamic resistance is not included in the I - V characteristic. In the present paper, a clear relationship between the capacitance and the voltage is investigated. Moreover, the dynamic resistance R_{d} is introduced in the current-voltage characteristic, which leads to a new expression of this current-voltage characteristic. The value of R_{d} is evaluated at the maximum power point in order to know whether it has an effect or not on the maximum output power of the PV cell.

2. Method

The capacitance and its variation with time as well as the voltage have an importance effect on the dynamic behaviour of a PV cell. Hence, the investigation undertaken in this paper will, in a first time, focus on the capacitance, to have a clear relation of its variation with the voltage. Then, in a second step, the dynamic resistance will be introduced in the dynamic model of the solar cell, leading to its introduction in the current-voltage characteristic.

2.1. The Capacitance of the Solar Cell

A variation of the voltage applied to a p-n junction is linked to a variation of charge, which yields a capacitance. Two different capacitances have to be taken into account: the transition (or space charge) capacitance *C _{t}* and the diffusion capacitance

$C={C}_{d}+{C}_{t}$ (1)

In dynamic regime, the junction schematic representation is given in Figure 1.

Therefore, the effect of the transition and diffusion capacitance on the diode is sketched in Figure 2.

In several papers [12] [18], the junction or transition capacitance, which is due to the charges stored in the depletion region, is expressed as:

${C}_{t}=\frac{{C}_{0}}{\sqrt{1-\frac{{V}_{d}}{\Phi}}}$ (2)

In Equation (2), Φ is the diffusion potential; *V _{d}* =

The diffusion capacitance *C _{d}* is due to the charges stored outside the depletion region, in the neutral regions, when the solar cell is under illumination. In some works [12], it is expressed as:

${C}_{d}=\frac{\tau {I}_{d}}{n{V}_{T}}$ (3)

In Equation (3), *τ* is the mean carrier lifetime. *I _{d}* is the current conducted by the forward biased diode.

Figure 1. Schematic representation of the junction in dynamic regime.

Figure 2. Effect of the transition and diffusion capacitances on the diode.

${V}_{T}=\frac{KT}{q}$ (4)

*C _{d}* dominates for forward biased p-n junctions.

In several works, the Negative Capacitance (NC) is studied. The admittance *Y*(*ω*), for a frequency *ω*, is expressed as [15]:

$Y\left(\omega \right)=\frac{\delta I\left(\omega \right)}{\delta V\left(\omega \right)}$ (5)

for *δI* and
$\delta V~{\text{e}}^{i\omega t}$.

A small *δV *voltage impulse yields a transient harmonic current *δI*. The capacitance is expressed as:

$C\left(\omega \right)=\frac{1}{\omega}\mathrm{Im}\left[Y\left(\omega \right)\right]$ (6)

The conductance *g*(*ω*) is expressed as:

$g\left(\omega \right)=\mathrm{Re}\left[Y\left(\omega \right)\right]$ (7)

In all cases, the capacitance is linked to the voltage. In several works [19] [20], the capacitance is linked to the excess minority carrier’s density under illumination. Let a p-n junction under illumination be considered. In the p base, the excess minority carriers density generated *δ*(*x*) is gotten from the continuity equation. The resolution of this equation yields the photocurrent density *J _{ph}*, the photovoltage

${V}_{ph}={V}_{T}\mathrm{ln}\left[1+\frac{{N}_{b}\delta \left(0\right)}{{n}_{i}^{2}}\right]$ (8)

where
$\delta \left(0\right)=\delta \left(x=0\right)$ ; *N _{b}* is the base doping density;

${C}_{t}=\frac{q}{{V}_{T}}\frac{{n}_{i}^{2}}{{N}_{b}}$ (9)

${C}_{d}=\frac{q}{{V}_{T}}\delta \left(0\right)$ (10)

Hence, from Equation (1), one gets:

$C=\frac{q}{{V}_{T}}\frac{{n}_{i}^{2}}{{N}_{b}}+\frac{q}{{V}_{T}}\delta \left(0\right)$ (11)

Moreover, from Equation (8), for *V* = *V _{ph}*, one gets:

$1+\frac{{N}_{b}\delta \left(0\right)}{{n}_{i}^{2}}={\text{e}}^{\frac{V}{{V}_{T}}}$ (12)

Consequently, one gets:

$\delta \left(0\right)=\frac{{n}_{i}^{2}}{{N}_{b}}\left({\text{e}}^{\frac{V}{{V}_{T}}}-1\right)$ (13)

Then the capacitance, given by Equation (11), becomes:

$C=\frac{q}{{V}_{T}}\frac{{n}_{i}^{2}}{{N}_{b}}+\frac{q}{{V}_{T}}\frac{{n}_{i}^{2}}{{N}_{b}}\left({\text{e}}^{\frac{V}{{V}_{T}}}-1\right)$ (14)

Finally, one gets:

$C=\frac{q}{{V}_{T}}\frac{{n}_{i}^{2}}{{N}_{b}}{\text{e}}^{\frac{V}{{V}_{T}}}$ (15)

Equation (15) shows explicitly, for the first time, the voltage dependence of the capacitance.

2.2. Dynamic Model and Dynamic Resistance Estimate

Very often, the single diode dynamic model is used for the dynamic study of solar cells. In some of the studies, only the forward-bias conducting diode *D _{f}*, the parallel capacitance

In such model, the impedance *Z* of the circuit is given by:

$Z=\mathrm{Re}\left(Z\right)+j\mathrm{Im}\left(Z\right)$ (16)

where:

$\mathrm{Re}\left(Z\right)={R}_{s}+\frac{{R}_{sh}}{1+{\omega}^{2}{C}_{p}^{2}{R}_{sh}^{2}}$ (17)

$\mathrm{Im}\left(Z\right)=-\frac{\omega {C}_{p}{R}_{sh}^{2}}{1+{\omega}^{2}{C}_{p}^{2}{R}_{sh}^{2}}$ (18)

In Equations (17) and (18), *ω* is the frequency. In other works [12], a more extended model is used, which, in addition, includes a series inductance *L _{s} *and a reverse-bias diode

Figure 3. Single diode solar cell dynamic model.

Figure 4. PV cell dynamic model including the reverse-bias diode and the series inductance.

For such model, the impedance Z is expressed as:

$Z=\left({R}_{s}+\frac{{R}_{sh}}{\alpha +1}\right)+j\left({L}_{s}\omega -\frac{\omega {C}_{p}{R}_{sh}^{2}}{\alpha +1}\right)$ (19)

where $\alpha ={\omega}^{2}{C}_{p}^{2}{R}_{sh}^{2}$.

The study of the variation of *Z* with the frequency *ω* is a mean to determine the cell parameters such as *R _{sh}*,

In the present paper, by taking into account Figure 1 and Figure 2, the dynamic resistance *R _{d}* is included in the model, without the reverse-bias diode and the series inductance, as shown in Figure 5.

The introduction of *R _{d}* in the model is justified by the following calculations. Let the capacitance current be considered. The diode voltage

${V}_{d}=V+{R}_{s}I$ (20)

The charge *Q* is:

$Q=C{V}_{d}$ (21)

Therefore, the current *i* is given by:

$i=\frac{\partial Q}{\partial t}=\frac{\partial \left(C{V}_{d}\right)}{\partial t}=C\frac{\partial {V}_{d}}{\partial t}+{V}_{d}\frac{\partial C}{\partial t}$ (22)

$i=C\frac{\partial \left(V+{R}_{s}I\right)}{\partial t}+\left(V+{R}_{s}I\right)\frac{\partial C}{\partial t}$ (23)

Then the current-voltage characteristic is given by [10]:

$I={I}_{ph}-{I}_{s}\left[\mathrm{exp}\left(\frac{V+{R}_{s}I}{n{V}_{T}}\right)-1\right]-\frac{V+{R}_{s}I}{{R}_{sh}}-\left(V+{R}_{s}I\right)\frac{\partial C}{\partial t}-C\frac{\partial \left(V+{R}_{s}I\right)}{\partial t}$ (24)

In Equation (24), *I _{ph}* is the photocurrent,

In the present paper, *R _{d}* is introduced as indicated below:

For a voltage *u* over the capacitance, *C* is expressed as:

$C=\frac{\delta Q}{\delta u}=\frac{\delta it}{\delta u}=\frac{t}{\frac{\delta u}{\delta i}}$ (25)

Since $\delta u={R}_{d}\delta i$

One gets:

Figure 5. Single diode solar cell dynamic model including the dynamic resistance.

$\frac{\partial C}{\partial t}=\frac{1}{{R}_{d}}$ (26)

Moreover:

$\frac{\partial \left(V+{R}_{s}I\right)}{\partial t}=\frac{\partial V}{\partial t}+{R}_{s}\frac{\partial I}{\partial t}=\frac{\partial V}{\partial t}+{R}_{s}\frac{\partial I}{\partial V}\frac{\partial V}{\partial t}=\frac{\partial V}{\partial t}+g{R}_{s}\frac{\partial V}{\partial t}$ (27)

In Equation (27), *g* is the conductance. Hence:

$\frac{\partial \left(V+{R}_{s}I\right)}{\partial t}=\frac{\partial V}{\partial t}\left(1+g{R}_{s}\right)$ (28)

Therefore, Equation (24) becomes:

$I={I}_{ph}-{I}_{s}\left[\mathrm{exp}\left(\frac{V+{R}_{s}I}{n{V}_{T}}\right)-1\right]-\frac{V+{R}_{s}I}{{R}_{sh}}-\frac{V+{R}_{s}I}{{R}_{d}}-C\frac{\partial V}{\partial t}\left(1+g{R}_{s}\right)$ (29)

The solution of the nonlinear Equation (29) yields the output of the solar cell under illumination:

If the term $g{R}_{s}\approx 0$, Equation (29) becomes:

$I={I}_{ph}-{I}_{s}\left[\mathrm{exp}\left(\frac{V+{R}_{s}I}{n{V}_{T}}\right)-1\right]-\frac{V+{R}_{s}I}{{R}_{sh}}-\frac{V+{R}_{s}I}{{R}_{d}}-C\frac{\partial V}{\partial t}$ (30)

3. Results and Discussion

3.1. The Variation of the Capacitance

Let a solar cell be considered with the following parameters:
${V}_{oc}=0.5925\text{\hspace{0.17em}}\text{V}$ ;
${n}_{i}=1.45\times {10}^{10}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\text{cm}}^{-3}$ ;
${N}_{b}={10}^{17}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\text{cm}}^{-3}$ ;
$q=1.602\times {10}^{-19}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{C}$ ;
$T=299.16\text{\hspace{0.17em}}\text{K}$. A plotting of log*C* against *V* is sketched on Figure 6. It is a straight line with slope

$a=\frac{1}{\mathrm{ln}10}\frac{1}{{V}_{T}}=16847$ and intercept $b=\mathrm{log}{C}_{t}=-1.8839$. For the chosen solar

cell, the voltage increases quickly from 0.013 pF (for the short-circuit) to 125.3 μF (for open voltage ${V}_{oc}=0.5925\text{\hspace{0.17em}}\text{V}$ ).

Hence; for any value of the voltage, the capacitance can be calculated and that is useful for the evaluation of the current as shown by Equation (29) and the impedance.

Figure 6. Plotting of logC against the voltage *V*.

3.2. Evaluation of the Dynamic Resistance

The evaluation of the dynamic resistance *R _{d}* was made by using the following reduced equation:

$I={I}_{ph}-{I}_{s}\left[\mathrm{exp}\left(\frac{V+{R}_{s}I}{n{V}_{T}}\right)-1\right]-\frac{V+{R}_{s}I}{{R}_{sh}}-\frac{V+{R}_{s}I}{{R}_{d}}$ (31)

The calculation was made for the following PV cell described in [10]:

${I}_{ph}=1.035\text{\hspace{0.17em}}\text{A}$ ; ${I}_{s}=1.05\times {10}^{-10}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{A}$ ; $n=1.2$ ; ${R}_{s}=0.08\text{\hspace{0.17em}}\Omega $ ; ${R}_{sh}=3.15\text{\hspace{0.17em}}\Omega $

At the maximum power point ( ${V}_{m}=0.5\text{\hspace{0.17em}}\text{V}$ ; ${I}_{m}=0.8\text{\hspace{0.17em}}\text{A}$ ), it is found: ${R}_{d}=23.53\text{\hspace{0.17em}}\Omega $.

The present work gives, for the first time, an estimation of the dynamic resistance.

This estimation of *R _{d}* at the maximum power point shows that, like

Equation (31) shows that the variation Δ*I _{m}* of the current, due to the introduction of

$\Delta {I}_{m}=\frac{{V}_{m}+{R}_{s}{I}_{m}}{{R}_{d}}=0.024\text{\hspace{0.17em}}\text{A}$ (32)

Hence the relative variation of *I _{m}* is
$\frac{\Delta {I}_{m}}{{I}_{m}}=0.03$. It shows that

It should be noted that *R _{d}* is intensity dependent, as shown by Equations (30) and (31). For the short-circuit (

${I}_{sc}\approx \frac{{I}_{ph}}{1+\frac{{R}_{s}}{{R}_{sh}}}$ (33)

As for the open voltage *V _{oc}*, it can be expressed by the semi-empiric relation [21]:

${V}_{oc}\approx {V}_{T}\mathrm{ln}\left(\frac{1}{1.5}\frac{{J}_{ph}}{{10}^{8}\left[\frac{\text{mA}}{{\text{cm}}^{\text{2}}}\right]}\right)+\frac{{E}_{g}}{q}$ (34)

In Equation (34), *J _{ph}* is the photocurrent density,

4. Conclusion

The capacitance of the solar cell was investigated and its variation with the voltage was analyzed. The dynamic resistance *R _{d}*, linked to the parasitic capacitance has been introduced in the dynamic modelling of the solar cell.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

[1] | Fezzani, A., Mohamed, I.H., Drid, S., Zaghba, L., Bouchakour, A. and Benbitour, M.K. (2017) Experimental Investigation of Effects of Partial Shading and Faults on Photovoltaic Modules Performances. Leonardo Electronic Journal of Practices and Technologies, No. 31, 183-200. |

[2] |
Izadian, A., Pourtaherian, A. and Motahari, S. (2012) Basic Model and Governing Equation of Solar Cells Used Bin Power Control Applications. 2012 IEEE Energy Conversion Congress and Exposition, Raleigh, 15-20 September 2012, 1483-1488. https://doi.org/10.1109/ECCE.2012.6342639 |

[3] | Kebir, S.T. (2020) A Summary of Methods to Get Parameters Values of Photovoltaic Cells/Panels. Journal of Renewable Energies, 23, 20-40. |

[4] |
Ishaque, K., Salam, Z. and Taheri, H. (2011) Accurate MATLAB Simulink PV System Simulator Based on a Two-Diode Model. Journal of Power Electronics Korea Science, 11, 179-187. https://doi.org/10.6113/JPE.2011.11.2.179 |

[5] |
Dhar, S., Sridhar, R. and Avasthy, V. (2012) Modeling and Simulation of Photovoltaic Arrays. https://www.iitk.ac.in/npsc/Papers/NPSC2012/papers/12118.pdf |

[6] |
Ishaque, K., Salam, Z., Taheri, H. and Syafaruddin (2011) Modeling and Simulation of Photovoltaic (PV) System during Partial Shading Based on a Two-Diode Model. Simulation Modelling Practice and Theory, 19, 1613-1625. https://doi.org/10.1016/j.simpat.2011.04.005 |

[7] |
Babu, C. and Gurjar, S. (2014) A Novel Simplified Two-Diode Model of Photovoltaic (PV). IEEE Journal of Photovoltaic, 4, 1156-1161. https://doi.org/10.1109/JPHOTOV.2014.2316371 |

[8] |
De Rosa, M., Tagliafico, L.A., et al. (2015) Dynamic Thermal Model for Hybrid Photovoltaic Panels. Energy Procedia, 81, 345-353. https://doi.org/10.1016/j.egypro.2015.12.104 |

[9] |
Zondag, H.A., De Vries, D.W. and Van Halder, W.G.J. (2003) The Yield of Different Combined PV-Thermal Collector Designs. Solar Energy, 74, 253-269. https://doi.org/10.1016/S0038-092X(03)00121-X |

[10] |
Merhej, P., Dallago, E. and Finarelli, D. (2006) Effect of Capacitance on the Output Characteristics of Solar Cells. https://www.ims.unipv.it/FIRB2006/pub/Merhej10.pdf |

[11] |
Oprea, M.I., Sopataru, S.V., Sera, D., Poulsen, P.B. and Thorsteinsson, S. (2017) Degradation in c-Si PV Panels Using Electrical Impedance Spectroscopy. https://doi.org/10.1109/PVSC.2016.7749885 |

[12] |
Kim, K.A., Xu, C.Y., Jin, L. and Krein, P.T. (2013) A Dynamic Photovoltaic Model Incorporating Capacitive and Reverse-Bias Characteristics. IEEE Journal of Photovoltaics, 3, 1334-1341. https://doi.org/10.1109/JPHOTOV.2013.2276483 |

[13] |
Jumar, R.A., Suresh, M.S. and Nagaraju, J. (2006) Effect of Solar Array Capacitance on the Performance of Switching Shunt Voltage Regulator. IEEE Transactions on Power Electronics, 21, 543-548. https://doi.org/10.1109/TPEL.2005.869779 |

[14] |
Mora-Sero, I., Garcia-Belmonte, G., Boix, P.P., Vazquez, M.A. and Bisquert, J. (2009) Impedance Spectroscopy Characterisation of Highly Efficient Silicon Solar Cells under Different Light Illumination Intensities. Energy & Environmental Science, 2, 678. https://doi.org/10.1039/b812468j |

[15] |
Ershov, M., Liu, H.C., Li, L., Buchanan, M., Wasilevski, Z.R. and Jonscher, A.K. (1998) Negative Capacitance Effect in Semiconductor Devices. IEEE Transactions on Electron Devices, 45, 2196-2206. https://doi.org/10.1109/16.725254 |

[16] |
Niu, Q., Wetzelaer, A.H., Blom, P.W.M., et al. (2018) Origin of Negative Capacitance in Bipàolar Organic Diodes. Physical Review Letters, 120, Article ID: 116602. https://doi.org/10.1103/PhysRevLett.120.116602 |

[17] |
Guan, M., Niu, L.T., Zhang, Y., Liu, X.F., Li, Y.Y. and Zeng, Y.P. (2017) Space Charges and Negative Capacitance Effect in Organic Light-Emitting Diodes by Transient Current Response Analysis. RSC Advances, 7, 50598-50602. https://doi.org/10.1039/C7RA07311A |

[18] | Grehant, B. (1990) Cours de Physique des Semiconducteurs. EYROLLES, Paris. |

[19] |
Sahin, G. (2016) Effect of Wavelength on the Electrical Parameters of a Vertical Parallel Junction Silicon Solar Cell Illuminated by Its Rear Side in Frequency Domain. Results in Physics, 6, 107-111. https://doi.org/10.1016/j.rinp.2016.02.003 |

[20] | Selma, M.S., Diatta, I., Traoré, U., Diouf, M.S., Habiboulahh, L., Wade, M. and Sissoko, G. (2018) Diffusion Capacitance in a Silicon Solar Cell under Frequency Modulated Illumination: Magnetic Field and Temperature Effects. Journal of Scientific and Engineering Research, 5, 317-324. |

[21] |
Shah, A.V., Platz, R. and Keppner, H. (1995) Thin-Film Silicon Solar Cells: A Review and Selected Trends. Solar Energy Materials and Solar Cells, 38, 501-520. https://doi.org/10.1016/0927-0248(94)00241-X |

Journals Menu

Contact us

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2023 by authors and Scientific Research Publishing Inc.

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