# Optimal Placement and Sizing Problem for Power Loss Minimization and Voltage Profile Improvement of Distribution Networks under Seasonal Loads Using Harris Hawks Optimizer

Authors: Habib Ur Rahman Habib; Asad Waqar; Sharoze Sohail; Abdul Khalique Junejo; Mahmoud F. Elmorshedy; Sheheryar Khan; Yun-Su Kim; Moustafa Magdi Ismail
Journal: International Transactions on Electrical Energy Systems
 (2022)

---

## Abstract

Improving efficiency with sustainable radial distribution networks (RDNs) is challenging for larger systems and small grid-connected RDNs. In this paper, the optimal placement of DGs with the Harris hawks optimizer (HHO) under seasonal load demands is proposed to simultaneously reduce total active and reactive power losses and minimize bus voltage drops with the consideration of operational constraints of RDNs. HHO is a newly inspired metaheuristic optimization algorithm primarily based on the Harris hawks’ intelligent behaviors during the chasing of the prey. Furthermore, the authors have investigated four stages of DGs. The first stage involves the optimal allocation of one DG. The second stage includes an investigation with two DGs, the third stage considers three DGs, and the fourth stage investigates the integration of four DGs. The effectiveness of the applied HHO is validated on IEEE 33 and 69 bus RDNs, and results are analyzed by comparing with the standard optimization methods. The Big-O test is also executed for statistical analysis with standard algorithms. The simulation results reveal the better performance of the applied HHO under different circumstances than other algorithms. Furthermore, the total active and reactive power losses and bus voltage drops are improved by adding more DGs into IEEE 33 and 69 bus RDNs.


---

## Body

## 1. Introduction

### 1.1. Motivation

DG penetration in the RDNs has changed its structural behavior from passive to active under power flow in various directions. DGs are small power generation units integrated into RDNs to enhance reliable power delivery, reduce the totalP and Q losses, and improve the bus voltage levels [1]. These DGs can be categorized into conventional power generation sources such as DEs and renewable energy resources (RERs) such as PV and wind power. Nevertheless, future expansion with high integration of RER-based DGs will have both positive and negative outcomes. The adverse effects are power flow reversal, undesired bus voltages level, and power loss [2]. Therefore, these adverse outcomes can be tackled with the help of carefully optimizing the best possible locations and sizing of these DGs.


### 1.2. Literature Review

Different optimization techniques have been proposed in the literature studies to allocate DGs with their optimal sizes [3, 4]. These optimization issues have been tackled with single and multiple objectives. Minimizing the power loss has been taken as the prime objective in single-objective optimization problems. Meanwhile, both single and multi-objective issues have been optimized using metaheuristic optimization methods during the optimal allocation of DGs.During single-objective optimization problems, GA has been used for optimal DG allocation with minimization of the totalP and Q losses [5]. In [6, 7], PSO has been employed to minimize the active power losses with different load types. AI (artificial intelligence)-based optimization techniques have been introduced in [8, 9] for determining the optimal allocation of multiple DG units. The fuzzy-based technique has been used in [10] for the optimal allocation of DGs. Various nature-inspired optimization methods have been employed in the recent literature for the optimal placement of DGs. These methods include BFOA [11], SKHA (stud krill herd algorithm) [12], WOA [13], and CSCA [14].Two schemes have been used in the literature studies to handle multi-objective problems. The first method includes the weighting sum of each objective function. Many researchers have used the same method for optimizing different objective functions such as power losses, VDI, and VSI. The optimization algorithms which have been used in this methodology include GA [15], PSO [15], GA-PSO [15], TLBO [16], QOTLBO [16], SIMBO-Q (swine influenza model-based optimization with quarantine) [17], QOSIMBO-Q [17], and ICA-GA [18]. This first method faces a few challenging problems while selecting the weighting factors. The second method incorporates multi-objective methodologies with trade-off concepts among different objectives using Pareto dominance (PD) criteria. In PD, all the feasible solution sets are categorized into two sets such as dominant and non-dominant. The decision-making expert can select the best possible solution from the non-dominant solution set [4]. Various algorithms have been used based on this second methodology. These algorithms include PAES, NSGA-II, SPEA (strength Pareto evolutionary algorithm), SPEA-II, and MOPSO [20]. MOPSO with fuzzy has been used in [19] to minimize power losses and VDI improvement. MOWOA has been applied in [20] to enhance VSI and reduce VDI and power losses. MOSBA has been employed in [21] for analyzing the DG influence under different load scenarios. TM and MOTA have been implemented for the optimal integration of DGs [22].In the literature studies, researchers have also analyzed DG allocation using IEEE 33 [10] and 69 bus RDNs [23] to minimize power losses and VDI.In [24], the authors have presented ALO (ant-lion optimizer) for optimal sizing and DG allocation with two cases for DGs such as one and two DGs. The authors in [25] applied one DG-based PSAT (power system analysis toolbox) method. In [26], the authors have employed ALO to allocate sites for PV-based one and two DGs. In [27], CSFS (chaotic stochastic fractal search) optimizer has been used to find many DGs with their sizing and site selection under multiple DGs scenarios. GWO (grey wolf optimizer) has been used in [28] for the DG allocation with one or two DGs for enhancing VSI and reducing power losses. In [29], an EMA has been proposed for optimal DG sizes and site selection of one or two DGs. The authors in [1] introduced the optimal DG allocation for accommodating the power flow, VSI, power factor, and lines loss. CDE (chaotic differential evolution) optimizer has been implemented in [30] for the optimal placement of DGs. Hybrid GA-GSA (genetic-gravitational search algorithm) has been employed [31] for DG allocation. In [32], a combined WIPSO-GSA (mixed weight improved particle swarm optimization-gravitational search) has also been used for the installation of two DGs along with capacitor banks.An HSA has been introduced to scale three DGs [33]. The authors in [34] have implemented a BFA (bacterial foraging algorithm) to scale three DGs. In [35], a combined ACO-ABC method has been proposed for the optimal allocation of three DGs. In [36], the BA (bat algorithm) method has been employed to scale PV sizes. The HA (hybrid algorithm) method has been introduced in [7] to handle power loss and VSI for optimal allocation of three DGs. The HA method utilized an analytical variance of PSO. In [37], HGWO (hybrid grey wolf optimizer) has been employed on the Indian network for the optimal allocation of three DGs. In [12], the integration of three DGs in the Portuguese 94-bus grid has been implemented with SKHA. The authors in [29] have introduced an EMA for deploying three DGs in the chosen network. In [38], SPEA-II has been used for the optimal integration of three DGs into the testing network. In [39], WCA (water cycle algorithm) has been suggested for optimal sizing and site selection with three DGs parallel to capacitor banks. An SSA (salp swarm algorithm) has been chosen in [40] for the optimal allocation of three DGs with capacitor banks. The authors in [41] applied an ASFLA (adaptive shuffled frog leaping algorithm) for optimally allocating three DGs. A combined GSA-GAMS (gravitational search algorithm-general algebraic modeling system) method has been implemented in [31] for optimal DG allocation. In [42], the QOCSOS (quasi-oppositional chaotic symbiotic organisms search) algorithm has been applied to search optimal sites of three DGs. The authors in [14] have introduced the CSCA technique for the optimal placement of three DGs. A combined TLCHS optimizer has been validated in [43] for optimal scaling and siting of four DGs. In [44], GA has been implemented for optimal deployment of one, two, and three DG units using IEEE 33-bus RDN. The authors in [45] have proposed an AIS (artificial immune system) for optimal deployment of DG units in IEEE 33 and 15 node RDNs. In [46], the authors have employed EHO [47] for optimal sizing and placement of DGs in IEEE 15, 33, and 69 bus systems. In [48], IHSA has been introduced for optimal sizing and site selection of three DGs in the IEEE 33-bus system.HHO algorithm has been employed in recent literature studies. Optimal DN topology with network reconfiguration is proposed in [49] using the HHO algorithm to minimize power losses and bus voltage deviation. In [50], renewable DGs are integrated into the DNs for power loss reduction with improved bus voltage profiles by considering the seasonal uncertainties of load demands. The authors in [51] investigated the optimal sizing and deployment of WTs in DNs, considering the correlation between load demands and wind power generation.


### 1.3. Main Contributions and Paper Organization

This research study proposes the optimal sizing and allocation of DGs in radial distribution networks (RDNs). HHO [52] is applied based on the hunting method of Harris hawks. The main benefit of HHO is a simple implementation technique with few scenarios of exploration and exploitation. HHO algorithm has been introduced in literature studies for solving different optimization problems, namely, identifying parameters of FC [53] PV cell [54] modules. Nevertheless, HHO is presented in this research study for optimally allocating DGs with single-objective function-based optimization issues. In HHO, the exploration phase represents a vast expansion of a searching space, while the exploitation phase relates to finding the best local solution. That is why the HHO algorithm has enhanced diversity. However, the significant contributions of this research study are summed up as follows:(i)
HHO algorithm is applied for optimally allocating DGs in the radial distribution networks (RDNs) for simultaneously reducing the total active and reactive power losses and bus voltage drops.(ii)
Four seasons, namely, summer, autumn, winter, and spring, integrate one, two, three, and four DGs, missing in the literature studies. This study’s main objectives include reducing total active and reactive power losses and minimizing the bus voltage drop.(iii)
Detailed analysis under multiple deployment scenarios of four types of DGs integration with four seasonal load demands by employing HHO and two IEEE RDNs is also unattended in the past literature studies.(iv)
The performance of the applied HHO algorithm is validated under various operating conditions.(v)
The statistical analysis uses Big-O under different DG integration and weather conditions.The remainder of the paper is organized as follows. In Section2, the formulation of the problem with methodology is presented. It includes the formulation of objectives with system constraints. In Section 3, the mathematical modeling of the HHO optimization algorithm is explained under different searching scenarios. Section 4 presents the results and analysis with a detailed analysis of the IEEE 33 and 69 bus RDNs under various scenarios. The conclusion is drawn in Section 5.


## 1.1. Motivation

DG penetration in the RDNs has changed its structural behavior from passive to active under power flow in various directions. DGs are small power generation units integrated into RDNs to enhance reliable power delivery, reduce the totalP and Q losses, and improve the bus voltage levels [1]. These DGs can be categorized into conventional power generation sources such as DEs and renewable energy resources (RERs) such as PV and wind power. Nevertheless, future expansion with high integration of RER-based DGs will have both positive and negative outcomes. The adverse effects are power flow reversal, undesired bus voltages level, and power loss [2]. Therefore, these adverse outcomes can be tackled with the help of carefully optimizing the best possible locations and sizing of these DGs.


## 1.2. Literature Review

Different optimization techniques have been proposed in the literature studies to allocate DGs with their optimal sizes [3, 4]. These optimization issues have been tackled with single and multiple objectives. Minimizing the power loss has been taken as the prime objective in single-objective optimization problems. Meanwhile, both single and multi-objective issues have been optimized using metaheuristic optimization methods during the optimal allocation of DGs.During single-objective optimization problems, GA has been used for optimal DG allocation with minimization of the totalP and Q losses [5]. In [6, 7], PSO has been employed to minimize the active power losses with different load types. AI (artificial intelligence)-based optimization techniques have been introduced in [8, 9] for determining the optimal allocation of multiple DG units. The fuzzy-based technique has been used in [10] for the optimal allocation of DGs. Various nature-inspired optimization methods have been employed in the recent literature for the optimal placement of DGs. These methods include BFOA [11], SKHA (stud krill herd algorithm) [12], WOA [13], and CSCA [14].Two schemes have been used in the literature studies to handle multi-objective problems. The first method includes the weighting sum of each objective function. Many researchers have used the same method for optimizing different objective functions such as power losses, VDI, and VSI. The optimization algorithms which have been used in this methodology include GA [15], PSO [15], GA-PSO [15], TLBO [16], QOTLBO [16], SIMBO-Q (swine influenza model-based optimization with quarantine) [17], QOSIMBO-Q [17], and ICA-GA [18]. This first method faces a few challenging problems while selecting the weighting factors. The second method incorporates multi-objective methodologies with trade-off concepts among different objectives using Pareto dominance (PD) criteria. In PD, all the feasible solution sets are categorized into two sets such as dominant and non-dominant. The decision-making expert can select the best possible solution from the non-dominant solution set [4]. Various algorithms have been used based on this second methodology. These algorithms include PAES, NSGA-II, SPEA (strength Pareto evolutionary algorithm), SPEA-II, and MOPSO [20]. MOPSO with fuzzy has been used in [19] to minimize power losses and VDI improvement. MOWOA has been applied in [20] to enhance VSI and reduce VDI and power losses. MOSBA has been employed in [21] for analyzing the DG influence under different load scenarios. TM and MOTA have been implemented for the optimal integration of DGs [22].In the literature studies, researchers have also analyzed DG allocation using IEEE 33 [10] and 69 bus RDNs [23] to minimize power losses and VDI.In [24], the authors have presented ALO (ant-lion optimizer) for optimal sizing and DG allocation with two cases for DGs such as one and two DGs. The authors in [25] applied one DG-based PSAT (power system analysis toolbox) method. In [26], the authors have employed ALO to allocate sites for PV-based one and two DGs. In [27], CSFS (chaotic stochastic fractal search) optimizer has been used to find many DGs with their sizing and site selection under multiple DGs scenarios. GWO (grey wolf optimizer) has been used in [28] for the DG allocation with one or two DGs for enhancing VSI and reducing power losses. In [29], an EMA has been proposed for optimal DG sizes and site selection of one or two DGs. The authors in [1] introduced the optimal DG allocation for accommodating the power flow, VSI, power factor, and lines loss. CDE (chaotic differential evolution) optimizer has been implemented in [30] for the optimal placement of DGs. Hybrid GA-GSA (genetic-gravitational search algorithm) has been employed [31] for DG allocation. In [32], a combined WIPSO-GSA (mixed weight improved particle swarm optimization-gravitational search) has also been used for the installation of two DGs along with capacitor banks.An HSA has been introduced to scale three DGs [33]. The authors in [34] have implemented a BFA (bacterial foraging algorithm) to scale three DGs. In [35], a combined ACO-ABC method has been proposed for the optimal allocation of three DGs. In [36], the BA (bat algorithm) method has been employed to scale PV sizes. The HA (hybrid algorithm) method has been introduced in [7] to handle power loss and VSI for optimal allocation of three DGs. The HA method utilized an analytical variance of PSO. In [37], HGWO (hybrid grey wolf optimizer) has been employed on the Indian network for the optimal allocation of three DGs. In [12], the integration of three DGs in the Portuguese 94-bus grid has been implemented with SKHA. The authors in [29] have introduced an EMA for deploying three DGs in the chosen network. In [38], SPEA-II has been used for the optimal integration of three DGs into the testing network. In [39], WCA (water cycle algorithm) has been suggested for optimal sizing and site selection with three DGs parallel to capacitor banks. An SSA (salp swarm algorithm) has been chosen in [40] for the optimal allocation of three DGs with capacitor banks. The authors in [41] applied an ASFLA (adaptive shuffled frog leaping algorithm) for optimally allocating three DGs. A combined GSA-GAMS (gravitational search algorithm-general algebraic modeling system) method has been implemented in [31] for optimal DG allocation. In [42], the QOCSOS (quasi-oppositional chaotic symbiotic organisms search) algorithm has been applied to search optimal sites of three DGs. The authors in [14] have introduced the CSCA technique for the optimal placement of three DGs. A combined TLCHS optimizer has been validated in [43] for optimal scaling and siting of four DGs. In [44], GA has been implemented for optimal deployment of one, two, and three DG units using IEEE 33-bus RDN. The authors in [45] have proposed an AIS (artificial immune system) for optimal deployment of DG units in IEEE 33 and 15 node RDNs. In [46], the authors have employed EHO [47] for optimal sizing and placement of DGs in IEEE 15, 33, and 69 bus systems. In [48], IHSA has been introduced for optimal sizing and site selection of three DGs in the IEEE 33-bus system.HHO algorithm has been employed in recent literature studies. Optimal DN topology with network reconfiguration is proposed in [49] using the HHO algorithm to minimize power losses and bus voltage deviation. In [50], renewable DGs are integrated into the DNs for power loss reduction with improved bus voltage profiles by considering the seasonal uncertainties of load demands. The authors in [51] investigated the optimal sizing and deployment of WTs in DNs, considering the correlation between load demands and wind power generation.


## 1.3. Main Contributions and Paper Organization

This research study proposes the optimal sizing and allocation of DGs in radial distribution networks (RDNs). HHO [52] is applied based on the hunting method of Harris hawks. The main benefit of HHO is a simple implementation technique with few scenarios of exploration and exploitation. HHO algorithm has been introduced in literature studies for solving different optimization problems, namely, identifying parameters of FC [53] PV cell [54] modules. Nevertheless, HHO is presented in this research study for optimally allocating DGs with single-objective function-based optimization issues. In HHO, the exploration phase represents a vast expansion of a searching space, while the exploitation phase relates to finding the best local solution. That is why the HHO algorithm has enhanced diversity. However, the significant contributions of this research study are summed up as follows:(i)
HHO algorithm is applied for optimally allocating DGs in the radial distribution networks (RDNs) for simultaneously reducing the total active and reactive power losses and bus voltage drops.(ii)
Four seasons, namely, summer, autumn, winter, and spring, integrate one, two, three, and four DGs, missing in the literature studies. This study’s main objectives include reducing total active and reactive power losses and minimizing the bus voltage drop.(iii)
Detailed analysis under multiple deployment scenarios of four types of DGs integration with four seasonal load demands by employing HHO and two IEEE RDNs is also unattended in the past literature studies.(iv)
The performance of the applied HHO algorithm is validated under various operating conditions.(v)
The statistical analysis uses Big-O under different DG integration and weather conditions.The remainder of the paper is organized as follows. In Section2, the formulation of the problem with methodology is presented. It includes the formulation of objectives with system constraints. In Section 3, the mathematical modeling of the HHO optimization algorithm is explained under different searching scenarios. Section 4 presents the results and analysis with a detailed analysis of the IEEE 33 and 69 bus RDNs under various scenarios. The conclusion is drawn in Section 5.


## 2. Problem Formulation

In this section, the objectives for optimal DG placement are included. In this research, the analysis is performed for the rural town (Shah Allah Ditta) in Islamabad, Pakistan (33.7209642 °N 72.9143201 °E). This rural town is 700 years old and was recognized as an essential route from Kabul (Afghanistan) to the Gandharan city of Taxila (Hindustan) by Alexander the Great and Sher Shah Suri. Other emperors, including Mughal rulers, frequently moved through this route during their travel from Afghanistan to Hindustan.


### 2.1. Objective Functions

The main objective of DG allocation in the distributed power system is to minimize total active and reactive power losses and bus voltage drop. The mathematical modeling of all these objectives is discussed in this section as follows.


#### 2.1.1. Minimizing Active Power Loss

Due to the radial distribution network, active power losses are more prominent. Therefore, these power losses (Ploss) should be reduced as follows:(1)fObj1=minPloss.The totalPloss can be calculated using the branch current losses relationship as follows [55]:(2)Ploss=∑b=1MbIb2Rb,where b denotes the branch number, Mb represents the total branches, Ib shows the current flowing, and Rb means branch resistance.


#### 2.1.2. Minimizing Total Reactive Power Loss

Due to the radial distribution network, reactive power losses are more prominent. Therefore, these power losses (Qloss) should be reduced as follows:(3)fObj2=minQloss.The totalP loss can be calculated using the branch current loss relationship as follows [55]:(4)Qloss=∑b=1MbIb2Xb,where Xb means branch reactance.


#### 2.1.3. Minimizing Voltage Deviation Index

The total deviated voltageVD indicated the voltage level of the RDNs and how far VD is from the targeted value Vt. Therefore, VD of the RDN can be found by taking voltage magnitude vj at bus j based on the targeted voltage as follows [56]:(5)VD=∑j=1NjVt−Vj,where Vt is considered as 1.0 per unit.


#### 2.1.4. Maximizing Voltage Stability Index

The voltage stability index (VSI) is the ability of the RDN to maintain the voltage value within the specified range. The main objective is to maximize VSI as follows [57]:(6)VSIR=VS4−4PRRSR+QRXSRVS2−4PRXSR−QRRSR,where SR denotes the sending and receiving end; PR, QR represent real and imaginary power at the receiving end; and RSR, XSR show the resistance and reactance between sending and receiving ends.


