# Selection of Design Parameters to Reduce the Zero-Sequence Circulating Current Flow in Parallel Operation of DC Linked Multiple Shunt APF Units

Authors: Shafiuzzaman K. Khadem; Malabika Basu; Michael F. Conlon
Journal: Advances in Power Electronics
 (2013)

---

## Abstract

Capacity enhancement and operational flexibility are two of the important limitations of the centralized shunt APF (APFsh) unit. These limitations can be conquered by the operation of multiple APF units in parallel and connected back to back by a common DC link capacitor. In that case, a circulating current (CC) flows within the APFsh units. This CC flow becomes out of control when the APFsh units operate in hysteresis based current controlled mode. One of the difficulties of this CC flow control or reduction is the variable switching frequency of the APFsh units. In this paper, the model for CC flow is derived by the switching dynamics study of the APFsh units. It is found that the selection of design parameters plays an important role in the amount of CC flow. Detailed simulation, analysis, and real-time performance show how the selection of design parameters affects the CC flow and the reduction of CC flow can also be achieved at an acceptable level by the proper selection of design parameters.


---

## Body

## 1. Introduction

The power quality, at all time, is a matter of concern where a number of nonlinear, harmonics producing, and sophisticated loads are connected in an electrical distribution network. In that case, for power quality (PQ) improvement, shunt active power filters (APFsh) are finding greater applications as interfacing and compensating devices in the distributed network. The power rating and switching frequency of the APFsh converters are determined by the magnitude of harmonic currents and required filter bandwidth. In high power applications, the filtering task cannot be performed for the whole spectrum of harmonics by using a single converter due to the limitations on switching frequency and power rating of the semiconductor devices. Therefore, compensating the reactive harmonic components to improve the power quality of the distributed network as well as to avoid the large capacity centralised APFsh, parallel operation of multiple low power APFsh units is increasing. Different controlling mechanisms and topologies are available in handling the difficulties of parallel operation of APFsh either in active load sharing or distributed mode. Although the harmonic current compensation is the primary function of parallel APF, it can also be used as a compensator for voltage harmonics. A detailed technical review on parallel operation of APFsh for current and voltage harmonic compensation in a distributed generation (DG) network has been done in [1, 2] where the pros and cons of the different control methods have been discussed. For these cases, there is no physical/electrical link between the APFsh units. Another kind of configuration is also available where the system components and cost are reduced by maintaining a common DC link back-to-back connection between the APFsh units [3]. But it then raises the control complexity by introducing the zero-sequence circulating current (ZSCC) flow within the APFsh units which can overload the inverters of the APF units and thus increases the power losses. Other problems associated with the ZSCC flow is that it can reduce the ampacity (capacity of carrying current) of the cable, it can increase the DC link voltage, and it is harmful for the DG system like DFIG application [4–8].The control methods for the reduction of ZSCC discussed in the literature are mainly for PWM based control system. In most of the cases, the methods are discussed for the parallel operation of inverter or rectifier [6–10]. Very few of the articles show the ZSCC reduction for the parallel operation of APFsh units [11]. None of these discuss the ZSCC flow issues and control method for the parallel operation of converter/APFsh based on hysteresis control. Therefore, an attempt has been made here to discuss the issues related to hysteresis control.


## 2. Working Principle

According to the working principle of anAPFsh unit in parallel to the load to compensate the reactive and harmonic current, there are two possible modes of operation: capacitive mode (where the current flows from the capacitor) and inductive mode (the current flows towards the capacitor). When multiple APF units work in a current sharing mode, there could be four possible modes of operation: (i) capacitive-capacitive, (ii) inductive-inductive, (iii) capacitive-inductive, and (iv) inductive-capacitive. As an example, Figure 1(a) shows a single line diagram of a 3-phase system, where icc represents the circulating current flow between the APF units. In the case of a 3-phase system, shown in Figure 1(b), this circulating current flow exists as a zero-sequence harmonics in the zero-sequence current flow when a circulating loop is created within the APFsh units and hence it is termed zero-sequence circulating current (ZSCC) flow. In general these harmonics are (3+n*6) order, where n=0,1,2….(a) Single line diagram of a 3ph 2 unitsAPFsh and (b) VSI switching with common DC link presenting icc=ZSCC flow.
(a)
(b)


