JREE1994601719948600
JREE1994601719948600
JREE1994601719948600
Journal of Renewable
Energy and Environment
Journal Homepage: www.jree.ir MERC
Research Article
Comparative Analysis of Dynamic and Steady State Performances of Hill Climbing and
Incremental Conductance MPPT Controllers for PV Systems
Noussaiba Mennaia*, Ammar Medoueda, Youcef Soufib
a Department of Electrical Engineering, LES Laboratory, University of 20 August 1955, Skikda, Algeria.
b Deportment of Electrical Engineering, LABGET Laboratory, Larbi Tebessi University, Tebessa, Algeria.
PAPER INFO
A B S T R A C T
Paper history:
Received 25 December 2023 The integration of photovoltaic (PV) solar energy into the utility grid is expanding progressively to meet
Revised 06 May 2024 increasing energy demand. A crucial aspect of optimizing the output of PV systems involves implementing
Accepted 23 June 2024 efficient maximum power point tracking (MPPT) controllers, necessary due to the nonlinear characteristics of
these systems. This study conducts a simulation-based comparative analysis of two prominent MPPT techniques:
hill climbing (HC) and incremental conductance (INC) methods. The emphasis is on the dynamic response and
Keywords: steady-state efficiency of these controllers. Using a modeled 500 kW PV array alongside both MPPT techniques
Photovoltaic Energy
MPPT
and other DC stage components, simulation tests were conducted in Matlab/Simulink under standard test
Hill Climbing conditions (STC) and varying meteorological conditions. The simulation results indicate that both techniques
Incremental Conductance successfully tracked the MPP. However, the INC algorithm exhibits superior speed, precision, and efficiency,
PV System Modeling particularly in scenarios involving sudden fluctuations in irradiance and temperature. Furthermore, investigating
the effect of perturbation step size on dynamic response and steady-state efficiency provided valuable insights
for enhancing MPPT controller performance.
https://doi.org/10.30501/jree.2024.432251.1788
adoption of a sophisticated maximum power point tracking
1. INTRODUCTION1
(MPPT) controller is essential. This controller continuously
Utilizing renewable energy sources (RES) is becoming a more adjusts the operating point to track the PV panel's maximum
appealing option to meet rising electricity demand (Geravandi power point (MPP) by controlling the DC-DC converter, even
& Moradi CheshmehBeigi, 2022). In this context, photovoltaic in the face of changing irradiance and temperature conditions
(PV) systems—which transform solar energy into electrical (Jately et al., 2021; Jha, 2022). However, as reviewed by (Al
power—have become quite well-known. The output power of Garni et al., 2024), their effectiveness in power point tracking
these systems is quite versatile; it can range from a few under varying conditions remains a concern, particularly for
milliwatts in applications such as scientific calculators to large-scale PV power plants. On the other hand, the design of
megawatts in solar farms (Jately et al., 2021), where PV arrays an MPPT controller must begin first by designing the boost
are formed by combining numerous modules. Solar farms, converter components, namely the inductance, the input, and
along with distributed generation (DG) and microgrid (MG) the output capacitance, then choosing the MPPT algorithm.
systems, play important roles in harnessing RES to satisfy the Many works have been published in this regard. The authors
growing demand for power. Indeed, despite the advantages of (Srdic et al., 2012) implemented the incremental conductance
reducing dependence on fossil fuels and the low-cost (INC) MPPT algorithm on a simplified PV string model and a
maintenance of PV systems, there are challenges related to 1 kW boost converter prototype; however, they did not provide
efficiency and the impact of external factors. The practical any comparison with other MPPT techniques, and the effect of
efficiency of PV modules is approximately 20%, and the temperature variation is missing. While S.V. Rajani and V.J.
production of energy is directly influenced by temperature and Pandya (Rajani & Pandya, 2015) focused their study only on
solar irradiation variations, significantly affecting the output comparing INC and perturb and observe (P&O) MPPT
power of PV arrays (Kumar et al., 2023; Priambodo et al., techniques via simulations in Matlab/Simulink, but did not
2022). Moreover, the operational point of a PV panel is located accurately model PV arrays or boost converter designs. On the
in the nonlinear region of current-voltage (I-V) characteristics. other hand, a model for DC components of a grid-connected PV
To address these challenges and optimize energy output, the system was developed in (Motan et al., 2018), including the PV
Finally, both parameters of I0 and IPV are used to obtain Eq. (6)
of the model current I m
𝑁𝑠𝑠
𝑉+𝐼.𝑅𝑠 ( )
𝑁𝑝𝑝
𝐼𝑚 = 𝐼𝑝𝑣 . 𝑁𝑝𝑝 − 𝐼0 . 𝑁𝑝𝑝 [𝑒𝑥𝑝 ( ) − 1] (6)
𝑎.𝑉𝑡 .𝑁𝑠𝑠
(𝑅𝑝 +𝑅𝑠 )
𝐼𝑝𝑣,𝑛 = 𝐼𝑠𝑐,𝑛 (10)
𝑅𝑝
(b)
Figure 4. (a) I-V and P-V curves at different solar irradiation; (b) I-V
and P-V curves at different temperatures.
Pulse width modulation (PWM) is employed to control the
DC-DC converter by modulating the on/off switches at a high
frequency. The input signals of PWM are generated using the
MPPT technique. For continuous current conduction, the
following formulas are utilized to calculate the boost
parameters: the inductance (L), the input capacitor (Ci), and the
output capacitor (Co) (Ayop & Tan, 2018; Kumar et al., 2023).
The results are summarized in Table 2.
𝑉𝑝𝑣 𝑉𝑝𝑣
𝑉𝑜 = →𝐷 =1− (11)
(1−𝐷) 𝑉𝑜
𝑉𝑜 2
𝑅= (12)
𝑃𝑜
𝑉𝑜
Figure 3. Iterative method flowchart for identifying Rs and Rp. 𝐼𝐿 = (13)
𝑅(1−𝐷)
𝐷.𝑉𝑝𝑣
The Simulink model of the PV array was designed using the 𝐿= (14)
𝑓𝑠𝑤 .𝛥𝐼𝐿
previous equations along with the data from the manufacturer 𝐷
𝐶𝑖 = ∆𝑉𝑚𝑝 2
(15)
SCHOTT Solar of the Perform Poly 230 PV module ("Schott 8𝐿( )𝑓𝑠𝑤
𝑉𝑚𝑝
𝐷
Datasheet parameters at 25 Adjsuted model parameters 𝐶𝑜 = ∆𝑉𝑜 (16)
𝑅( )𝑓
𝑉𝑜 𝑠𝑤
°C, 1.5AM, 1000 W/m2 Table 2. Parameters of DC-DC converter.
Imp 7.66 A
Rs 170 mΩ
V mp 30 V V pv (V mp.Nss) 660 V
P max,e 229.8 W Rp 186.690053 Ω Vo 1000 V
Isc 8.33 A P max,m 229.8 W P pv ≈ P o 500 kW
Voc,n 36.9 V I0 8.394565.10 -8 A fsw 10 kHz
Kv -0.32 %/K Ipv 8.337632 A ∆IL/IL 15 %
Ki 0.04 %/K ∆V mp/V mp = ∆Vo/Vo 1%
a 1.3
Ns 60 D 0.34
PV array parameters R 1.978 Ω
Nss 22 IL 766 A
Npp 100 L 196 µH
Perform Polycrystalline Solar Modules ", 2011) and the values Ci 218 µF
of Rs and Rp. All parameters of the 500 kw PV array used in Co 1719 µF
this work are shown in Table 1. The I-V and P-V characteristics
of this PV array are shown in Figure 4. 2.2. MPPT Techniques
Table 1. Parameters of the Schott Perform Poly 230 based PV array. In Figure 4, the output power of the PV array varies with
irradiation level and temperature, displaying a distinct
Maximum Power Point (MPP) characterized by V mp and Imp.
While the exact position of this point is not specified, it must
be identified by implementing an MPPT algorithm to ensure
optimal operation of the PV array at its MPP. There are various
methods available in the literature to track the MPP (Mohamed
& Abd El Sattar, 2019). This study compares two MPPT
strategies: Hill Climbing (HC) and Incremental Conductance
(INC).
(a) 2.2.1. Hill Climbing (HC) MPPT
The conventional HC algorithm follows a procedure similar to
the Perturb and Observe (P&O) algorithm. However, in HC, the
perturbation is applied to the duty ratio D, whereas in P&O, it
is applied to the voltage of the PV array. As depicted in Figure
5, the changes in PV power (dPpv) and PV voltage (dVpv)
?. ???????? et al. / JREE: Vol. ?, No. ?, (???????? 2024) 1-??
𝑑𝑖 𝑖 𝑑𝑝
determine the duty cycle perturbation in the HC algorithm 𝑖𝑓 = − = 0 𝑎𝑡 𝑀𝑃𝑃
𝑑𝑣 𝑣 𝑑𝑣
(Mohamed & Abd El Sattar, 2019). 𝑑𝑖 𝑖 𝑑𝑝
𝑖𝑓 < − < 0 𝑟𝑖𝑔ℎ𝑡 𝑜𝑓 𝑀𝑃𝑃 (19)
𝑑𝑣 𝑣 𝑑𝑣
The duty cycle D is determined by the formula provided in 𝑑𝑖 𝑖 𝑑𝑝
Eq. (17), where k is the iteration number and ∆D is the {𝑖𝑓 𝑑𝑣
> −
𝑣 𝑑𝑣
> 0 𝑙𝑒𝑓𝑡 𝑜𝑓 𝑀𝑃𝑃
perturbation step size. As a result of comparing di/dv and i/v, the incremental
𝐷(𝑘) = 𝐷(𝑘 − 1) ± ∆𝐷 (17) conductance, and instantaneous conductance, repectively, the
duty cycle D is also governed by Eq. (17), allowing for the
tracking of the MPP (Jately et al., 2021). Figure 6 shows the
flowchart of the INC approach.
Figure 8. Simulation results at 1000 W/m2 & 25°C & step size : 5e-4.
Table 3. Comparaison of HC and INC MPPT results at 1000 W/m2 & 25°C and different step size.
?. ???????? et al. / JREE: Vol. ?, No. ?, (???????? 2024) 1-??
MPPT HC INC
Tr (ms) 48 76 17 24.73
76.1
∆P (kW) 16 18.8 13.9
At STC 62.4
1000 ∆V (V) 16.1 18.9 13.4
W/m2 &
25 °C 40.2
∆I (A) 8.1 9.5 7
Based on the results obtained, both the Hill Climbing (HC) Additionally, another test was performed at Standard Test
and Incremental Conductance (INC) methods nearly achieve Conditions (STC) by increasing the perturbation step size from
Maximum Power Point (MPP). However, the INC method 5e-4 to 1e-3. The findings presented in Table 3 clearly delineate
achieves MPP more rapidly and accurately compared to HC. the impact of a larger step size on MPPT controllers: the
The boost action exhibits a delay of 50 ms relative to INC, as dynamic response quickens but introduces significant ripple
illustrated in Figures 8b and 8c. Table 3 indicates that INC around the MPP in steady state, thereby reducing MPPT
induces minimal power oscillations, resulting in higher controller efficiency.
efficiency at 99.7% compared to 97.57% with HC.
3.2. Test at variable irradiance and temperature
Figure 9. Simulation results at variable irradiance, temperature, and step size : 5e-4.
MPPT HC INC
Tr (ms) 54.14 20
∆P (kW) 18.8 8.1
700 W/m2 & ∆V (V) 22.8 9.9
25 °C
∆I (A) 11.5 5.1
Efficiency (%) 98.29 99.31
step size 5e-4