Next Article in Journal
Design Thinking as a Framework for the Design of a Sustainable Waste Sterilization System: The Case of Piedmont Region, Italy
Next Article in Special Issue
A Novel Opposition-Based Arithmetic Optimization Algorithm for Parameter Extraction of PEM Fuel Cell
Previous Article in Journal
Neural Network Explainable AI Based on Paraconsistent Analysis: An Extension
Previous Article in Special Issue
An Ultra-Low-Power CMOS Supercapacitor Storage Unit for Energy Harvesting Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Three-Phase Harmonic Power Flow Algorithm for Unbalanced Radial Distribution Networks with the Presence of D-STATCOM Devices

1
Department of Electrical & Electronics Engineering, Anil Neerukonda Institute of Technology and Sciences (A), Visakhapatnam 531162, India
2
Department of Electrical Engineering, Andhra University, Visakhapatnam 530003, India
3
Department of Electrical Power & Machines, Faculty of Engineering, Ain Shams University, Cairo 11517, Egypt
4
Energy Technology Program, School of Engineering Technology, Purdue University, West Lafayette, IN 47907, USA
*
Author to whom correspondence should be addressed.
Electronics 2021, 10(21), 2663; https://doi.org/10.3390/electronics10212663
Submission received: 6 October 2021 / Revised: 24 October 2021 / Accepted: 28 October 2021 / Published: 30 October 2021
(This article belongs to the Special Issue Energy Harvesting and Energy Storage Systems)

Abstract

:
Due to the rapid advancement in power electronic devices in recent years, there is a fast growth of non-linear loads in distribution networks (DNs). These non-linear loads can cause harmonic pollution in the networks. The harmonic pollution is low, and the resonance problem is absent in distribution static synchronous compensators (D-STATCOM), which is the not case in traditional compensating devices such as capacitors. The power quality issue can be enhanced in DNs with the interfacing of D-STATCOM devices. A novel three-phase harmonic power flow algorithm (HPFA) for unbalanced radial distribution networks (URDN) with the existence of linear and non-linear loads and the integration of a D-STATCOM device is presented in this paper. The bus number matrix (BNM) and branch number matrix (BRNM) are developed in this paper by exploiting the radial topology in DNs. These matrices make the development of HPFA simple. Without D-STATCOM integration, the accuracy of the fundamental power flow solution and harmonic power flow solution are tested on IEEE−13 bus URDN, and the results are found to be precise with the existing work. Test studies are conducted on the IEEE−13 bus and the IEEE−34 bus URDN with interfacing D-STATCOM devices, and the results show that the fundamental r.m.s voltage profile is improved and the fundamental harmonic power loss and total harmonic distortion (THD) are reduced.

1. Introduction

In terms of harmonics, the loads are classified into two types, linear loads and non-linear loads. A linear load [1] is one which, when supplied by an AC source at fundamental frequency, can produce only fundamental AC currents. Non-linear loads, however, generate harmonic currents. The use of non-linear loads can inject harmonic currents into URDN. These harmonic injections can cause overheating of the equipment, insulation stress on winding in electric machines, added power loss in the equipment, and interference with the communication. Therefore, HPFAs are essential for finding the harmonic distortion level on URDN. In [2], based on current injection, graph theory, and the sparse matrix technique, a three-phase HPFA is proposed. The authors of [3] utilized the decoupled harmonic power flow (DHPF) algorithm to present the results of harmonic power flow calculations. In [4,5], a forward/backward-based HPFA for DN is proposed that considered the special topology of radial distribution networks (RDN). The authors of [6] developed an iterative time-dependent, computer-aided HPFA by combining the time-dependent cross-coupled harmonic model. To obtain this model, large data are received from the practical DNs. Tracing THD in secondary RDN with photovoltaic uncertainties by multiphase HPFA is discussed in [7]. The authors of [8] propose a new combined analytical technique (CAT) for HPFA in the presence of correlated input uncertainties from photovoltaic (PV) systems in RDN. In [9], static capacitors are allocated in shunt along RDN using a Cuckoo search optimization method. For allocating and sizing of capacitors optimally, a flower pollination algorithm is proposed in [10].In [11,12], a novel three-phase power flow algorithm for URDN with multiple integrations of distributed generations (DGs) and a D-STATCOM device is presented. In [13], an electrical energy management in unbalanced distribution networks using virtual power plant concept is presented. In [14], an efficient multi-objective optimization approach based on the supervised big bang–big crunch method for the optimal planning of a dispatchable distributed generator is presented. This approach aims to enhance the system performance indices by the optimal sizing and placement of distributed generators connected to balanced/unbalanced distribution networks. The optimal planning of distributed generators in unbalanced distribution networks using a modified firefly method is presented in [15].
The authors of [16] examine the utilization of D-STATCOM without a capacitor to compensate for power quality in DNs. The optimal D-STATCOM allocation in DNs is discussed in [17,18]. In [19], an optimal algorithm to control a three-phase D-STATCOM is proposed. This algorithm can give harmonic compensation as well as reactive power compensation in linear and non-linear loads, which are connected in three-phase. In [20], for minimizing the total real power loss in DNs with the interfacing of DGs, plug-in-hybrid electric vehicles (PHEVs), and D-STATCOM, a genetic algorithm is proposed. A control technique is developed in [21] for D-STATCOM for extracting the fundamental weight components from non-sinusoidal load currents to produce grid reference currents. For harmonics elimination, the injection of reactive power and balancing of load, this D-STATCOM is developed. The D-STATCOM’s performance is examined in different working modes. The combination of two problems such as the reconfiguration and interfacing of D-STATCOM can be solved by using the grey wolf optimization (GWO) method proposed in [22].
The proposed power flow algorithm (PFA) can give both fundamental and harmonic solutions. The solution of the fundamental power flow algorithm (FPFA) discussed in this paper is used in modelling the linear and non-linear loads for HPFA. The BNM and BRNM developed in this paper make the implementation of the PFA simple. The bus numbers and branch numbers of newly created sections of RDN are stored in BNM and BRNM, respectively. This paper is arranged in the following order. The network components’ modelling is addressed in Section 2. The algorithm to develop BNM and BRNM is discussed in Section 3. In Section 4, the three-phase HPFA with the integration of the D-STATCOM device is discussed. Section 5 presents the test studies and discussions on the IEEE−13 bus and IEEE−34 bus URDN. Section 6 discusses the concluding remarks.

2. Network Components and Their Modeling

The URDN includes the main components such as lines, three-phase transformers, three-phase capacitor banks, and loads. These components are briefly modeled in this section.

2.1. Overhead or Underground Distribution Lines