## 3. Model for ZSCC Flow

For simplicity, derivation of ZSCC flow has been carried out on a per phase basis. Figure2 shows the possible mode of operation between the APF units (1 phase), where sh-i and sh-c represent the corresponding inductive and capacitive mode, respectively. If there is a difference in any of the parameters, such as switching frequency, current sharing, interfacing inductor, and hysteresis band (in case of hysteresis current controller), then the APFsh units can operate in inductive-capacitive or capacitive-inductive mode. In these cases, ZSCC will flow and these are reflected in Figures 2(f) and 2(h).Figure 2
Possible mode of operation between two APF units to calculate circulating current.From Figures2(a) and 2(b), during the capacitive mode, the current flow can be obtained as
(1)dish-cdt=vsh-c-vpccLsh1∥Lsh2,
while vsh-c1=δ1vdc=vsh-c2=(1-δ2)vdc=vsh-c; here δ is the ducy cycle; that is, the Ac voltages of the APFsh units are equal.Similarly, Figures2(c) and 2(d) show that when both the APF units work in inductive mode during its switching cycle, the resulting current flows can be found as
(2)dish-idt=-vsh-i+vpccLsh1∥Lsh2,
while vsh-i1=(1-δ1)vdc=vsh-i2=δ2vdc=vsh-i.From Figures2(e) and 2(f), the circulating current flow in capacitive-inductive mode can be obtained as
(3)dicc-cdt=vcc-cLsh1+Lsh2=vsh-c1-vsh-i2Lsh1+Lsh2=(δ1-δ2)vdcLsh1+Lsh2,
while vsh-c1>vsh-i2.Similarly, the circulating current flow for inductive-capacitive mode is(4)dicc-idt=vcc-iLsh1+Lsh2=vsh-i1-vsh-c2Lsh1+Lsh2=-(δ1-δ2)vdcLsh1+Lsh2,
while vsh-i1>vsh-c2.Therefore, in general, the equation for the ZSCC flow can be written as(5)diccdt=δccvdcLsh1+Lsh2,
where δcc=δ1-δ2 for capacitive mode and δcc=δ2-δ1 for inductive mode.The possible maximum value for CC flow (Icc-max<Vdc/(Lsh1+Lsh2)) in one leg can be obtained as
(6)Icc-max=(S11ON~S12ON)·Vdc(Lsha1+Lsha2)=Δt·Vdc(Lsha1+Lsha2),
where S, S11ON, and S12ON are the switch and switching time, as shown in Figure 1, respectively. The possible value of Δt can be obtained as Δt=k·Tsw. Here, k=0,…,1 and Tsw is the time period.For 3-phase system, the ZSCC can be found as(7)ZSCC=ish1-zs~ish2-zs,
where ish1-zs and ish2-zs are the zero-sequence components of the APFsh units.This is also confirmed in [11] for APFsh topology analysis based on a parallel interleaved inverter. It is also similar to the rectifier analysis in [6]. In most of the research articles, derivation or ZSCC cancellation has been carried out for two inverters or APFsh units [7–11]. Moreover, the control methods discussed in the literature are mainly for PWM based control system. None of them discuss the CC flow issues for the parallel operation of APFs based on hysteresis control. Therefore, an attempt has been made here to discuss the issues related to hysteresis control.


## 4. Control Issues for the ZSCC Flow

There are mainly two ways to eliminate or reduce the ZSCC flow: (i) breaking the route of CC flow by introducing physical devices such as an isolation transformer or common mode inductor or (ii) by proper control methods. From the derivation of ZSCC flow it is found that the CC value depends on the duty cycle (δ) or switching frequency (fsw), DC link voltage (Vdc), and the interfacing inductor (Lsh) of the APF units. Therefore, before selecting the control method for multiple APFsh units with a common Dc link, some design consideration should be made to reduce the CC flow. The design options that can be implemented to reduce the CC flow are as follows.(i)
Decrease theVdc. According to the design criteria of single/3-phase APFsh system there is a minimum value of Vdc that has to be maintained to compensate the reactive and harmonic current.(ii)
As theδ=τ/T=τfsw, reducing the fsw and decreasing the difference between the fsw or δ will reduce CC flow.(iii)
Reducedfsw will further increase the Lsh (in the case of hysteresis current control) and it will then in turn reduce the circulating current flow. But increased Lsh will reduce the VAR rating of APFsh.It is also clear from the derivation that the CC flow does not directly depend on the compensating current rating of theAPFsh units. Therefore, once the current or VAR rating of APFsh is fixed, the design components (fsw and Lsh) should be chosen in such a way so that the Lsh is possibly high and the fsw difference between the APFsh units is as low as possible. In that case, introducing the common mode inductor in the APFsh units can be a better choice. Implementation of a CC control method in the APFsh control strategy will further reduce the CC flow. But it will then increase the control complexity.


