# Analytical Approach to Circulating Current Mitigation in Hexagram Converter-Based Grid-Connected Photovoltaic Systems Using Multiwinding Coupled Inductors.

1. IntroductionRecently, photovoltaic (PV) energy systems have gained more attention as distributed generation units, as they offer low cost of generation closer to that of conventional plants, as well as less maintenance and no grid noise [1, 2]. Moreover, PV systems can solve multiple typical problems present in conventional AC power systems. However, PV systems present frequency and voltage fluctuations when islanding operation occurs. Therefore, PV plants should be integrated with the power system in order to maintain the overall frequency and voltage at a stable condition. Several studies suggest interconnection methods of PV systems to the grid through voltage source converters (VSCs), because they provide versatile functions enhancing capabilities of the system [3-5]. The main purpose of the VSC is to connect the PV plant to the grid while guaranteeing power quality (PQ) standards. However, the high-frequency switching of VSCs introduces additional harmonic components to the system, hence creating PQ problems if not implemented accurately [6].

Most of the available VSCs for PV systems are traditional two-level three-phase VSCs with low power capacity [7]. In the literature, researchers have suggested numerous multilevel topologies for grid-connected PV plants [8-10]. Multilevel VSCs are considered more attractive than traditional two-level VSCs, as they improve the output voltage quality and reduce electromagnetic interference, voltage stress on IGBTs, and common-mode voltage. Moreover, multilevel VSCs operate at a low switching frequency and hence increase system efficiency [10]. Consequently, multilevel VSCs have been widely used in chemical, oil, and different kinds of plants, as well as in power plants, transmission systems, and PQ compensators [11].

Current research is mainly focused on three specific multilevel VSC topologies, namely, the neutral-point clamped (NPC) [12], flying capacitor (FC) [13] and the Cascaded H-bridge (CHB) [14] topologies. The NPC VSC requires a high number of clamping diodes to increase the number of voltage levels, which can cause problems related to switches with different ratings, high-voltage rating in blocking diodes, and capacitor voltage imbalance. However, the FC VSC consists of a large number of capacitors, which leads to complications in regulating capacitor voltages. Finally, the CHB VSC provides isolated DC sources, hence being appropriate for use in PV systems [15]. In addition, other advantages of this topology include its modular structure and the reduced number of components compared to the other multilevel converters (i.e., NPC and FC). Therefore, the CHB VSC can reach the same number of voltage levels with a simpler assembly and maintenance. However, this topology presents three main drawbacks: (i) high overall component count; (ii) high energy storage requirement because the instantaneous power related to each H-bridge varies at twice the fundamental frequency, given its single-phase modular structure; and (iii) difficulty to control the voltage across DC-link capacitors.

The concept of interconnecting three traditional three-phase VSCs to produce a multilevel converter was first proposed in [16] and applied to medium-voltage variable-speed drives. In [17], three three-phase two-level VSCs are interconnected using three single-phase transformers with a 1 : 1 turn ratio. These interconnected transformers increase the output voltage and suppress the circulating current inside the converter. The power capacity of the overall converter is three times the capacity of each interconnected converter, whereas the volt-ampere rating of each intermediate transformer is equal to that of each interconnected converter. In [18], another topology known as the hexagram multilevel converter (HMC) is proposed, which combines six two-level converters by using six inductors. This topology shares many advantages of the CHB but uses fewer switches and reduces the size of the DC-link capacitor [19]. The advantages of the proposed topology can be summarized as follows: (1) only six standard three-phase VSCs are necessary to generate multilevel output voltage; (2) each VSC module is balanced in operation, equally loaded, and supplies 1/6 output power; (3) modular construction which facilitates system maintenance and spare management; (4) only six isolated DC links with no voltage unbalance problem; (5) low number of power electronics switches and low DC energy storage requirement; and (6) the output transformer contributes to higher output voltage. Figures 1 and 2 are schematic diagrams illustrating multilevel topologies based on cascaded two-level VSCs according to previous research.

However, selecting an inductor with suitable inductance is no simple process, where the inductance value should be adequate to minimize the circulating current as well as to reduce an extra voltage drop on the inductors that affects the terminal voltage. This difficulty can be circumvented by using a multiwinding coupled inductor. Therefore, this paper analyzes the utilization of a multiwinding (e.g., two, three, and six windings) coupled inductor to interconnect the six VSCs instead of six single inductors. Consequently, both goals minimize the circulating current and the minimal effects to the output voltage can be accomplished instantaneously. Then, an analytical model to calculate the total inductance imposed to the circulating current path is derived. The equivalent circuit model of the HMC is detailed in the abc reference frame and then transformed into the orthogonal dq0 reference frame. Moreover, to extract the maximum power from the PV arrays, a control algorithm for maximum power point tracking is also presented.

The paper is organized as follows: Section 2 describes the HMC. Section 3 presents the mathematical model for a grid-connected PV system. Section 4 demonstrates the benefits of the multiwinding coupled inductor inside the HMC. Section 5 describes the control system and modulation strategy. Section 6 shows the simulation results followed by the corresponding discussion. Finally, the conclusions achieved from the present work are summarized in Section 7.

2. System Description

The proposed configuration of the three-phase grid-connected PV plants is shown in Figure 3. This configuration is composed of six traditional three-phase two-level VSC modules, which have a hexagonal interconnection to produce higher voltage levels as shown in Figure 2. The circulating current in the obtained loop can be suppressed using six inductors, which present a small impedance at the switching frequency. Each VSC module supplies 1/6 of the converter output power. In addition, each module has one of its AC terminals (i.e., a, b, or c) designated as converter output. The three-phase AC output terminals of the HMC are labeled as follows: A (AC terminal a1 of module 1), B (AC terminal b3 of module 3), and C (AC terminal c5 of module 5). The remaining three-phase AC output terminals are labeled as follows: A' (AC terminal a4 of module 4), B' (AC terminal b6 of module 6), and C' (AC terminal c2 of module 2). The other two AC terminals of each module are, respectively, connected to an adjacent module through an inductor. For instance, AC terminal b1 of module 1 is coupled to AC terminal b2 of module 2 through an inductor [L.sub.d]; AC terminal c1 of module 1 is coupled to AC terminal c6 of module 6 through another inductor. The remaining modules consist of similar connections, as shown in Figure 2.