With the Carson’s equations presented in [23], the primitive impedance matrices for three-phase overhead and underground lines can be formed. For a grounded neural system, these matrices are reduced to phase impedance matrices of 3 × 3 size using Kron reduction. Figure 1 shows the three-phase distribution line model, and its shunt admittance is neglected due to its small effect. The phase impedance matrix for the line section ‘jk’ is given in Equation (1).
Z a b c j k = Z a a Z a b Z a c Z b a Z b b Z b c Z c a Z c b Z c c j k
From Figure 1, the relationship between the phase voltage matrices of bus-j and bus-k is given in Equation (2):
V a V b V c k = V a V b V c j Z a a Z a b Z a c Z b a Z b b Z b c Z c a Z c b Z c c j k I a I b I c j k
The reactance of line is regarded as proportionate to the harmonic order for HPFA. For h-order harmonic frequency, the self-impedance of phase ‘a’ is given in Equation (3),
Z a a h = R a a + j h X a a

2.2. Loads

The phase current matrix and line current matrix serving the different types of three-phase loads are outlined in Table 1. Detailed discussion on Table 1 is provided in [23].

2.2.1. Linear Loads

These loads produce only fundamental sinusoidal response upon supplied by sinusoidal source. The liner loads can be modelled in several ways [1]. Each model will show a different impact on harmonic analysis. The impedance modelling of these loads is taken as series combination of R and X.

2.2.2. Non-Linear Loads

With the harmonic spectrum of non-linear loads and their load current obtained from the fundamental power flow, these loads are modelled as constant current sources [24]. The magnitude of the current source is obtained with Equation (4), and its phase angle is obtained with Equation (5):
I h = I r a t e d I h _ s p e c t r u m I 1 _ s p e c t r u m
The phase angle of the current source is obtained as:
θ h = θ h _ s p e c t r u m + h θ 1 θ h _ s p e t r u m
where:
θ 1 : Phase angle of the rated current at fundamental frequency;
θ h _ s p e c t r u m : Phase angle of the harmonic source current spectrum.

2.3. Capacitor Banks

Modelling of the capacitor banks is presented in Table 2.
For HPFA, the capacitive susceptance (B) is to be multiplied with ‘h’ for ‘h’ order frequency.

2.4. Tree-Phase Transformer

The fundamental voltage and current relationships between the primary and secondary sides for different transformer connections are presented in [25]. The modelling of the three-phase transformers for HPFA is given in [24,26].

2.5. D-STATCOM

The D-STATCOM is commonly regarded as a shunt compensator which supplies reactive power in PFAs. The voltage magnitude at the D-STATCOM bus can be controlled by adjusting the reactive power injection of D-STATCOM.
The interface of the D-STATCOM at ith bus shown in Figure 2a, and its traditional modelling for PFAs is shown in Figure 2b. The specified reactive power of the load is combined with the reactive power output of D-STATCOM, so that reactive power varies as magnitude of Vi varies. This is absolutely a PV bus modelling with the real power output of the D-STATCOM set to zero [27,28]. The hypothesis in this model is that losses in the D-STATCOM and its connection are ignored. The D-STATCOMs have low harmonic content, so the harmonic current injected by the D-STATCOM is considered as zero for HPFA.

3. Algorithm for Developing BNM and BRNM

The performance of the HPFA of URDN is enhanced by the systematic numbering of buses and branches. From [29], the numbering scheme to buses and branches is taken. The following steps are to be followed to write a Software Code in order to split the URDN into different sections, as shown in Figure 3.
  • From the distribution network shown in Figure 3, form Table 3.
  • Start with BN = 1.Read the RE of BN, i.e., 2. Then, check how many times this 2 appears in the SE row. Inthe above table, it appears one time. That means bus 2 is the sending end for only one branch. Fill these RE 2 and BN 1 in two different matrices (BNM and BRNM) as the first row and first column elements. Then, increase the column number by one.
  • Increase the BN (i.e., BN = 2), and read the RE of BN, i.e., 3. Then, as in step 1, check for the appearance of 3 in SE row. The bus 3 appears two times. That means that, from the bus 3, two branches are leaving. Then, fill these RE 3 and BN 2 into the same variables as the first row and present the column elements. Name this row elements as section-I. Now increase the row number by one and set the column number to one.
  • Similarly, increase the BN, and read the RE of BN. Then, check for the appearance of this RE in the SE row. If it appears one time, then fill these RE and BN values as the present row and present column elements of the variables BNM and BRNM. Then increase the column number by one and repeat step 4. If it does not appears or appears more than one time in the SE row, then fill the corresponding RE and BN values as present row and present column elements. Then identify this row as a section. Then increase the row number by one and set the column number to one and repeat the step 4.
The above steps are repeated until the BN value reaches the last branch number. At the end, the BNM and BRNM are obtained as follows: Electronics 10 02663 i001

4. Three-Phase HPFA with Non-Linear Loads and D-STATCOM Devices