### 2.2. System Constraints

The main constraints of DG allocation in the distributed power system are defined as follows.


#### 2.2.1. Equality Constraints

The generation-demand balance with the consideration of total active and reactive power losses can be expressed as follows [58]:(7)∑j=1NdgPgj=Ploss+Pload,(8)∑j=1NdgQgj=Qloss+Qload,where Ndg indicates the total number of installed DGs, Pgi is the generated power of the jth DG, and Pload denotes the total demanded power of the loads.


#### 2.2.2. Inequality Constraints

The operating limits of the distributed power system (such as active/reactive power and voltage limits) can be considered as follows [59]:(9)Pgjmin≤Pgj≤Pgjmax,(10)Qgjmin≤Qgj≤Qgjmax,(11)0.95≤Vj≤1.05.


## 2.1. Objective Functions

The main objective of DG allocation in the distributed power system is to minimize total active and reactive power losses and bus voltage drop. The mathematical modeling of all these objectives is discussed in this section as follows.


### 2.1.1. Minimizing Active Power Loss

Due to the radial distribution network, active power losses are more prominent. Therefore, these power losses (Ploss) should be reduced as follows:(1)fObj1=minPloss.The totalPloss can be calculated using the branch current losses relationship as follows [55]:(2)Ploss=∑b=1MbIb2Rb,where b denotes the branch number, Mb represents the total branches, Ib shows the current flowing, and Rb means branch resistance.


### 2.1.2. Minimizing Total Reactive Power Loss

Due to the radial distribution network, reactive power losses are more prominent. Therefore, these power losses (Qloss) should be reduced as follows:(3)fObj2=minQloss.The totalP loss can be calculated using the branch current loss relationship as follows [55]:(4)Qloss=∑b=1MbIb2Xb,where Xb means branch reactance.


### 2.1.3. Minimizing Voltage Deviation Index

The total deviated voltageVD indicated the voltage level of the RDNs and how far VD is from the targeted value Vt. Therefore, VD of the RDN can be found by taking voltage magnitude vj at bus j based on the targeted voltage as follows [56]:(5)VD=∑j=1NjVt−Vj,where Vt is considered as 1.0 per unit.


### 2.1.4. Maximizing Voltage Stability Index

The voltage stability index (VSI) is the ability of the RDN to maintain the voltage value within the specified range. The main objective is to maximize VSI as follows [57]:(6)VSIR=VS4−4PRRSR+QRXSRVS2−4PRXSR−QRRSR,where SR denotes the sending and receiving end; PR, QR represent real and imaginary power at the receiving end; and RSR, XSR show the resistance and reactance between sending and receiving ends.


## 2.1.1. Minimizing Active Power Loss

Due to the radial distribution network, active power losses are more prominent. Therefore, these power losses (Ploss) should be reduced as follows:(1)fObj1=minPloss.The totalPloss can be calculated using the branch current losses relationship as follows [55]:(2)Ploss=∑b=1MbIb2Rb,where b denotes the branch number, Mb represents the total branches, Ib shows the current flowing, and Rb means branch resistance.


## 2.1.2. Minimizing Total Reactive Power Loss

Due to the radial distribution network, reactive power losses are more prominent. Therefore, these power losses (Qloss) should be reduced as follows:(3)fObj2=minQloss.The totalP loss can be calculated using the branch current loss relationship as follows [55]:(4)Qloss=∑b=1MbIb2Xb,where Xb means branch reactance.


## 2.1.3. Minimizing Voltage Deviation Index

The total deviated voltageVD indicated the voltage level of the RDNs and how far VD is from the targeted value Vt. Therefore, VD of the RDN can be found by taking voltage magnitude vj at bus j based on the targeted voltage as follows [56]:(5)VD=∑j=1NjVt−Vj,where Vt is considered as 1.0 per unit.


## 2.1.4. Maximizing Voltage Stability Index

The voltage stability index (VSI) is the ability of the RDN to maintain the voltage value within the specified range. The main objective is to maximize VSI as follows [57]:(6)VSIR=VS4−4PRRSR+QRXSRVS2−4PRXSR−QRRSR,where SR denotes the sending and receiving end; PR, QR represent real and imaginary power at the receiving end; and RSR, XSR show the resistance and reactance between sending and receiving ends.


## 2.2. System Constraints

The main constraints of DG allocation in the distributed power system are defined as follows.


### 2.2.1. Equality Constraints

The generation-demand balance with the consideration of total active and reactive power losses can be expressed as follows [58]:(7)∑j=1NdgPgj=Ploss+Pload,(8)∑j=1NdgQgj=Qloss+Qload,where Ndg indicates the total number of installed DGs, Pgi is the generated power of the jth DG, and Pload denotes the total demanded power of the loads.


### 2.2.2. Inequality Constraints

The operating limits of the distributed power system (such as active/reactive power and voltage limits) can be considered as follows [59]:(9)Pgjmin≤Pgj≤Pgjmax,(10)Qgjmin≤Qgj≤Qgjmax,(11)0.95≤Vj≤1.05.


## 2.2.1. Equality Constraints

The generation-demand balance with the consideration of total active and reactive power losses can be expressed as follows [58]:(7)∑j=1NdgPgj=Ploss+Pload,(8)∑j=1NdgQgj=Qloss+Qload,where Ndg indicates the total number of installed DGs, Pgi is the generated power of the jth DG, and Pload denotes the total demanded power of the loads.


## 2.2.2. Inequality Constraints

The operating limits of the distributed power system (such as active/reactive power and voltage limits) can be considered as follows [59]:(9)Pgjmin≤Pgj≤Pgjmax,(10)Qgjmin≤Qgj≤Qgjmax,(11)0.95≤Vj≤1.05.


## 3. Harris Hawks Optimizer (HHO)

HHO is a population-based method with two phases for implementation: exploration and exploitation. The mathematical derivation of both phases is explained as follows.


### 3.1. Exploration Stage

The prime purpose of the Harris hawks is hunting the prey, which is primarily a rabbit. Therefore, the hawks search for the rabbit. This exploring process can be categorized into two scenarios. The first scenario supposes that the hawks’ locations are near the family members and the prey. At the same time, the following scenario considers the hawks’ locations at random trees. The mathematical formulation of these two scenarios can be written as follows [60]:(12)Xi+1=Xrandi−rand1Xrandi−2rand2Xi,p≥0.5,Xpreyi−Xmi−rand3SL+rand4SU−SL,p≥0.5,where i indicates the current iteration, X(i + 1) is the hawks’ position at iteration i + 1, X(i) represents the hawks’ position at the current iteration i, Xprey(i) represents the prey’s position at the current iteration i, Xrand(i) represents the random position selection at the current iteration i, Rand1, Rand2, Rand3, and Rand4 represent the random position selection within [0, 1], and SL and SU denote the searching space with the lower and upper limits, respectively. Both exploration scenarios can be activated using a random variable p between 0 and 1.Xm(i) represents the mean hawks’ position at the current iteration i and can be expressed as follows [61]:(13)Xmi=1h∑y=1hXyi,where Xy(i) denotes the hawk’s position y and h represents the total number of hawks.


### 3.2. Transformation from Exploration to Exploitation

The running away of the energy of the rabbitER during the chasing process is considered for the transformation from exploration to exploitation in the HHO algorithm. It is defined as follows [61]:(14)ER=2Eo1−iI,where I denotes the maximum number of iterations and Eo denotes the randomly generated initial energy of the rabbit within [−1, 1]. If ER ≥ 1, it means the hawks still carry on prey’s tendency to escape and explore. If ER ≤ 1, the hawks will begin the process of exploitation close to the prey position, and the target will shift to avoid running.


### 3.3. Exploitation Stage

The exploitation stage is subjected to the chances of the prey escaping (e) and running away with the energy ER. The target successfully runs away when e < 0.5 and is unsuccessful during e > 0.5. The hawks can have two options based on running away energy: soft surround during |ER| ≥ 0.5 and hard surround during |ER| < 0.5.As a result, the exploitation method is categorized in four steps as follows.


#### 3.3.1. Easy surround

The easy surround indicates the prey’s efforts to escape with the help of haphazard jumps; nevertheless, the hawks quickly surround it. The mathematical derivation of this easy surround is defined as follows [52]:(15)Xi+1=ΔXi−EKXpreyi−Xi,(16)K=21−rand5,(17)ΔXi=Xpreyi−Xi,where ΔX(i) represents the distance of prey from hawks’ location, K represents the randomly taken jumps of the target, and rand5 denotes the random number within [0, 1].


#### 3.3.2. Difficult Surround

The difficult surround occurs whene ≥ 0.5 and ER < 0.5. In this step, the prey is very tired, and it is difficult for the hawks to surround the rabbit (target). The mathematical formulation is written as follows [62]:(18)Xi+1=Xpreyi−EΔXi.


#### 3.3.3. Easy Surround with Progressive Speedy Dives

This surround is assumed to be an intelligent scheme that differentiates the HHO algorithm from the other swarm optimization techniques. Whene < 0.5 and |ER| ≥ 0.5, the prey has the energy to run away, and the hawks can quickly surround the prey. A Lévy flight (LF) idea is used for the formulation of this surround step as follows [52]:(19)W=Xpreyi−EKXpreyi−Xi,where W represents the easy surround location. Based on LF, the hawks’ speedy dives are formulated as follows [62]:(20)V=W+SR∗LFDP,where DP represents the problem dimension and SR denotes the vector with randomly generated values with the matrix size of 1 ∗ DP. The LF is defined as follows [62]:(21)LFDP=0.01×α×σβ1/c,(22)σ=Γ1+c×sinπc/2Γ1+c/2×c×2c−1/21/c,where c represents a constant value taken as 1.5 and α and β denote randomly generated values [0, 1]. Therefore, the hawks’ location at the next iteration is defined as follows [62]:(23)Xi+1=W,FW<FXi,V,FZ<FXi.


#### 3.3.4. Difficult surround with Progressive Speedy Dives

The surround is considered whene < 0.5 and |ER| < 0.5, the prey is very tired, and the hawks have difficulty surrounding the prey. A similar Lévy flight (LF) idea is used for the formulation in equations (18) to (21), while the estimation of W is as follows [62]:(24)W=Xpreyi−EKXpreyi−Xmi.The implementation steps of the HHO algorithm for DGs allocation are defined as follows:Step 1: collecting the system data for the input of the HHO algorithm, such as line/load data. The objective functions are also defined along with the constraints.Step 2: initializing the set of hawks’ searches with randomly generated values for upper and lower limits of DGs sizing with the location. It also includes the initialization of HHO parameters with maximum iterationsimax.Step 3: simulating the model with power flow analysis to calculate the objective functions such asPloss for every searching hawk.Step 4: saving the best value of solutionXprey(i).Step 5: updating the HHO parameters such asER, Eo, and K.Step 6: updating the sizing and location of the best possible solution sets based on the two stages such as exploration and exploitation.Step 7: updating the limitations on sizing and location of DGs while updating positionX(i + 1).Step 8: observing ifi < imax. Go to step 3 again.Step 9: saving the finalized best solution, i.e., DGs location and sizing.Step 10: simulating the power analysis for obtaining the bus voltage profiles.


## 3.1. Exploration Stage

The prime purpose of the Harris hawks is hunting the prey, which is primarily a rabbit. Therefore, the hawks search for the rabbit. This exploring process can be categorized into two scenarios. The first scenario supposes that the hawks’ locations are near the family members and the prey. At the same time, the following scenario considers the hawks’ locations at random trees. The mathematical formulation of these two scenarios can be written as follows [60]:(12)Xi+1=Xrandi−rand1Xrandi−2rand2Xi,p≥0.5,Xpreyi−Xmi−rand3SL+rand4SU−SL,p≥0.5,where i indicates the current iteration, X(i + 1) is the hawks’ position at iteration i + 1, X(i) represents the hawks’ position at the current iteration i, Xprey(i) represents the prey’s position at the current iteration i, Xrand(i) represents the random position selection at the current iteration i, Rand1, Rand2, Rand3, and Rand4 represent the random position selection within [0, 1], and SL and SU denote the searching space with the lower and upper limits, respectively. Both exploration scenarios can be activated using a random variable p between 0 and 1.Xm(i) represents the mean hawks’ position at the current iteration i and can be expressed as follows [61]:(13)Xmi=1h∑y=1hXyi,where Xy(i) denotes the hawk’s position y and h represents the total number of hawks.


## 3.2. Transformation from Exploration to Exploitation

The running away of the energy of the rabbitER during the chasing process is considered for the transformation from exploration to exploitation in the HHO algorithm. It is defined as follows [61]:(14)ER=2Eo1−iI,where I denotes the maximum number of iterations and Eo denotes the randomly generated initial energy of the rabbit within [−1, 1]. If ER ≥ 1, it means the hawks still carry on prey’s tendency to escape and explore. If ER ≤ 1, the hawks will begin the process of exploitation close to the prey position, and the target will shift to avoid running.


## 3.3. Exploitation Stage

The exploitation stage is subjected to the chances of the prey escaping (e) and running away with the energy ER. The target successfully runs away when e < 0.5 and is unsuccessful during e > 0.5. The hawks can have two options based on running away energy: soft surround during |ER| ≥ 0.5 and hard surround during |ER| < 0.5.As a result, the exploitation method is categorized in four steps as follows.


### 3.3.1. Easy surround

The easy surround indicates the prey’s efforts to escape with the help of haphazard jumps; nevertheless, the hawks quickly surround it. The mathematical derivation of this easy surround is defined as follows [52]:(15)Xi+1=ΔXi−EKXpreyi−Xi,(16)K=21−rand5,(17)ΔXi=Xpreyi−Xi,where ΔX(i) represents the distance of prey from hawks’ location, K represents the randomly taken jumps of the target, and rand5 denotes the random number within [0, 1].


### 3.3.2. Difficult Surround

The difficult surround occurs whene ≥ 0.5 and ER < 0.5. In this step, the prey is very tired, and it is difficult for the hawks to surround the rabbit (target). The mathematical formulation is written as follows [62]:(18)Xi+1=Xpreyi−EΔXi.


### 3.3.3. Easy Surround with Progressive Speedy Dives

This surround is assumed to be an intelligent scheme that differentiates the HHO algorithm from the other swarm optimization techniques. Whene < 0.5 and |ER| ≥ 0.5, the prey has the energy to run away, and the hawks can quickly surround the prey. A Lévy flight (LF) idea is used for the formulation of this surround step as follows [52]:(19)W=Xpreyi−EKXpreyi−Xi,where W represents the easy surround location. Based on LF, the hawks’ speedy dives are formulated as follows [62]:(20)V=W+SR∗LFDP,where DP represents the problem dimension and SR denotes the vector with randomly generated values with the matrix size of 1 ∗ DP. The LF is defined as follows [62]:(21)LFDP=0.01×α×σβ1/c,(22)σ=Γ1+c×sinπc/2Γ1+c/2×c×2c−1/21/c,where c represents a constant value taken as 1.5 and α and β denote randomly generated values [0, 1]. Therefore, the hawks’ location at the next iteration is defined as follows [62]:(23)Xi+1=W,FW<FXi,V,FZ<FXi.


### 3.3.4. Difficult surround with Progressive Speedy Dives

The surround is considered whene < 0.5 and |ER| < 0.5, the prey is very tired, and the hawks have difficulty surrounding the prey. A similar Lévy flight (LF) idea is used for the formulation in equations (18) to (21), while the estimation of W is as follows [62]:(24)W=Xpreyi−EKXpreyi−Xmi.The implementation steps of the HHO algorithm for DGs allocation are defined as follows:Step 1: collecting the system data for the input of the HHO algorithm, such as line/load data. The objective functions are also defined along with the constraints.Step 2: initializing the set of hawks’ searches with randomly generated values for upper and lower limits of DGs sizing with the location. It also includes the initialization of HHO parameters with maximum iterationsimax.Step 3: simulating the model with power flow analysis to calculate the objective functions such asPloss for every searching hawk.Step 4: saving the best value of solutionXprey(i).Step 5: updating the HHO parameters such asER, Eo, and K.Step 6: updating the sizing and location of the best possible solution sets based on the two stages such as exploration and exploitation.Step 7: updating the limitations on sizing and location of DGs while updating positionX(i + 1).Step 8: observing ifi < imax. Go to step 3 again.Step 9: saving the finalized best solution, i.e., DGs location and sizing.Step 10: simulating the power analysis for obtaining the bus voltage profiles.


## 3.3.1. Easy surround

The easy surround indicates the prey’s efforts to escape with the help of haphazard jumps; nevertheless, the hawks quickly surround it. The mathematical derivation of this easy surround is defined as follows [52]:(15)Xi+1=ΔXi−EKXpreyi−Xi,(16)K=21−rand5,(17)ΔXi=Xpreyi−Xi,where ΔX(i) represents the distance of prey from hawks’ location, K represents the randomly taken jumps of the target, and rand5 denotes the random number within [0, 1].


## 3.3.2. Difficult Surround

The difficult surround occurs whene ≥ 0.5 and ER < 0.5. In this step, the prey is very tired, and it is difficult for the hawks to surround the rabbit (target). The mathematical formulation is written as follows [62]:(18)Xi+1=Xpreyi−EΔXi.


## 3.3.3. Easy Surround with Progressive Speedy Dives

This surround is assumed to be an intelligent scheme that differentiates the HHO algorithm from the other swarm optimization techniques. Whene < 0.5 and |ER| ≥ 0.5, the prey has the energy to run away, and the hawks can quickly surround the prey. A Lévy flight (LF) idea is used for the formulation of this surround step as follows [52]:(19)W=Xpreyi−EKXpreyi−Xi,where W represents the easy surround location. Based on LF, the hawks’ speedy dives are formulated as follows [62]:(20)V=W+SR∗LFDP,where DP represents the problem dimension and SR denotes the vector with randomly generated values with the matrix size of 1 ∗ DP. The LF is defined as follows [62]:(21)LFDP=0.01×α×σβ1/c,(22)σ=Γ1+c×sinπc/2Γ1+c/2×c×2c−1/21/c,where c represents a constant value taken as 1.5 and α and β denote randomly generated values [0, 1]. Therefore, the hawks’ location at the next iteration is defined as follows [62]:(23)Xi+1=W,FW<FXi,V,FZ<FXi.


## 3.3.4. Difficult surround with Progressive Speedy Dives

The surround is considered whene < 0.5 and |ER| < 0.5, the prey is very tired, and the hawks have difficulty surrounding the prey. A similar Lévy flight (LF) idea is used for the formulation in equations (18) to (21), while the estimation of W is as follows [62]:(24)W=Xpreyi−EKXpreyi−Xmi.The implementation steps of the HHO algorithm for DGs allocation are defined as follows:Step 1: collecting the system data for the input of the HHO algorithm, such as line/load data. The objective functions are also defined along with the constraints.Step 2: initializing the set of hawks’ searches with randomly generated values for upper and lower limits of DGs sizing with the location. It also includes the initialization of HHO parameters with maximum iterationsimax.Step 3: simulating the model with power flow analysis to calculate the objective functions such asPloss for every searching hawk.Step 4: saving the best value of solutionXprey(i).Step 5: updating the HHO parameters such asER, Eo, and K.Step 6: updating the sizing and location of the best possible solution sets based on the two stages such as exploration and exploitation.Step 7: updating the limitations on sizing and location of DGs while updating positionX(i + 1).Step 8: observing ifi < imax. Go to step 3 again.Step 9: saving the finalized best solution, i.e., DGs location and sizing.Step 10: simulating the power analysis for obtaining the bus voltage profiles.


## 4. Results and Discussion