As shown in Figure 3, the PV system based on the HMC has six output terminals. Therefore, when linked to the three-phase electrical grid, an open-end winding (OEW) transformer is necessary to provide these six terminals. The secondary windings of the OEW transformer are connected differentially. In the proposed converter, AC terminals A-A', B-B', and C-C' are used to provide phases A, B, and C, respectively. High-voltage windings are arranged in a star configuration and coupled directly to the three-phase grid. Each VSC module is supplied with a separate PV string to produce two-level individual outputs. The proposed converter topology presents several advantages, such as the use of fewer power electronic switches, diodes, and storage capacitors when compared with other topologies, making it suitable for renewable energy and other applications.

3. System Modeling

3.1. PV System Modeling and MPPT Control. Figure 4 shows the equivalent model of a PV array. The PV array is composed of several series-parallel connected solar cells. The basic equation of the PV array is given by [20].

I = [I.sub.pv] - [I.sub.0] [exp (V + [R.sub.s]I/a[V.sub.t])] - V + [R.sub.s]I/[R.sub.p], (1)

where [I.sub.pv] = [Np.sub.][I.sub.pv,cell] is the PV current of the array, [I.sub.pv,cell] is the current generated by incident light (directly proportional to Sun irradiation), [N.sub.p] is the number of cells connected in parallel, [I.sub.0] = [N.sub.p][I.sub.0,cell] is the saturation current of the array, [I.sub.0,cell] is the reverse saturation or leakage current of the diode, [V.sub.t] = [N.sub.s]kT/q is the thermal voltage of the array, [N.sub.s] is the number of cells connected in series, q is the electron charge, k is the Boltzmann constant, T is the temperature of the p-n junction, a is the diode ideality constant, [R.sub.s] is the equivalent series resistance of the array, and [R.sub.p] is the equivalent parallel resistance. Table 1 shows the nominal parameters of the KC200GT PV array. Figure 5 shows the I-V and P-V curves obtained using (1) of the PV array under variable solar radiation.

3.2. Maximum Power Point Tracking. The proposed system consists of 12 PV modules distributed as six pairs of series-connected modules and coupled directly to a single DC link (see Figure 3). Therefore, the power rating of the proposed converter is 2.4 kW. The controller is aimed to regulate the voltage of the PV-module pair at 52.6 V, to guarantee maximum power transfer to the grid.

Next, the aim of the MPPT control algorithm used in this paper is to ensure that under any solar radiation and temperature conditions, the maximum power is extracted from the PV modules. This is achieved by matching the PV-array maximum power point to the corresponding operating voltage and current of the converters. The perturb and observe (P&O) algorithm is a commonly used MPPT technique because it is easy to implement [21]. The operating principle of the P&O algorithm is shown in Figure 6. This algorithm measures the PV plant voltage and current, then it varies the operating voltage and compares the power received between the two voltage values. After each perturbation, the algorithm compares the output power from the PV before and after the perturbation. The direction of a new perturbation depends upon the output power: the perturbation will follow the same direction if higher power is measured when the output voltage varies and the opposite direction otherwise. These procedures are repeated continuously, and the reference voltage is generated and fed to the converter controller.

The numerical illustration of the P&O algorithm [2] is given as follows:

[dP.sub.pv](k)/[dV.sub.pv](k) = [P.sub.pv](k) - [P.sub.pv](k - 1)/[V.sub.pv](k)

- [V.sub.pv](k - 1). (1)

Here, [P.sub.pv](k) and [P.sub.pv](k - 1) stand for current power and previous measured power, while [V.sub.pv](k) and [V.sub.pv](k - 1) stand for current PV voltage and previous one.

3.3. Voltage and Current Analysis. Under the symmetrical operation condition, the fundamental components of the six VSC modules are the same. The corresponding open-circuit voltage of the HMC is presented in Figure 7.

3.3.1. Current Relations. The phase currents of each VSC modules fulfill the following expressions:

[mathematical expression not reproducible]. (3)

Assuming that the intermediate coupled inductors have large magnetizing inductance, the circulating currents are suppressed to a low value and can be ignored, hence

[i.sub.b1] + [i.sub.a2] + [i.sub.c3] + [i.sub.b4] + [i.sub.a5] + [i.sub.c6] = 0 (4)

Because every pair of the six VSC modules is coupled, the current within the proposed converter has the following expressions:

[mathematical expression not reproducible]. (5)

Next, suppose that the HMC is linked to a three-phase grid. Then, the output currents will satisfy the following equation:

[mathematical expression not reproducible]. (6)

Using (2), (3), and (4), the output current of each VSC module can be shown to be

[mathematical expression not reproducible], (7)

where [[[I.sub.A], [I.sub.B], [I.sub.C]].sup.T] and I are the converter output phase currents and their RMS values, respectively. It can be concluded from (5) that each VSC module within the HMC will have an identical current under the symmetrical operation condition.

3.3.2. Voltage Relations. Under the symmetrical operation condition, the fundamental component (RMS) of the phase voltages at the VSC modules can be described as

[mathematical expression not reproducible], (8)

V = 1/2[square root of 2] [V.sub,dc][m.sub.a], (9)

where V, [V.sub.dc], and ma are the RMS values of each VSC-module phase voltage, the DC-link voltage of each VSC module, and the amplitude modulation index, respectively; subscripts o1, o2, o3, o4, o5, and o6 represent the virtual neutral points of each VSC.

The output voltages of the proposed converter are written as

[mathematical expression not reproducible]. (10)

where impedance [Z.sub.d] = j[omega][L.sub.d].

Thus, using (5), (6), (7), and (8), the net output voltages of the HMC under the symmetrical operation condition are given by

[mathematical expression not reproducible]. (11)

Consequently, (9) demonstrates that the three-phase voltage of the HMC is approximately six times the phase voltage of each VSC module. In other words, the voltage stress is reduced by a factor of six.