For modeling the linear and non-linear loads for HPFA, the fundamental power flow solution is required. Hence, the algorithm consists of two parts. PartA illustrates the iterative procedure for FPFA with the D-STATCOM device and PartB illustrates the iterative procedure for HPFA with linear loads, non-linear loads, and D-STATCOM devices.
After developing the BNM and BRNM for the URDN, the iterative procedure is explained with the following steps.
PartA: FPFA with D-STATCOM
  • The voltages at all busses are assigned as substation bus voltage.
    V a V b V c = 1 0 o 1 120 o 1 120 o
  • Find the line current matrix serving the load at all buses.
  • Start with collecting line current matrix at bus−23 (the tail bus in section-VII in BNM), and thereby find the line current matrix for branch−22 (the tail branch in section-VII in BRNM). Then, continue to the bus−22 and branch−21 to find the line current matrix at the bus and line current matrix in branch, respectively. From Figure 4, the following equations are obtained by applying the KCL at every bus:
    I a b c k = I L a b c k + I s h a b c k + I C a b c k
    I a b c j k = I a b c k
    I a b c j = I a b c j k + I L a b c j + I s h a b c j + I C a b c j
    I a b c i j = I a b c j
    I a b c i = I a b c i j
    where:
    I a b c k : Line current matrix at bus-k;
    I L a b c j k : Line current in branch-jk;
    I L a b c k : Load current matrix at bus-k;
    I s h a b c k : Line current matrix drawn by shunt admittance at bus-k;
    I C a b c k : Line current matrix drawn by capacitor bank at bus-k, if any.
  • Now go to section-VI and repeat procedure as in step 5 to find the line current matrix at the head bus and line current matrix for head branch. Similarly, proceed up to section-I and find the line current matrix up to bus−1 and line current matrix up to branch−1.
  • Now start with head bus in section-I and continue to the tail bus in section-I by finding the phase voltage matrix at all buses with Equation (2). Then, go to the next section and repeat the same procedure.
  • Steps 4 to 6 are to be repeated until the convergence criterion as given in Equation (14) is satisfied:
    V a b c i r V a b c i r 1 ε a b c
    where ‘r’ is the iteration number.
  • D-STATCOM location is selected and model as PV bus for the outside γthiteration.
  • The mismatches in voltages at D-STATCOM buses are obtained with Equation (13):
    Δ V a Δ V b Δ V c γ = V a s p V b s p V c s p V a c a l V b c a l V c c a l γ
    Δ V a b c γ ε a b c
    where Δ V γ is the mismatch matrix for the voltage and its size is 3 n × 1 , and ‘n’ is the total number of PV buses.
  • If the Equation (16) is not satisfied, then the incremental current injection matrix at D-STATCOM bus is calculated with Equation (17) to maintain the specified voltages:
    Δ I γ = Z P V 1 Δ V γ
    where Z P V is the sensitivity matrix for the PV bus with its size 3 n × 3 n . The formation of this matrix is presented in [30].
  • The incremental reactive current injection matrix at D-STATCOM bus is obtained with Equation (18):
    Δ I D , a Δ I D , b Δ I D , c j γ = Δ I a c o s 9 0 0 + δ v , a + j * s i n 9 0 0 + δ v , a Δ I b c o s 9 0 0 + δ v , b + j * s i n 9 0 0 + δ v , b Δ I c c o s 9 0 0 + δ v , c + j * s i n 9 0 0 + δ v , c j γ
  • In Figure 5, by applying the KCL at bus-j, the line current matrix in branch-ij is obtained as:
    I a I b I c i j γ = I L a I L b I L c j γ Δ I D , a Δ I D , b Δ I D , c j γ
    With V a b c j γ and I a b c i j γ , the reactive power flow in the line Q a b c i j γ is evaluated. Then, the incremental reactive current injection matrix is obtained with Equation (20):
    Δ Q D , a Δ Q D , b Δ Q D , c j γ = Q L a Q L b Q L c j γ Q a Q b Q c i j γ
    The reactive power generation matrix needed at D-STATCOM bus-j is obtained with Equation (21):
    Q D , a Q D , b Q D , c j γ = Q D , a Q D , b Q D , c j γ 1 + Δ Q D , a Δ Q D , b Δ Q D , c j γ
  • If the D-STATCOM device is able to generate limited reactive power, then find the total reactive power generation of D-STATCOM device with Equation (22). The total reactive power generation of D-STATCOM is now compared with the maximum and minimum limits of reactive power generation of D-STATCOM device limits. Equation (22) is calculated as follows:
    Q D j γ = Q D , a j γ + Q D , b j γ + Q D , c j γ
    If Q j , min Q D j γ Q j , max
    Then set complex power generation is as in Equation (21)
    If Q D j γ Q j , min
    Then set Q D j γ = Q j , min and Q D , a j γ = Q D , b j γ = Q D , c j γ = Q j , m i n / 3
    If Q D j γ Q j , max
    Then set Q D j γ = Q j , max and Q D , a j γ = Q D , b j γ = Q D , c j γ = Q j , m a x / 3
  • Now, find the complex power generation matrix at D-STATCM bus with Equation (23):
    S D , a S D , b S D , c j γ = P D , a P D , b P D , c j + j Q D , a Q D , b Q D , c j γ
    where P D , a b c j is the specified real power generation matrix of the D-STATCOM device and its value is set to zero.
  • The line current matrix injected by the D-STATCOM is obtained with the complex power generation matrix obtained in Equation (23) and bus voltage matrix as:
    I D , a b c j γ = S D , a / V a S D , b / V b S D , c / V c j γ
  • Using the current injection matrix at the D-STATCOM buses, repeat from step 7 by setting γ = γ+1.
  • If Equation (16) is satisfied at all D-STATCOM buses, then stop the FPFA algorithm.
  • With the complex power loss in branch-ij in Equation (25), find the total power loss in the network by summing up the losses in all branches:
    S L o s s a S L o s s b S L o s s c i j = V a i I a i j V b i I b i j V c i I c i j V a j I a j i V b j I b j i V c j I c j i
PartB: HPFA with non-linear loads and D-STATCOM device.
18.
With the converged bus voltages and specified load, the impedances of the linear loads are calculated for the harmonic order-h of interest.
19.
Find the harmonic current injection matrix for the non-linear loads for the selected h-order harmonic of interest. The harmonic current injection matrix of D-STATCOM is taken as zero.
20.
The harmonic voltage at the substation bus is taken as zero since the supply voltage is assumed to bea pure sinusoidal voltage waveform.
21.
The harmonic voltages at all other buses for the first iteration are assumed to be zeros as that of the substation bus:
V a V b V c h = 0 0 0
22.
Find the net harmonic current matrix at all the buses with the harmonic current matrix drawn by the linear loads and the harmonic current injection matrix of non-linear loads and the D-STATCOM device. The current matrix drawn by the linear loads at all the buses is zero for the first iteration as the harmonic voltage at all the buses is zero for the first iteration. This is illustrated with the sample section as shown in Figure 6. The net harmonic current matrix at bus-j is given by Equation (27), and the harmonic current matrix in branch-ij is given by Equation (28):
I a b c j h = I S a b c j h I D a b c j h + I L a b c j h
I a b c i j h = I a b c j h
where:
I a b c j h : Harmonic current matrix at bus-j for harmonic order-h;
I a b c i j h : Harmonic current matrix in branch-ij for harmonic order-h;
I L a b c j h : Harmonic current matrix drawn by linear load at bus-j for harmonic order-h;
I S a b c j h : Harmonic current injection matrix by non-linear load at bus-j for harmonic order-h;
I D a b c j h : Harmonic current injection matrix by D-STATCOM device at bus-j for harmonic order-h.
Likewise, the harmonic currents in all branches are to be obtained by moving up to the substation as explained in step 3 to step 4 in PartA for FPFA.
23.
Then, start finding the harmonic voltages at all buses downstream from the substation bus with Equation (29) as explained in step 5 in PartA:
V a V b V c j h = V a V b V c i h Z a a Z a b Z a c Z b a Z b b Z b c Z c a Z c b Z c c i j h I a I b I c i j h
24.
Repeat the steps 22 to 23 until the magnitude mismatch of harmonic voltages of h-order at all the busses is within the tolerance limit.
25.
Find the harmonic power loss in all branches using Equation (30). Then find the total harmonic power loss in the network for the selected harmonic order-h using Equation (31):
S L o s s a S L o s s b S L o s s c i j h = V a i I a i j V b i I b i j V c i I c i j h V a j I a j i V b j I b j i V c j I c j i h
T S _ l o s s a b c h = b r = 1 N b r S L o s s a b c b r h
26.
Likewise, repeat the steps from 10 to 16 for all the harmonics of selected harmonic orders (h = 3, 5, 7, 9, 11, 13, and 15).
27.
Find the total harmonic loss of the network using Equation (32):
T o t a l _ l o s s = h = h o h m b r = 1 N b r S L o s s a b c b r h
28.
The total r.m.s voltage at bus-i, say, phase ‘a’, is calculated as:
V a i = V a i 1 2 + h = h o h m V a i h 2
29.
The total harmonic distortion at every bus is calculated using Equation (34):
T H D i a = h = h o h m V a i h 2 V a i 1
where:
h o : Minimum harmonic order;
h m : Maximum harmonic order;
b r : Branch number;
N b r : Total number of branches.