This section analyzes the simulation results of different scenarios using two RDNs, namely, IEEE 33 [63] and 69 bus RDNs [64]. The optimal sizing and placement of DGs are simulated using the HHO algorithm. The main objectives are to minimize the total active and reactive power losses and voltage drop and maximize VSI. The MATLAB simulation environment is used, and the simulations are run fifteen times for different scenarios to ensure result accuracy and algorithm robustness. Four strategies are analyzed for each RDN under four seasonal uncertainties of load demand. Case 1 includes one DG, case 2 considers deployment of two DGs, case 3 has three DGs, and case 4 incorporates four DGs. The result comparison with feasible solutions of four cases is analyzed with the standard optimization techniques under four seasonal uncertainties of load demands. Table 1 shows the simulation parameters. Table 2 shows exponent values for different load types.Table 1
Simulation parameters.
Parameters33-bus RDNIEEE 69-bus RDNUnitsPopulation size200200—Maximum iterations100100—Size limit of DGs10∼600010∼6000kWKPR [38]350350$/kWGT [65]0.190.19$/kWhIR [38]12.512.5%IF [38]99%No. of Iterations5050—Table 2
Different values of exponent values and load categories.
LoadαβConstant [66]00Residential [67]Spring0.7202.960Summer0.9204.040Autumn1.0404.190Winter1.3004.380Commercial [68]Spring1.2503.500Summer0.9903.950Autumn1.5003.150Winter1.5103.400


### 4.1. IEEE 33-Bus RDN