## 5. Selection of Control Method forAPFsh Units

With the control issues of the CC flow in mind, selection of a control strategy for theAPFsh units should be made. Different control strategies for parallel operation of multiple APF units were reviewed [1, 2]. In general, control can be based on active current sharing where the compensating power/current can be divided equally or up to the compensating capacity of the APF units. Another prominent control is droop control where information exchange between the APF units is required. For simplicity, a current sharing (power splitting) method based on hysteresis current control has been chosen for the remainder of the analysis. The control of switches mainly depends on the control method for the APF units. In the case of the current sharing method, the switching can be determined by current sharing (equally distributed) or a capacity limitation technique. As an example, if Ish, Ish-ref, IAPFsh-rated, and Isw-ref are the Rms values of APF compensating current, reference compensating current, rated current, and the switching reference current for ON/OFF condition, thenIsh≤Isw-ref for the current sharing technique where Isw-ref=Ish-ref/N andIsh≤IAPFsh-rated for the capacity limitation technique.


## 6. Simulation Study

In-depth information on parallel operation of multipleAPFsh based on hysteresis current controller is not available yet, especially in the case of ZSCC flow. Therefore, in the present work, the objective is to reduce the ZSCC flow based on the design parameter selection. Two units of APFsh have been modelled here to operate in power sharing mode. In terms of sharing the following proportions have been chosen:(A)
IAPFsh1=IAPFsh2=0.5*Ilqh-max and Ilqh-max=100A; each APFsh compensates half of the load reactive and harmonic components,(B)
IAPFsh1=0.75*Ilqh-max; IAPFsh2=0.25*Ilqh-max and Ilqh-max=100A,(C)
IAPFsh1=0.6*Ilqh-max; IAPFsh2=0.4*Ilqh-max; and Ilqh-max=100A,(D)
IAPFsh1=0.6*Ilqh-max; IAPFsh2=0.4*Ilqh-max; and Ilqh-max=200A.In the case of a hysteresis current control basedAPFsh, according to the selection of design parameters for APFsh, described in [12–14] and the control issues of CC flow, the CC flow can be reduced by the proper selection of the design components including fsw, Lsh, h, and Vdc. As the Vdc is fixed, variation of other components has been made here to analyze the reduction of CC flow. Results and discussion are given below.Table1 shows the parameters that have been considered for the proportion (A), as mentioned above. Corresponding results have been given in Figure 3 and discussed below.Table 1
A:IAPFsh1 = IAPFsh2=0.5*Ilqh-max and Ilqh-max=100 A.
Frequency (kHz)
Cases
L
sh1
/
L
sh2 (mH)
h
1
/
h
2 (A)
THD@Ipcc (%)
f
sw1
=
17.3 fsw2=17.3
A1
0.625
/
0.625
5
/
5
2.57
A2
1.25
/
1.25
2.5
/
2.5
2.02
A3
1.25
/
0.625
2.5
/
5
1.74
f
sw1
=
17.3 fsw2=11.6
A4
1.25
/
1.25
2.5
/
3.75
1.67
A5
1.25
/
1.875
2.5
/
2.5
2.35
f
sw1
=
13.0 fsw2=11.6
A6
1.67
/
0.625
2.5
/
5
2.88
A7
1.67
/
1.25
2.5
/
3.75
2.82
A8
1.67
/
1.875
2.5
/
2.5
3.76Performance ofAPFsh for the cases A1–A5 and the corresponding CC flow.
(a)
(b)Figure3 shows the performance study of the two APFsh units in power sharing mode for the cases A1a–A5a where each of the units compensates half of the load reactive and harmonic components. APFsh1 starts working at the beginning, whereas APFsh2 starts at 0.5 sec. As there is no ZSCC that occurs before 0.5 sec, therefore the ZSCC waveforms for the corresponding simulations have been shown in A1b–A5b from 0.55 sec to 0.75 sec. For the cases A1–A3, the switching frequency of both units was equal and the variations have been made for Lsh and h. Cases A1 and A2 show that the possible CC flow can be zero when Lsh1=Lsh2 and h1=h2; that is, both the APFsh units are identical. Although the rating and fsw of each unit are the same, ZSCC can flow as shown in case A3. This happens due to the difference in design parameters Lsh and h. Cases A4, and A5 show the results for different fsw, Lsh, and h. For both the cases, ZSCC flows. In a comparison with A3, A4 and A5, A5 shows better results in the reduction of ZSCC flow. It is due to the higher values of Lsh and lower values of h, though the parameter fsw is different.Table2 shows the parameters that have been considered for the proportion (B). Corresponding results have been given in Figure 4 and discussed below.Table 2
B:IAPFsh1=0.75*Ilqh-max;IAPFsh2=0.25*Ilqh-max; and Ilqh-max=100 A.
Frequency (kHz)
Cases
L
sh1
/
L
sh2 (mH)
H
1
/
h
2 (A)
THD@Ipcc (%)
f
sw1
=
17.3 fsw2=17.3
B1
0.835
/
2.5
3.75
/
1.25
1.74
B2
0.418
/
1.25
7.5
/
2.5
1.81
B3
0.835
/
1.25
3.75
/
2.5
2.25
f
sw1
=
17.3 fsw2=11.6
B4
0.835
/
1.87
3.75
/
2.5
2.25
B5
0.835
/
3.74
3.75
/
1.25
2.33Performance ofAPFsh for the cases B1–B5 and the corresponding CC flow.
(a)
(b)Figure4 shows that, for all the cases, ZSCC flows. Even though the fsw values are different, the reduced ZSCC flow occurs in the case of B5 where the values for Lsh are high and h are low compared to other cases. The worst case occurs for B2 where Lsh is low and h is high.Tables3 and 4 show the parameters for the cases C and D where the sharing proportion is different than that for A and B. For both the cases, values of Lsh are kept low, whereas the values for h are high. In both cases, fsw are equal. And the same components have been chosen for D where the compensating current is two times higher than that of case C.Table 3
C:IAPFsh1=0.6*Ilqh-max;IAPFsh2=0.4*Ilqh-max; and Ilqh-max=100 A.
Frequency (kHz)
Cases
L
sh1
/
L
sh2 (mH)
H
1
/
h
2 (A)
THD@Ipcc (%)
f
sw1
=
17.3 fsw2=17.3
C1
0.625
/
0.625
5
/
5
2.65Table 4
D:IAPFsh1=0.6*Ilqh-max;IAPFsh2=0.4*Ilqh-max; and Ilqh-max=200 A.
Frequency (kHz)
Cases
L
sh1
/
L
sh2 (mH)
H
1
/
h
2 (A)
THD@Ipcc (%)
f
sw1
=
17.3 fsw2=17.3
D1
0.625
/
0.625
5
/
5
2.12Figure5 shows the performances for the cases C and D in terms of ZSCC flow. In both cases, ZSCC flow is the same even though for case D, APFsh units are compensating higher current than that of C. These results indicate that the ZSCC flow does not depend on the amount of compensating current. These ZSCC flows are higher than that of A and B and this is due to the lower values of Lsh.CC flow for the cases C and D.Finally, a comparison has been made for the ZSCC between the A3–A5 and A6–A8 cases for one cycle, as shown in Figure6. For A3, fsw is the same for both of the APFsh units, whereas for A4 and A5 it is different. But the values of Lsh are comparatively higher for A5. Therefore, within A3 to A5, the ZSCC flow is comparatively lower in the case of A5, as shown in Figure 6(a).ZSCC comparison between the cases (a) A3, A4, and A5 and cases (b) A6, A7, and A8.
(a)
(b)For the cases A6–A8, the switching frequency (fsw1) of APFsh1 is reduced and therefore the value of Lsh increases, as shown in Table 1. At the same time, the switching frequency difference between the APFsh units is also reduced. The performance of the APFsh units in terms of ZSCC flow for these cases is reflected in Figure 6(b). This shows that ZSCC flow can further be reduced by reducing the fsw (to increase the value of Lsh) and the difference between the fsw of the APFsh units. The reduction of zero-sequence harmonics content for the cases A6–A8 is also shown in Figure 7, where it shows that the magnitude of zero-sequence harmonics components in the compensating current is reduced gradually for the cases from A6 to A8. This also confirms the outcome as shown in Figure 6.FFT analysis of the ZSCC for the cases (a) A6, (b) A7, and (c) A8.
(a)
(b)
(c)But it also has to be kept in mind that increasingLsh could increase the THD at the PCC, as given in Table 1. Therefore, in terms of design criteria for the selection of APFsh units to place in multiple units with common DC link system, one should consider the higher values of Lsh and lower the difference between fsw of the APFsh units with care about THD.