5. Results and Discussions

5.1. IEEE−13 Bus URDN

5.1.1. Fundamental Power Flow Solution for Accuracy Test

The proposed three-phase FPFA is examined on IEEE−13 bus unbalanced test feeder without interfacing of D-STATCOM device. Figure 7 shows the IEEE−13 bus feeder and its data is collected from [31]. 5000 kVA and 4.16 kV are the chosen base values for this network. The FPFA is taken 5 iterations for its convergence with tolerance for convergence is 10−4. The comparison of obtained power flow solution with IEEE solution and errors in voltage magnitudes and phase angles at every bus are presented in Table 4. Table 5 presents the comparison of obtained power loss with the IEEE losses. Insignificant values of maximum errors of 0.0005 p.u and 0.010o for voltage magnitudes and phase angles are observed in Table 6. So that, in terms of accuracy the test results are consistent with IEEE results.

5.1.2. Fundamental and Harmonic Power Flow Solutions without D-STATCOM

The regulator between buses 650 and 632 is removed and the capacitor banks at bus 675 and 611 are removed from the network. The data for the harmonic load composition and current spectra of harmonic loads aretaken from [32]. The convergence tolerance is taken as 10−4. Table 6 presents the harmonic power losses and total power loss of the network including fundamental and harmonic loss. The harmonic voltage solutions for the selected range of harmonics of order 3, 5, 7, 9, 11, 13, and 15 are presented in Table 7. Table 8 presents the fundamental r.m.s profile, the total harmonic voltage profile, and the THD %. It is observed from Table 8 that the maximum THD % on the network is 5.2263 at bus−611 for c-phase, and in [2], it was reported that the maximum THD % at bus−611 is 5.23. Therefore, the results of the proposed HPFA are almost matches the literature in terms of accuracy. To see the impact of the D-STATCOM on the fundamental r.m.s voltage profile, the total r.m.s voltage profile, the fundamental and harmonic power loss, and the THD %, the results of this case study are taken as benchmarks.

5.1.3. IEEE−13 Bus URDN with D-STATCOM

In this case, a three-phase D-STATCOM is integrated at bus 680. The D-STATCOM is modelled as a PV model with its real power generation set to zero and the lower limit and upper limit for the three-phase reactive power generation are 100 kVAR and 1000 kVAR, respectively. The phase voltages specified at this bus are 1 p.u. Table 8 presents the harmonic power loss and total power loss (including fundamental and harmonic power loss). In comparison with Table 6, it is observed that the integration of the D-STATCOM into the network reduces both the fundamental and harmonic power losses, thereby the total power loss in the network is also reduced. Table 9 presents the fundamental r.m.s voltage profile, the total r.m.s voltage profile, and the THD %.In comparison with Table 7, it is observed that there is an improvement in fundamental r.m.s voltage profile. The minimum fundamental r.m.s voltage in the network without D-STATCOM is 0.8651 p.uat bus−611 for c-phase, whereas its value is 0.8763 p.u at bus−611 for c-phase with integration of D-STATCOM. The maximum THD % in the network is reduced from 5.2263 to 5.1133 with integration of D-STATCOM. Figure 8 shows the comparison of THD % with and without integration of D-STATCOM. Figure 9 presents the comparison of fundamental r.m.s voltages on the network for the two case studies.

5.2. IEEE−34 Bus URDN

The date for the IEEE−34 bus URDN is taken from [31]. The base values selected for the system are 2500 kVA and 24.9 kV. The load composition at spot loads for harmonic analysis is presented in Table 10. The data for the current spectra of harmonic loads aretaken from [32]. The convergence tolerance for both fundamental and harmonic power flows is 10−4. The case studies on the network are presented in Table 11. The rating and location of the D-STATCOM device for Case2 is presented in Table 11. Table 12 presents the fundamental r.m.s voltage profile, the total r.m.s voltage profile, and the THD % for Case1. The test results of Case1 are used as a benchmark to see the fundamental and harmonic impacts of D-STATCOM on the network. The summary of results for the case studies is presented in Table 13. In Case2, which has integrations of two D-STATCOM devices, the maximum THD% is observed to be 5.2567 which is less than in Case 1. The number of phases effected with a THD% more than fiveis reduced from fourto twofrom Case 1 to Case 2. From Case 1 to Case 2, it is found that the minimum fundamental voltage on the network is improved from 0.7641 p.u to 0.8137 p.u at bus 890 for the a-phase.The fundamental power loss and the total power loss including harmonic loss of the network reduced in Case 2 in comparison with Case 1. Figure 10 shows the comparison of THD % with and without the integration of the D-STATCOM. Figure 11 shows the comparison of the fundamental r.m.s voltages on the network for the two case studies.

6. Conclusions

This paper proposes new three-phase PFAs for URDN with the presence of linear and non-linear loads and D-STATCOM devices. These PFAs can give both fundamental and harmonic power flow solutions with/without the presence of D-STATCOM devices. The developed BNM and BRNM make both the FPFA and HPFA simple. These matrices are developed by exploiting the radial topology in distribution networks. This method uses the basic concepts of circuit theory, and they can be easily understood. In this paper, the linear loads are modeled as a series combination of resistance and reactance, and non-linear loads are modeled as constant current sources with its magnitude and angle obtained from the current spectra. The harmonic current injections from the D-STATCOM are assumed as zero. The proposed FPFA and HPFA are tested on the IEEE−13 bus URDN, and the results are found to be accurate with the literature. The test studies are carried on the IEEE−13 bus and IEEE−34 bus URDN, and the results of the case studies show thatthere is an improvement in the fundamental voltage profile, a reduction in the fundamental and harmonic power loss, and a reduction in THD% with the integration of D-STATCOM devices.

Author Contributions