Figure1 [63] shows IEEE 33-bus RDN [63], which comprises thirty-seven lines connected with 33 buses. The total active load demand is 3720 kW, while the total reactive load demand is 2300 kVAR. The total P loss is 208.4592 kW, and the bus voltage magnitude is 0.929 (per unit) [69].Figure 1
IEEE 33-bus RDN.Table3 compares the minimum and average bus voltage values of literature with the results of proposed HHO for 33-bus RDN for exponent values of α = 0.720, β = 2.960. Table 4 compares the minimum and average bus voltage values of literature with the results of proposed HHO for 33-bus RDN for exponent values of α = 0.920, β = 4.040. Table 5 compares the minimum and average bus voltage values of literature with the results of proposed HHO for 33-bus RDN for exponent values of α = 1.040, β = 4.190. Table 6 compares the minimum and average bus voltage values of literature with the results of proposed HHO for 33-bus RDN for exponent values of α = 1.300, β = 4.380.Table 3
Simulation results of system I (α = 0.720, β = 2.960).
AlgorithmVave (per unit)vmin (per unit)FWA [70]0.96940.9459MBFOA [71]0.96980.9430ITS [72]0.96770.9397SLR [73]0.96750.9397PSO [68]0.95570.9221MPSO [68]0.96750.9432Proposed HHO0.97860.9391Table 4
Simulation results of system I (α = 0.920, β = 4.040).
AlgorithmVave (per unit)vmin (per unit)FWA [70]0.97000.9474MBFOA [71]0.97050.9446ITS [72]0.96840.9416SLR [73]0.96820.9416PSO [68]0.95520.9246MPSO [68]0.96820.9449Proposed HHO0.97580.9472Table 5
Simulation results of system I (α = 1.040, β = 4.190).
AlgorithmVave (per unit)vmin (per unit)FWA [70]0.97020.9477MBFOA [71]0.97060.9450ITS [72]0.96860.9420SLR [73]0.96840.9420PSO [68]0.95560.9254MPSO [68]0.96840.9453Proposed HHO0.98430.9417Table 6
Simulation results of system I (α = 1.300, β = 4.380).
AlgorithmVave (per unit)vmin (per unit)FWA [70]0.97050.9484MBFOA [71]0.97090.9457ITS [72]0.96890.9428SLR [73]0.96880.9428PSO [68]0.95640.9268MPSO [68]0.96880.9460Proposed HHO0.98710.9526Table7 compares the literature simulation results with the proposed HHO based on four DGs (IEEE 33-bus RDN). It is observed that the maximum active power loss reduction of 66.8% is obtained compared to the other literature results except for one case with the ABC algorithm. The total Q loss reduction of 56.4% is obtained, higher than different algorithm results except for two cases with BFOA and ABC methods. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total active and reactive power loss reduction is 52.98 kW and 42.23 kVAR. At the same time, the bus voltage magnitude is upgraded to 0.9471 (per unit). The objective function for this case is 0.36579.Table 7
Comparison of literature results with four DGs (33-bus RDN).
AlgorithmRef.Yearvmin(per unit)P loss(kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DGsizes (kW)DG’s location(bus no.)Proposed HHO20210.9471452.984342.22866.804756.4191318.67,42.94,997.65,322.883, 6, 8, 4PPA [69]20200.9790058.36—72—1058.0,2201.0,171.70,74.8003,6,4,5GA [74]20200.981064.19019.933047.6050.701500.0,422.80,1071.4,5781.011,29,30,8PSO [74]20200.980070.32016.542048.2050.401176.0,981.60,829.70,959.008,13,32,6GA/PSO [74]20200.967069.21021.320049.2043.70925.00,863.00,1200.0,1400.011,16,32,5HSA [65]20200.967077.140—52.30—572.40,107.00,1046.2,1034.017,18,33,15IWD [14]20200.969668.720—57.70—600.30,300.00,1011.2,1350.09,16,30,13BFOA [75]20200.9640——51.5057.30633.00,90.000,947.00,1000.017,18,33,11IWO [14]20200.971668.710—57.7054.30624.00,104.90,1056.0,879.4014,18,32,22PSO & analytical [76]2020—68.43016.540064.1052.30790.00,1070.0,1010.0,950.0013,24,30,20SKHA [77]20200.968773.450—64.7052.10801.80,1091.0,1053.6,950.0013,24,30,43ABC [65]20200.977079.50044.854067.3058.40514.00,948.00,635.16,479.5014,28,23,3SSA [40]20190.968667.430—64.80—753.60,1100.4,1070.0,570.0013,23,29,5CDE [30]20190.9702067.12—67.8—926.69, 646.70, 967.30,679.306,14,24,31WCA [39]20180.973074.530—64.00—854.6,1101.7,1180.0,750.0014,24,29,11BSOA [75]20180.955465.350—56.1053.4032.000,487.00,550.00,1000.013,28,31,7BA [78]20180.980066.430—63.0051.00816.00,952.50,952.35,598.0015,25,30,2TLCHS [43]20180.9770067.28—67.72—941.20,684.70,966.40,710.706,14,24,31ACO-ABC [68]20170.973571.670—62.8055.40754.00,1099.9,1071.0,589.4014,24,30,9HGWO [79]2013—77.430—64.4053.30802.00,1090.0,1054.0,800.0013,24,30,21Table8 compares the literature simulation results with the proposed HHO based on three DGs (IEEE 33-bus RDN). It is observed that the maximum active power loss reduction of 69.65% is obtained as compared to the other literature results. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total P loss reduction is 107.93 kW. At the same time, the bus voltage magnitude is upgraded to 0.9706 (per unit). The objective function for this case is 0.39752.Table 8
Comparison of literature results with three DGs (33-bus RDN).
Algorithm Ref.Yearvmin (per unit)FOP loss (kW)PLOSS (%)DG size (kW)DG’s locationACO-ABC [35]20150.97350.8007577.54062.800754.700,1099.90,1071.4014,24,30BA [36]20160.98000.8050077.12063.000816.300,952.350,952.35015,25,30HA [7]2016—0.3205074.83064.100790.000,1070.00,1010.0013,24,30HGWO [37]2017—0.3220074.20064.400802.000,1090.00,1054.0013,24,30EMA [29]20180.96840.8058074.37064.320976.600, 1169.09, 943.54030,24,12SPEA2 [38]20180.96160.8363060.24071.100691.000,733.400,742.90018,29,8SSA [40]20190.96860.8083073.37064.800753.600, 1100.40, 1070.6013,23,29GSA-GAMS [31]20190.96860.8125071.62065.640801.220, 1091.30, 1053.5913,24,30GA [44]20190.98600.8280068.78067.000761.000,1170.00,1082.0014,24,30QOCSOS [42]2020—0.3275071.91065.500801.700, 1091.30, 1053.6013,24,30CSCA [14]20200.96900.8070073.99064.500871.000, 1091.47, 954.08013,24,30IHSA [48]2020—0.3260072.54065.200800.800, 1087.60, 1050.7013,24,30PPA [69]20200.97000.8445058.41071.9001141.80, 161.71, 2214.3013,23,28Proposed HHO20210.97060.39752107.92569.6471304.8,751.02,1735.53,6,8


#### 4.1.1. Case 1 (33-Bus RDN) for Winter

Figure2(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.90998 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.39695. Figure 2(b) shows the active power loss, reduced from 159.61 kW to 121.21 kW, while the decrease of total P loss is 24.06%. Figure 2(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 96.89 kVAR to 76.82 kVAR, while the decrease of total Q loss is 20.72%. The optimal DG size is found as 2578.5 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 2
Case 1 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 9
Results of case 1 (33-bus RDN with HHO).
Seasonvmin (per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.909980.39695121.210776.81724.059920.72212578.53Without DGs0.89587—159.613696.8958————Spring0.840890.52997368.0215234.536425.592522.40614888.13Without DGs0.81372—494.6027302.2616————Summer0.814840.59875495.3813316.116125.942922.90145805.53Without DGs0.7819—668.9179410.0154————Autumn0.88820.43087184.8168117.382324.692121.36223291.63Without DGs0.87025—245.4149149.2695————


#### 4.1.2. Case 1 (33-Bus RDN) for Spring

Figure3(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.84089 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.52997. Figure 3(b) shows the total P loss, which is reduced from 494.60 kW to 368.02 kW, while the decrease of total P loss is 25.59%. Figure 3(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 302.26 kVAR to 234.53 kVAR, while the decrease of total Q loss is 22.40%. The optimal DG size is found as 4888.1 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 3
Case 1 (33-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.3. Case 1 (33-Bus RDN) for Summer

Figure4(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.81484 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.59875. Figure 4(b) shows the total P loss reduced from 668.91 kW to 495.38 kW, while the decrease of total P loss is 25.94%. Figure 4(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 410.01 kVAR to 316.11 kVAR, while the decrease of total Q loss is 22.90%. The optimal DG size is found as 5805.5 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 4
Case 1 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.4. Case 1 (33-Bus RDN) for Autumn

Figure5(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.8882 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.43087. Figure 5(b) shows the total P loss, which is reduced from 245.41 kW to 184.81 kW, while the decrease of total P loss is 24.69%. Figure 5(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 149.26 kVAR to 117.38 kVAR, while the decrease of total Q loss is 21.36%. The optimal DG size is found as 3291.6 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 5
Case 1 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.5. Case 2 (33-Bus RDN) for Winter

Figure6(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.9131 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.4016. Figure 6(b) shows the total P loss, which increased from 159.61 kW to 245.46 kW, while the percentage increase of total P loss is 53.78%. Figure 6(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 96.89 kVAR to 142.43 kVAR, while the percentage increase of total Q loss is 46.9967%. Optimal DG sizes are found as 347.21 kW and 989.34 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 6
Case 2 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 10
Results of case 2 (33-bus RDN with HHO).
Seasonvmin(per unit)FOP loss (kW)Q loss (kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’ssize (kW)DG’s location(bus no.)Winter0.91310.4016245.4671142.4336−53.7883−46.9967347.21, 989.343,6Without DGs0.89587—159.613696.8958————Spring0.847380.54649736.4943428.922−48.9062−41.9042715.63, 1844.43,6Without DGs0.81372—494.6027302.2616————Summer0.82290.6219980.4337571.8549−46.5701−39.4716873.16, 2181.43,6Without DGs0.7819—668.9179410.0154————Autumn0.892290.43844374.0595217.289−52.4192−45.5682457.26, 1254.13,6Without DGs0.87025—245.4149149.2695————


#### 4.1.6. Case 2 (33-Bus RDN) for Spring

Figure7(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.84738 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.54649. Figure 7(b) shows the total P loss, which is increased from 494.60 kW to 736.49 kW, while the percentage increase of total P loss is 48.90%. Figure 7(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 302.26 kVAR to 428.9 kVAR, while the percentage increase of total Q loss is 41.90%. Optimal DG sizes are found as 715.63 and 1844.4 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 7
Case 2 (33-bus RDN) for spring. (a) Bus voltage profile. (b) Total P loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.7. Case 2 (33-Bus RDN) for Summer

Figure8(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.8229 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.6219. Figure 8(b) shows the total P loss, which increased from 668.91 kW to 980.43 kW, while the percentage increase of total P loss is 46.57%. Figure 8(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 410.01 kVAR to 571.85 kVAR, while the percentage increase of total Q loss is 39.47%. Optimal DG sizes are found as 873.16 and 2181.4 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 8
Case 2 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.8. Case 2 (33-Bus RDN) for Autumn

Figure9(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.89229 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.43844. Figure 9(b) shows the total P loss, which increased from 245.41 kW to 374.05 kW, while the percentage increase of total P loss is 52.41%. Figure 9(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 149.26 kVAR to 217.28 kVAR, while the percentage increase of total Q loss is 45.56%. Optimal DG sizes are found as 457.26 and 1254.1 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 9
Case 2 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.9. Case 3 (33-Bus RDN) for Winter

Figure10(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.95257 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.36559. Figure 10(b) shows the total P loss, reduced from 159.61 kW to 52.52 kW, while the percentage decrease in P loss is 67.09%. Figure 10(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 96.89 kVAR to 41.63 kVAR, while the percentage decrease of total Q loss is 57.03%. Optimal DG sizes are found as 728.58, 612.37, and 1092.6 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 10
Case 3 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 11
Results of case 3 (33-bus RDN with HHO).
Seasonvmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.952570.3655952.52841.631567.090557.0347728.58, 612.37, 1092.63,6,8Without DGs0.89587—159.613696.8958——Spring0.919070.4185144.1232115.724570.860861.71381161,1177.6, 1966.13,6,8Without DGs0.81372—494.6027302.2616————Summer0.907250.44344187.6912150.771571.941163.22781419.6, 1288.6, 2287.93,6,8Without DGs0.7819—668.9179410.0154————Autumn0.94170.3797777.198861.623968.543558.7164804.21,826.6,1369.43,6,8Without DGs0.87025—245.4149149.2695——


#### 4.1.10. Case 3 (33-Bus RDN) for Spring

Figure11(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.91907 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.4185. Figure 11(b) shows the total P loss, reduced from 494.60 kW to 144.12 kW, while the percentage decrease of total P loss is 70.86%. Figure 11(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 302.26 kVAR to 115.72 kVAR, while the percentage decrease of total Q loss is 61.71%. Optimal DG sizes are found as 1161, 1177.6, and 1966.1 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 11
Case 3 (33-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.11. Case 3 (33-Bus RDN) for Summer

Figure12(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.90725 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.44344. Figure 12(b) shows the total P loss reduced from 668.91 kW to 187.69 kW, while the percentage decrease of total P loss is 71.94%. Figure 12(c) depicts the entire Q loss profile. For this case, the whole Q loss is reduced from 410.01 kVAR to 150.77 kVAR, while the total Q loss percentage decrease is 63.22%. Optimal DG sizes are found as 1419.6, 1288.6, and 2287.9 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 12
Case 3 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.12. Case 3 (33-Bus RDN) for Autumn

Figure13(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.9417 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.37977. Figure 13(b) shows the total P loss, which is reduced from 245.41 kW to 77.19 kW, while the percentage decrease of total P loss is 68.54%. Figure 13(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 149.26 kVAR to 61.62 kVAR, while the percentage decrease of total Q loss is 58.71%. Optimal DG sizes are found as 804.21, 826.6, and 1369.4 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 13
Case 3 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.13. Case 4 (33-Bus RDN) for Winter

Figure14(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.95257 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.37277. Figure 14(b) shows the total P loss, reduced from 159.61 kW to 53.74 kW, while the percentage decrease of total P loss is 66.32%. Figure 14(c) depicts the entire Q loss profile. For this case, the whole Q loss is reduced from 96.89 kVAR to 42.04 kVAR, while the total Q loss percentage decrease is 56.60%. Optimal DG sizes are 579.36, 587.25, 1102.9, and 174.16 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 14
Case 4 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 12
Results of case 4 (33-bus RDN with HHO).
Seasonvmin(per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.952570.3727753.748142.045966.326156.607579.36, 587.25, 1102.9, 174.163,6,8,4Without DGs0.89587—159.613696.8958————Spring0.919070.44263151.6158117.466169.345961.1376987.93, 1116.1, 1997.2, 290.953,6,8,4Without DGs0.81372—494.6027302.2616————Summer0.907250.47614199.115153.377470.233362.59231246, 1376.5, 2300.6, 188.483,6,8,4Without DGs0.7819—668.9179410.0154———Autumn0.94170.391379.703962.253567.522858.2945749.27, 753.49, 1389.8, 184.743,6,8,4Without DGs0.87025—245.4149149.2695———


#### 4.1.14. Case 4 (33-Bus RDN) for Spring

Figure15(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.91907 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.44263. Figure 15(b) shows the total P loss, reduced from 494.60 kW to 151.61 kW, while the percentage decrease in P loss is 69.34%. Figure 15(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 302.26 kVAR to 117.46 kVAR, while the total Q loss percentage decrease is 61.1376%. Optimal DG sizes are 987.93, 1116.1, 1997.2, and 290.95 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 15
Case 4 (33-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.15. Case 4 (33-Bus RDN) for Summer

Figure16(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.90725 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.47614. Figure 16(b) shows the total P loss, reduced from 668.91 kW to 199.11 kW, while the percentage decrease of total P loss is 70.23%. Figure 16(c) illustrates the entire Q loss profile. For this case, the complete Q loss is reduced from 410.01 kVAR to 153.37 kVAR, while the total Q loss percentage decrease is 62.59%. Optimal DG sizes are 1246, 1376.5, 2300.6, and 188.48 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 16
Case 4 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.1.16. Case 4 (33-Bus RDN) for Autumn

Figure17(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.9417 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.3913. Figure 17(b) shows the total P loss, which is reduced from 245.41 kW to 79.70 kW, while the percentage decrease of total P loss is 67.52%. Figure 17(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 149.26 kVAR to 62.25 kVAR, while the total Q loss percentage decrease is 58.29%. Optimal DG sizes are 749.27, 753.49, 1389.8, and 184.74 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 17
Case 4 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2. IEEE 69-Bus RDN

Figure18 shows IEEE 69-bus RDN [80]. The total active load demand capacity of this IEEE 69-bus RDN is 3802 kW, while this system’s total reactive load demand is 2696 kVAR. The total P loss is 225.007 kW, and the bus voltage is 0.9091 (per unit) [69].Figure 18
IEEE 69-bus RDN.Table13 compares the literature simulation results with the proposed HHO based on four DGs (IEEE 33-bus RDN). It is observed that the maximum total P loss reduction of 72.75% is obtained as compared to the other literature results. The total Q loss reduction of 77.86% is obtained, which is higher than different algorithms results. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total active and Q loss reduction is 28.48 kW and 22.62 kVAR. In comparison, the bus voltage magnitude is upgraded to 0.97827 (per unit). The objective function for this case is 0.35001.Table 13
Results of literature with four DGs (69-bus RDN).
Algorithm Ref.Yearvmin(per unit)P loss (kW)PLOSSreduction (%)Q loss (kVAR)QLOSSreduction (%)DG sizes (kW)DG’s location(bus no.)Proposed HHO20210.9782728.481872.753722.6277.86329.88,229.42,112.03,924.1357,7,6,58PPA [69]20200.9782917.2392.3415——1014,47.22,145,1586.37,57,58,61GA [74]20200.97018.4360.4014.3259.30929.00,1075.0,984.80,1321.321,62,64,45PSO [74]2020—20.3163.0215.4554.301199.0,795.60,992.50,889.2061,63,17,32GA/PSO [74]2020—19.5663.9014.3260.20910.50,1192.0,884.90,978.6021,61,63,57HSA [65]20200.96734.2161.40—58.301302.0,369.00,101.80,978.5063,64,65,56BFOA [75]2020——66.56—61.50295.40,447.00,1345.1,1000.027,65,61,32ABC [65]20200.96032.0485.8019.9876.401000.0,200.00,338.20,1100.061,51,62,32SSA [40]20190.97023.4069.10—62.40380.00,527.00,1718.0,467.0017,10,60,32WCA [39]20180.97022.4368.20—61.20775.00,1105.0,438.00,890.0061,62,23,21HGWO [79]20130.96020.4169.14—64.30527.00,380.00,1718.0,398.5011,17,61,45Table14 compares the literature simulation results with the proposed HHO based on three DGs (IEEE 33-bus RDN). It is observed that the maximum total P loss reduction of 68.55% is obtained compared to the other literature results. The total Q loss reduction of 69.86% is obtained, which is higher than different algorithms results. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total active and Q loss reduction is 32.87 kW and 30.79 kVAR. At the same time, the bus voltage magnitude is upgraded to 0.9710 (per unit). The objective function for this case is 0.35373.Table 14
Results of literature with three DGs (69-bus RDN).
Algorithm Ref.Yearvmin (per unit)P loss (kW)PLOSS reduction (%)Q loss (kVAR)QLOSS reduction (%)DG sizes (kW)DG’s location (bus no.)Proposed HHO20210.9710032.871968.55430.788369.861764.3,617.55,399.6657,7,6PPA [69]20200.98240.5082——27.919,1108.3,1558.351,61,62GA [74]20200.970019.87560.407.11252.2929.70,1075.0,984.8021,62,64PSO [74]2020—19.43063.0212.32051.21199.0,795.60,992.5061,63,17GA/PSO [74]2020—28.72063.9013.57053.0910.50,1192.0,884.9021,61,63HSA [65]20200.967035.43061.40—51.01302.4,369.00,101.8063,64,65BFOA [75]20200.9608—66.56—56.5295.40,447.60,1345.027,65,61ABC [65]20200.960032.05686.8019.98077.31000.0,200.00,338.2061,51,62QOCSOS [42]2020—69.4369.14——526.8,380.4,1719.011,18,61CSCA [14]20200.9870.2068.8——365.9,1675.85,65.5217,61,67SSA [40]20190.978969.5269.1——380,527,171817,10,60FWA [41]20190.97477.8765.39——480.5,1198.6,225.865,61,27SFLA [41]20190.975277.78465.43——1088.7,167.3,980.957,63,26SSA [40]20190.978927.54069.10—59.2380.00,527.00,1718.017,10,60WCA [39]20180.98771.5568.2——775,1105,43861,62,23WCA [39]20180.967021.45068.20—51.1775.00,1105.0,438.0061,62,23HGWO [37]20170.9869.4369.14——527,380,171811,17,61BFOA [34]20140.980875.2466.56——295.4,447.6,1345.127,65,61HSA [33]20130.96786.8561.4——1302.4,369,101.863,64,65HGWO [79]20130.960024.23069.14—59.1527.00,380.00,1718.011,17,61Table15 compares the literature simulation results with the proposed HHO based on two DGs (IEEE 33-bus RDN). It is observed that the maximum total P loss reduction of 67.16% is obtained as compared to the other literature results. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total P loss reduction is 34.32 kW, while the bus voltage magnitude is upgraded to 0.9710 (per unit). The objective function for this case is 0.35361.Table 15
Comparison of literature results with two DGs (case 2 (69-bus RDN)).
Algorithm Ref.Yearvmin (per unit)P loss (kW)PLOSS (%)DG’s size (kW)DG’s locationFOALO [24]20170.9801070.7768.5470538.70,1700.017,610.83200ALO [26]20180.9789071.6868.1400531.48,1781.517,610.83000EMA [29]20180.9794071.8668.06001886.9,649.3061,690.83000CSFS-III [27]2018—71.6868.1400531.00,1781.017,610.34000PPA [69]20200.9710569.1169.24001040.0,1856.07,570.83500Proposed HHO20210.9710034.324567.16441232.0,414.9957,70.35361Table16 compares the literature simulation results with the proposed HHO based on one DG (IEEE 33-bus RDN). It is observed that the maximum total P loss reduction of 65.56% is obtained compared to the other literature results. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total P loss reduction is 36.01 kW. In comparison, the bus voltage magnitude is upgraded to 0.9710 (per unit). The objective function for this case is 0.35442.Table 16
Comparison of literature results with one DG (case 1 (69-bus RDN)).
Algorithm Ref.Yearvmin (per unit)P loss (kW)PLOSS (%)DG’s size (kW)DG’s locationFOALO [24]20170.9679081.801063.64501800.0610.80200ALO [26]20180.9682083.250063.00001872.7610.79900EMA [29]20180.9689082.870063.17001910.3570.80000GA [1]20190.9682063.000072.00001850.0600.84400CSFS-III [27]2018—83.230063.01001873.0610.31500CDE [30]20190.9683083.250063.00001872.4610.79900PPA [69]20200.9710574.250067.00002088.0570.82900Proposed HHO20210.9710036.006965.55511330.2570.35442


#### 4.2.1. Case 1 (69-Bus RDN) for Winter

Figure19(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for winter load. The bus voltage magnitude is upgraded from 0.93825 (per unit) to 0.97991 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.36388. Figure 19(b) shows the total P loss reduced from 104.53 kW to 35.74 kW while the percentage decrease of total P loss is 65.80%. Figure 19(c) depicts the entire Q loss profile. For this case, the whole Q loss is reduced from 47.61 kVAR to 15.93 kVAR, while the total Q loss percentage decrease is 66.53%. The optimal DG size is found as 1346.9 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 19
Case 1 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 17
Results with one DG (case 1 (69-bus RDN)) with proposed HHO.
Seasonvmin (per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.979910.3638835.748615.935565.802266.53081346.957Without DGs0.93825104.534747.6124———Spring0.966560.4103897.976443.653568.00368.55192357.257Without DGs0.89392306.2051138.8112———Summer0.962040.43178126.690856.419468.585469.08752710.857Without DGs0.87809403.2863182.5134———Autumn0.97550.3765152.724823.505166.670967.32841675.657Without DGs0.92395158.194671.9434———


#### 4.2.2. Case 1 (69-Bus RDN) for Spring

Figure20(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.96656 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.41038. Figure 20(b) shows the total P loss reduced from 306.20 kW to 97.97 kW while the percentage decrease of total P loss is 68.00%. Figure 20(c) illustrates the entire Q loss profile. For this case, the complete Q loss is reduced from 138.81 kVAR to 43.65 kVAR, while the total Q loss percentage decrease is 68.55%. The optimal DG size is found as 2357.2 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 20
Case 1 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.3. Case 1 (69-Bus RDN) for Summer

Figure21(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.96204 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.43178. Figure 21(b) shows the total P loss reduced from 403.28 kW to 126.69 kW while the percentage decrease of total P loss is 68.58%. Figure 21(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 182.51 kVAR to 56.41 kVAR, while the total Q loss percentage decrease is 69.08%. The optimal DG size is found as 2710.8 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 21
Case 1 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.4. Case 1 (69-Bus RDN) for Autumn

Figure22(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.9755 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.37651. Figure 22(b) shows the total P loss reduced from 158.1946 kW to 52.7248 kW, while the percentage decrease of total P loss is 66.67%. Figure 22(c) depicts the entire Q loss profile. For this case, the total Q loss is reduced from 71.94 kVAR to 23.50 kVAR, while the total Q loss percentage decrease is 67.32%. The optimal DG size is found as 1675.6 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 22
Case 1 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.5. Case 2 (69-Bus RDN) for Winter

Figure23(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for winter load. The bus voltage magnitude is upgraded from 0.93825 (per unit) to 0.97991 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.37332. Figure 23(b) shows the total P loss, which increased from 104.5347 kW to 346.82 kW, while the percentage increase of total P loss is 231.77%. Figure 23(c) depicts the entire Q loss profile. The total Q loss for this case is increased from 47.61 kVAR to 155.20 kVAR, while the percentage increase of total Q loss is 225.97%. Optimal DG sizes are found as 868.32 and 419.28 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 23
Case 2 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 18
Results with two DGs (case 2 (69-bus RDN)) with proposed HHO.
Seasonvmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)DG’s size (kW)DG’s location (bus no.)Winter0.979910.37332346.8215155.2039−231.7765−225.9733868.32, 419.2857, 7Without DGs0.93825104.534747.6124———Spring0.966560.44221122.0975504.2769−266.4529−263.28261484.7, 788.6657, 7Without DGs0.89392306.2051138.8112———Summer0.962040.47521542.1811694.2981−282.4035−280.40941699.7, 923.4257, 7Without DGs0.87809403.2863182.5134———Autumn0.97550.39162540.5331242.1968−241.6888−236.64921068.5, 540.7757, 7Without DGs0.92395—158.194671.9434————


#### 4.2.6. Case 2 (69-Bus RDN) for Spring

Figure24(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.96656 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.4422. Figure 24(b) shows the total P loss, which is increased from 306.2051 kW to 1122.09 kW, while the percentage increase of total P loss is 266.45%. Figure 24(c) shows the entire Q loss profile. The total Q loss for this case is increased from 138.81 kVAR to 504.27 kVAR, while the percentage increase of total Q loss is 263.28%. Optimal DG sizes are found as 1484.7 and 788.66 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 24
Case 2 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.7. Case 2 (69-Bus RDN) for Summer

Figure25(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.96204 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.4752. Figure 25(b) shows the total P loss, which is increased from 403.28 kW to 1542.18 kW, while the percentage increase of total P loss is 282.40%. Figure 25(c) shows the entire Q loss profile. The total Q loss for this case is increased from 182.51 kVAR to 694.29 kVAR, while the percentage increase of total Q loss is 280.40%. Optimal DG sizes are found as 1699.7 and 923.42 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 25
Case 2 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.8. Case 2 (69-Bus RDN) for Autumn

Figure26(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.9755 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.39162. Figure 26(b) shows the total P loss, which increased from 158.19 kW to 540.53 kW, while the percentage increase of total P loss is 241.68%. Figure 26(c) depicts the entire Q loss profile. The total Q loss for this case is increased from 71.94 kVAR to 242.19 kVAR, while the percentage increase of total Q loss is 236.64%. Optimal DG sizes are found as 1068.5 and 540.77 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 26
Case 2 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.9. Case 3 (69-Bus RDN) for Winter

Figure27(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for winter load. The bus voltage magnitude is upgraded from 0.9091 (per unit) to 0.97991 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.36328. Figure 27(b) shows the total P loss, reduced from 225.00 kW to 33.27 kW, while the percentage decrease of total P loss is 33.27%. Figure 27(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 33.27 kVAR to 14.74 kVAR, while the total Q loss percentage decrease is 14.74%. Optimal DG sizes are found as 1212.5, 552.87, and 19.166 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 27
Case 3 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 19
Results with three DGs (case 3 (69-bus RDN)) with proposed HHO.
Seasonvmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.979910.3632833.276314.749768.167269.02141212.5, 552.87, 19.16657, 7, 6Without DGs0.9091225.00747.6124———Spring0.966560.379992.490140.904469.794770.53242108.7, 706.76, 247.2857, 7, 6Without DGs0.89392306.2051138.8112———Summer0.962040.39179117.109651.634270.961271.70942361.5, 1059.3, 128.5557, 7, 6Without DGs0.87809403.2863182.5134———Autumn0.97550.3607549.354221.853568.801669.6241483.8, 674.15, 23.57457, 7, 6Without DGs0.92395158.194671.9434———


#### 4.2.10. Case 3 (69-Bus RDN) for Spring

Figure28(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.96656 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.3799. Figure 28(b) shows the total P loss reduced from 306.20 kW to 92.49 kW while the percentage decrease of total P loss is 69.79%. Figure 28(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 138.81 kVAR to 40.90 kVAR, while the percentage decrease in Q loss is 70.53%. Optimal DG sizes are found as 2108.7, 706.76, and 247.28 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 28
Case 3 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.11. Case 3 (69-Bus RDN) for Summer

Figure29(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.96204 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.39179. Figure 29(b) shows the total P loss, reduced from 403.28 kW to 117.10 kW, while the percentage decrease of total P loss is 70.96%. Figure 29(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 70.96 kVAR to 51.63 kVAR, while the percentage decrease is 71.70%. Optimal DG sizes are found as 2361.5, 1059.3, and 128.55 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 29
Case 3 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.12. Case 3 (69-Bus RDN) for Autumn

Figure30(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.9755 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.36075. Figure 30(b) shows the total P loss reduced from 158.19 kW to 49.35 kW while the percentage decrease of total P loss is 68.80%. Figure 30(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 71.94 kVAR to 21.85 kVAR, while the total Q loss percentage decrease is 69.62%. Optimal DG sizes are found as 1483.8, 674.15, and 23.574 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 30
Case 3 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.13. Case 4 (69-Bus RDN) for Winter

Figure31(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for winter load. The bus voltage magnitude is upgraded from 0.93825 (per unit) to 0.9849 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.36025. Figure 31(b) shows the total P loss reduced from 104.53 kW to 22.76 kW while the percentage decrease of total P loss is 78.22%. Figure 31(c) depicts the entire Q loss profile. For this case, the total Q loss is reduced from 47.61 kVAR to 11.21 kVAR, while the total Q loss percentage decrease is 76.45%. Optimal DG sizes are 40.66, 553.91, 19.998, and 1156.8 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 31
Case 4 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 20
Results with four DGs (case 4 (69-bus RDN)) with proposed HHO.
Seasonvmin(per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no).Winter0.98490.3602522.762611.212578.224876.450540.66, 553.91, 19.998, 1156.857, 7, 6, 58Without DGs0.93825104.534747.6124———Spring0.974960.3992161.156230.042480.027778.3574108.65, 1032.2, 31.762, 1957.557, 7, 6, 58Without DGs0.89392306.2051138.8112———Summer0.971610.41778.855838.683680.446778.805123.63, 1198.9, 41.268, 2240.557, 7, 6, 58Without DGs0.87809403.2863182.5134———Autumn0.981610.3708933.246816.367978.983677.248954.404, 700.2, 39.572, 1427.157, 7, 6, 58Without DGs0.92395158.194671.9434———


#### 4.2.14. Case 4 (69-Bus RDN) for Spring

Figure32(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.97496 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.39921. Figure 32(b) shows the total P loss, reduced from 306.20 kW to 61.15 kW, while the percentage decrease of total P loss is 80.02%. Figure 32(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 138.81 kVAR to 30.04 kVAR, while the percentage decrease in Q loss is 78.35%. Optimal DG sizes are 108.65, 1032.2, 31.762, and 1957.5 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 32
Case 4 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.15. Case 4 (69-Bus RDN) for Summer

Figure33(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.97161 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.417. Figure 33(b) shows the total P loss, which is reduced from 403.28 kW to 78.85 kW, while the percentage decrease of total P loss is 80.44%. Figure 33(c) shows the total Q loss profile. For this case, the total Q loss is reduced from 182.51 kVAR to 38.68 kVAR, while the percentage decrease in Q loss is 78.80%. Optimal DG sizes are 123.63, 1198.9, 41.268, and 2240.5 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 33
Case 4 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


#### 4.2.16. Case 4 (69-Bus RDN) for Autumn

Figure34(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.98161 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.37089. Figure 34(b) shows the total P loss, reduced from 158.19 kW to 33.24 kW, while the percentage decrease of total P loss is 33.24%. Figure 34(c) shows the total Q loss profile. For this case, the total Q loss is reduced from 71.94 kVAR to 16.36 kVAR, while the percentage decrease in Q loss is 77.24%. Optimal DG sizes are 54.404, 700.2, 39.572, and 1427.1 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 34
Case 4 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.3. Comparison of Cases

Figures35(a) and 35(b) show the bus voltage and total P loss profiles of 33-bus RDN with the deployment of one, two, three, and four DGs for all four cases, respectively.Figure 35
Comparison of (a) bus voltage profiles and (b) totalP loss profiles in IEEE 33-bus RDN.
(a)(b)For winter load, the bus voltage drop is minimized in case 3 of 33-bus RDN, which is 0.95257 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 3 of 33-bus RDN is 52.52 kW, while the percentage decrease of total P loss, in this case, is 67.09% which shows that deployment of three DGs in a 33-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 3 of 33-bus RDN is 41.63 kVAR, while the percentage decrease of total Q loss is 57.03% which shows that deployment of three DGs in a 33-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of 33-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 245.46 kW and 142.43 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 53.78% and 46.99%, respectively, showing that deploying two DGs in 33-bus RDS is not feasible to minimize total P and Q losses with improved bus voltage. Table 21 presents the result summary for all four cases of 33-bus RDN with seasonal load data of winter.Table 21
Result comparison of 33-bus RDN for winter.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.909980.39695121.210776.81724.059920.7221Case 20.91310.4016245.4671142.4336−53.7883−46.9967Case 30.952570.3655952.52841.631567.090557.0347Case 40.952570.3727753.748142.045966.326156.607For spring load, the bus voltage drop is minimized in case 3 of 33-bus RDN, 0.91907 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 3 of 33-bus RDN is 144.12 kW, while the percentage decrease of total P loss, in this case, is 144.12% which shows that deployment of three DGs in a 33-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 3 of 33-bus RDN is 144.12 kVAR, while the percentage decrease of total Q loss is 144.12% which shows that deploying three DGs in a 33-bus system is still more feasible, minimizing total Q loss with improved bus voltage. However, case 2 of 33-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 736.49 kW and 428.92 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 48.90% and 41.90%, respectively, which shows that deploying two DGs in 33-bus RDS is not feasible in minimizing total P and Q losses with improved bus voltage. Table 22 presents the result summary for all four cases of 33-bus RDN with seasonal load data of spring.Table 22
Result comparison of 33-bus RDN for spring.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.840890.52997368.0215234.536425.592522.4061Case 20.847380.54649736.4943428.922−48.9062−41.9042Case 30.919070.4185144.1232115.724570.860861.7138Case 40.919070.44263151.6158117.466169.345961.1376For summer load, the bus voltage drop is minimized in case 3 of 33-bus RDN, which is 0.90725 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 3 of 33-bus RDN is 187.69 kW, while the percentage decrease of total P loss, in this case, is 187.69% which shows that deployment of three DGs in a 33-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 3 of 33-bus RDN is 187.69 kVAR, while the percentage decrease of total Q loss is 63.22% which shows that deployment of three DGs in a 33-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of 33-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 63.22 kW and 571.85 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 46.57% and 39.47%, respectively, which shows that deploying two DGs in 33-bus RDS is not feasible in minimizing total P and Q losses with improved bus voltage. Table 23 presents the result summary for all four cases of 33-bus RDN with seasonal load data of summer.Table 23
Result comparison of 33-bus RDN for summer.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.814840.59875495.3813316.116125.942922.9014Case 20.82290.6219980.4337571.8549−46.5701−39.4716Case 30.907250.44344187.6912150.771571.941163.2278Case 40.907250.47614199.115153.377470.233362.5923For autumn load, the bus voltage drop is minimized in case 3 of 33-bus RDN, which is 0.9417 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 3 of 33-bus RDN is 77.19 kW, while the percentage decrease of total P loss, in this case, is 68.54% which shows that deployment of three DGs in a 33-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 3 of 33-bus RDN is 61.62 kVAR, while the percentage decrease of total Q loss is 58.71% which shows that deployment of three DGs in a 33-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of 33-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 374.05 kW and 217.28 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 52.41% and 45.56%, respectively, which shows that deploying two DGs in 33-bus RDS is not feasible in minimizing total P and Q losses with improved bus voltage. Table 24 presents the result summary for all four cases of 33-bus RDN with seasonal load data of autumn.Table 24
Result comparison of 33-bus RDN for autumn.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.88820.43087184.8168117.382324.692121.3622Case 20.892290.43844374.0595217.289−52.4192−45.5682Case 30.94170.3797777.198861.623968.543558.7164Case 40.94170.391379.703962.253567.522858.2945Figures36(a) and 36(b) show the bus voltage and total P loss profiles of 33-bus RDN with the deployment of one, two, three, and four DGs for all four cases, respectively.Figure 36
Comparison of (a) bus voltage profiles and (b) totalP loss profiles in IEEE 69-bus RDN.
(a)(b)For winter load, the bus voltage drop is minimized in case 4 of IEEE 69-bus RDN, which is 0.9849 (per unit), which shows more improved value than other cases. The total P loss reduction for case 4 of IEEE 69-bus RDN is 22.76 kW, while the percentage decrease of totalP loss, in this case, is 78.22% which shows that deployment of three DGs in a 69-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 4 of IEEE 69-bus RDN is 11.21 kVAR, while the percentage decrease of total Q loss is 76.45% which shows that deployment of three DGs in a 69-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of IEEE 69-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 346.82 kW and 155.20 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 231.77% and 225.97%, respectively, showing that deploying two DGs in 69-bus RDS is not feasible to minimize total P and Q losses with improved bus voltage. Table 25 presents the result summary for all four cases of IEEE 69-bus RDN with seasonal load data of winter.Table 25
Result comparison of IEEE 69-bus RDN for winter.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.979910.3638835.748615.935565.802266.5308Case 20.979910.37332346.8215155.2039−231.7765−225.9733Case 30.979910.3632833.276314.749768.167269.0214Case 40.98490.3602522.762611.212578.224876.4505For spring load, the bus voltage drop is minimized in case 4 of IEEE 69-bus RDN, 0.97496 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 4 of IEEE 69-bus RDN is 61.15 kW, while the percentage decrease of total P loss, in this case, is 80.02% which shows that deployment of three DGs in a 69-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 4 of IEEE 69-bus RDN is 30.04 kVAR, while the percentage decrease of total Q loss is 78.35% which shows that deployment of three DGs in a 69-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of IEEE 69-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 1122.09 kW and 504.27 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 266.45% and 263.28%, respectively, showing that deploying two DGs in 69-bus RDS is not feasible in minimizing total P and Q losses with improved bus voltage. Table 26 presents the result summary for all four cases of IEEE 69-bus RDN with seasonal load data of spring.Table 26
Result comparison of IEEE 69-bus RDN for spring.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.966560.4103897.976443.653568.00368.5519Case 20.966560.44221122.0975504.2769−266.4529−263.2826Case 30.966560.379992.490140.904469.794770.5324Case 40.974960.3992161.156230.042480.027778.3574For summer load, the bus voltage drop is minimized in case 4 of IEEE 69-bus RDN, which is 0.97161 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 4 of IEEE 69-bus RDN is 78.85 kW, while the percentage decrease of total P loss, in this case, is 80.44% which shows that deployment of three DGs in a 69-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 4 of IEEE 69-bus RDN is 38.68 kVAR, while the percentage decrease of total Q loss is 78.80% which shows that deployment of three DGs in a 69-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of IEEE 69-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 1542.18 kW and 694.29 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 282.40% and 280.40%, respectively, showing that deploying two DGs in 69-bus RDS is not feasible to minimize total P and Q losses with improved bus voltage. Table 27 presents the result summary for all four cases of IEEE 69-bus RDN with seasonal load data of summer.Table 27
Result comparison of IEEE 69-bus RDN for summer.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.962040.43178126.690856.419468.585469.0875Case 20.962040.47521542.1811694.2981−282.4035−280.4094Case 30.962040.39179117.109651.634270.961271.7094Case 40.971610.41778.855838.683680.446778.805For autumn load, the bus voltage drop is minimized in case 4 of IEEE 69-bus RDN, which is 0.98161 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 4 of IEEE 69-bus RDN is 33.24 kW, while the percentage decrease of total P loss, in this case, is 78.98% which shows that deployment of three DGs in a 69-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 4 of IEEE 69-bus RDN is 16.36 kVAR, while the percentage decrease of total Q loss is 77.24% which shows that deployment of three DGs in a 69-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of IEEE 69-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 540.53 kW and 242.19 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 241.68% and 236.64%, respectively, showing that deploying two DGs in 69-bus RDS is not feasible in minimizing total P and Q losses with improved bus voltage. Table 28 presents the result summary for all four cases of IEEE 69-bus RDN with seasonal load data of autumn.Table 28
Result comparison of IEEE 69-bus RDN for autumn.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.97550.3765152.724823.505166.670967.3284Case 20.97550.39162540.5331242.1968−241.6888−236.6492Case 30.97550.3607549.354221.853568.801669.624Case 40.981610.3708933.246816.367978.983677.2489Moreover, Figures37 and 38 illustrate the graphical comparison of literature results with four and three DGs for IEEE 33-bus RDN, respectively. Figures 39 and 40 illustrate the graphical comparison of literature results with four and three DGs for IEEE 69-bus RDN, respectively. For four DGs in 33-bus RDN, the voltage profile improvement of the proposed HHO is less than that of other literature-based algorithms such as GA, BA, PSO, PPA, ABC and TLCHS. Active power loss reduction is better with the proposed HHO. For three DGs in 33-bus RDN, the voltage profile improvement of the proposed HHO is better than that of other literature-based algorithms such as GA, BA, GSA-GAMS, SSA, EMA, SPEA2, CSCA, PPA, and ACO-ABC. Active power loss reduction is higher with the proposed HHO. For four DGs in 69-bus RDN, the voltage profile improvement of the proposed HHO is better than that of other literature-based algorithms such as GA, SSA, PPA, WCA, HAS, HGWO, and ABC. Active power loss reduction is higher with the proposed HHO, PPA, and ABC. For three DGs in 69-bus RDN, the voltage profile improvement of the proposed HHO is better than that of other literature-based algorithms. Active power loss reduction is higher with the proposed HHO, PPA, and ABC.Figure 37
Graphical comparison of literature results with four DGs (33-bus RDN).Figure 38
Graphical comparison of literature results with three DGs (33-bus RDN).Figure 39
Graphical comparison of literature results with four DGs (69-bus RDN).Figure 40
Graphical comparison of literature results with three DG (69-bus RDN).


### 4.4. Statistical Analysis

This section investigates the statistical analysis for the objective functions taken using the HHO algorithm. The statistical analysis is carried out on the objective function in which voltage and PQ loss reduction (LR) is part of the objective function. The formula for real and reactive power loss is taken from [55]. The technique to calculate the losses is applied from existing literature [55, 81, 82]. The algorithm is simulated fifteen times during the load flow analysis for both IEEE RDNs, and results from the Big-O analysis are shown in Table 29 [54]. It is observed that the changes in the objective functions for the maximum and minimum values taken from the HHO have no significant variations. It validates the efficient and robust behavior of the applied HHO method. However, case 2 of IEEE 33-bus RDN shows a significant difference among best, worst, and mean values of total P loss. It shows that case 2 of 33-bus RDN cannot be taken under consideration during the planning stages. Moreover, case 4 of IEEE 69-bus RDN shows a significant difference among best, worst, and mean values of minimum bus voltages. It also shows that case 4 of IEEE 69-bus RDN cannot be considered during the planning stages.Table 29
Statistical analysis of HHO.
CasesOFBest caseWorst caseMean caseIEEE 33-bus RDNCase 1vmin (per unit)0.865240.865240.86524Case 1PLOSS0.2520120.2190340.245402Case 2vmin (per unit)0.90270.880340.897694Case 2PLOSS0.5804640.2777590.519263Case 3vmin (per unit)0.930540.930540.93054Case 3PLOSS0.6971690.6786740.688568Case 4vmin (per unit)0.930540.930540.93054Case 4PLOSS0.6958930.6784370.683718IEEE 69-bus RDNCase 1vmin (per unit)0.971050.971050.97105Case 1PLOSS0.6732650.6730480.673112Case 2vmin (per unit)0.971050.971050.97105Case 2PLOSS0.6964950.6853160.691088Case 3vmin (per unit)0.971050.971050.97105Case 3PLOSS0.6955840.6804380.690053Case 4vmin (per unit)0.978360.978280.978318Case 4PLOSS0.7748490.7076720.73695Table30 illustrates the analysis of variance (ANOVA). ANOVA tests help determine the variance of the objective function under various algorithms [83]. For IEEE 33-bus RDN, the ANOVA test has been carried out between fourteen (14) algorithms, while for IEEE 69-bus RDN, the test has been performed between twenty (20) algorithms. The value of F for both RDNs is less than the critical level. It proves that the variation found while computing the objection function is significant and not by chance [84].Table 30
ANOVA test.
Source of varianceSum of square (SS)Degree of freedom (df)Mean Square (MS)F ratioP valueF critical [83]Between groups0.008810.00880.30070.58734.1596Within groups0.9116140.0294IEEE 69-bus RDNBetween groups0.073030.02430.71770.54712.8327Within groups1.3903200.0339Moreover, Figures41 and 42 illustrate the box plots for IEEE 33-bus and 69-bus RDNs, respectively. It is apparent that no outliers are present and significant data are available for statistical analysis.Figure 41
Box plot to compare objective function for four cases in IEEE 33-bus RDN.Figure 42
Box plot to compare objective function for four cases in IEEE 69-bus RDN.


## 4.1. IEEE 33-Bus RDN

Figure1 [63] shows IEEE 33-bus RDN [63], which comprises thirty-seven lines connected with 33 buses. The total active load demand is 3720 kW, while the total reactive load demand is 2300 kVAR. The total P loss is 208.4592 kW, and the bus voltage magnitude is 0.929 (per unit) [69].Figure 1
IEEE 33-bus RDN.Table3 compares the minimum and average bus voltage values of literature with the results of proposed HHO for 33-bus RDN for exponent values of α = 0.720, β = 2.960. Table 4 compares the minimum and average bus voltage values of literature with the results of proposed HHO for 33-bus RDN for exponent values of α = 0.920, β = 4.040. Table 5 compares the minimum and average bus voltage values of literature with the results of proposed HHO for 33-bus RDN for exponent values of α = 1.040, β = 4.190. Table 6 compares the minimum and average bus voltage values of literature with the results of proposed HHO for 33-bus RDN for exponent values of α = 1.300, β = 4.380.Table 3
Simulation results of system I (α = 0.720, β = 2.960).
AlgorithmVave (per unit)vmin (per unit)FWA [70]0.96940.9459MBFOA [71]0.96980.9430ITS [72]0.96770.9397SLR [73]0.96750.9397PSO [68]0.95570.9221MPSO [68]0.96750.9432Proposed HHO0.97860.9391Table 4
Simulation results of system I (α = 0.920, β = 4.040).
AlgorithmVave (per unit)vmin (per unit)FWA [70]0.97000.9474MBFOA [71]0.97050.9446ITS [72]0.96840.9416SLR [73]0.96820.9416PSO [68]0.95520.9246MPSO [68]0.96820.9449Proposed HHO0.97580.9472Table 5
Simulation results of system I (α = 1.040, β = 4.190).
AlgorithmVave (per unit)vmin (per unit)FWA [70]0.97020.9477MBFOA [71]0.97060.9450ITS [72]0.96860.9420SLR [73]0.96840.9420PSO [68]0.95560.9254MPSO [68]0.96840.9453Proposed HHO0.98430.9417Table 6
Simulation results of system I (α = 1.300, β = 4.380).
AlgorithmVave (per unit)vmin (per unit)FWA [70]0.97050.9484MBFOA [71]0.97090.9457ITS [72]0.96890.9428SLR [73]0.96880.9428PSO [68]0.95640.9268MPSO [68]0.96880.9460Proposed HHO0.98710.9526Table7 compares the literature simulation results with the proposed HHO based on four DGs (IEEE 33-bus RDN). It is observed that the maximum active power loss reduction of 66.8% is obtained compared to the other literature results except for one case with the ABC algorithm. The total Q loss reduction of 56.4% is obtained, higher than different algorithm results except for two cases with BFOA and ABC methods. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total active and reactive power loss reduction is 52.98 kW and 42.23 kVAR. At the same time, the bus voltage magnitude is upgraded to 0.9471 (per unit). The objective function for this case is 0.36579.Table 7
Comparison of literature results with four DGs (33-bus RDN).
AlgorithmRef.Yearvmin(per unit)P loss(kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DGsizes (kW)DG’s location(bus no.)Proposed HHO20210.9471452.984342.22866.804756.4191318.67,42.94,997.65,322.883, 6, 8, 4PPA [69]20200.9790058.36—72—1058.0,2201.0,171.70,74.8003,6,4,5GA [74]20200.981064.19019.933047.6050.701500.0,422.80,1071.4,5781.011,29,30,8PSO [74]20200.980070.32016.542048.2050.401176.0,981.60,829.70,959.008,13,32,6GA/PSO [74]20200.967069.21021.320049.2043.70925.00,863.00,1200.0,1400.011,16,32,5HSA [65]20200.967077.140—52.30—572.40,107.00,1046.2,1034.017,18,33,15IWD [14]20200.969668.720—57.70—600.30,300.00,1011.2,1350.09,16,30,13BFOA [75]20200.9640——51.5057.30633.00,90.000,947.00,1000.017,18,33,11IWO [14]20200.971668.710—57.7054.30624.00,104.90,1056.0,879.4014,18,32,22PSO & analytical [76]2020—68.43016.540064.1052.30790.00,1070.0,1010.0,950.0013,24,30,20SKHA [77]20200.968773.450—64.7052.10801.80,1091.0,1053.6,950.0013,24,30,43ABC [65]20200.977079.50044.854067.3058.40514.00,948.00,635.16,479.5014,28,23,3SSA [40]20190.968667.430—64.80—753.60,1100.4,1070.0,570.0013,23,29,5CDE [30]20190.9702067.12—67.8—926.69, 646.70, 967.30,679.306,14,24,31WCA [39]20180.973074.530—64.00—854.6,1101.7,1180.0,750.0014,24,29,11BSOA [75]20180.955465.350—56.1053.4032.000,487.00,550.00,1000.013,28,31,7BA [78]20180.980066.430—63.0051.00816.00,952.50,952.35,598.0015,25,30,2TLCHS [43]20180.9770067.28—67.72—941.20,684.70,966.40,710.706,14,24,31ACO-ABC [68]20170.973571.670—62.8055.40754.00,1099.9,1071.0,589.4014,24,30,9HGWO [79]2013—77.430—64.4053.30802.00,1090.0,1054.0,800.0013,24,30,21Table8 compares the literature simulation results with the proposed HHO based on three DGs (IEEE 33-bus RDN). It is observed that the maximum active power loss reduction of 69.65% is obtained as compared to the other literature results. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total P loss reduction is 107.93 kW. At the same time, the bus voltage magnitude is upgraded to 0.9706 (per unit). The objective function for this case is 0.39752.Table 8
Comparison of literature results with three DGs (33-bus RDN).
Algorithm Ref.Yearvmin (per unit)FOP loss (kW)PLOSS (%)DG size (kW)DG’s locationACO-ABC [35]20150.97350.8007577.54062.800754.700,1099.90,1071.4014,24,30BA [36]20160.98000.8050077.12063.000816.300,952.350,952.35015,25,30HA [7]2016—0.3205074.83064.100790.000,1070.00,1010.0013,24,30HGWO [37]2017—0.3220074.20064.400802.000,1090.00,1054.0013,24,30EMA [29]20180.96840.8058074.37064.320976.600, 1169.09, 943.54030,24,12SPEA2 [38]20180.96160.8363060.24071.100691.000,733.400,742.90018,29,8SSA [40]20190.96860.8083073.37064.800753.600, 1100.40, 1070.6013,23,29GSA-GAMS [31]20190.96860.8125071.62065.640801.220, 1091.30, 1053.5913,24,30GA [44]20190.98600.8280068.78067.000761.000,1170.00,1082.0014,24,30QOCSOS [42]2020—0.3275071.91065.500801.700, 1091.30, 1053.6013,24,30CSCA [14]20200.96900.8070073.99064.500871.000, 1091.47, 954.08013,24,30IHSA [48]2020—0.3260072.54065.200800.800, 1087.60, 1050.7013,24,30PPA [69]20200.97000.8445058.41071.9001141.80, 161.71, 2214.3013,23,28Proposed HHO20210.97060.39752107.92569.6471304.8,751.02,1735.53,6,8


### 4.1.1. Case 1 (33-Bus RDN) for Winter

Figure2(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.90998 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.39695. Figure 2(b) shows the active power loss, reduced from 159.61 kW to 121.21 kW, while the decrease of total P loss is 24.06%. Figure 2(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 96.89 kVAR to 76.82 kVAR, while the decrease of total Q loss is 20.72%. The optimal DG size is found as 2578.5 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 2
Case 1 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 9
Results of case 1 (33-bus RDN with HHO).
Seasonvmin (per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.909980.39695121.210776.81724.059920.72212578.53Without DGs0.89587—159.613696.8958————Spring0.840890.52997368.0215234.536425.592522.40614888.13Without DGs0.81372—494.6027302.2616————Summer0.814840.59875495.3813316.116125.942922.90145805.53Without DGs0.7819—668.9179410.0154————Autumn0.88820.43087184.8168117.382324.692121.36223291.63Without DGs0.87025—245.4149149.2695————


### 4.1.2. Case 1 (33-Bus RDN) for Spring

Figure3(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.84089 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.52997. Figure 3(b) shows the total P loss, which is reduced from 494.60 kW to 368.02 kW, while the decrease of total P loss is 25.59%. Figure 3(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 302.26 kVAR to 234.53 kVAR, while the decrease of total Q loss is 22.40%. The optimal DG size is found as 4888.1 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 3
Case 1 (33-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.3. Case 1 (33-Bus RDN) for Summer

Figure4(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.81484 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.59875. Figure 4(b) shows the total P loss reduced from 668.91 kW to 495.38 kW, while the decrease of total P loss is 25.94%. Figure 4(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 410.01 kVAR to 316.11 kVAR, while the decrease of total Q loss is 22.90%. The optimal DG size is found as 5805.5 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 4
Case 1 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.4. Case 1 (33-Bus RDN) for Autumn

Figure5(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.8882 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.43087. Figure 5(b) shows the total P loss, which is reduced from 245.41 kW to 184.81 kW, while the decrease of total P loss is 24.69%. Figure 5(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 149.26 kVAR to 117.38 kVAR, while the decrease of total Q loss is 21.36%. The optimal DG size is found as 3291.6 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 5
Case 1 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.5. Case 2 (33-Bus RDN) for Winter

Figure6(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.9131 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.4016. Figure 6(b) shows the total P loss, which increased from 159.61 kW to 245.46 kW, while the percentage increase of total P loss is 53.78%. Figure 6(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 96.89 kVAR to 142.43 kVAR, while the percentage increase of total Q loss is 46.9967%. Optimal DG sizes are found as 347.21 kW and 989.34 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 6
Case 2 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 10
Results of case 2 (33-bus RDN with HHO).
Seasonvmin(per unit)FOP loss (kW)Q loss (kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’ssize (kW)DG’s location(bus no.)Winter0.91310.4016245.4671142.4336−53.7883−46.9967347.21, 989.343,6Without DGs0.89587—159.613696.8958————Spring0.847380.54649736.4943428.922−48.9062−41.9042715.63, 1844.43,6Without DGs0.81372—494.6027302.2616————Summer0.82290.6219980.4337571.8549−46.5701−39.4716873.16, 2181.43,6Without DGs0.7819—668.9179410.0154————Autumn0.892290.43844374.0595217.289−52.4192−45.5682457.26, 1254.13,6Without DGs0.87025—245.4149149.2695————


### 4.1.6. Case 2 (33-Bus RDN) for Spring

Figure7(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.84738 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.54649. Figure 7(b) shows the total P loss, which is increased from 494.60 kW to 736.49 kW, while the percentage increase of total P loss is 48.90%. Figure 7(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 302.26 kVAR to 428.9 kVAR, while the percentage increase of total Q loss is 41.90%. Optimal DG sizes are found as 715.63 and 1844.4 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 7
Case 2 (33-bus RDN) for spring. (a) Bus voltage profile. (b) Total P loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.7. Case 2 (33-Bus RDN) for Summer

Figure8(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.8229 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.6219. Figure 8(b) shows the total P loss, which increased from 668.91 kW to 980.43 kW, while the percentage increase of total P loss is 46.57%. Figure 8(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 410.01 kVAR to 571.85 kVAR, while the percentage increase of total Q loss is 39.47%. Optimal DG sizes are found as 873.16 and 2181.4 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 8
Case 2 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.8. Case 2 (33-Bus RDN) for Autumn

Figure9(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.89229 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.43844. Figure 9(b) shows the total P loss, which increased from 245.41 kW to 374.05 kW, while the percentage increase of total P loss is 52.41%. Figure 9(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 149.26 kVAR to 217.28 kVAR, while the percentage increase of total Q loss is 45.56%. Optimal DG sizes are found as 457.26 and 1254.1 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 9
Case 2 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.9. Case 3 (33-Bus RDN) for Winter

Figure10(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.95257 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.36559. Figure 10(b) shows the total P loss, reduced from 159.61 kW to 52.52 kW, while the percentage decrease in P loss is 67.09%. Figure 10(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 96.89 kVAR to 41.63 kVAR, while the percentage decrease of total Q loss is 57.03%. Optimal DG sizes are found as 728.58, 612.37, and 1092.6 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 10
Case 3 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 11
Results of case 3 (33-bus RDN with HHO).
Seasonvmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.952570.3655952.52841.631567.090557.0347728.58, 612.37, 1092.63,6,8Without DGs0.89587—159.613696.8958——Spring0.919070.4185144.1232115.724570.860861.71381161,1177.6, 1966.13,6,8Without DGs0.81372—494.6027302.2616————Summer0.907250.44344187.6912150.771571.941163.22781419.6, 1288.6, 2287.93,6,8Without DGs0.7819—668.9179410.0154————Autumn0.94170.3797777.198861.623968.543558.7164804.21,826.6,1369.43,6,8Without DGs0.87025—245.4149149.2695——


### 4.1.10. Case 3 (33-Bus RDN) for Spring

Figure11(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.91907 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.4185. Figure 11(b) shows the total P loss, reduced from 494.60 kW to 144.12 kW, while the percentage decrease of total P loss is 70.86%. Figure 11(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 302.26 kVAR to 115.72 kVAR, while the percentage decrease of total Q loss is 61.71%. Optimal DG sizes are found as 1161, 1177.6, and 1966.1 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 11
Case 3 (33-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.11. Case 3 (33-Bus RDN) for Summer

Figure12(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.90725 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.44344. Figure 12(b) shows the total P loss reduced from 668.91 kW to 187.69 kW, while the percentage decrease of total P loss is 71.94%. Figure 12(c) depicts the entire Q loss profile. For this case, the whole Q loss is reduced from 410.01 kVAR to 150.77 kVAR, while the total Q loss percentage decrease is 63.22%. Optimal DG sizes are found as 1419.6, 1288.6, and 2287.9 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 12
Case 3 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.12. Case 3 (33-Bus RDN) for Autumn

Figure13(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.9417 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.37977. Figure 13(b) shows the total P loss, which is reduced from 245.41 kW to 77.19 kW, while the percentage decrease of total P loss is 68.54%. Figure 13(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 149.26 kVAR to 61.62 kVAR, while the percentage decrease of total Q loss is 58.71%. Optimal DG sizes are found as 804.21, 826.6, and 1369.4 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 13
Case 3 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.13. Case 4 (33-Bus RDN) for Winter

Figure14(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.95257 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.37277. Figure 14(b) shows the total P loss, reduced from 159.61 kW to 53.74 kW, while the percentage decrease of total P loss is 66.32%. Figure 14(c) depicts the entire Q loss profile. For this case, the whole Q loss is reduced from 96.89 kVAR to 42.04 kVAR, while the total Q loss percentage decrease is 56.60%. Optimal DG sizes are 579.36, 587.25, 1102.9, and 174.16 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 14
Case 4 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 12
Results of case 4 (33-bus RDN with HHO).
Seasonvmin(per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.952570.3727753.748142.045966.326156.607579.36, 587.25, 1102.9, 174.163,6,8,4Without DGs0.89587—159.613696.8958————Spring0.919070.44263151.6158117.466169.345961.1376987.93, 1116.1, 1997.2, 290.953,6,8,4Without DGs0.81372—494.6027302.2616————Summer0.907250.47614199.115153.377470.233362.59231246, 1376.5, 2300.6, 188.483,6,8,4Without DGs0.7819—668.9179410.0154———Autumn0.94170.391379.703962.253567.522858.2945749.27, 753.49, 1389.8, 184.743,6,8,4Without DGs0.87025—245.4149149.2695———


### 4.1.14. Case 4 (33-Bus RDN) for Spring

Figure15(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.91907 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.44263. Figure 15(b) shows the total P loss, reduced from 494.60 kW to 151.61 kW, while the percentage decrease in P loss is 69.34%. Figure 15(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 302.26 kVAR to 117.46 kVAR, while the total Q loss percentage decrease is 61.1376%. Optimal DG sizes are 987.93, 1116.1, 1997.2, and 290.95 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 15
Case 4 (33-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.15. Case 4 (33-Bus RDN) for Summer

Figure16(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.90725 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.47614. Figure 16(b) shows the total P loss, reduced from 668.91 kW to 199.11 kW, while the percentage decrease of total P loss is 70.23%. Figure 16(c) illustrates the entire Q loss profile. For this case, the complete Q loss is reduced from 410.01 kVAR to 153.37 kVAR, while the total Q loss percentage decrease is 62.59%. Optimal DG sizes are 1246, 1376.5, 2300.6, and 188.48 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 16
Case 4 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.1.16. Case 4 (33-Bus RDN) for Autumn

Figure17(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.9417 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.3913. Figure 17(b) shows the total P loss, which is reduced from 245.41 kW to 79.70 kW, while the percentage decrease of total P loss is 67.52%. Figure 17(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 149.26 kVAR to 62.25 kVAR, while the total Q loss percentage decrease is 58.29%. Optimal DG sizes are 749.27, 753.49, 1389.8, and 184.74 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 17
Case 4 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.1. Case 1 (33-Bus RDN) for Winter

Figure2(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.90998 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.39695. Figure 2(b) shows the active power loss, reduced from 159.61 kW to 121.21 kW, while the decrease of total P loss is 24.06%. Figure 2(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 96.89 kVAR to 76.82 kVAR, while the decrease of total Q loss is 20.72%. The optimal DG size is found as 2578.5 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 2
Case 1 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 9
Results of case 1 (33-bus RDN with HHO).
Seasonvmin (per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.909980.39695121.210776.81724.059920.72212578.53Without DGs0.89587—159.613696.8958————Spring0.840890.52997368.0215234.536425.592522.40614888.13Without DGs0.81372—494.6027302.2616————Summer0.814840.59875495.3813316.116125.942922.90145805.53Without DGs0.7819—668.9179410.0154————Autumn0.88820.43087184.8168117.382324.692121.36223291.63Without DGs0.87025—245.4149149.2695————


## 4.1.2. Case 1 (33-Bus RDN) for Spring

Figure3(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.84089 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.52997. Figure 3(b) shows the total P loss, which is reduced from 494.60 kW to 368.02 kW, while the decrease of total P loss is 25.59%. Figure 3(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 302.26 kVAR to 234.53 kVAR, while the decrease of total Q loss is 22.40%. The optimal DG size is found as 4888.1 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 3
Case 1 (33-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.3. Case 1 (33-Bus RDN) for Summer

Figure4(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.81484 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.59875. Figure 4(b) shows the total P loss reduced from 668.91 kW to 495.38 kW, while the decrease of total P loss is 25.94%. Figure 4(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 410.01 kVAR to 316.11 kVAR, while the decrease of total Q loss is 22.90%. The optimal DG size is found as 5805.5 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 4
Case 1 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.4. Case 1 (33-Bus RDN) for Autumn

Figure5(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of one DG for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.8882 (per unit). The feasible bus selection is bus 3. The objective function for this case is obtained as 0.43087. Figure 5(b) shows the total P loss, which is reduced from 245.41 kW to 184.81 kW, while the decrease of total P loss is 24.69%. Figure 5(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 149.26 kVAR to 117.38 kVAR, while the decrease of total Q loss is 21.36%. The optimal DG size is found as 3291.6 kW. Table 9 presents the result summary for case 1 of 33-bus RDN for all four seasons’ data.Figure 5
Case 1 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.5. Case 2 (33-Bus RDN) for Winter

Figure6(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.9131 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.4016. Figure 6(b) shows the total P loss, which increased from 159.61 kW to 245.46 kW, while the percentage increase of total P loss is 53.78%. Figure 6(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 96.89 kVAR to 142.43 kVAR, while the percentage increase of total Q loss is 46.9967%. Optimal DG sizes are found as 347.21 kW and 989.34 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 6
Case 2 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 10
Results of case 2 (33-bus RDN with HHO).
Seasonvmin(per unit)FOP loss (kW)Q loss (kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’ssize (kW)DG’s location(bus no.)Winter0.91310.4016245.4671142.4336−53.7883−46.9967347.21, 989.343,6Without DGs0.89587—159.613696.8958————Spring0.847380.54649736.4943428.922−48.9062−41.9042715.63, 1844.43,6Without DGs0.81372—494.6027302.2616————Summer0.82290.6219980.4337571.8549−46.5701−39.4716873.16, 2181.43,6Without DGs0.7819—668.9179410.0154————Autumn0.892290.43844374.0595217.289−52.4192−45.5682457.26, 1254.13,6Without DGs0.87025—245.4149149.2695————


## 4.1.6. Case 2 (33-Bus RDN) for Spring

Figure7(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.84738 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.54649. Figure 7(b) shows the total P loss, which is increased from 494.60 kW to 736.49 kW, while the percentage increase of total P loss is 48.90%. Figure 7(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 302.26 kVAR to 428.9 kVAR, while the percentage increase of total Q loss is 41.90%. Optimal DG sizes are found as 715.63 and 1844.4 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 7
Case 2 (33-bus RDN) for spring. (a) Bus voltage profile. (b) Total P loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.7. Case 2 (33-Bus RDN) for Summer

Figure8(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.8229 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.6219. Figure 8(b) shows the total P loss, which increased from 668.91 kW to 980.43 kW, while the percentage increase of total P loss is 46.57%. Figure 8(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 410.01 kVAR to 571.85 kVAR, while the percentage increase of total Q loss is 39.47%. Optimal DG sizes are found as 873.16 and 2181.4 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 8
Case 2 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.8. Case 2 (33-Bus RDN) for Autumn

Figure9(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of two DGs for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.89229 (per unit). The feasible buses are buses 3 and 6. The objective function for this case is obtained as 0.43844. Figure 9(b) shows the total P loss, which increased from 245.41 kW to 374.05 kW, while the percentage increase of total P loss is 52.41%. Figure 9(c) shows the entire Q loss profile. For this case, the whole Q loss increases from 149.26 kVAR to 217.28 kVAR, while the percentage increase of total Q loss is 45.56%. Optimal DG sizes are found as 457.26 and 1254.1 kW. Table 10 presents the result summary for case 2 of 33-bus RDN for all four seasons’ data.Figure 9
Case 2 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.9. Case 3 (33-Bus RDN) for Winter

Figure10(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.95257 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.36559. Figure 10(b) shows the total P loss, reduced from 159.61 kW to 52.52 kW, while the percentage decrease in P loss is 67.09%. Figure 10(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 96.89 kVAR to 41.63 kVAR, while the percentage decrease of total Q loss is 57.03%. Optimal DG sizes are found as 728.58, 612.37, and 1092.6 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 10
Case 3 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 11
Results of case 3 (33-bus RDN with HHO).
Seasonvmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.952570.3655952.52841.631567.090557.0347728.58, 612.37, 1092.63,6,8Without DGs0.89587—159.613696.8958——Spring0.919070.4185144.1232115.724570.860861.71381161,1177.6, 1966.13,6,8Without DGs0.81372—494.6027302.2616————Summer0.907250.44344187.6912150.771571.941163.22781419.6, 1288.6, 2287.93,6,8Without DGs0.7819—668.9179410.0154————Autumn0.94170.3797777.198861.623968.543558.7164804.21,826.6,1369.43,6,8Without DGs0.87025—245.4149149.2695——


## 4.1.10. Case 3 (33-Bus RDN) for Spring

Figure11(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.91907 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.4185. Figure 11(b) shows the total P loss, reduced from 494.60 kW to 144.12 kW, while the percentage decrease of total P loss is 70.86%. Figure 11(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 302.26 kVAR to 115.72 kVAR, while the percentage decrease of total Q loss is 61.71%. Optimal DG sizes are found as 1161, 1177.6, and 1966.1 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 11
Case 3 (33-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.11. Case 3 (33-Bus RDN) for Summer

Figure12(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.90725 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.44344. Figure 12(b) shows the total P loss reduced from 668.91 kW to 187.69 kW, while the percentage decrease of total P loss is 71.94%. Figure 12(c) depicts the entire Q loss profile. For this case, the whole Q loss is reduced from 410.01 kVAR to 150.77 kVAR, while the total Q loss percentage decrease is 63.22%. Optimal DG sizes are found as 1419.6, 1288.6, and 2287.9 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 12
Case 3 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.12. Case 3 (33-Bus RDN) for Autumn

Figure13(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of three DGs for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.9417 (per unit). The feasible buses are buses 3, 6, and 8. The objective function for this case is obtained as 0.37977. Figure 13(b) shows the total P loss, which is reduced from 245.41 kW to 77.19 kW, while the percentage decrease of total P loss is 68.54%. Figure 13(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 149.26 kVAR to 61.62 kVAR, while the percentage decrease of total Q loss is 58.71%. Optimal DG sizes are found as 804.21, 826.6, and 1369.4 kW. Table 11 presents the result summary for case 3 of 33-bus RDN for all four seasons’ data.Figure 13
Case 3 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.13. Case 4 (33-Bus RDN) for Winter

Figure14(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for winter load. The bus voltage magnitude is upgraded from 0.89587 (per unit) to 0.95257 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.37277. Figure 14(b) shows the total P loss, reduced from 159.61 kW to 53.74 kW, while the percentage decrease of total P loss is 66.32%. Figure 14(c) depicts the entire Q loss profile. For this case, the whole Q loss is reduced from 96.89 kVAR to 42.04 kVAR, while the total Q loss percentage decrease is 56.60%. Optimal DG sizes are 579.36, 587.25, 1102.9, and 174.16 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 14
Case 4 (33-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 12
Results of case 4 (33-bus RDN with HHO).
Seasonvmin(per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.952570.3727753.748142.045966.326156.607579.36, 587.25, 1102.9, 174.163,6,8,4Without DGs0.89587—159.613696.8958————Spring0.919070.44263151.6158117.466169.345961.1376987.93, 1116.1, 1997.2, 290.953,6,8,4Without DGs0.81372—494.6027302.2616————Summer0.907250.47614199.115153.377470.233362.59231246, 1376.5, 2300.6, 188.483,6,8,4Without DGs0.7819—668.9179410.0154———Autumn0.94170.391379.703962.253567.522858.2945749.27, 753.49, 1389.8, 184.743,6,8,4Without DGs0.87025—245.4149149.2695———


## 4.1.14. Case 4 (33-Bus RDN) for Spring

Figure15(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for spring load. The bus voltage magnitude is upgraded from 0.81372 (per unit) to 0.91907 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.44263. Figure 15(b) shows the total P loss, reduced from 494.60 kW to 151.61 kW, while the percentage decrease in P loss is 69.34%. Figure 15(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 302.26 kVAR to 117.46 kVAR, while the total Q loss percentage decrease is 61.1376%. Optimal DG sizes are 987.93, 1116.1, 1997.2, and 290.95 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 15
Case 4 (33-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.15. Case 4 (33-Bus RDN) for Summer

Figure16(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for summer load. The bus voltage magnitude is upgraded from 0.7819 (per unit) to 0.90725 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.47614. Figure 16(b) shows the total P loss, reduced from 668.91 kW to 199.11 kW, while the percentage decrease of total P loss is 70.23%. Figure 16(c) illustrates the entire Q loss profile. For this case, the complete Q loss is reduced from 410.01 kVAR to 153.37 kVAR, while the total Q loss percentage decrease is 62.59%. Optimal DG sizes are 1246, 1376.5, 2300.6, and 188.48 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 16
Case 4 (33-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.1.16. Case 4 (33-Bus RDN) for Autumn

Figure17(a) illustrates the bus voltage magnitude of 33-bus RDN with and without considering the deployment of four DGs for the autumn load. The bus voltage magnitude is upgraded from 0.87025 (per unit) to 0.9417 (per unit). The feasible buses are buses 3, 6, 8, and 4. The objective function for this case is obtained as 0.3913. Figure 17(b) shows the total P loss, which is reduced from 245.41 kW to 79.70 kW, while the percentage decrease of total P loss is 67.52%. Figure 17(c) shows the entire Q loss profile. For this case, the whole Q loss is reduced from 149.26 kVAR to 62.25 kVAR, while the total Q loss percentage decrease is 58.29%. Optimal DG sizes are 749.27, 753.49, 1389.8, and 184.74 kW. Table 12 presents the result summary for case 4 of 33-bus RDN for all four seasons’ data.Figure 17
Case 4 (33-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2. IEEE 69-Bus RDN

Figure18 shows IEEE 69-bus RDN [80]. The total active load demand capacity of this IEEE 69-bus RDN is 3802 kW, while this system’s total reactive load demand is 2696 kVAR. The total P loss is 225.007 kW, and the bus voltage is 0.9091 (per unit) [69].Figure 18
IEEE 69-bus RDN.Table13 compares the literature simulation results with the proposed HHO based on four DGs (IEEE 33-bus RDN). It is observed that the maximum total P loss reduction of 72.75% is obtained as compared to the other literature results. The total Q loss reduction of 77.86% is obtained, which is higher than different algorithms results. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total active and Q loss reduction is 28.48 kW and 22.62 kVAR. In comparison, the bus voltage magnitude is upgraded to 0.97827 (per unit). The objective function for this case is 0.35001.Table 13
Results of literature with four DGs (69-bus RDN).
Algorithm Ref.Yearvmin(per unit)P loss (kW)PLOSSreduction (%)Q loss (kVAR)QLOSSreduction (%)DG sizes (kW)DG’s location(bus no.)Proposed HHO20210.9782728.481872.753722.6277.86329.88,229.42,112.03,924.1357,7,6,58PPA [69]20200.9782917.2392.3415——1014,47.22,145,1586.37,57,58,61GA [74]20200.97018.4360.4014.3259.30929.00,1075.0,984.80,1321.321,62,64,45PSO [74]2020—20.3163.0215.4554.301199.0,795.60,992.50,889.2061,63,17,32GA/PSO [74]2020—19.5663.9014.3260.20910.50,1192.0,884.90,978.6021,61,63,57HSA [65]20200.96734.2161.40—58.301302.0,369.00,101.80,978.5063,64,65,56BFOA [75]2020——66.56—61.50295.40,447.00,1345.1,1000.027,65,61,32ABC [65]20200.96032.0485.8019.9876.401000.0,200.00,338.20,1100.061,51,62,32SSA [40]20190.97023.4069.10—62.40380.00,527.00,1718.0,467.0017,10,60,32WCA [39]20180.97022.4368.20—61.20775.00,1105.0,438.00,890.0061,62,23,21HGWO [79]20130.96020.4169.14—64.30527.00,380.00,1718.0,398.5011,17,61,45Table14 compares the literature simulation results with the proposed HHO based on three DGs (IEEE 33-bus RDN). It is observed that the maximum total P loss reduction of 68.55% is obtained compared to the other literature results. The total Q loss reduction of 69.86% is obtained, which is higher than different algorithms results. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total active and Q loss reduction is 32.87 kW and 30.79 kVAR. At the same time, the bus voltage magnitude is upgraded to 0.9710 (per unit). The objective function for this case is 0.35373.Table 14
Results of literature with three DGs (69-bus RDN).
Algorithm Ref.Yearvmin (per unit)P loss (kW)PLOSS reduction (%)Q loss (kVAR)QLOSS reduction (%)DG sizes (kW)DG’s location (bus no.)Proposed HHO20210.9710032.871968.55430.788369.861764.3,617.55,399.6657,7,6PPA [69]20200.98240.5082——27.919,1108.3,1558.351,61,62GA [74]20200.970019.87560.407.11252.2929.70,1075.0,984.8021,62,64PSO [74]2020—19.43063.0212.32051.21199.0,795.60,992.5061,63,17GA/PSO [74]2020—28.72063.9013.57053.0910.50,1192.0,884.9021,61,63HSA [65]20200.967035.43061.40—51.01302.4,369.00,101.8063,64,65BFOA [75]20200.9608—66.56—56.5295.40,447.60,1345.027,65,61ABC [65]20200.960032.05686.8019.98077.31000.0,200.00,338.2061,51,62QOCSOS [42]2020—69.4369.14——526.8,380.4,1719.011,18,61CSCA [14]20200.9870.2068.8——365.9,1675.85,65.5217,61,67SSA [40]20190.978969.5269.1——380,527,171817,10,60FWA [41]20190.97477.8765.39——480.5,1198.6,225.865,61,27SFLA [41]20190.975277.78465.43——1088.7,167.3,980.957,63,26SSA [40]20190.978927.54069.10—59.2380.00,527.00,1718.017,10,60WCA [39]20180.98771.5568.2——775,1105,43861,62,23WCA [39]20180.967021.45068.20—51.1775.00,1105.0,438.0061,62,23HGWO [37]20170.9869.4369.14——527,380,171811,17,61BFOA [34]20140.980875.2466.56——295.4,447.6,1345.127,65,61HSA [33]20130.96786.8561.4——1302.4,369,101.863,64,65HGWO [79]20130.960024.23069.14—59.1527.00,380.00,1718.011,17,61Table15 compares the literature simulation results with the proposed HHO based on two DGs (IEEE 33-bus RDN). It is observed that the maximum total P loss reduction of 67.16% is obtained as compared to the other literature results. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total P loss reduction is 34.32 kW, while the bus voltage magnitude is upgraded to 0.9710 (per unit). The objective function for this case is 0.35361.Table 15
Comparison of literature results with two DGs (case 2 (69-bus RDN)).
Algorithm Ref.Yearvmin (per unit)P loss (kW)PLOSS (%)DG’s size (kW)DG’s locationFOALO [24]20170.9801070.7768.5470538.70,1700.017,610.83200ALO [26]20180.9789071.6868.1400531.48,1781.517,610.83000EMA [29]20180.9794071.8668.06001886.9,649.3061,690.83000CSFS-III [27]2018—71.6868.1400531.00,1781.017,610.34000PPA [69]20200.9710569.1169.24001040.0,1856.07,570.83500Proposed HHO20210.9710034.324567.16441232.0,414.9957,70.35361Table16 compares the literature simulation results with the proposed HHO based on one DG (IEEE 33-bus RDN). It is observed that the maximum total P loss reduction of 65.56% is obtained compared to the other literature results. Moreover, the bus voltage deviation is also improved with the applied HHO algorithm. The total P loss reduction is 36.01 kW. In comparison, the bus voltage magnitude is upgraded to 0.9710 (per unit). The objective function for this case is 0.35442.Table 16
Comparison of literature results with one DG (case 1 (69-bus RDN)).
Algorithm Ref.Yearvmin (per unit)P loss (kW)PLOSS (%)DG’s size (kW)DG’s locationFOALO [24]20170.9679081.801063.64501800.0610.80200ALO [26]20180.9682083.250063.00001872.7610.79900EMA [29]20180.9689082.870063.17001910.3570.80000GA [1]20190.9682063.000072.00001850.0600.84400CSFS-III [27]2018—83.230063.01001873.0610.31500CDE [30]20190.9683083.250063.00001872.4610.79900PPA [69]20200.9710574.250067.00002088.0570.82900Proposed HHO20210.9710036.006965.55511330.2570.35442


### 4.2.1. Case 1 (69-Bus RDN) for Winter

Figure19(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for winter load. The bus voltage magnitude is upgraded from 0.93825 (per unit) to 0.97991 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.36388. Figure 19(b) shows the total P loss reduced from 104.53 kW to 35.74 kW while the percentage decrease of total P loss is 65.80%. Figure 19(c) depicts the entire Q loss profile. For this case, the whole Q loss is reduced from 47.61 kVAR to 15.93 kVAR, while the total Q loss percentage decrease is 66.53%. The optimal DG size is found as 1346.9 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 19
Case 1 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 17
Results with one DG (case 1 (69-bus RDN)) with proposed HHO.
Seasonvmin (per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.979910.3638835.748615.935565.802266.53081346.957Without DGs0.93825104.534747.6124———Spring0.966560.4103897.976443.653568.00368.55192357.257Without DGs0.89392306.2051138.8112———Summer0.962040.43178126.690856.419468.585469.08752710.857Without DGs0.87809403.2863182.5134———Autumn0.97550.3765152.724823.505166.670967.32841675.657Without DGs0.92395158.194671.9434———


### 4.2.2. Case 1 (69-Bus RDN) for Spring

Figure20(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.96656 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.41038. Figure 20(b) shows the total P loss reduced from 306.20 kW to 97.97 kW while the percentage decrease of total P loss is 68.00%. Figure 20(c) illustrates the entire Q loss profile. For this case, the complete Q loss is reduced from 138.81 kVAR to 43.65 kVAR, while the total Q loss percentage decrease is 68.55%. The optimal DG size is found as 2357.2 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 20
Case 1 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.3. Case 1 (69-Bus RDN) for Summer

Figure21(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.96204 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.43178. Figure 21(b) shows the total P loss reduced from 403.28 kW to 126.69 kW while the percentage decrease of total P loss is 68.58%. Figure 21(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 182.51 kVAR to 56.41 kVAR, while the total Q loss percentage decrease is 69.08%. The optimal DG size is found as 2710.8 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 21
Case 1 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.4. Case 1 (69-Bus RDN) for Autumn

Figure22(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.9755 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.37651. Figure 22(b) shows the total P loss reduced from 158.1946 kW to 52.7248 kW, while the percentage decrease of total P loss is 66.67%. Figure 22(c) depicts the entire Q loss profile. For this case, the total Q loss is reduced from 71.94 kVAR to 23.50 kVAR, while the total Q loss percentage decrease is 67.32%. The optimal DG size is found as 1675.6 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 22
Case 1 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.5. Case 2 (69-Bus RDN) for Winter

Figure23(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for winter load. The bus voltage magnitude is upgraded from 0.93825 (per unit) to 0.97991 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.37332. Figure 23(b) shows the total P loss, which increased from 104.5347 kW to 346.82 kW, while the percentage increase of total P loss is 231.77%. Figure 23(c) depicts the entire Q loss profile. The total Q loss for this case is increased from 47.61 kVAR to 155.20 kVAR, while the percentage increase of total Q loss is 225.97%. Optimal DG sizes are found as 868.32 and 419.28 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 23
Case 2 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 18
Results with two DGs (case 2 (69-bus RDN)) with proposed HHO.
Seasonvmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)DG’s size (kW)DG’s location (bus no.)Winter0.979910.37332346.8215155.2039−231.7765−225.9733868.32, 419.2857, 7Without DGs0.93825104.534747.6124———Spring0.966560.44221122.0975504.2769−266.4529−263.28261484.7, 788.6657, 7Without DGs0.89392306.2051138.8112———Summer0.962040.47521542.1811694.2981−282.4035−280.40941699.7, 923.4257, 7Without DGs0.87809403.2863182.5134———Autumn0.97550.39162540.5331242.1968−241.6888−236.64921068.5, 540.7757, 7Without DGs0.92395—158.194671.9434————


### 4.2.6. Case 2 (69-Bus RDN) for Spring

Figure24(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.96656 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.4422. Figure 24(b) shows the total P loss, which is increased from 306.2051 kW to 1122.09 kW, while the percentage increase of total P loss is 266.45%. Figure 24(c) shows the entire Q loss profile. The total Q loss for this case is increased from 138.81 kVAR to 504.27 kVAR, while the percentage increase of total Q loss is 263.28%. Optimal DG sizes are found as 1484.7 and 788.66 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 24
Case 2 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.7. Case 2 (69-Bus RDN) for Summer

Figure25(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.96204 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.4752. Figure 25(b) shows the total P loss, which is increased from 403.28 kW to 1542.18 kW, while the percentage increase of total P loss is 282.40%. Figure 25(c) shows the entire Q loss profile. The total Q loss for this case is increased from 182.51 kVAR to 694.29 kVAR, while the percentage increase of total Q loss is 280.40%. Optimal DG sizes are found as 1699.7 and 923.42 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 25
Case 2 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.8. Case 2 (69-Bus RDN) for Autumn

Figure26(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.9755 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.39162. Figure 26(b) shows the total P loss, which increased from 158.19 kW to 540.53 kW, while the percentage increase of total P loss is 241.68%. Figure 26(c) depicts the entire Q loss profile. The total Q loss for this case is increased from 71.94 kVAR to 242.19 kVAR, while the percentage increase of total Q loss is 236.64%. Optimal DG sizes are found as 1068.5 and 540.77 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 26
Case 2 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.9. Case 3 (69-Bus RDN) for Winter

Figure27(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for winter load. The bus voltage magnitude is upgraded from 0.9091 (per unit) to 0.97991 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.36328. Figure 27(b) shows the total P loss, reduced from 225.00 kW to 33.27 kW, while the percentage decrease of total P loss is 33.27%. Figure 27(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 33.27 kVAR to 14.74 kVAR, while the total Q loss percentage decrease is 14.74%. Optimal DG sizes are found as 1212.5, 552.87, and 19.166 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 27
Case 3 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 19
Results with three DGs (case 3 (69-bus RDN)) with proposed HHO.
Seasonvmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.979910.3632833.276314.749768.167269.02141212.5, 552.87, 19.16657, 7, 6Without DGs0.9091225.00747.6124———Spring0.966560.379992.490140.904469.794770.53242108.7, 706.76, 247.2857, 7, 6Without DGs0.89392306.2051138.8112———Summer0.962040.39179117.109651.634270.961271.70942361.5, 1059.3, 128.5557, 7, 6Without DGs0.87809403.2863182.5134———Autumn0.97550.3607549.354221.853568.801669.6241483.8, 674.15, 23.57457, 7, 6Without DGs0.92395158.194671.9434———


### 4.2.10. Case 3 (69-Bus RDN) for Spring

Figure28(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.96656 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.3799. Figure 28(b) shows the total P loss reduced from 306.20 kW to 92.49 kW while the percentage decrease of total P loss is 69.79%. Figure 28(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 138.81 kVAR to 40.90 kVAR, while the percentage decrease in Q loss is 70.53%. Optimal DG sizes are found as 2108.7, 706.76, and 247.28 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 28
Case 3 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.11. Case 3 (69-Bus RDN) for Summer

Figure29(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.96204 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.39179. Figure 29(b) shows the total P loss, reduced from 403.28 kW to 117.10 kW, while the percentage decrease of total P loss is 70.96%. Figure 29(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 70.96 kVAR to 51.63 kVAR, while the percentage decrease is 71.70%. Optimal DG sizes are found as 2361.5, 1059.3, and 128.55 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 29
Case 3 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.12. Case 3 (69-Bus RDN) for Autumn

Figure30(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.9755 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.36075. Figure 30(b) shows the total P loss reduced from 158.19 kW to 49.35 kW while the percentage decrease of total P loss is 68.80%. Figure 30(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 71.94 kVAR to 21.85 kVAR, while the total Q loss percentage decrease is 69.62%. Optimal DG sizes are found as 1483.8, 674.15, and 23.574 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 30
Case 3 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.13. Case 4 (69-Bus RDN) for Winter

Figure31(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for winter load. The bus voltage magnitude is upgraded from 0.93825 (per unit) to 0.9849 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.36025. Figure 31(b) shows the total P loss reduced from 104.53 kW to 22.76 kW while the percentage decrease of total P loss is 78.22%. Figure 31(c) depicts the entire Q loss profile. For this case, the total Q loss is reduced from 47.61 kVAR to 11.21 kVAR, while the total Q loss percentage decrease is 76.45%. Optimal DG sizes are 40.66, 553.91, 19.998, and 1156.8 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 31
Case 4 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 20
Results with four DGs (case 4 (69-bus RDN)) with proposed HHO.
Seasonvmin(per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no).Winter0.98490.3602522.762611.212578.224876.450540.66, 553.91, 19.998, 1156.857, 7, 6, 58Without DGs0.93825104.534747.6124———Spring0.974960.3992161.156230.042480.027778.3574108.65, 1032.2, 31.762, 1957.557, 7, 6, 58Without DGs0.89392306.2051138.8112———Summer0.971610.41778.855838.683680.446778.805123.63, 1198.9, 41.268, 2240.557, 7, 6, 58Without DGs0.87809403.2863182.5134———Autumn0.981610.3708933.246816.367978.983677.248954.404, 700.2, 39.572, 1427.157, 7, 6, 58Without DGs0.92395158.194671.9434———


### 4.2.14. Case 4 (69-Bus RDN) for Spring

Figure32(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.97496 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.39921. Figure 32(b) shows the total P loss, reduced from 306.20 kW to 61.15 kW, while the percentage decrease of total P loss is 80.02%. Figure 32(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 138.81 kVAR to 30.04 kVAR, while the percentage decrease in Q loss is 78.35%. Optimal DG sizes are 108.65, 1032.2, 31.762, and 1957.5 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 32
Case 4 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.15. Case 4 (69-Bus RDN) for Summer

Figure33(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.97161 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.417. Figure 33(b) shows the total P loss, which is reduced from 403.28 kW to 78.85 kW, while the percentage decrease of total P loss is 80.44%. Figure 33(c) shows the total Q loss profile. For this case, the total Q loss is reduced from 182.51 kVAR to 38.68 kVAR, while the percentage decrease in Q loss is 78.80%. Optimal DG sizes are 123.63, 1198.9, 41.268, and 2240.5 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 33
Case 4 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


### 4.2.16. Case 4 (69-Bus RDN) for Autumn

Figure34(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.98161 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.37089. Figure 34(b) shows the total P loss, reduced from 158.19 kW to 33.24 kW, while the percentage decrease of total P loss is 33.24%. Figure 34(c) shows the total Q loss profile. For this case, the total Q loss is reduced from 71.94 kVAR to 16.36 kVAR, while the percentage decrease in Q loss is 77.24%. Optimal DG sizes are 54.404, 700.2, 39.572, and 1427.1 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 34
Case 4 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.1. Case 1 (69-Bus RDN) for Winter

Figure19(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for winter load. The bus voltage magnitude is upgraded from 0.93825 (per unit) to 0.97991 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.36388. Figure 19(b) shows the total P loss reduced from 104.53 kW to 35.74 kW while the percentage decrease of total P loss is 65.80%. Figure 19(c) depicts the entire Q loss profile. For this case, the whole Q loss is reduced from 47.61 kVAR to 15.93 kVAR, while the total Q loss percentage decrease is 66.53%. The optimal DG size is found as 1346.9 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 19
Case 1 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 17
Results with one DG (case 1 (69-bus RDN)) with proposed HHO.
Seasonvmin (per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.979910.3638835.748615.935565.802266.53081346.957Without DGs0.93825104.534747.6124———Spring0.966560.4103897.976443.653568.00368.55192357.257Without DGs0.89392306.2051138.8112———Summer0.962040.43178126.690856.419468.585469.08752710.857Without DGs0.87809403.2863182.5134———Autumn0.97550.3765152.724823.505166.670967.32841675.657Without DGs0.92395158.194671.9434———


## 4.2.2. Case 1 (69-Bus RDN) for Spring

Figure20(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.96656 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.41038. Figure 20(b) shows the total P loss reduced from 306.20 kW to 97.97 kW while the percentage decrease of total P loss is 68.00%. Figure 20(c) illustrates the entire Q loss profile. For this case, the complete Q loss is reduced from 138.81 kVAR to 43.65 kVAR, while the total Q loss percentage decrease is 68.55%. The optimal DG size is found as 2357.2 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 20
Case 1 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.3. Case 1 (69-Bus RDN) for Summer

Figure21(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.96204 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.43178. Figure 21(b) shows the total P loss reduced from 403.28 kW to 126.69 kW while the percentage decrease of total P loss is 68.58%. Figure 21(c) illustrates the entire Q loss profile. For this case, the whole Q loss is reduced from 182.51 kVAR to 56.41 kVAR, while the total Q loss percentage decrease is 69.08%. The optimal DG size is found as 2710.8 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 21
Case 1 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.4. Case 1 (69-Bus RDN) for Autumn

Figure22(a) illustrates the bus voltage of 69-bus RDN considering the deployment of one DG for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.9755 (per unit). The feasible bus selection is bus 57. The objective function for this case is obtained as 0.37651. Figure 22(b) shows the total P loss reduced from 158.1946 kW to 52.7248 kW, while the percentage decrease of total P loss is 66.67%. Figure 22(c) depicts the entire Q loss profile. For this case, the total Q loss is reduced from 71.94 kVAR to 23.50 kVAR, while the total Q loss percentage decrease is 67.32%. The optimal DG size is found as 1675.6 kW. Table 17 presents the result summary for case 1 of IEEE 69-bus RDN for all four seasons’ data.Figure 22
Case 1 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.5. Case 2 (69-Bus RDN) for Winter

Figure23(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for winter load. The bus voltage magnitude is upgraded from 0.93825 (per unit) to 0.97991 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.37332. Figure 23(b) shows the total P loss, which increased from 104.5347 kW to 346.82 kW, while the percentage increase of total P loss is 231.77%. Figure 23(c) depicts the entire Q loss profile. The total Q loss for this case is increased from 47.61 kVAR to 155.20 kVAR, while the percentage increase of total Q loss is 225.97%. Optimal DG sizes are found as 868.32 and 419.28 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 23
Case 2 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 18
Results with two DGs (case 2 (69-bus RDN)) with proposed HHO.
Seasonvmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)DG’s size (kW)DG’s location (bus no.)Winter0.979910.37332346.8215155.2039−231.7765−225.9733868.32, 419.2857, 7Without DGs0.93825104.534747.6124———Spring0.966560.44221122.0975504.2769−266.4529−263.28261484.7, 788.6657, 7Without DGs0.89392306.2051138.8112———Summer0.962040.47521542.1811694.2981−282.4035−280.40941699.7, 923.4257, 7Without DGs0.87809403.2863182.5134———Autumn0.97550.39162540.5331242.1968−241.6888−236.64921068.5, 540.7757, 7Without DGs0.92395—158.194671.9434————


## 4.2.6. Case 2 (69-Bus RDN) for Spring

Figure24(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.96656 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.4422. Figure 24(b) shows the total P loss, which is increased from 306.2051 kW to 1122.09 kW, while the percentage increase of total P loss is 266.45%. Figure 24(c) shows the entire Q loss profile. The total Q loss for this case is increased from 138.81 kVAR to 504.27 kVAR, while the percentage increase of total Q loss is 263.28%. Optimal DG sizes are found as 1484.7 and 788.66 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 24
Case 2 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.7. Case 2 (69-Bus RDN) for Summer

Figure25(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.96204 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.4752. Figure 25(b) shows the total P loss, which is increased from 403.28 kW to 1542.18 kW, while the percentage increase of total P loss is 282.40%. Figure 25(c) shows the entire Q loss profile. The total Q loss for this case is increased from 182.51 kVAR to 694.29 kVAR, while the percentage increase of total Q loss is 280.40%. Optimal DG sizes are found as 1699.7 and 923.42 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 25
Case 2 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.8. Case 2 (69-Bus RDN) for Autumn

Figure26(a) illustrates the bus voltage of 69-bus RDN considering the deployment of two DGs for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.9755 (per unit). The feasible buses are 57 and 7. The objective function for this case is obtained as 0.39162. Figure 26(b) shows the total P loss, which increased from 158.19 kW to 540.53 kW, while the percentage increase of total P loss is 241.68%. Figure 26(c) depicts the entire Q loss profile. The total Q loss for this case is increased from 71.94 kVAR to 242.19 kVAR, while the percentage increase of total Q loss is 236.64%. Optimal DG sizes are found as 1068.5 and 540.77 kW. Table 18 presents case 2 of IEEE 69-bus RDN for all four seasons’ data.Figure 26
Case 2 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.9. Case 3 (69-Bus RDN) for Winter

Figure27(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for winter load. The bus voltage magnitude is upgraded from 0.9091 (per unit) to 0.97991 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.36328. Figure 27(b) shows the total P loss, reduced from 225.00 kW to 33.27 kW, while the percentage decrease of total P loss is 33.27%. Figure 27(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 33.27 kVAR to 14.74 kVAR, while the total Q loss percentage decrease is 14.74%. Optimal DG sizes are found as 1212.5, 552.87, and 19.166 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 27
Case 3 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 19
Results with three DGs (case 3 (69-bus RDN)) with proposed HHO.
Seasonvmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no.)Winter0.979910.3632833.276314.749768.167269.02141212.5, 552.87, 19.16657, 7, 6Without DGs0.9091225.00747.6124———Spring0.966560.379992.490140.904469.794770.53242108.7, 706.76, 247.2857, 7, 6Without DGs0.89392306.2051138.8112———Summer0.962040.39179117.109651.634270.961271.70942361.5, 1059.3, 128.5557, 7, 6Without DGs0.87809403.2863182.5134———Autumn0.97550.3607549.354221.853568.801669.6241483.8, 674.15, 23.57457, 7, 6Without DGs0.92395158.194671.9434———


## 4.2.10. Case 3 (69-Bus RDN) for Spring

Figure28(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.96656 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.3799. Figure 28(b) shows the total P loss reduced from 306.20 kW to 92.49 kW while the percentage decrease of total P loss is 69.79%. Figure 28(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 138.81 kVAR to 40.90 kVAR, while the percentage decrease in Q loss is 70.53%. Optimal DG sizes are found as 2108.7, 706.76, and 247.28 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 28
Case 3 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.11. Case 3 (69-Bus RDN) for Summer

Figure29(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.96204 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.39179. Figure 29(b) shows the total P loss, reduced from 403.28 kW to 117.10 kW, while the percentage decrease of total P loss is 70.96%. Figure 29(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 70.96 kVAR to 51.63 kVAR, while the percentage decrease is 71.70%. Optimal DG sizes are found as 2361.5, 1059.3, and 128.55 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 29
Case 3 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.12. Case 3 (69-Bus RDN) for Autumn

Figure30(a) illustrates the bus voltage of 69-bus RDN considering the deployment of three DGs for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.9755 (per unit). The feasible buses are 57, 7, and 6. The objective function for this case is obtained as 0.36075. Figure 30(b) shows the total P loss reduced from 158.19 kW to 49.35 kW while the percentage decrease of total P loss is 68.80%. Figure 30(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 71.94 kVAR to 21.85 kVAR, while the total Q loss percentage decrease is 69.62%. Optimal DG sizes are found as 1483.8, 674.15, and 23.574 kW. Table 19 presents case 3 of IEEE 69-bus RDN for all four seasons’ data.Figure 30
Case 3 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.13. Case 4 (69-Bus RDN) for Winter

Figure31(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for winter load. The bus voltage magnitude is upgraded from 0.93825 (per unit) to 0.9849 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.36025. Figure 31(b) shows the total P loss reduced from 104.53 kW to 22.76 kW while the percentage decrease of total P loss is 78.22%. Figure 31(c) depicts the entire Q loss profile. For this case, the total Q loss is reduced from 47.61 kVAR to 11.21 kVAR, while the total Q loss percentage decrease is 76.45%. Optimal DG sizes are 40.66, 553.91, 19.998, and 1156.8 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 31
Case 4 (69-bus RDN) for winter. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)Table 20
Results with four DGs (case 4 (69-bus RDN)) with proposed HHO.
Seasonvmin(per unit)FOP loss (kW)Q loss(kVAR)PLOSSreduction (%)QLOSSreduction (%)DG’s size (kW)DG’s location(bus no).Winter0.98490.3602522.762611.212578.224876.450540.66, 553.91, 19.998, 1156.857, 7, 6, 58Without DGs0.93825104.534747.6124———Spring0.974960.3992161.156230.042480.027778.3574108.65, 1032.2, 31.762, 1957.557, 7, 6, 58Without DGs0.89392306.2051138.8112———Summer0.971610.41778.855838.683680.446778.805123.63, 1198.9, 41.268, 2240.557, 7, 6, 58Without DGs0.87809403.2863182.5134———Autumn0.981610.3708933.246816.367978.983677.248954.404, 700.2, 39.572, 1427.157, 7, 6, 58Without DGs0.92395158.194671.9434———


## 4.2.14. Case 4 (69-Bus RDN) for Spring

Figure32(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for spring load. The bus voltage magnitude is upgraded from 0.89392 (per unit) to 0.97496 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.39921. Figure 32(b) shows the total P loss, reduced from 306.20 kW to 61.15 kW, while the percentage decrease of total P loss is 80.02%. Figure 32(c) illustrates the entire Q loss profile. For this case, the total Q loss is reduced from 138.81 kVAR to 30.04 kVAR, while the percentage decrease in Q loss is 78.35%. Optimal DG sizes are 108.65, 1032.2, 31.762, and 1957.5 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 32
Case 4 (69-bus RDN) for spring. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.15. Case 4 (69-Bus RDN) for Summer

Figure33(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for summer load. The bus voltage magnitude is upgraded from 0.87809 (per unit) to 0.97161 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.417. Figure 33(b) shows the total P loss, which is reduced from 403.28 kW to 78.85 kW, while the percentage decrease of total P loss is 80.44%. Figure 33(c) shows the total Q loss profile. For this case, the total Q loss is reduced from 182.51 kVAR to 38.68 kVAR, while the percentage decrease in Q loss is 78.80%. Optimal DG sizes are 123.63, 1198.9, 41.268, and 2240.5 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 33
Case 4 (69-bus RDN) for summer. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.2.16. Case 4 (69-Bus RDN) for Autumn

Figure34(a) illustrates the bus voltage of 69-bus RDN considering the deployment of four DGs for the autumn load. The bus voltage magnitude is upgraded from 0.92395 (per unit) to 0.98161 (per unit). The feasible buses are 57, 7, 6, and 58. The value of the objective function for this case is 0.37089. Figure 34(b) shows the total P loss, reduced from 158.19 kW to 33.24 kW, while the percentage decrease of total P loss is 33.24%. Figure 34(c) shows the total Q loss profile. For this case, the total Q loss is reduced from 71.94 kVAR to 16.36 kVAR, while the percentage decrease in Q loss is 77.24%. Optimal DG sizes are 54.404, 700.2, 39.572, and 1427.1 kW. Table 20 presents the result summary for case 4 of IEEE 69-bus RDN for all four seasons’ data.Figure 34
Case 4 (69-bus RDN) for autumn. (a) Bus voltage profile. (b) TotalP loss. (c) Total Q loss.
(a)(b)(c)


## 4.3. Comparison of Cases

Figures35(a) and 35(b) show the bus voltage and total P loss profiles of 33-bus RDN with the deployment of one, two, three, and four DGs for all four cases, respectively.Figure 35
Comparison of (a) bus voltage profiles and (b) totalP loss profiles in IEEE 33-bus RDN.
(a)(b)For winter load, the bus voltage drop is minimized in case 3 of 33-bus RDN, which is 0.95257 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 3 of 33-bus RDN is 52.52 kW, while the percentage decrease of total P loss, in this case, is 67.09% which shows that deployment of three DGs in a 33-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 3 of 33-bus RDN is 41.63 kVAR, while the percentage decrease of total Q loss is 57.03% which shows that deployment of three DGs in a 33-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of 33-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 245.46 kW and 142.43 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 53.78% and 46.99%, respectively, showing that deploying two DGs in 33-bus RDS is not feasible to minimize total P and Q losses with improved bus voltage. Table 21 presents the result summary for all four cases of 33-bus RDN with seasonal load data of winter.Table 21
Result comparison of 33-bus RDN for winter.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.909980.39695121.210776.81724.059920.7221Case 20.91310.4016245.4671142.4336−53.7883−46.9967Case 30.952570.3655952.52841.631567.090557.0347Case 40.952570.3727753.748142.045966.326156.607For spring load, the bus voltage drop is minimized in case 3 of 33-bus RDN, 0.91907 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 3 of 33-bus RDN is 144.12 kW, while the percentage decrease of total P loss, in this case, is 144.12% which shows that deployment of three DGs in a 33-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 3 of 33-bus RDN is 144.12 kVAR, while the percentage decrease of total Q loss is 144.12% which shows that deploying three DGs in a 33-bus system is still more feasible, minimizing total Q loss with improved bus voltage. However, case 2 of 33-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 736.49 kW and 428.92 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 48.90% and 41.90%, respectively, which shows that deploying two DGs in 33-bus RDS is not feasible in minimizing total P and Q losses with improved bus voltage. Table 22 presents the result summary for all four cases of 33-bus RDN with seasonal load data of spring.Table 22
Result comparison of 33-bus RDN for spring.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.840890.52997368.0215234.536425.592522.4061Case 20.847380.54649736.4943428.922−48.9062−41.9042Case 30.919070.4185144.1232115.724570.860861.7138Case 40.919070.44263151.6158117.466169.345961.1376For summer load, the bus voltage drop is minimized in case 3 of 33-bus RDN, which is 0.90725 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 3 of 33-bus RDN is 187.69 kW, while the percentage decrease of total P loss, in this case, is 187.69% which shows that deployment of three DGs in a 33-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 3 of 33-bus RDN is 187.69 kVAR, while the percentage decrease of total Q loss is 63.22% which shows that deployment of three DGs in a 33-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of 33-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 63.22 kW and 571.85 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 46.57% and 39.47%, respectively, which shows that deploying two DGs in 33-bus RDS is not feasible in minimizing total P and Q losses with improved bus voltage. Table 23 presents the result summary for all four cases of 33-bus RDN with seasonal load data of summer.Table 23
Result comparison of 33-bus RDN for summer.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.814840.59875495.3813316.116125.942922.9014Case 20.82290.6219980.4337571.8549−46.5701−39.4716Case 30.907250.44344187.6912150.771571.941163.2278Case 40.907250.47614199.115153.377470.233362.5923For autumn load, the bus voltage drop is minimized in case 3 of 33-bus RDN, which is 0.9417 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 3 of 33-bus RDN is 77.19 kW, while the percentage decrease of total P loss, in this case, is 68.54% which shows that deployment of three DGs in a 33-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 3 of 33-bus RDN is 61.62 kVAR, while the percentage decrease of total Q loss is 58.71% which shows that deployment of three DGs in a 33-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of 33-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 374.05 kW and 217.28 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 52.41% and 45.56%, respectively, which shows that deploying two DGs in 33-bus RDS is not feasible in minimizing total P and Q losses with improved bus voltage. Table 24 presents the result summary for all four cases of 33-bus RDN with seasonal load data of autumn.Table 24
Result comparison of 33-bus RDN for autumn.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.88820.43087184.8168117.382324.692121.3622Case 20.892290.43844374.0595217.289−52.4192−45.5682Case 30.94170.3797777.198861.623968.543558.7164Case 40.94170.391379.703962.253567.522858.2945Figures36(a) and 36(b) show the bus voltage and total P loss profiles of 33-bus RDN with the deployment of one, two, three, and four DGs for all four cases, respectively.Figure 36
Comparison of (a) bus voltage profiles and (b) totalP loss profiles in IEEE 69-bus RDN.
(a)(b)For winter load, the bus voltage drop is minimized in case 4 of IEEE 69-bus RDN, which is 0.9849 (per unit), which shows more improved value than other cases. The total P loss reduction for case 4 of IEEE 69-bus RDN is 22.76 kW, while the percentage decrease of totalP loss, in this case, is 78.22% which shows that deployment of three DGs in a 69-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 4 of IEEE 69-bus RDN is 11.21 kVAR, while the percentage decrease of total Q loss is 76.45% which shows that deployment of three DGs in a 69-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of IEEE 69-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 346.82 kW and 155.20 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 231.77% and 225.97%, respectively, showing that deploying two DGs in 69-bus RDS is not feasible to minimize total P and Q losses with improved bus voltage. Table 25 presents the result summary for all four cases of IEEE 69-bus RDN with seasonal load data of winter.Table 25
Result comparison of IEEE 69-bus RDN for winter.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.979910.3638835.748615.935565.802266.5308Case 20.979910.37332346.8215155.2039−231.7765−225.9733Case 30.979910.3632833.276314.749768.167269.0214Case 40.98490.3602522.762611.212578.224876.4505For spring load, the bus voltage drop is minimized in case 4 of IEEE 69-bus RDN, 0.97496 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 4 of IEEE 69-bus RDN is 61.15 kW, while the percentage decrease of total P loss, in this case, is 80.02% which shows that deployment of three DGs in a 69-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 4 of IEEE 69-bus RDN is 30.04 kVAR, while the percentage decrease of total Q loss is 78.35% which shows that deployment of three DGs in a 69-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of IEEE 69-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 1122.09 kW and 504.27 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 266.45% and 263.28%, respectively, showing that deploying two DGs in 69-bus RDS is not feasible in minimizing total P and Q losses with improved bus voltage. Table 26 presents the result summary for all four cases of IEEE 69-bus RDN with seasonal load data of spring.Table 26
Result comparison of IEEE 69-bus RDN for spring.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.966560.4103897.976443.653568.00368.5519Case 20.966560.44221122.0975504.2769−266.4529−263.2826Case 30.966560.379992.490140.904469.794770.5324Case 40.974960.3992161.156230.042480.027778.3574For summer load, the bus voltage drop is minimized in case 4 of IEEE 69-bus RDN, which is 0.97161 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 4 of IEEE 69-bus RDN is 78.85 kW, while the percentage decrease of total P loss, in this case, is 80.44% which shows that deployment of three DGs in a 69-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 4 of IEEE 69-bus RDN is 38.68 kVAR, while the percentage decrease of total Q loss is 78.80% which shows that deployment of three DGs in a 69-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of IEEE 69-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 1542.18 kW and 694.29 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 282.40% and 280.40%, respectively, showing that deploying two DGs in 69-bus RDS is not feasible to minimize total P and Q losses with improved bus voltage. Table 27 presents the result summary for all four cases of IEEE 69-bus RDN with seasonal load data of summer.Table 27
Result comparison of IEEE 69-bus RDN for summer.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.962040.43178126.690856.419468.585469.0875Case 20.962040.47521542.1811694.2981−282.4035−280.4094Case 30.962040.39179117.109651.634270.961271.7094Case 40.971610.41778.855838.683680.446778.805For autumn load, the bus voltage drop is minimized in case 4 of IEEE 69-bus RDN, which is 0.98161 (per unit), which shows more improved value than other cases. The totalP loss reduction for case 4 of IEEE 69-bus RDN is 33.24 kW, while the percentage decrease of total P loss, in this case, is 78.98% which shows that deployment of three DGs in a 69-bus system is more feasible in terms of minimized total P loss with improved bus voltage. The total Q loss for case 4 of IEEE 69-bus RDN is 16.36 kVAR, while the percentage decrease of total Q loss is 77.24% which shows that deployment of three DGs in a 69-bus system is still more feasible in terms of minimized total Q loss with improved bus voltage. However, case 2 of IEEE 69-bus RDN is not feasible due to increasing total P and Q losses. The increase of total P and Q losses of case 2 is 540.53 kW and 242.19 kVAR, respectively. The percentage increase of total P and Q losses of case 2 is 241.68% and 236.64%, respectively, showing that deploying two DGs in 69-bus RDS is not feasible in minimizing total P and Q losses with improved bus voltage. Table 28 presents the result summary for all four cases of IEEE 69-bus RDN with seasonal load data of autumn.Table 28
Result comparison of IEEE 69-bus RDN for autumn.
Casevmin (per unit)FOP loss (kW)Q loss (kVAR)PLOSS reduction (%)QLOSS reduction (%)Case 10.97550.3765152.724823.505166.670967.3284Case 20.97550.39162540.5331242.1968−241.6888−236.6492Case 30.97550.3607549.354221.853568.801669.624Case 40.981610.3708933.246816.367978.983677.2489Moreover, Figures37 and 38 illustrate the graphical comparison of literature results with four and three DGs for IEEE 33-bus RDN, respectively. Figures 39 and 40 illustrate the graphical comparison of literature results with four and three DGs for IEEE 69-bus RDN, respectively. For four DGs in 33-bus RDN, the voltage profile improvement of the proposed HHO is less than that of other literature-based algorithms such as GA, BA, PSO, PPA, ABC and TLCHS. Active power loss reduction is better with the proposed HHO. For three DGs in 33-bus RDN, the voltage profile improvement of the proposed HHO is better than that of other literature-based algorithms such as GA, BA, GSA-GAMS, SSA, EMA, SPEA2, CSCA, PPA, and ACO-ABC. Active power loss reduction is higher with the proposed HHO. For four DGs in 69-bus RDN, the voltage profile improvement of the proposed HHO is better than that of other literature-based algorithms such as GA, SSA, PPA, WCA, HAS, HGWO, and ABC. Active power loss reduction is higher with the proposed HHO, PPA, and ABC. For three DGs in 69-bus RDN, the voltage profile improvement of the proposed HHO is better than that of other literature-based algorithms. Active power loss reduction is higher with the proposed HHO, PPA, and ABC.Figure 37
Graphical comparison of literature results with four DGs (33-bus RDN).Figure 38
Graphical comparison of literature results with three DGs (33-bus RDN).Figure 39
Graphical comparison of literature results with four DGs (69-bus RDN).Figure 40
Graphical comparison of literature results with three DG (69-bus RDN).


## 4.4. Statistical Analysis

This section investigates the statistical analysis for the objective functions taken using the HHO algorithm. The statistical analysis is carried out on the objective function in which voltage and PQ loss reduction (LR) is part of the objective function. The formula for real and reactive power loss is taken from [55]. The technique to calculate the losses is applied from existing literature [55, 81, 82]. The algorithm is simulated fifteen times during the load flow analysis for both IEEE RDNs, and results from the Big-O analysis are shown in Table 29 [54]. It is observed that the changes in the objective functions for the maximum and minimum values taken from the HHO have no significant variations. It validates the efficient and robust behavior of the applied HHO method. However, case 2 of IEEE 33-bus RDN shows a significant difference among best, worst, and mean values of total P loss. It shows that case 2 of 33-bus RDN cannot be taken under consideration during the planning stages. Moreover, case 4 of IEEE 69-bus RDN shows a significant difference among best, worst, and mean values of minimum bus voltages. It also shows that case 4 of IEEE 69-bus RDN cannot be considered during the planning stages.Table 29
Statistical analysis of HHO.
CasesOFBest caseWorst caseMean caseIEEE 33-bus RDNCase 1vmin (per unit)0.865240.865240.86524Case 1PLOSS0.2520120.2190340.245402Case 2vmin (per unit)0.90270.880340.897694Case 2PLOSS0.5804640.2777590.519263Case 3vmin (per unit)0.930540.930540.93054Case 3PLOSS0.6971690.6786740.688568Case 4vmin (per unit)0.930540.930540.93054Case 4PLOSS0.6958930.6784370.683718IEEE 69-bus RDNCase 1vmin (per unit)0.971050.971050.97105Case 1PLOSS0.6732650.6730480.673112Case 2vmin (per unit)0.971050.971050.97105Case 2PLOSS0.6964950.6853160.691088Case 3vmin (per unit)0.971050.971050.97105Case 3PLOSS0.6955840.6804380.690053Case 4vmin (per unit)0.978360.978280.978318Case 4PLOSS0.7748490.7076720.73695Table30 illustrates the analysis of variance (ANOVA). ANOVA tests help determine the variance of the objective function under various algorithms [83]. For IEEE 33-bus RDN, the ANOVA test has been carried out between fourteen (14) algorithms, while for IEEE 69-bus RDN, the test has been performed between twenty (20) algorithms. The value of F for both RDNs is less than the critical level. It proves that the variation found while computing the objection function is significant and not by chance [84].Table 30
ANOVA test.
Source of varianceSum of square (SS)Degree of freedom (df)Mean Square (MS)F ratioP valueF critical [83]Between groups0.008810.00880.30070.58734.1596Within groups0.9116140.0294IEEE 69-bus RDNBetween groups0.073030.02430.71770.54712.8327Within groups1.3903200.0339Moreover, Figures41 and 42 illustrate the box plots for IEEE 33-bus and 69-bus RDNs, respectively. It is apparent that no outliers are present and significant data are available for statistical analysis.Figure 41
Box plot to compare objective function for four cases in IEEE 33-bus RDN.Figure 42
Box plot to compare objective function for four cases in IEEE 69-bus RDN.


## 5. Conclusion and Future Directions

This paper investigates reduction in totalP and Q losses and minimization of bus voltage drop with the optimal DG allocation using the HHO. Four stages of DG integration were investigated in two IEEE RDNs. The objective functions were simultaneous reduction of the total P and Q losses and minimization of bus voltage drop. The summary of the conclusion is as follows:(i)
With the applied HHO, the overall result outcomes have been improved compared to other optimization algorithms for all four cases in both IEEE RDNs.(ii)
The comparison among the four cases for both IEEE RDNs was carried out. The best results were obtained when the given number of DGs could be integrated into the RDNs, as shown by Case 3 of a 33-bus RDN and Case 4 of a 69-bus RDN. Hence, the higher deployment rate of DGs in RDNs has significantly reduced the totalP and Q losses and bus voltage drops with improved results with a robust trend of HHO.



---

*Source: 8640423-2022-06-30.xml*