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


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.

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 Rd is introduced in the current-voltage characteristic, which leads to a new expression of this current-voltage characteristic. The value of Rd 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 Ct and the diffusion capacitance Cd. Cd and Ct are parallel, so that the cell capacitance is expressed as:

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 = C 0 1 V d Φ (2)

In Equation (2), Φ is the diffusion potential; Vd = U is the voltage over the diode; C0 is the voltage for U = 0. The voltage V = Φ − U is the potential between the limits of the depletion region. The transition capacitance dominates for small voltages.

The diffusion capacitance Cd 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 = τ I d n V T (3)

In Equation (3), τ is the mean carrier lifetime. Id is the current conducted by the forward biased diode. VT is the thermal voltage and is expressed as:

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

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

V T = K T q (4)

Cd 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 ( ω ) = δ I ( ω ) δ V ( ω ) (5)

for δI and δ V ~ e i ω t .

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

C ( ω ) = 1 ω Im [ Y ( ω ) ] (6)

The conductance g(ω) is expressed as:

g ( ω ) = Re [ Y ( ω ) ] (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 Jph, the photovoltage Vph and the capacitance. The photovoltage Vph is expressed as:

V p h = V T ln [ 1 + N b δ ( 0 ) n i 2 ] (8)

where δ ( 0 ) = δ ( x = 0 ) ; Nb is the base doping density; ni is the intrinsic carriers density. The transition capacitance Ct and the diffusion capacitance Cd are expressed as [19] [20]:

C t = q V T n i 2 N b (9)

C d = q V T δ ( 0 ) (10)

Hence, from Equation (1), one gets:

C = q V T n i 2 N b + q V T δ ( 0 ) (11)

Moreover, from Equation (8), for V = Vph, one gets:

1 + N b δ ( 0 ) n i 2 = e V V T (12)

Consequently, one gets:

δ ( 0 ) = n i 2 N b ( e V V T 1 ) (13)

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

C = q V T n i 2 N b + q V T n i 2 N b ( e V V T 1 ) (14)

Finally, one gets:

C = q V T n i 2 N b e 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 Df, the parallel capacitance Cp (which is voltage dependant), the shunt resistance Rsh and the series resistance Rs are taken into account, as shown in Figure 3 [10] [11].

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

Z = Re ( Z ) + j Im ( Z ) (16)


Re ( Z ) = R s + R s h 1 + ω 2 C p 2 R s h 2 (17)

Im ( Z ) = ω C p R s h 2 1 + ω 2 C p 2 R s h 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 Ls and a reverse-bias diode Dr, with a breakdown voltage Vbd, as shown in Figure 4.

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 = ( R s + R s h α + 1 ) + j ( L s ω ω C p R s h 2 α + 1 ) (19)

where α = ω 2 C p 2 R s h 2 .

The study of the variation of Z with the frequency ω is a mean to determine the cell parameters such as Rsh, Rs, Cp and Ls.

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

The introduction of Rd in the model is justified by the following calculations. Let the capacitance current be considered. The diode voltage Vd is expressed as:

V d = V + R s I (20)

The charge Q is:

Q = C V d (21)

Therefore, the current i is given by:

i = Q t = ( C V d ) t = C V d t + V d C t (22)

i = C ( V + R s I ) t + ( V + R s I ) C t (23)

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

I = I p h I s [ exp ( V + R s I n V T ) 1 ] V + R s I R s h ( V + R s I ) C t C ( V + R s I ) t (24)

In Equation (24), Iph is the photocurrent, Is is the diode saturation current and n is the diode ideality factor.

In the present paper, Rd is introduced as indicated below:

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

C = δ Q δ u = δ i t δ u = t δ u δ i (25)

Since δ u = R d δ i

One gets:

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

C t = 1 R d (26)


( V + R s I ) t = V t + R s I t = V t + R s I V V t = V t + g R s V t (27)

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

( V + R s I ) t = V t ( 1 + g R s ) (28)

Therefore, Equation (24) becomes:

I = I p h I s [ exp ( V + R s I n V T ) 1 ] V + R s I R s h V + R s I R d C V t ( 1 + g R s ) (29)

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

If the term g R s 0 , Equation (29) becomes:

I = I p h I s [ exp ( V + R s I n V T ) 1 ] V + R s I R s h V + R s I R d C V t (30)

3. Results and Discussion

3.1. The Variation of the Capacitance

Let a solar cell be considered with the following parameters: V o c = 0.5925 V ; n i = 1.45 × 10 10 cm 3 ; N b = 10 17 cm 3 ; q = 1.602 × 10 19 C ; T = 299.16 K . A plotting of logC against V is sketched on Figure 6. It is a straight line with slope

a = 1 ln 10 1 V T = 16847 and intercept b = 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 o c = 0.5925 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 Rd was made by using the following reduced equation:

I = I p h I s [ exp ( V + R s I n V T ) 1 ] V + R s I R s h V + R s I R d (31)

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

I p h = 1.035 A ; I s = 1.05 × 10 10 A ; n = 1.2 ; R s = 0.08 Ω ; R s h = 3.15 Ω

At the maximum power point ( V m = 0.5 V ; I m = 0.8 A ), it is found: R d = 23.53 Ω .

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

This estimation of Rd at the maximum power point shows that, like Rsh, R d R s .

Equation (31) shows that the variation ΔIm of the current, due to the introduction of Rd in the solar cell characteristics, at the maximum power point is:

Δ I m = V m + R s I m R d = 0.024 A (32)

Hence the relative variation of Im is Δ I m I m = 0.03 . It shows that Rd leads to a 3% reduction of the current, at the Maximum Power Point (MPP).

It should be noted that Rd is intensity dependent, as shown by Equations (30) and (31). For the short-circuit (V = 0), it is higher, so that its influence on the short-circuit intensity Isc is very weak. The short-circuit intensity Isc is linked to the photocurrent, Iph by the following relation:

I s c I p h 1 + R s R s h (33)

As for the open voltage Voc, it can be expressed by the semi-empiric relation [21]:

V o c V T ln ( 1 1.5 J p h 10 8 [ mA cm 2 ] ) + E g q (34)

In Equation (34), Jph is the photocurrent density, Eg is the band gap and q is the electron charge. Hence the introduction of Rd which leads to a very small decrease of Isc results in a weak decrease of Voc. But as shown above, its influence on Im at the maximum power point is not negligible. Consequently, the dynamic resistance Rd has an influence on the maximum output power P m = I m V m , while its influence on Isc and Voc is negligible.

4. Conclusion

The capacitance of the solar cell was investigated and its variation with the voltage was analyzed. The dynamic resistance Rd, linked to the parasitic capacitance has been introduced in the dynamic modelling of the solar cell. Rd has been introduced in the current-voltage characteristics and its value has been estimated at the maximum power point. It has been shown that it has a significant influence on the maximum power current Im and consequently on the maximum power Pm. Hence, the introduction of the dynamic resistance Rd in the I - V characteristics is quite relevant.

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.
[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.
[5] Dhar, S., Sridhar, R. and Avasthy, V. (2012) Modeling and Simulation of Photovoltaic Arrays.
[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.
[7] Babu, C. and Gurjar, S. (2014) A Novel Simplified Two-Diode Model of Photovoltaic (PV). IEEE Journal of Photovoltaic, 4, 1156-1161.
[8] De Rosa, M., Tagliafico, L.A., et al. (2015) Dynamic Thermal Model for Hybrid Photovoltaic Panels. Energy Procedia, 81, 345-353.
[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.
[10] Merhej, P., Dallago, E. and Finarelli, D. (2006) Effect of Capacitance on the Output Characteristics of Solar Cells.
[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.
[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.
[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.
[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.
[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.
[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.
[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.
[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.
[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.

Copyright © 2023 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.