R.S. and K.V. designed the problem under study, performed the simulations and obtained the results. R.S. and K.V. wrote the paper, which was further reviewed by A.Y.A. and A.E.-S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burch, R.; Chang, G.K.; Grady, M.; Hatziadoniu, C.; Liu, Y.; Marz, M.; Ortmeyer, T.; Xu, W.; Ranade, S.; Ribeiro, P. Impact of aggregate linear load modeling on harmonic analysis: A comparison of common practice and analytical models. IEEETrans. Power Deliv. 2003, 18, 625–630. [Google Scholar] [CrossRef]
  2. Yang, N.-C.; Le, M.-D. Three-phase harmonic power flow by direct ZBUS method for unbalanced radial distribution systems with passive power filters. IET Gener. Transm. Distrib. 2016, 10, 3211–3219. [Google Scholar] [CrossRef]
  3. Milovanović, M.; Radosavljević, J.; Perović, B.; Dragičević, M. Power flow in radial distribution systems in the presence of harmonics. Int. J. Electr. Eng. Comput. 2018, 2, 11–19. [Google Scholar] [CrossRef]
  4. Amini, M.A.; Jalilian, A.; Behbahani, M.R.P. Fast network reconfiguration in harmonic polluted distribution network based on developed backward/forward sweep harmonic load flow. Electr. Power Syst. Res. 2019, 168, 295–304. [Google Scholar] [CrossRef]
  5. Milovanović, M.; Radosavljević, J.; Perović, B. A backward/forward sweep power flow method for harmonic polluted radial distribution systems with distributed generation units. Int. Trans. Electr. Energy Syst. 2019, 30, e12310. [Google Scholar] [CrossRef]
  6. Nduka, O.S.; Ahmadi, A.R. Data-driven robust extended computer-aided harmonic power flow analysis. IET Gener. Transm. Distrib. 2020, 14, 4398–4409. [Google Scholar] [CrossRef]
  7. Hernandez, J.C.; Ruiz-Rodriguez, F.J.; Jurado, F.; Sanchez-Sutil, F. Tracing harmonic distribution and voltage unbalance in secondary radial distribution networks with photovoltaic uncertainties by a multiphase harmonic load flow. Electr. Power Syst. Res. 2020, 185, 1–18. [Google Scholar] [CrossRef]
  8. Ruiz-Rodriguez, F.J.; Hernandez, J.C.; Jurado, F. Iterative harmonic load flow by using the point-estimate method and complex affine arithmetic for radial distribution systems with photovoltaic uncertainties. Electr. Power Energy Syst. 2020, 118, 1–16. [Google Scholar] [CrossRef]
  9. El-Fergany, A.; Abdelaziz, A.Y. Cuckoo Search-based Algorithm for Optimal Shunt Capacitors Allocations in Distribution Networks. Electr. Power Compon. Syst. 2013, 41, 1567–1581. [Google Scholar] [CrossRef]
  10. Abdelaziz, A.Y.; Ali, E.S.; Abd Elazim, S.M. Flower Pollination Algorithm for Optimal Capacitor Placement and Sizing in Distribution Systems. Electr. Power Compon. Syst. J. 2016, 44, 544–555. [Google Scholar] [CrossRef]
  11. Satish, R.; Kantarao, P.; Vaisakh, K. A new algorithm for impacts of multiple DGs and D-STATCOM in unbalanced radial distribution networks. Int. J. Renew. Energy Technol. 2021, 12, 221–242. [Google Scholar] [CrossRef]
  12. Satish, R.; Vaisakh, K.; Abdelaziz, A.Y.; El-Shahat, A. A Novel Three-phase Power Flow Algorithm for Evaluation of Impact of Renewable Energy Sources and D-STATCOM Device in Unbalanced Radial Distribution Networks. Energies 2021, 14, 6152. [Google Scholar] [CrossRef]
  13. Othman, M.M.; Hegazy, Y.G.; Abdelaziz, A.Y. Electrical Energy Management in Unbalanced Distribution Networks using Virtual Power Plant Concept. Electr. Power Syst. Res. 2017, 145, 157–165. [Google Scholar] [CrossRef]
  14. Abdelaziz, A.Y.; Hegazy, Y.G.; El-Khattam, W.; Othman, M.M. A Multi objective Optimization for Sizing and Placement of Voltage Controlled Distributed Generation Using Supervised Big Bang Big Crunch Method. Electr. Power Compon. Syst. 2015, 43, 105–117. [Google Scholar] [CrossRef]
  15. Abdelaziz, A.Y.; Hegazy, Y.G.; El-Khattam, W.; Othman, M.M. Optimal Planning of Distributed Generators in Distribution Networks Using Modified Firefly Method. Electr. Power Compon. Syst. 2015, 43, 320–333. [Google Scholar] [CrossRef]
  16. Rohouma, W.; Balog, R.S.; Peerzada, A.A.; Begovic, M.M. D-STATCOM for harmonic mitigation in low voltage distribution network with high penetration of nonlinear loads. Renew. Energy 2020, 145, 1449–1464. [Google Scholar] [CrossRef]
  17. Sirjani, R.; Jordehi, A.R. Optimal placement and sizing of distribution static compensator (D-STATCOM) in electric distribution networks: A review. Renew. Sustain. Energy Rev. 2017, 77, 688–694. [Google Scholar] [CrossRef]
  18. Rezaeian-Marjani, S.; Galvani, S.; Talavat, V.; Farhadi-Kangarlu, M. Optimal allocation of D-STATCOM in distribution networks including correlated renewable energy sources. Electr. Power Energy Syst. 2020, 122, 1–14. [Google Scholar] [CrossRef]
  19. Patel, S.K.; Arya, S.R.; Maurya, R. Optimal Step LMS-Based Control Algorithm for DSTATCOM in Distribution System. Electr. Power Compon. Syst. 2019, 47, 675–691. [Google Scholar] [CrossRef]
  20. Singh, B.; Singh, S. GA-based optimization for integration of DGs, STATCOM and PHEVs in distribution systems. Energy Rep. 2019, 5, 84–103. [Google Scholar] [CrossRef]
  21. Badoni, M.; Singh, A.; Singh, B.; Saxena, H. Real-time implementation of active shunt compensator with adaptive SRLMMN control technique for power quality improvement in the distribution system. Int. Trans. Electr. Energy Syst. 2020, 14, 1598–1606. [Google Scholar] [CrossRef]
  22. Selvaraj, G.; Rajangam, K. Multi-objective grey wolf optimizer algorithm for combination of network reconfiguration and D-STATCOM allocation in distribution system. Int. Trans. Electr. Energy Syst. 2019, 29, e12100. [Google Scholar] [CrossRef]
  23. Kersting, W.H. Distribution System Modeling and Analysis, 4th ed.; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  24. Task Force on Harmonics Modeling and Simulation. Modeling and simulation of the propagation of harmonics in electric power networks I. Concepts, models, and simulation techniques. IEEE Trans. Power Deliv. 1996, 11, 452–465. [Google Scholar] [CrossRef]
  25. Chen, T.H.; Chang, J.D.; Chang, Y.L. Models of grounded mid-tap open-wye and open-delta connected transformers for rigorous analysis of a distribution system. IEEE Proc. Gener. Transm. Distrib. 1996, 143, 82–88. [Google Scholar] [CrossRef]
  26. Arrillaga, J.; Bradley, D.A.; Bodger, P.S. Power System Harmonics, 1st ed.; Wiley: Hoboken, NJ, USA, 1985. [Google Scholar]
  27. Yang, Z.; Shen, C.; Crow, M.L.; Zhang, L. An improved STATCOM model for power flow analysis. In Proceedings of the 2000 Power Engineering Society Summer Meeting, Seattle, WA, USA, 16–20 July 2000; pp. 1121–1126. [Google Scholar]
  28. Jazebi, S.; Hosseinian, S.H.; Vahidi, B. DSTATCOM allocation in distribution networks considering reconfiguration using differential evolution algorithm. Energy Convers. Manag. 2011, 52, 2777–2783. [Google Scholar] [CrossRef]
  29. Das, D.; Nagi, H.S.; Kothari, D.P. Novel method for solving radial distribution networks. IEEE Proc. Gener. Transm. Distrib. 1994, 141, 291–298. [Google Scholar] [CrossRef]
  30. Shirmohammadi, D.; Carol, S.; Cheng, A. A Three phase power flow method for real time distribution system analysis. IEEE Trans. Power Syst. 1995, 10, 671–679. [Google Scholar]
  31. Radial Distribution Test Feeders. Available online: http://sites.ieee.org/pes-testfeeders/resources (accessed on 28 September 2021).
  32. Abu-Hashim, R.; Burch, R.; Chang, G.; Grady, M.; Gunther, E.; Halpin, M.; Xu, W.; Marz, M.; Sim, T.; Liu, Y.; et al. Test systems for harmonic modeling and simulation. IEEE Trans. Power Deliv. 1999, 14, 579–587. [Google Scholar] [CrossRef]