3.4. Equivalent Circuit Model. Based on Figure 3, the line voltages of the HMC are expressed by

[mathematical expression not reproducible], (12)

where [[[E.sub.AB], [E.sub.BC], [E.sub.CA]].sup.T] are the RMS values of the line-to-line voltages in the grid [[[U.sub.AB] [U.sub.BC] [U.sub.CA]].sup.T] and [[[I.sub.A], [I.sub.B], [I.sub.C]].sup.T] are the RMS values of line voltages and the RMS values of the output current, respectively, and [L.sub.f] is the AC-side filtering inductance.

When the conversion powers between the VSC modules within the HMC are balanced, the six DC-link voltages are equal to [V.sub.DCav], hence

[mathematical expression not reproducible]. (13)

Moreover, the equivalent DC-link voltage can be expressed as

[U.sub.DC_eq] = [6V.sub.DC_av]. (14)

If each VSC is driven using the same PWM, the relationship between the DC-link voltages and the phase voltages can be expressed as

V = 1/2[square root of 2]] [U.sub.DC-eq] [m.sub.a]. (15)

From (12) and (14), the HMC can be modeled by

[mathematical expression not reproducible]. (16)

where

[mathematical expression not reproducible]. (17)

Using (14) and (17), the HMC can be modeled as a conventional three-phase two-level VSC, as shown in Figure 8.

Figure 8 shows the equivalent circuit model of the HMC. The equivalent circuit is composed of an open-circuit voltage source linked in series to the HMC output impedance. The open-circuit voltage of the HMC is shown in Figure 7. As shown in Figure 8, the equivalent interface inductor of the HMC is 2[L.sub.d] + [L.sub.f]. This feature is advantageous, because the filtering inductor Lf could be minimalized or even removed. However, this can be inopportune when the line transformer has sufficient inductance for filtering the output current. In this case, essentially if the HMC supplies reactive power, the output impedance causes a voltage drop, decreasing the power capability. Thus, it might be better to minimize the interconnected total inductances. However, the inductances are essential to suppress the circulating current within the HMC. The aforementioned decision shows that the interconnected inductance value selection is not a simple process, where the inductance should be large enough to minimize the circulating current but small enough to avoid an extra voltage drop. In the following section, this feature is investigated using multiwinding coupled inductors.

4. The Role of the Coupled Inductors

The function of the multiwinding coupled inductors within the HMC is investigated in this section. The HMC performance analysis can be demonstrated based on the function of the coupled inductors. The six VSC modules are interconnected to each other, creating a hexagon as explained in Section 3.4. However, instead of using six inductors, multi-winding coupled inductors are used to interconnect the VSC modules.

4.1. HMC with Three Two-Winding Coupled Inductors. Two inductors with an equal number of turns are coupled together; that is, the input to one side will produce an output on both, as shown in Figure 9. Since the turn ratio of the coupled inductor is approximately 1 : 1, the self-inductances of the primary and secondary windings are the same (e.g., [L.sub.p] = [L.sub.s] = [L.sub.B]). Applying the voltage relationships of coupled inductors, the following equations are given:

[V.sub.b1b2] = [L.sub.11] [di.sub.b1]/dt + [L.sub.12] [di.sub.b5]/dt, (18)

and

[V.sub.b5b4] = [L.sub.22] [di.sub.b5]/dt + [L.sub.21] [di.sub.b1]/dt, (19)

where [L.sub.22] = [L.sub.11] = [L.sub.B] + [L.sub.m], [L.sub.m] = [L.sub.21] = [L.sub.12] = [i.sub.b1] + [i.sub.cir] and [i.sub.b5] = -[i.sub.b1] + [i.sub.cir], and where [i.sub.x] (x = a, b, c) is the line current, Lb is the leakage inductance of the inductor windings, which is expected to be identical for all windings, and [L.sub.m] is the magnetizing inductance of the coupled inductor.

It is obvious in Figure 10 that the voltage drop on one of the coupled inductors is

v = [v.sub.b1b2] + [v.sub.b5b4]. (20)

Using (18) and (19), the voltage across the coupled inductor can be expressed as

[mathematical expression not reproducible], (21)

Therefore, the total voltage drop across the three coupled inductors can be expressed as

[V.sub.total] - 3(2[L.sub.B] + 4[L.sub.m]) d([i.sub.cir])/dt = (6[L.sub.b] + 12[L.sub.m]) d([i.sub.cir])/dt. (22)

4.2. HMC with Two Three-Winding Coupled Inductors. Three inductors with an equal number of turns are coupled together as shown in Figure 11. Writing the voltage equations of one coupled inductor with winding inductance Z1, the following equations are given:

[mathematical expression not reproducible], (23)

where

[L.sub.a] - [L.sub.b] - [L.sub.c] = L + [L.sub.m], [L.sub.m] = [L.sub.ab] = [L.sub.ba] = [L.sub.bc] = [L.sub.cb] = [L.sub.ca] = [L.sub.ac], (24)

and

[i.sub.a] - [i.sub.a5] + [i.sub.cir],

[i.sub.b] - [i.sub.b1] + [i.sub.cir],

[i.sub.c] - [i.sub.c3] + [i.sub.cir], (25)

where [i.sub.x] (x - a, b, c) is the line current, L is the leakage inductance of the inductor windings, and [L.sub.m] is the magnetizing inductance. It is obvious in Figure 12 that the voltage drop on the coupled inductor is

[mathematical expression not reproducible]. (26)

Using (23), the voltage drop across one of the coupled inductor can be expressed as

[mathematical expression not reproducible]. (27)

The voltage drop across the coupled inductor windings is