## 7. Real-Time Performance Study

With the advancement of technology, real-time performance of any system can be observed using a real-time simulator. Instead of developing the complete actual system at full capacity, the controller/system can either be modelled in software or can be built in hardware or can be a combination of both. In real-time simulation, the accuracy of the computations depends upon the precise dynamic representation of the system and the processing time to produce the results [15].An electrical power system including 2 units ofAPFsh has been modelled in MATLAB using RT-LAB (real-time simulation) tools to observe the performance in the real-time environment. The system is then tested in SIL (Software-in-Loop; both the controller and the plant are simulated in real-time in the same simulator) with the hardware synchronization mode (data communication is in real-time through the ADC/DAC) which is similar to hardware-in-loop (HIL) test. Figure 8 shows the real-time simulation structure in a SIL configuration used to develop the real-time environment by OPAL-RT. With the combination of MATLAB SPS (Sim Power System) from MATHWORKS and the RT-LAB toolbox from OPAL-RT, the real-time model of the power system and controller is developed for the real-time simulator.Figure 8
Real-time simulation structure in SIL configuration with hardware synchronization mode.Figure9 shows the results for the cases A6, A7, and A8, where both the APFsh units are compensating harmonic current in a power sharing mode. The waveform of the utility current (is) contains fewer harmonics for A7 compared to A6 and A8. This happens because of the moderate value of Lsh and h for A7 compared to that of A6 (Lsh is too low) and A8 (Lsh is too high). Therefore, the THD value is also reduced.Performance of APFsh units for the cases (a) A6, (b) A7, and (c) A8.
(a)
(b)
(c)For the case of A8, although the THD increases, ZSCC is reduced. Figures10(a), 10(c), and 10(e) show the ZSCC flow for these (A6–A8) cases and it clearly depicts the gradual reduction of ZSCC from A6 to A8. It is also clear from the zoom-in for one cycle study in Figures 10(b), 10(d), and 10(f) that a considerable amount of ZSCC can be reduced by proper selection of design parameters for the APFsh units to operate in load sharing mode. Thus, it validates the simulation, analysis of the performance, and ZSCC flow study.ZSCC flow for the cases (a) A6, (c) A7, and (e) A8; the zoom-in for one cycle (b) A6, (d) A7, and (f) A8.
(a)
(b)
(c)
(d)
(e)
(f)


## 8. Conclusion

Capacity enhancement and operation flexibility are two of the important limitations of the centralizedAPFsh unit when placed in a DG integrated network, especially in the low voltage distribution level where the compensating current could be high. Therefore multiple APFsh units can operate in parallel. This can be controlled in an active load sharing mode (intercommunication is required) or in droop control mode. The most important part is the common DC link between the APF units. Therefore, a circulating current could flow. In the case of a hysteresis current controller based multiple APFsh units in load sharing mode, this CC control is not yet achieved. One of the difficulties of this CC flow control or reduction is the variable switching frequency of the APFsh units. By proper selection of design parameters, this CC flow can be reduced in an acceptable level.



---

*Source: 381581-2013-12-08.xml*