Figure 1. A sample three-phase distribution line.
Figure 1. A sample three-phase distribution line.
Electronics 10 02663 g001
Figure 2. (a) D-STATCOM interface at ith bus 2 (b) Traditional modelling of D-STATCOM as PV bus.
Figure 2. (a) D-STATCOM interface at ith bus 2 (b) Traditional modelling of D-STATCOM as PV bus.
Electronics 10 02663 g002
Figure 3. Divided sections for sample distribution network.
Figure 3. Divided sections for sample distribution network.
Electronics 10 02663 g003
Figure 4. A simple URDN three busses.
Figure 4. A simple URDN three busses.
Electronics 10 02663 g004
Figure 5. A simple URDN with two buses with D-STATCOM placed at bus-j.
Figure 5. A simple URDN with two buses with D-STATCOM placed at bus-j.
Electronics 10 02663 g005
Figure 6. Sample section of two buses for HPFA.
Figure 6. Sample section of two buses for HPFA.
Electronics 10 02663 g006
Figure 7. IEEE 13 Bus URDN.
Figure 7. IEEE 13 Bus URDN.
Electronics 10 02663 g007
Figure 8. Comparison of THD % for case studies on IEEE−13 URDN.
Figure 8. Comparison of THD % for case studies on IEEE−13 URDN.
Electronics 10 02663 g008
Figure 9. Comparison of fundamental r.m.svoltages for case studies on IEEE−13 URDN.
Figure 9. Comparison of fundamental r.m.svoltages for case studies on IEEE−13 URDN.
Electronics 10 02663 g009
Figure 10. Comparison of THD% for case studies on IEEE−34 URDN.
Figure 10. Comparison of THD% for case studies on IEEE−34 URDN.
Electronics 10 02663 g010
Figure 11. Comparison of fundamental r.m.svoltages for case studies on IEEE−34 URDN.
Figure 11. Comparison of fundamental r.m.svoltages for case studies on IEEE−34 URDN.
Electronics 10 02663 g011
Table 1. Load modelling.
Table 1. Load modelling.
Wye ConnectionDelta Connection
Phase voltage matrix and specified load matrix at bus. V a n δ a V b n δ b V c n δ c , S a θ a S b θ b S c θ c V a b δ a b V b c δ b c V c a δ c a , S a b θ a b S b c θ b c S c a θ c a
Phase current matrix serving constant power load I L a I L b I L c = S a V a n S b V b n S c V c n I L a b I L b c I L c a = S a b V a b S b c V b c S c a V c a
Phase current matrix serving the constant impedance load Z L a Z L b Z L c = V a n 2 S a V b n 2 S b V c n 2 S c
I L a I L b I L c = V a n Z L a V b n Z L b V c n Z L c
Z L a b Z L b c Z L c a = V a b 2 S a b V b c 2 S b c V c a 2 S c a
I L a b I L b c I L c a = V a b Z L a b V b c Z L b c V c a Z L c a
Phase current matrix serving the constant current load I L a I L b I L c = I L a δ a θ a I L b δ b θ b I L c δ c θ c I L a b I L b c I L c a = I L a b δ a b θ a b I L b c δ b c θ b c I L c a δ c a θ c a
Line current matrix entering the load I L a I L b I L c I L a I L b I L c = 1 0 1 1 1 0 0 1 1 I L a b I L b c I L c a
Distributed LoadsCreate a duplicate node at a distance of one-fourth the length from the sending end and connect a two-third of lode. At the receiving end one-third of load is connected.
Table 2. Capacitor banks modeling.
Table 2. Capacitor banks modeling.
Wye ConnectedDelta Connected
Phase voltage matrix and specified reactive power matrix at bus. V a n θ a V b n θ b V c n θ c , Q a Q b Q c V a b θ a b V b c θ b c V c a θ c a , Q a b Q b c Q c a
Phase current matrix serving the capacitor bank B a b c = Q a V a 2 Q b V b 2 Q c V c 2
I C a I C b I C c = j B a V a n j B b V b n j B c V c n
B a b c = Q a b V a b 2 Q b c V b c 2 Q c a V c a 2
I C a b I C b c I C c a = j B a b V a b j B b c V b c j B c a V c a
Line current matrix serving the capacitor bank I C a I C b I C c I C a I C b I C c = 1 0 1 1 1 0 0 1 1 I C a b I C b c I C c a
Table 3. Branch numbering of distribution network in Figure 2.
Table 3. Branch numbering of distribution network in Figure 2.
Branch Number(BN)Sending Bus(SE)Receiving Bus(RE)
112
223
334
445
556
667
778
889
9910
101011
11312
121213
131314
141415
151516
16617
171718
181819
191920
202021
211822
222223
Table 4. Fundamental voltage solution for IEEE−13 bus URDN.
Table 4. Fundamental voltage solution for IEEE−13 bus URDN.
BusPhaseObtained SolutionIEEE Solution [31]Error in Voltage Mag.Error in Voltage Ang.
650a1 0o1 0o0.00000.00
b1 −120o1 −120o0.00000.00
c1 120o1 120o0.00000.00
RGa1.0625 0o1.0625 0o0.00000.00
b1.0500 −120o1.0500 −120o0.00000.00
c1.0687 120o1.0687 120o0.00000.00
632a1.0210 −2.49o1.0210 −2.49o0.00000.00
b1.0420 −121.72o1.0420 −121.72o0.00000.00
c1.0175 117.83o1.0170 117.83o−0.00050.00
671a0.9900 −5.30o0.9900 −5.30o0.00000.00
b1.0529 −122.34o1.0529 −122.34o0.00000.00
c0.977 116.03o0.9778 116.02o0.0001−0.01
680a0.9900 −5.30o0.9900 −5.30o0.00000.00
b1.0529 −122.34o1.0529 −122.34o0.00000.00
c0.9778 116.03o0.977 116.02o0.0001−0.01
633a1.0180 −2.55o1.0180 −2.56o0.00000.01
b1.0401 −121.77o1.0401 −121.77o0.00000.00
c1.0148 117.82o1.0148 117.82o0.00000.00
634a0.9940 −3.23o0.9940 −3.23o0.00000.00
b1.0218 −122.22o1.0218 −122.22o0.00000.00
c0.9960 117.35o0.9960 117.34o0.0000−0.01
645b1.0328 −121.90o1.0329 −121.90o0.00010.00
c1.0155 117.86o1.0155 117.86o0.00010.00
646b1.0311 −121.98o1.0311 −121.98o0.00000.00
c1.0134 117.90o1.0134 117.90o0.00000.01
692a0.9900 −5.30o0.9900 −5.31o0.00000.01
b1.0529 −122.34o1.0529 −122.34o0.00000.00
c0.9778 116.03o0.9777 116.02o−0.0001−0.01
675a0.9835 −5.55o0.9835 −5.56o0.00000.01
b1.0553 −122.52o1.0553 −122.52o0.00000.00
c0.9759 116.04o0.9758 116.03o−0.0001−0.01
684a0.9881 −5.32o0.9881 −5.32o0.00000.00
c0.9758 115.92o0.9758 115.92o0.00000.00
611c0.9738 115.78o0.9738 115.78o0.00000.00
652a0.9825 −5.24o0.9825 −5.25o0.00000.01
Table 5. Power loss in IEEE−13 bus URDN.
Table 5. Power loss in IEEE−13 bus URDN.
PhaseObtained Power LossIEEE Loss [31]
Active
(kW)
Reactive
(kVAR)
Active
(kW)
Reactive
(kVAR)
a39.13152.6239.11152.59
b−4.7442.27−4.7042.22
c76.59129.6976.65129.85
Total110.98324.57111.13324.66
Table 6. Power loss in IEEE−13 bus URDN.
Table 6. Power loss in IEEE−13 bus URDN.
Harmonic OrderHarmonic Power Loss
APL (kW)RPL (kVAR)
30.79586.5165
50.08561.1483
70.00720.1183
90.00430.0902
110.00080.0164
130.00080.0226
150.00100.0340
Total harmonic loss0.89837.9464
Fundamental loss147.33 433.54
Total power loss 148.23441.49
Table 7. Fundamental r.m.s voltage, total r.m.s voltages, and THD % in IEEE−13 bus URDN.
Table 7. Fundamental r.m.s voltage, total r.m.s voltages, and THD % in IEEE−13 bus URDN.
BusPhaseS. NoFundamental r.m.s VoltageTotal r.m.s VoltageTHD %
650a11 0o10
b21 −120o10
c31 120o10
632a40.9498 −2.7462o0.95001.9173
b50.9839 −121.6817o0.98390.4974
c60.9300 117.8000o0.93022.2737
671a70.9109 −5.8987o0.91174.0623
b80.9875 −122.2091o0.98751.0363
c90.8717 115.9500o0.87284.9409
680a100.9109 −5.8987o0.91174.0623
b110.9875 −122.2091o0.98751.0363
c120.8717 115.9500o0.87284.9409
633a130.9466 −2.8223o0.94681.9098
b140.9819 −121.7315o0.98190.4919
c150.9271 117.7946o0.92732.2648
634a160.9207 −3.6073o0.92091.8801
b170.9624 −122.2445o0.96240.4873
c180.9064 117.2178o0.90662.2406
645b190.9745 −121.8646o0.97450.4991
c200.9283 117.8225o0.92862.2769
646b210.9729 −121.9382o0.97290.5000
c220.9264 117.8696o0.92672.2815
692a230.9109 −5.8987o0.91174.0623
b240.9875 −122.2091o0.98751.0363
c250.8717 115.9500o0.87284.9409
675a260.9025 −6.0795o0.90344.3128
b270.9887 −122.3037o0.98871.0491
c280.8678 116.0660o0.86895.0687
684a290.9093 −5.9502o0.91004.0765
c300.8684 115.9163o0.86955.0741
611c310.8651 115.8365o0.86635.2263
652a320.9041 −5.8755o0.90494.0900
Table 8. Fundamental and Harmonic power loss for IEEE−13 URDN with D-STATCOM.
Table 8. Fundamental and Harmonic power loss for IEEE−13 URDN with D-STATCOM.
Harmonic OrderHarmonic Power Loss
Active
(kW)
Reactive
(kVAR)
30.78366.3935
50.08411.1282
70.00710.1168
90.00430.0886
110.00080.0160
130.00080.0223
150.00090.0336
Total harmonic loss 0.88167.7991
Fundamental Loss135.34 396.63
Total power loss 136.22 404.43
Table 9. Fundamental r.m.s voltages, total r.m.s voltages, and THD % in IEEE−13 bus URDN with D-STATCOM.
Table 9. Fundamental r.m.s voltages, total r.m.s voltages, and THD % in IEEE−13 bus URDN with D-STATCOM.
BusPhaseFundamental r.m.s VoltageTotal r.m.s Voltage in p.uTHD %
650a1.0000 0o10
b1.0000 −120o10
c1.0000 120o10
632a0.9545 −2.8433o0.95471.8920
b0.9883 −121.7039o0.98830.4923
c0.9355 117.6975o0.93572.2404
671a0.9204 −6.0689o0.92113.9865
b0.9961 −122.2496o0.99621.0211
c0.8829 115.7512o0.88394.8346
680a0.9227 −6.1274o0.92343.9764
b0.9983 −122.2605o0.99831.0189
c0.8855 115.6856o0.88654.8205
633a0.9512 −2.9186o0.95141.8846
b0.9863 −121.7535o0.98630.4869
c0.9326 117.6923o0.93292.2317
634a0.9255 −3.6957o0.92561.8556
b0.9669 −122.2618o0.96690.4823
c0.9121 117.1225o0.91232.2082
645b0.9789 −121.8864o0.97890.4940
c0.9338 117.7205o0.93412.2436
646b0.9773 −121.9600o0.97730.4948
c0.9319 117.7674o0.93222.2482
692a0.9204 −6.0689o0.92113.9865
b0.9961 −122.2496o0.99621.0211
c0.8829 115.7512o0.88394.8346
675a0.9121 −6.2456o0.91294.2312
b0.9973 −122.3429o0.99731.0335
c0.8790 115.8645o0.88014.9586
684a0.9187 −6.1203o0.91944.0007
c0.8796 115.7184o0.88074.9647
611c0.8763 115.6396o0.87755.1133
652a0.9135 −6.0456o0.91424.0140
Table 10. Load composition of spot loads in IEEE−34 bus URDN.
Table 10. Load composition of spot loads in IEEE−34 bus URDN.
Bus No.Load Composition
Non-Linear LoadsLinear Loads
Fluorescent Light BanksAdjustable Speed DrivesComposite Residential Loads
830NoneNone80%20%
84430%30%30%10%
84830%30%30%10%
89030%None60%10%
86030%30%30%10%
84030%30%30%10%
Table 11. Case studies on IEEE−34 bus URDN.
Table 11. Case studies on IEEE−34 bus URDN.
Case StudyDescription
Case 1
(Without D-STATCOM)
  • Voltage regulator is taken off between bus−614 and bus−650 and bus−852 and bus−832.
  • Capacitor banks are taken off at bus−844 and bus−848.
  • Load composition of spot loads as in Table 9