v = (L + 3[L.sub.m] d([i.sub.ab] + [i.sub.b1] + [i.sub.c3])/d + (3L + 9[L.sub.m] d([i.sub.cir])/dt. (28)

Under balanced conditions,

[i.sub.a5] + [i.sub.b1] + [i.sub.c3] - [i.sub.a] + [i.sub.b] + [i.sub.c] = 0. (29)

Accordingly, by applying (28) and (29), the voltage drop on the coupled inductor windings is

v = (3L + 9[L.sub.m]d([i.sub.cir])/dt. (30)

The total voltage drop across the two coupled inductors can be expressed as

[v.sub.total] = (6[L.sub.B] + 18[L.sub.m] d([i.sub.cir])/dt. (31)

4.3. HMC with One Six-Winding Coupled Inductor. Instead of using six inductances, one six-winding coupled inductor can be used to interconnect the VSC modules. Since the turn ratio of the coupled inductor is approximately 1 : 1, the self-inductances of all the windings are the same as shown in Figure 13. The HMC with one six-winding coupled inductor is shown in Figure 14. Applying the voltage equations of the coupled inductor, the following equations are given:

[mathematical expression not reproducible]. (32)

Since,

[mathematical expression not reproducible]. (33)

The voltage drop equations across each coupled inductor can be expressed as

[mathematical expression not reproducible], (34)

where

[L.sub.a] = [L.sub.b] = [L.sub.c] = L + [L.sub.m], (35)

and

[i.sub.a] = [i.sub.a] + [i.sub.cir],

[i.sub.b] = [i.sub.b] + [i.sub.cir],

[i.sub.c] = [i.sub.c] + [i.sub.cir]. (36)

Using (34), the voltage drop across the coupled inductor can be expressed as

v = [v.sub.a5a6] + [v.sub.b1b2] + [v.sub.c3c4] + [v.sub.a3a2] + [v.sub.b5b4] + [v.sub.c1c6], (37)

v = 2[v.sub.a5a6] + 2[v.sub.b1b2] + 2[v.sub.c3c4]. (38)

v = (6L + 36[L.sub.m]) [di.sub.cir]/dt + (2L + 12[L.sub.m]) d([i.sub.a] + [i.sub.b] + [i.sub.c])/dt. (39)

Under balanced conditions,

[i.sub.a] + [i.sub.b] + [i.sub.c] = 0. (40)

Consequently, by applying (39) and (40), the voltage drop on the coupled inductor windings is

v = (6L + 36[L.sub.m]) [di.sub.cir]/dt. (41)

It is considered that the magnetizing inductance is much higher than the leakage inductance. Therefore, neglecting the leakage inductance, the two-winding coupled inductor imposes twelve times the magnetizing inductance for the circulating current, while the impedance to the circulating current using three-winding coupled inductors is eighteen as much as the magnetizing inductance. Thanks to the six-winding coupled inductor, the impedance on the circulating current path is thirty six times the magnetizing inductance.

In a general way if k is the number of the coupled inductor windings, the impedance to the circulating current can be expressed as

[Z.sub.cir] = 2[pi]f(6k[L.sub.m]). (42)

5. Control Scheme and Modulation Strategy

5.1. Control Scheme. The main function of the controller is to generate reference currents such that the proposed converter only provides available active power from the DC links to the grid at the point of common coupling (PCC) [22, 23]. Using the equivalent circuit model presented in Figure 7 and applying Kirchhoff s voltage and current laws at the PCC, the following two equations in the abc frame can be obtained:

[mathematical expression not reproducible]. (43)

[mathematical expression not reproducible]. (44)

where [L.sub.fd] = [L.sub.f] + [L.sub.d] is the equivalent interface inductor and [S.sub.A], [S.sub.B], and [S.sub.C] represent the switching states of the equivalent circuit model under balanced conditions.

Assuming that the voltages are balanced and the zero-sequence component is zero, the voltage between the neutral virtual point of equivalent circuit model (O) and the grid neutral point (O') is given by

[U.sub.O'O] = -([bar.[U.sub.AO]] + [bar.[U.sub.BO]] + [bar.[U.sub.CO]])/3, (45)

[mathematical expression not reproducible] (46)

Substituting (46) into (43) and (44), the following relation is obtained: [mathematical expression not reproducible]. (47)

The dynamic model in the abc frame of the HMC equivalent circuit is represented by (47). The switching state functions, [d.sub.i] (i = A, B, C), are defined as

[mathematical expression not reproducible]. (48)

The dynamic model of the equivalent circuit model in the abc frame is achieved by combining (47) and (48) in the following equation:

[mathematical expression not reproducible]. (49)

The DC side differential equation can be written as

[mathematical expression not reproducible]. (50)

d[U.sub.DC_aw]/dt = 1[C.sub.eq] (2[d.sub.A] + [d.sub.B]) [I.sub.A] + 1/[C.sub.eq] ([d.sub.A] + 2[d.sub.B]) [I.sub.B]. (51)

It can be seen that the model represented by (49) and (51) is time varying. Thus, to facilitate the control algorithm implementation, the model can be expressed in the synchronous reference frame rotating at constant frequency w. The corresponding conversion matrix is

[mathematical expression not reproducible], (52)

where [theta] - [omega]t.

Applying the coordinate transformation to (49) we obtain

[mathematical expression not reproducible]. (53)

Similarly, applying this transformation to (50) we obtain

[C.sub.eq] = d[U.sub.DC_av]/dt = [d.sub.d][I.sub.d] + [d.sub.q][I.sub.q]). (54)

The obtained model, represented by (53) and (54), is nonlinear owing to the product between the state variables (i.e., [d.sub.d][I.sub.d], and [U.sub.DC_av]) and the inputs (i.e., [d.sub.d] and [d.sub.q]).

Figure 11 shows the control principle of the HMC. Because the proposed configuration is a three-wire system, only two phase currents are required to be measured.

Currents [I.sub.A] and [I.sub.B] are measured and converted to the dq0 frame to obtain the corresponding currents [I.sub.d] and [I.sub.q]. Accordingly, (53) is rewritten as follows:

[mathematical expression not reproducible]. (55)

Demonstrated that currents [I.sub.d] and [I.sub.q] can be controlled separately by acting on inputs [u.sub.d] and [u.sub.q], respectively. Hence, the controller is designed using the following expressions:

[mathematical expression not reproducible], (56)

where [[??].sub.d] = [I.sub.d_REF] - [I.sub.d] and [[??].sub.q] = [I.sub.q_REF] - [I.sub.q] are the current error signals, with [I.sub.d_REF] and [I.sub.q_REF] being the reference values for currents [I.sub.d] and [I.sub.q], respectively.

Using (55), the current control law is given by the following expression:

[mathematical expression not reproducible]. (57)

The d-axis reference current ([I.sub.d_REF]) is produced using the DC-link voltage controller, while the q-axis reference current ([I.sub.q_REF]) is taken as the load q-axis current ([I.sub.Lq]). To aim for a unity power factor, the [I.sub.q_REF] is set to zero (i.e., [i.sub.q_REF] = 0). The active power exchange between the DC links and the grid is proportional to the direct-axis current Id and can be expressed as

[P.sub.dc] = 3/2 ([E.sub.d][I.sub.d] + [E.sub.q][I.sub.q]) = 3/2 [E.sub.d][I.sub.d]. (58)

(58) it is shown that direct-axis current Id is responsible for maintaining the DC-link voltages at a desired value. Thus, using (54), one deduces that

[C.sub.eq] d[U.sub.DC_av]/dt = [d.sub.d][I.sub.d] = [u.sub.dc]. (59)

Consequently, the active current is

[I.sub.d] = [u.sub.dc]/[d.sub.d] = [u.sub.dc][V.sub.DC_av]/[d.sub.d][V.sub.DC_av]. (60)

Assuming that the current loop is perfect and the HMC works under balanced conditions, the following expressions hold:

[mathematical expression not reproducible], (61)

where [??] and [E.sub.d] are the RMS voltage and the direct-axis phase voltage at the PCC, respectively. Thus, the control effort of the DC link voltage loop is given by

[I.sub.d_REF] = [u.sub.dc]/[d.sub.d] = [square root of 2/3] [V.sub.DC_av]/[??] [u.sub.dc]. (62)

To control the DC-link voltage, a PI controller is used, which is expressed as

[mathematical expression not reproducible], (63)

where [[??].sub.dc] = [V.sub.DC_REF]-[V.sub.DC-av] is the DC-link voltage error, whereas [V.sub.DC_REF] and [V.sub.DC-av] are the DC-link reference and average voltages, respectively.

5.2. DC-Link Voltage Controller. With reference to Figure 3, the HMC is based on a symmetric configuration, having six converters with identical power capabilities that are supplied by six equal PV strings. The PV strings are directly connected to each converter, [V.sub.pv1] = [V.sub.pv2] = [V.sub.pv3] = [V.sub.pv4] = [V.sub.pv5] = [V.sub.pv6] = [V.sub.pv], the MPPT must be achieved by the converters, and the DC-link voltages continuously fluctuate. Because the PV strings are supposed to be identical, being created by a single PV string divided into six identical parts, a single MPPT algorithm can be considered. For this reason, the same DC-link voltage reference for the six converters has been considered. The DC-link voltage reference is compared to the sum of the actual six DC-link voltages, and the error is passed through a PI controller to determine the control parameter [u.sub.dc].

5.3. Modulation Strategy. To obtain an output voltage with low total harmonic distortion (THD), a multicarrier phase-shifted PWM (PS-PWM) switching strategy [24, 25] is implemented to drive each IGBT in the HMC. Optimum harmonic cancellation is accomplished by shifting each carrier cell by 2[pi][T.sub.s]/3T in sequence, where [T.sub.s] is the switching time and T is the cycle modulation time. A reference signal of 50 Hz is generated using the control algorithm represented in Figure 15, with the switching frequency fixed at 2500 Hz.

Figure 16 shows the relationship between the modulation waveforms and the three groups of carriers within the HMC. As shown in Figure 16, triangular carriers (i.e., carrier #1, carrier #2, and carrier #3) are phase-shifted 120[degrees] to each other and directly compared with the modulation signals to drive the IGBTs within module #1, module #2, and module #3. In order to generate the switching signals used to drive the IGBTs within module #2, module #4, and module #6, the modulation signals are inverted and then compared with the triangular carriers.

6. Simulation Results

In order to demonstrate the performance of the HMC and its control algorithm, the complete grid-connected PV system was simulated using the MATLAB/Simulink environment. The technical characteristics of the system parameters and the coupled inductors are provided in Table 2. In order to validate the performance of the HMC, a PS-PWM technique has been implemented, as shown in Figure 16. The output active and reactive power supplies in response to the fluctuations in solar radiation value are shown in the following subsections.

6.1. Performance Analysis. The seventeen-level phase voltages of the HMC are generated at the steady state and shown in Figure 17. It can be seen that the voltages are balanced, as shown in the voltage phasor diagram. Figure 18 demonstrates the currents of modules 1, 3, and 5, and the currents of modules 2, 4, and 6. Given the configuration of the three-phase HMC, the line currents of every VSC module within the HMC are symmetrical, hence verifying the relationship in (5).

From this figure, it can also be acknowledged that the direction of the currents of modules 2, 4, and 6 is reversed from those of modules 1, 3, and 5. Moreover, the currents inside every VSC module within the HMC are the same as the output currents.

The THD values of the output voltage and current are calculated using the following equation.

%[THD.sub.x] = [square root of [summation over (h[not equal to]1)] [([x.sub.sh]/[x.sub.s1]).sup.2], (64)

where the subscript x indicates the THD in the signal (voltage or current), [x.sub.s1] is the fundamental component, [x.sub.sh] is the component at the h harmonic frequency.

The harmonic spectra of the output voltage and current are shown in Figures 19 and 20, respectively. The THD of the HMC current is 0.68%, which is fewer than 5% and meets the power quality standard. Using the suggested modulation technique, the highest harmonic family of the phase voltage appears at the band of the 100th harmonic order. Consequently, the effective switching frequency of the phase voltage is two times higher than the switching frequency.

As a comparison, the HMC using two-, three-, and six-winding coupled inductors are simulated under the conditions. The magnetizing inductance of each coupled inductor is 3.5 mH. According to the current from the equations provided in Section 4, if the voltage of the DC links is unbalanced the circulating current will be introduced on the converter output currents. Therefore, to intentionally produce a circulating current, the DC link voltage of Module 3 is decreased from 52.6 to 26.3 V. The created loop voltage is computed using the following equation:

[V.sub.loop,rms] = [square root of 3]/2[square root of 2] (52.6 - 26.3)[m.sub.a]. (65)

Thus, using (42) and (65), the circulating current inside the HMC is expressed as

[I.sub.cir,peak] = [square root of 3]/24[pi] (52.6 - 26.3) [m.sub.a]/fk[L.sub.m]. (66)

Investigations with different magnetizing inductance levels have been carried. Figure 21 shows the theoretical values of the peak circulating currents with different magnetizing inductance levels. Figure 22 shows the simulation results with a relatively large magnetizing inductance (3.5 mH). The circulating currents are found by measuring the difference between the output currents of module 1 and module 3. As shown in the waveforms, the circulating currents inside the HMC with two-, three-, and six-winding coupled inductors are 2.3, 1.5, and 0.76 A peaks as calculated in (66). From the waveforms in Figure 22, using a six-winding coupled inductor, the circulating current is efficiently minimized.

In order to achieve the same circulating currents (e.g., 0.5 A), the magnetizing inductance of the two-winding, three-winding and six-winding coupled inductors should be increased to 17.5, 11.6, and 6mH, respectively. As shown in Figure 23, the waveform matches the theoretical analysis in Figure 21. The output current of module 1, module 3, and the circulating current inside the HMC with two-winding, three-winding, and six-winding coupled inductors under an unbalanced DC link voltage are shown in Figures 24, 25, and 26, respectively. The waveforms show the simulation result with an extremely small magnetizing inductance of 0.5 mH, which is 1/7 of that in Figure 22. The simulations validated the quantitative relationship that the circulating current will be too big to maintain a normal operation if the coupled inductors with low magnetizing inductance are used. As shown, the current ia3 nearly doubles the circulating current, which evidently proves that the HMC is sensitive to an unbalanced DC voltage. The waveforms match the theoretical analysis in Figure 21.

6.2. Active Power Variation. The system is tested for different solar radiation conditions and the results are shown in Table 3. Figure 27 shows the solar radiation variation which is considered in this study. The output current and grid voltage variation in response to the changes in solar radiation is shown in Figure 28. The transient behavior of the total DC-link voltage is presented in Figure 29. The parameters are successfully attuned by the controller to keep the DC-link voltage at the desired level of 52.6 V. Fluctuations are observed in the DC-links, due to the step variations in the solar radiation. Nevertheless, the controller brings the voltage to the reference level within 0.02 s.

The active power supplied by the HMC is directly proportional to the magnitude of direct axis, [i.sub.d], as indicated by (38). The solar radiation at all DC links is reduced by 40% at 0.25 s. The fluctuation in solar radiation is imitated through the reduction in the direct-axis current by approximately 40%, in a step. It decreased to 5.216 A from the initial value of 8.823 A, as shown in Figure 30. This reduction in id was to keep the DC-link voltage at the reference level by decreasing the output power drawn from the PV. Moreover, to guarantee maximum utilization of the PV system, the quadrature axis current was kept at zero. At 0.45 s, the solar radiation is further reduced to 400 W/[m.sup.2]. The direct-axis current is reduced to 3.267 A from the original value of 5.216 A. Later, the direct-axis current is increased to 8.823 A, because of the increment in the solar radiation.

6.3. Reactive Power Compensation. The described control algorithm permits the HMC to act as DSTATCOM in the absence of solar radiation. The output reactive power can be calculated as

[Q.sub.out] = - 3/2 [I.sub.q][V.sub.d]. (67)

In this condition, the reactive power is increased by 2700 VAR in a step, in the absence of solar radiation. The influence, for the DSTATCOM mode operation, is presented in Figure 31.

The direct-axis current is found to be zero, indicating the fact that no active power is being transferred to the grid in the absence of solar radiation. However, due to the step change in reactive power, the capacitors consume current from the grid. It is found that the direct-axis current takes 0.02 s to stabilize. The nature of the fluctuation in the DC-link voltage is depicted in Figure 29. The DC-link controller effectively keeps the DC-link voltage by regulating the power flow through the capacitor. Hence, it can be mentioned that the HMC effectively operates as DSTATCOM in the absence of solar radiation. Figure 31 shows the source voltage and converter output current in DSTATCOM mode. The waveforms show that the output current increases after the reactive power command comes at 0.8 s. Moreover, Figure 31 approves that the phase difference of the converter output current with grid voltage is 90[degrees] in this mode of operation.

7. Conclusion

The HMC for a grid-connected PV system shares many advantages of the CHB but uses fewer switches and reduces the size of the DC-link capacitor. However, the HMC is sensitive to the loop voltage produced by the probable DC-link voltage unbalance. The unbalanced conditions of the DC-link voltages will initiate a line-frequency circulating current. Therefore, to minimize the circulating current as well as to reduce an extra voltage drop on the inductors that affects the terminal voltage, the inductance value should be adequate. The multiwinding coupled inductors are the key to control the circulating current and ensure the proper operation of the HMC. The equivalent circuit model of the HMC configuration is derived to recognize the control scheme. The two-winding coupled inductor imposes twelve times the magnetizing inductance for the circulating current, while the impedance to the circulating current using three-winding coupled inductors is eighteen times as much as the magnetizing inductance. Owing to the six-winding coupled inductor, the impedance on the circulating current path is thirty-six times the magnetizing inductance. The results show the good performance of the control algorithm in both steady state and transient conditions. Moreover, it is interesting to note that in the absence of solar radiation, the controller acts in DSTATCOM mode to supply reactive power to the grid. The performance of the controller under different solar radiation conditions is found to be satisfactory.

https://doi.org/10.1155/2018/9164528

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to express their thanks to King Abdulaziz City for Science and Technology (KACST) for providing financial and technical support to this study.

References

[1] Z. M. Shahrbabak, A. Tabesh, and G. R. Yousefi, "Economical design of utility-scale photovoltaic power plants with optimum availability," IEEE Transactions on Industrial Electronics, vol. 61, no. 7, pp. 3399-3406, 2014.

[2] A. A. Al-Shamma'a and K. E. Addoweesh, "Techno-economic optimization of hybrid power system using genetic algorithm," International Journal of Energy Research, vol. 38, no. 12, pp. 1608-1623, 2014.

[3] J. H. R. Enslin and P. J. M. Heskes, "Harmonic interaction between a large number of distributed power inverters and the distribution network," IEEE Transactions on Power Electronics, vol. 19, no. 6, pp. 1586-1593, 2004.

[4] Y. Riffonneau, S. Bacha, F. Barruel, and S. Ploix, "Optimal power flow management for grid connected PV systems with batteries," IEEE Transactions on Sustainable Energy, vol. 2, no. 3, pp. 309-320, 2011.

[5] P. G. Arul, V. K. Ramachandaramurthy, and R. K. Rajkumar, "Control strategies for a hybrid renewable energy system: a review," Renewable and Sustainable Energy Reviews, vol. 42, pp. 597-608, 2015.

[6] J. M. Guerrero, P. C. Loh, T.-L. Lee, and M. Chandorkar, "Advanced control architectures for intelligent microgrids--part II: power quality, energy storage, and AC/DC microgrids," IEEE Transactions on Industrial Electronics, vol. 60, no. 4, pp. 1263-1270, 2013.

[7] M. A. Elsaharty, H. A. Ashour, E. Rakhshani, E. Pouresmaeil, and J. P. S. Catalao, "A novel DC-bus sensor-less MPPT technique for single-stage PV grid-connected inverters," Energies, vol. 9, no. 4, pp. 1-23, 2016.

[8] J. Rodriguez, J. S. Lai, and F. Z. Peng, "Multilevel inverters: a survey of topologies, controls, and applications," IEEE Transactions on Industrial Electronics, vol. 49, no. 4, pp. 724-738, 2002.

[9] J. A. Barrena, L. Marroyo, M. A. Rodriguez Vidal, and J. R. T. Apraiz, "Individual voltage balancing strategy for PWM cascaded H-bridge converter-based STATCOM," IEEE Transactions on Industrial Electronics, vol. 55, no. 1, pp. 21-29, 2008.

[10] K. K. Gupta, A. Ranjan, P. Bhatnagar, L. K. Sahu, and S. Jain, "Multilevel inverter topologies with reduced device count: a review," IEEE Transactions on Power Electronics, vol. 31, no. 1, pp. 135-151, 2016.

[11] L. Franquelo, J. Rodriguez, J. Leon, S. Kouro, R. Portillo, and M. Prats, "The age of multilevel converters arrives," IEEE Industrial Electronics Magazine, vol. 2, no. 2, pp. 28-39, 2008.

[12] A. Nabae, I. Takahashi, and H. Akagi, "A new neutral-point-clamped PWM inverter," IEEE Transactions on Industry Applications, vol. IA-17, no. 5, pp. 518-523, 1981.

[13] M. F. Escalante, J. C. Vannier, and A. Arzande, "Flying capacitor multilevel inverters and DTC motor drive applications," IEEE Transactions on Industrial Electronics, vol. 49, no. 4, pp. 809-815, 2002.

[14] M. A. Perez, S. Bernet, J. Rodriguez, S. Kouro, and R. Lizana, "Circuit topologies, modeling, control schemes, and applications of modular multilevel converters," IEEE Transactions on Power Electronics, vol. 30, no. 1, pp. 4-17, 2015.

[15] Y. Yu, G. Konstantinou, B. Hredzak, and V. G. Agelidis, "Power balance of cascaded H-bridge multilevel converters for large-scale photovoltaic integration," IEEE Transactions on Power Electronics, vol. 31, no. 1, pp. 292-303, 2016.

[16] E. Cengelci, P. Enjeti, C. Singh, F. Blaabjerg, and J. K. Pederson, "New medium voltage PWM inverter topologies for adjustable speed AC motor drive systems," in APEC '98. Conference Proceedings 1998, Thirteenth Annual Applied Power Electronics Conference and Exposition, pp. 565-571, Anaheim, CA, USA, 1998.

[17] R. Teodorescu, F. Blaabjerg, J. K. Pedersen, E. Cengelci, and P. N. Enjeti, "Multilevel inverter by cascading industrial VSI," IEEE Transactions on Industrial Electronics, vol. 49, no. 4, pp. 832-838, 2002.

[18] J. Wen and K. Ma Smedley, "Hexagram inverter for medium-voltage six-phase variable-speed drives," IEEE Transactions on Industrial Electronics, vol. 55, no. 6, pp. 2473-2481, 2008.

[19] J. Wen and K. M. Smedley, "Synthesis of multilevel converters based on single- and/or three-phase converter building blocks," IEEE Transactions on Power Electronics, vol. 23, no. 3, pp. 1247-1256, 2008.

[20] M. G. Vaillalva, J. R. Gazoli, and E. R. Filho, "Comprehensive approach to modeling and simulation of photovoltaic arrays," IEEE Transactions on Power Electronics, vol. 24, no. 5, pp. 1198-1208, 2009.

[21] L. Piegari, R. Rizzo, I. Spina, and P. Tricoli, "Optimized adaptive perturb and observe maximum power point tracking control for photovoltaic generation," Energies, vol. 8, no. 5, pp. 3418-3436, 2015.

[22] J. A. Barrena, L. Marroyo, M. A. R. Vidal, and J. R. T. Apraiz, "Individual voltage balancing strategy for PWM cascaded H-bridge converter-based STATCOM," IEEE Transactions on Industrial Electronics, vol. 55, no. 1, pp. 21-29, 2008.

[23] S. Debnath, J. Qin, B. Bahrani, M. Saeedifard, and P. Barbosa, "Operation, control, and applications of the modular multilevel converter: a review," IEEE Transactions on Power Electronics, vol. 30, no. 1, pp. 37-53, 2015.

[24] M. Moranchel, F. Huerta, I. Sanz, E. Bueno, and F. Rodriguez, "A comparison of modulation techniques for modular multilevel converters," Energies, vol. 9, no. 12, pp. 1-20, 2016.

[25] D. Holmes and T. Lipo, Pulse width modulation for power converters, IEEE Press, Piscataway, NJ, USA, 2003.

Abdullrahman A. Al-Shamma'a, (1,2) Abdullah M. Noman, (1,2) Khaled E. Addoweesh, (1) Ayman A. Alabduljabbar, (3) and A. I. Alolah (1)

(1) Department of Electrical Engineering, College of Engineering, KingSaud University, Riyadh 11421, Saudi Arabia

(2) Department of Mechatronics Engineering, College of Engineering, Taiz University, Taiz, Yemen

(3) King Abdulaziz City for Science and Technology (KACST), P.O. Box 6086, Riyadh 11442, Saudi Arabia

Correspondence should be addressed to Abdullrahman A. Al-Shamma'a; ashammaa@ksu.edu.sa

Received 31 December 2017; Accepted 2 April 2018; Published 7 May 2018

Academic Editor: Francesco Riganti-Fulginei

Caption: Figure 1: Topologies of cascade two-level converters.

Caption: Figure 2: HMC electrical diagram.

Caption: Figure 4: Equivalent model of a PV array.

Caption: Figure 5: I-V and P-V curves of the KC200GT PV array.

Caption: Figure 6: Flowchart for P&O MPPT method.

Caption: Figure 7: Phasor voltage of the HMC.

Caption: Figure 8: Equivalent circuit model of HMC.

Caption: Figure 9: Equivalent circuit of the two-winding coupled inductor.

Caption: Figure 10: The HMC with three two-winding coupled inductors.

Caption: Figure 11: Equivalent circuit of the three-winding coupled inductor.

Caption: Figure 12: The HMC with two three-winding coupled inductors.

Caption: Figure 13: Equivalent circuit of the six-winding coupled inductor.

Caption: Figure 14: The HMC with one six-winding coupled inductor.

Caption: Figure 15: HMC control block diagram.

Caption: Figure 16: Proposed PS-PWM diagram.

Caption: Figure 17: Three-phase output voltage of the HMC.

Caption: Figure 18: Currents inside the HMC.

Caption: Figure 19: THD of the HMC output voltage.

Caption: Figure 20: THD of the HMC output current.

Caption: Figure 21: Peak circulating current versus magnetizing inductance.

Caption: Figure 22: Circulating current inside the HMC with multiwinding coupled inductors.

Caption: Figure 23: Circulating current inside the HMC.

Caption: Figure 24: Circulating current with two-winding coupled inductors with a magnetizing inductance of 0.5 mH.

Caption: Figure 25: Circulating current with three-winding coupled inductors with a magnetizing inductance of 0.5 mH.

Caption: Figure 26: Circulating current with six-winding coupled inductors with a magnetizing inductance of 0.5 mH.

Caption: Figure 27: Step changes in solar radiation.

Caption: Figure 28: Response to changes in solar irradiance in the hexagram converter.

Caption: Figure 29: Variations in DC-link voltages in response to the changes in solar radiation.

Caption: Figure 30: Response to changes in solar radiation. (a) Output active and reactive power, and (b) direct and quadrature axis current.

Caption: Figure 31: System performance with the increase in the reactive power supply in DSTATCOM mode.

Table 1: Parameters of the KC200GT PV array at 25[degrees]C, AM1.5, and 1000W/[m.sup.2]. [I.sub.mp] [V.sub.mp] [P.sub.max] [I.sub.sc] (A) (V) (W) (A) 7.61 26.3 200.143 8.21 [V.sub.oc] [K.sub.v] [K.sub.I] [N.sub.ss] (V) (V/K) (A/K) (No.) 32.9 -0.123 0.0032 54 [R.sub.p] [R.sub.s] [I.sub.o,n] a ([OMEGA]) ([OMEGA]) (A) 415.405 0.221 9.825 x [10.sup.-8] 1.3 Table 2: General data of the proposed system. Parameter Value System parameters Nominal power P = 2.5 kVA Phase voltage and frequency [V.sub.ph] = 110 V (rms), [f.sub.sys] = 50Hz, [f.sub.swching] = 2500Hz Output current [I.sub.out] = 8.8 A (rms) DC bus voltages [V.sub.dc1] = [V.sub.dc2] = [V.sub.dc3] = [V.sub.dc4] = [V.sub.dc5] = [V.sub.dc6] = 52.6 V Open-end winding transformers [V.sub.pri]/[V.sub.sec] = 1, [f.sub.tran] = 50Hz Current controller parameters K.sub.p] = 55, K.sub.i] = 0.001 Voltage controller parameters K.sub.p] = 3, K.sub.i] = 20 Coupled inductor parameters RMS voltage, current 20 V (rms), 8.8 A (rms) Magnetizing inductance 3.2 mH Leakage inductance 150 [micro]H Table 3: Total harmonic distortion under different solar radiation conditions. THD (%) Solar radiation Voltage Current 1000 (W/[m.sup.2]) 13.64 0.68 600 (W/[m.sup.2]) 13.65 1.15 400 (W/[m.sup.2]) 12.78 2.80 0 (W/[m.sup.2]) 23.94 0.35 Fundamental value Solar radiation Voltage (V) Current (A) 1000 (W/[m.sup.2]) 184.0 8.80 600 (W/[m.sup.2]) 181.8 5.20 400 (W/[m.sup.2]) 180.9 3.44 0 (W/[m.sup.2]) 148.3 9.98

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Shammaa, Abdullrahman A. Al-; Noman, Abdullah M.; Addoweesh, Khaled E.; Alabduljabbar, Ayman A.; Alo |

Publication: | International Journal of Photoenergy |

Date: | Jan 1, 2018 |

Words: | 8062 |

Previous Article: | CFD Analysis and Electrical Efficiency Improvement of a Hybrid PV/T Panel Cooled by Forced Air Circulation. |

Next Article: | Research on Parameter Distribution Features of Photovoltaic Array under the Cover and Shadow Shading Conditions. |

Topics: |