Case 2
(With D-STATCOM)
  • Three-phase D-STATCOM placed at bus 890 is modeled as PV bus with P = 0 and three-phase reactive power limits: 100 ≤ Q ≤ 500 kVAR.
  • Three-phase D-STATCOM placed at bus 650 is modeled as PV bus with P = 0 and three-phase reactive power limits: 100 ≤ Q ≤ 500 kVAR.
Table 12. Fundamental r.m.s voltages, total r.m.s voltages, and THD % in the IEEE−34 bus URDN for Case 1.
Table 12. Fundamental r.m.s voltages, total r.m.s voltages, and THD % in the IEEE−34 bus URDN for Case 1.
Bus No.S. No.Ph.Total r.m.s
Voltages in p.u
THD %Bus NoS. No.Ph.Total r.m.sVoltages in p.uTHD %
8001a10 46a0.81694.8628
2b1083447b0.85314.0835
3c10 48c0.85704.5399
8024a0.99720.0541 49a0.81684.8691
5b0.99800.048184250b0.85304.0888
6c0.99810.0531 51c0.85694.5458
8067a0.99530.0908 52a0.81644.8990
8b0.99670.080884453b0.85244.1147
9c0.99690.0891 54c0.85654.5742
80810a0.96050.7955 55a0.81624.9094
11b0.97450.698384656b0.85194.1261
12c0.97420.7706 57c0.85624.5841
81213a0.92001.6825 58a0.81624.9109
14b0.94881.453184859b0.85184.1276
15c0.94791.6045 60c0.85624.5855
81416a0.88802.443581061b0.99750.0562
17b0.92842.081381862a0.88692.4466
18c0.92712.299482063a0.85992.5235
85019a0.88802.443882264a0.85642.5339
20b0.92842.081582665b0.91692.3333
21c0.92712.299785666b0.89542.8693
81622a0.88772.4526 67a0.78755.5150
23b0.92812.088888868b0.82484.8565
24c0.92682.3079 69c0.82755.1312
82425a0.87782.7468 70a0.78535.5791
26b0.91712.332889071b0.82284.9085
27c0.91732.5798 72c0.82525.1873
82828a0.87692.771386473a0.82034.7018
29b0.91622.352986074a0.81644.8734
30c0.91652.602575b0.85264.0921
83031a0.85743.382376c0.85664.5493
32b0.89612.8557 77a0.81624.8779
33c0.89743.167383678b0.85224.0967
85434a0.85693.3981 79c0.85654.5531
35b0.89562.8688 80a0.81624.8790
36c0.89693.182084081b0.85224.0976
85237a0.82334.5667 82c0.85644.5542
38b0.86033.8340 83a0.81624.8779
39c0.86344.262686284b0.85224.0968
83240a0.82334.5671 85c0.85654.5530
41b0.86033.834283886b0.85204.0977
42c0.86344.2629---
85843a0.82034.7018
44b0.85703.9475
45c0.86044.3890
Table 13. Summary of results for the case studies on IEEE−34 bus URDN.
Table 13. Summary of results for the case studies on IEEE−34 bus URDN.
Case StudyMin.
Fundamental Voltage, p.u
Min. Total r.m.sVoltage, p.uMax. THD%No. of Phases of Buses (THD > 5%)
Case 10.7841 at bus−890,
a-phase
0.7853 at bus−890, a-phase5.5791 at bus−890, a-phase4
Case 20.8137 at bus−890, a-phase0.8148 at bus−890,
a-phase
5.2567 at bus−890, a-phase2
Case studyTotal fundamental power lossTotal power loss including total harmonic loss
Active (kW)Reactive (kVAR)Active (kW)Reactive (kVAR)
Case 1260.89180.49264.56188.15
Case 2227.69 155.82231.23 163.21
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Satish, R.; Vaisakh, K.; Abdelaziz, A.Y.; El-Shahat, A. A Novel Three-Phase Harmonic Power Flow Algorithm for Unbalanced Radial Distribution Networks with the Presence of D-STATCOM Devices. Electronics 2021, 10, 2663. https://doi.org/10.3390/electronics10212663

AMA Style

Satish R, Vaisakh K, Abdelaziz AY, El-Shahat A. A Novel Three-Phase Harmonic Power Flow Algorithm for Unbalanced Radial Distribution Networks with the Presence of D-STATCOM Devices. Electronics. 2021; 10(21):2663. https://doi.org/10.3390/electronics10212663

Chicago/Turabian Style

Satish, Raavi, Kanchapogu Vaisakh, Almoataz Y. Abdelaziz, and Adel El-Shahat. 2021. "A Novel Three-Phase Harmonic Power Flow Algorithm for Unbalanced Radial Distribution Networks with the Presence of D-STATCOM Devices" Electronics 10, no. 21: 2663. https://doi.org/10.3390/electronics10212663

APA Style

Satish, R., Vaisakh, K., Abdelaziz, A. Y., & El-Shahat, A. (2021). A Novel Three-Phase Harmonic Power Flow Algorithm for Unbalanced Radial Distribution Networks with the Presence of D-STATCOM Devices. Electronics, 10(21), 2663. https://doi.org/10.3390/electronics10212663

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop