Next Article in Journal
Evaluation of Water Resources through Efficiency Index and Water Productivity in EU
Next Article in Special Issue
Solar UAVs—More Aerodynamic Efficiency or More Electrical Power?
Previous Article in Journal
System Performance and Pollutant Emissions of Micro Gas Turbine Combined Cycle in Variable Fuel Type Cases
Previous Article in Special Issue
Experimental Study on the Behavior of Aluminum Fuse Element Inside 24 kV, 50 kA High-Voltage Fuses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Investigation of the Properties of Transients in Unbalanced Three-Phase Four-Wire Networks

Department of Electronics, Information and Bioengineering, Politecnico di Milano, 20133 Milan, MI, Italy
Energies 2022, 15(23), 9122; https://doi.org/10.3390/en15239122
Submission received: 19 October 2022 / Revised: 20 November 2022 / Accepted: 29 November 2022 / Published: 1 December 2022

Abstract

:
In this paper, the theoretical properties of transients in three-phase networks, including the fourth wire and unbalanced source, are systematically investigated by resorting to several analytical tools. First, a Cartesian space vector is introduced to provide the geometrical three-dimensional representation of a three-phase voltage/current transient. It is shown that a voltage/current transient can be always represented as a three-dimensional trajectory leaving one plane (corresponding to the previous steady state) towards another plane (corresponding to the new steady state at the end of the transient). The quantity driving the trajectory from one plane to another is the Park/Clarke zero component. Second, the Clarke transformation is used to study the original three-phase circuit as a superposition of two independent transients: the transient of the Clarke space vectors and the transient of the Clarke zero components. Third, the Park transformation is used to evaluate the dq0 components of transient voltages/currents. Since the Park transformation is related to the Clarke transformation through a simple frequency shift, the effects on the frequency content of the steady-state and transient waveforms are put into evidence. Three case studies are presented and solved to prove the theoretical results.

1. Introduction

Modern power systems for the transmission and distribution of electrical power can be seen, from a circuit theory viewpoint, as interconnections of three-phase and single-phase networks [1]. Usually, electric power transmission is implemented by resorting to three-phase networks, whereas power distribution involves both three-phase and single-phase networks/components. A typical example of a single-phase network is the so-called fourth wire of a power system, which is connected to the pure three-phase system through wye (or star) connections. A single-phase load can also be connected from one phase of the three-phase network to the fourth wire. An important example of a huge single-phase load is a railway system.
A power system, as briefly described above, is a complex network requiring specific methodologies for an effective and meaningful solution within the framework of three-phase system analysis. In fact, in general terms, a power system could be solved by resorting to conventional techniques for circuit analysis (e.g., the nodal analysis), but this approach would not allow the exploitation of the specific properties and features of three-phase systems. For this reason, in the past decades, specific methodologies have been developed for three-phase systems in order to exploit the intrinsic configuration symmetry of many power systems. The most important approach for the solution of symmetrically configured three-phase systems (i.e., three-phase systems with three equal phases and with equal coupling between the phases) is based on the well-known Symmetrical Component Transformation (SCT) [2,3,4,5]. The main advantage of the SCT is its capability of resolving the three-phase system in three uncoupled circuits in the phasor domain (i.e., the so-called positive, negative, and zero-sequence circuits), whose solution is much simpler than the solution of the original coupled circuit. Moreover, since the SCT also operates on the voltages and currents, such transformed variables (called sequence voltages/currents) have a meaningful interpretation within the framework of three-phase circuit analysis. In fact, a positive-sequence phasor corresponds to three symmetrical (or balanced) phasors with the positive sequence abc in the original three-phase circuit. A negative-sequence phasor corresponds to three symmetrical (or balanced) phasors with the negative sequence acb. Finally, a zero-sequence phasor corresponds to three equal phasors (in magnitude and phase). Thus, a generic unbalanced set of three phasors in the original three-phase network is decomposed into the sum of two balanced sets and one set with three equal components.
The approach outlined above works in the phasor domain, thus it is suited for steady-state analysis. In power system analysis, however, transient analysis in the time domain is of paramount importance. In fact, it is well-known that connecting and disconnecting apparatus and loads to the main network can result in overvoltages/overcurrents, which can lead to the malfunctioning or damage of system components [6,7,8,9,10]. Therefore, the early-stage prediction and evaluation of network transients is a critical issue in modern power systems.
Until a few years ago, the conventional approach for transient analysis of power systems was the so-called quasi-static model [11,12,13]. In the quasi-static model, the key assumption is that phasors change slowly in comparison to the system frequency. Therefore, the system is modeled through time-varying phasors. The main advantage of this approach is that the frequency is assumed constant. Thus, conventional methods such as the nodal analysis can be used directly since the nodal admittance matrix is unchanged. The main drawback, however, is that the quasi-static approach can model only slow phenomena that occur in time frames of seconds to minutes, typical of machine stability studies.
In recent years, with the increasing penetration of small, distributed generators and power electronics devices working at a high switching frequency, the quasi-static phasor approach is becoming more and more inadequate to model the fast dynamical behavior of power networks. For this reason, today, the most common approach for the transient analysis of power systems is based on the so-called dynamic phasors [13,14,15,16]. The basic idea underlying the dynamic phasors is using the Fourier series expansion for each variable involved in a transient and retaining only a small number of significant harmonics. The time-domain equations, where the dynamic phasors are used, are written for the dq0 components. The dq0 components of time-domain voltages and currents can be obtained by using the Park transformation on three-phase time-domain voltages and currents. Notice that the Park transformation was originally introduced to analyze rotating machines, but since it can be used to model static components too and since it operates in the time domain, the Park transformation can be used to model a complete power system in the time domain for transient analysis. Within this framework, several authors pointed out that the dynamic phasors are used only to obtain a more efficient formulation for the numerical solution of the dq0 time-domain equations (e.g., see [13]).
A recent paper [17] proposed a thorough and comprehensive comparison between the performances of the quasi-static model, the dynamic phasor model, and the time-domain model consisting in the solution of the dq0 differential equations through conventional numerical methods. As expected, it was proven that the dynamic phasor model can capture faster transients than the quasi-static model, but it cannot deliver the same accuracy as the time-domain model. Moreover, when fast electromagnetic transients are considered, the computational performance of the dynamic phasor approach does not surpass the time-domain model. Thus, it was made clear that, as far as fast electromagnetic transients are considered, the time-domain approach is still unsurpassed by other approximate approaches.
Given its superiority in terms of accuracy in the analysis of fast electromagnetic transients, the time-domain approach to power system transients is investigated in this paper from a methodological and theoretical viewpoint. In particular, the theoretical properties of transients in three-phase networks, including the fourth wire and unbalanced source, are systematically investigated by resorting to analytical tools such as the Cartesian space vector and the Clarke and Park transformations. More specifically, the main findings of the proposed analytical investigation can be summarized as follows.
First, a Cartesian space vector is introduced to provide the geometrical three-dimensional representation of a three-phase voltage/current [18,19]. It is shown that, under steady-state conditions, the trajectory of the Cartesian space vector lies always on a plane. This plane is the α β plane of the well-known Clarke transformation (i.e., the Park transformation with fixed axes) in the case of a Park/Clarke zero component equal to zero. Otherwise, in the case of a zero component different from zero, the plane is different from α β but it can be still properly defined [20]. Thus, a voltage/current transient can be always represented as a three-dimensional trajectory leaving one plane (corresponding to the previous steady state) towards another plane (corresponding to the new steady state at the end of the transient). The quantity driving the trajectory from one plane to another is the Park/Clarke zero component.
Second, the Clarke transformation is used to study the original three-phase circuit as a superposition of two uncoupled circuits: the space-vector equivalent circuit (taking into account possible unbalanced sources) and the zero-component equivalent circuit (taking into account the fourth wire, if any) [21,22]. Thus, the original three-phase transient is split into the superposition of two independent transients: the transient of the Clarke space vectors and the transient of the Clarke zero components. This approach allows us to put into evidence the relationship between the Clarke space vector and the Cartesian space vector. In fact, in the case of null zero components, the trajectory of the Clarke space vector equals the trajectory of the Cartesian space vector. Otherwise, in the case of zero components different from zero, the trajectory of the Clarke space vector is the projection on the α β plane of the Cartesian space-vector trajectory.
Third, the Park transformation is used to evaluate the dq0 components of transient voltages/currents. Since the Park transformation is related to the Clarke transformation through a simple frequency shift (i.e., a negative frequency shift ω , where ω is the angular frequency of the source), the effects on the frequency content of the steady-state and transient waveforms are put into evidence. In particular, as far as the steady state is concerned, a positive-sequence component leads to constant values of dq variables, whereas a negative-sequence component leads to dq variables with double frequency. As far as transient components are considered, real eigenvalues lead to damped oscillating components, whereas complex conjugate eigenvalues lead to damped oscillations with shifted frequencies.
The paper is organized as follows. In Section 2, the analytical tools mentioned above and the related theoretical derivations are presented in detail. In Section 3, three case studies are presented and solved to prove the theoretical results obtained in Section 2. In particular, two three-phase networks are presented: the first leading to a first-order circuit in the space-vector domain; the second leading to a second-order circuit in the space-vector domain. Moreover, the solution of a three-bus power system is presented. In all the cases, the fourth wire and the effect of an unbalanced source are considered and discussed. Concluding remarks are drawn in Section 4.

2. Space-Vector Analysis of Three-Phase Transients

Let us consider a generic set of three-phase voltages v a ,   v b ,   v c in the time domain, under steady-state or transient conditions (the following derivations also hold for a generic set of three-phase currents i a ,   i b ,   i c ). By introducing an orthogonal abc space where each axis is associated to a phase voltage, the set of three-phase voltages can be represented by the Cartesian space vector [18].
v a b c = v a e ^ a + v b e ^ b + v c e ^ c
where e ^ a ,   e ^ b ,   e ^ c is the set of unit vectors of the orthogonal abc space. If no specific assumptions are made about the set of voltages v a ,   v b ,   v c , the trajectory (or locus) of v a b c is a non-planar trajectory (e.g., see Figure 1).
The analysis of three-phase quantities in the time domain is often performed by resorting to the Clarke transformation, defined as [10]:
v α v β v 0 = T 0 v a v b v c = 2 3 1 1 / 2 1 / 2 0 3 / 2 3 / 2 1 / 2 1 / 2 1 / 2 v a v b v c
where the transformation matrix T 0 is given in its power-invariant form, i.e., T 0 1 = T 0 T . Notice that, from the inverse of Equation (2), we obtain:
v a v b v c = 2 3 1 0 1 / 2 1 / 2 3 / 2 1 / 2 1 / 2 3 / 2 1 / 2 v α v β v 0 = v α e ^ α + v β e ^ β + v 0 e ^ 0
where
e ^ α = 2 3 1 1 / 2 1 / 2 ,   e ^ β = 2 3 0 3 / 2 3 / 2 ,   e ^ 0 = 2 3 1 / 2 1 / 2 1 / 2
is an orthonormal basis, i.e.,
e ^ α · e ^ α = e ^ β · e ^ β = e ^ 0 · e ^ 0 = 1   and   e ^ α · e ^ β = e ^ β · e ^ 0 = e ^ α · e ^ 0 = 0
In the special case where the zero component v 0 in Equation (2) is null, the set of three-phase voltages is such that v a + v b + v c = 0 , i.e., the trajectory of the Cartesian space vector v a b c lies on the plane α β with equation a + b + c = 0 and is defined by the unit vectors e ^ α ,   e ^ β (see Figure 2). Therefore, in this case, the three-phase voltages v a ,   v b ,   v c can be completely represented by the two voltages v α ,   v β , i.e., by the voltage space vector, defined as:
v ¯ = v α + j v β
In the general case, i.e., when v 0 0 , the Cartesian space vector v a b c does not lie on the plane α β since the e ^ 0 component is also present in Equation (3). In fact, it can be easily shown from Equation (3) that the general equations to recover the phase variables v a ,   v b ,   v c from the space vector Equation (6) and the zero component v 0 can be written:
v a = 2 3 R e v ¯ + 1 3 v 0 ,       v b = 2 3 R e a 2 v ¯ + 1 3 v 0 ,       v c = 2 3 R e a v ¯ + 1 3 v 0
where a = e j 2 π / 3 . Thus, the locus of the zero component v 0 is the line along the direction e ^ 0 perpendicular to the plane α β , moving the trajectory of the Cartesian space vector v a b c outside the plane α β . Notice that this property holds regardless of the time behavior of the variables, i.e., under steady-state or transient conditions.

2.1. Balanced and Unbalanced Sinusoidal Steady State

Let us consider a set of sinusoidal three-phase voltages v a ,   v b ,   v c with angular frequency ω and a zero component v 0 = 0 (i.e., the balanced case). In this case, the trajectory of the Cartesian space vector v a b c lies on the α β plane (Figure 2), and the corresponding space vector Equation (6) is given by:
v ¯ = V p e j ω t + V n e j ω t
where, according to the well-known Symmetrical Component Transformation (SCT) operating on the phasors V a , V b , V c , V p and V n are the phasors of the positive and negative (complex conjugate)-sequence components given by [1]:
V p V n V 0 = S V a V b V c = 1 3 1 a a 2 1 a 2 a 1 1 1 V a V b V c
where the transformation matrix S is in its power-invariant form, i.e., S 1 = S T .
The trajectory of the space vector v ¯ in Equation (8) is elliptical (Figure 2), with semi-major and semi-minor axes given by [22]:
M = V p + V n ,   m = V p V n
and the inclination angle by:
φ = a r g V p a r g V n / 2
As a special case, when the negative-sequence component V n = 0 , the trajectory of the space vector v ¯ becomes circular with a radius r = V p .
In the more general case, where the set of sinusoidal three-phase voltages v a ,   v b ,   v c , with angular frequency ω , has a zero component v 0 0 (i.e., the unbalanced case), the trajectory of the Cartesian space vector v a b c is still planar, but the plane is different from the α β plane. Thus, the effect of a sinusoidal steady-state zero component v 0 is the definition of a new plane where the trajectory of v a b c lies. This point can be proven by observing that, since each phase voltage (i.e., v a , v b , or v c ) is a sine wave with angular frequency ω , such voltage can be seen as a solution of a simple harmonic oscillator defined by a homogeneous second-order differential equation. Since this is true for each phase voltage, the second-order differential equation can be written for the Cartesian space vector [20]:
d 2 d t 2 v a b c + ω 2 v a b c = 0
Thus, the solution can be written as a linear combination of two sine waves:
v a b c = g 1 · cos ω t + g 2 · sin ω t
where g 1 and g 2 are constant vectors given by the initial conditions:
g 1 = v a b c t = 0 ,   g 2 = 1 ω d d t v a b c t = 0
Equation (13) is the vector representation of a plane defined by the vectors g 1 , g 2 and the axes origin. Therefore, the trajectory of the Cartesian space vector v a b c is always planar and lying on the plane g 1 g 2 . Such a plane becomes the plane α β in the special case of the null zero component (i.e., V 0 = 0 ), whereas the difference with the α β plane increases as the zero component V 0 increases, i.e., for increasing unbalance (see Figure 3). Notice that, in Figure 3, the vectors g 1 , g 2 terminate on the ellipse because of Equation (14), whereas the unit vectors e ^ α ,   e ^ β define the α β plane and, in general, they do not terminate on the ellipse on the α β plane.

2.2. Transient Conditions

The dynamics of a three-phase network can be analyzed through the Clarke transformation Equation (2) and the use of voltage/current space vectors defined as in Equation (6). In fact, under the assumption of a symmetrically configured three-phase network, it is well-known that the original network can be studied as two uncoupled networks [21]:
  • The space vector network, where all the voltages and currents are defined as space vectors, and the three-phase components are diagonalized by the Clarke transformation (e.g., a three-phase symmetrical mutual inductor with phase self-inductances L p h and mutual inductances L m , corresponds to the inductance L = L p h L m in the space-vector domain after diagonalization through the Clarke transformation);
  • The zero-component network, where all the circuit variables are the zero-component variables, the three-phase components are given by Clarke diagonalization (e.g., L 0 = L p h + 2 L m in the example mentioned above), and single-phase networks (e.g., the fourth wire) connected to the three-phase network are taken into account.
Once the two uncoupled equivalent circuits (i.e., the space-vector and the zero-component circuits) are defined, the conventional state-space equations can be written and solved. Therefore, for the space-vector circuit, we can write the dynamic model [21]:
d d t x ¯ = A x ¯ + B u ¯
where x ¯ is the vector of the space-vector state variables, and u ¯ is the vector of the space-vector inputs. The matrices A and B are given by the conventional matrices of the state-space approach.
The solution for each state variable (i.e., the generic component of the vector x ¯ ) can be written as:
x ¯ = k = 1 N C k e λ k t + X p e j ω t + X n e j ω t
where, for the sake of simplicity, the case of N distinct eigenvalues λ k k = 1 N was considered, and the complex coefficients C k k = 1 N can be determined by setting the initial conditions. Notice that the steady-state component of the solution Equation (16) has the form of Equation (8), whereas the transient component decreases to zero for a stable network.
As far as the zero components are considered, the solution of the corresponding state-space model for the generic state variable can be written in the following form:
x 0 = k = 1 M D k e j μ k t + 2 X 0 c o s ω t + ψ
where the eigenvalues μ k k = 1 M are different from the eigenvalues in Equation (16) since the zero-component circuit has different topology/parameters.
According to the properties previously derived, the space-vector solution Equation (16) lies on the α β plane, whereas the zero-component solution Equation (17) (if present) provides a component orthogonal to the α β plane, leading to a non-planar total solution for the corresponding Cartesian space vector. In general terms, by assuming a transient starting at t = 0 , the total transient solution Equations (16) and (17) describe the trajectory of the corresponding Cartesian space vector x a b c from the steady-state plane g 1 g 2 for t < 0 towards a new steady-state plane g 1 g 2 for t , where the steady-state plane g 1 g 2 can be the α β plane in the special case of the null zero component.

2.3. The Park Transformation

Let us consider an orthogonal reference frame dq (i.e., direct and quadrature axes) rotating at angular speed ω on the α β plane. The components of the Cartesian space vector on the rotating axes dq can be readily obtained by rotating the α β components by the angle ϑ = ω t . Thus, the Park transformation can be seen as a composite transformation, where the Clarke transformation is followed by axes rotation [18]:
v d v q v 0 = H ϑ T 0 v a v b v c = c o s ϑ s i n ϑ 0 s i n ϑ c o s ϑ 0 0 0 1 2 3 1 1 / 2 1 / 2 0 3 / 2 3 / 2 1 / 2 1 / 2 1 / 2 v a v b v c = T v a v b v c
Notice that, in Equation (18), the zero component v 0 remains unchanged. Moreover, for ϑ = 0 , the dq components equal the α β components.
The space vector defined on the rotating frame dq is given by:
v ¯ d q = v d + j v q
According to Equation (18), the relationship between the space-vector Equation (6), defined on the fixed axes α β , and the space-vector Equation (19), defined on the rotating axes dq, can be written:
v ¯ d q = v ¯ e j ϑ = v ¯ e j ω t
Using Equation (20) in Equation (16) (i.e., multiplying Equation (16) by e j ω t ) allows us to highlight the following important properties of transient space-vector solutions when the Park transformation is used:
  • The positive-sequence component X p e j ω t of the steady-state solution becomes the constant phasor X p . Thus, the dq components v d and v q are constant values.
  • The negative-sequence component X n e j ω t of the steady-state solution becomes the space vector X n e j 2 ω t rotating at double negative speed. Corresponding dq components are second-harmonic sine waves.
  • Transient components C k e j λ k t have a negative frequency shift equal to ω . In particular:
    • For real eigenvalues λ k = α k : each term C k e j λ k t becomes C k e α k t e j ω t , thus we obtain a damped space vector rotating at negative speed ω , and the corresponding dq components are damped sine waves at ω .
    • For a couple of complex conjugate eigenvalues α k   ±   j β k : we obtain the terms C k e α k t e j β k ω t and C k + 1 e α k t e j β k + ω t , corresponding to damped space vectors rotating at the shifted angular speeds β k ω and β k + ω . The related dq components contain damped sine waves with shifted angular frequencies β k ω and β k + ω .

3. Test Cases and Discussion

The theoretical properties and results derived in Section 2 are validated and highlighted in this section by resorting to three significant case studies. The first case study consists in a three-phase network leading to a first-order circuit in the space-vector domain, whereas the second case study leads to a second-order circuit in the space-vector domain. In both the cases, the fourth wire is also involved in the transients. The third case study is a three-bus power system. The proposed networks are solved analytically by resorting to the general approach based on the Clarke transformation and the space vectors. The analytical results, however, have been properly checked through numerical simulations with Matlab/Simulink.

3.1. First Case Study: First-Order Circuit in the Space-Vector Domain

Figure 4a shows a three-phase network including a three-phase mutual inductor and resistive elements. The three-phase resistor R 3 and the grounding resistor R 02 are both inserted at t = 0 by closing the corresponding switches. The equivalent circuits in the space-vector domain and in the zero-component domain are shown in Figure 4b,c, respectively. In particular, notice that, for the space vectors, the inductance L = L p h L m , whereas for the zero components, L 0 = L p h + 2 L m . Moreover, it is well-known that the single-phase ground network R 01 ,   R 02 has no effect on the space-vector equivalent circuit since each three-phase wye connection is equivalent to a short circuit for space vectors. On the contrary, in the zero-component equivalent circuit, the single-phase ground network holds its topology, but resistive parameters must be multiplied by the factor 3 (i.e., 3 R 01 ,   3 R 02 ).
The space-vector equivalent circuit in Figure 4b is a first order circuit whose analytical solution for the inductor current i ¯   t can be written in the usual form:
i ¯ t = i ¯ 0 i s ¯ 0 e t / τ + i s ¯ t
where i ¯ 0 is the initial condition, τ = L / R 2 + R 1 / / R 3 , and i s ¯ t is the steady-state solution given by:
i s ¯ t = e p ¯ t Z e q · R 3 R 3 + R 2 + Z L + e n ¯ t Z e q · R 3 R 3 + R 2 + Z L
where Z L = j ω L , Z e q = R 1 + R 3 // R 2 + Z L , and:
e ¯ t = E p e j ω t + E n e j ω t = e p ¯ t + e n ¯ t
The initial condition i ¯ 0 in Equation (21) can be readily calculated by evaluating the steady-state solution at t = 0 before the operation of the switch (i.e., by considering the switch as an open circuit):
i ¯ 0 = E p R 1 + R 2 + Z L + E n R 1 + R 2 + Z L
The zero-component equivalent circuit in Figure 4c is a conventional first-order circuit whose analytical solution can be written as:
i 0 t = i 0 0 i 0 s 0 e t / τ 0 + i 0 s t
where i 0 0 is the initial condition ( i 0 0 = 0 in this specific case), τ 0 = L 0 / R 2 + 3 R 02 + R 3 / / R 1 + 3 R 01 , and i 0 s t is the steady-state solution corresponding to the phasor solution:
I 0 = E 0 Z e q 0 · R 3 R 3 + R 2 + Z L 0 + 3 R 02
where Z L 0 = j ω L 0 , and Z e q 0 = R 1 + 3 R 01 + R 3 / / R 2 + Z L 0 + 3 R 02 .
Once the space vector i ¯ t and the zero component i 0 t are calculated, the phase variables abc can be recovered through Equation (7) when written in terms of currents. Therefore, the Cartesian space vector i a b c , the α and β components of the space vector i ¯ t , and the Park space vector i ¯ d q t with the corresponding d and q components can be readily calculated from the relevant relationships reported in Section 2.
A total of two different source conditions were considered for the three-phase network in Figure 4. First, a three-phase voltage source with a negative-sequence component E n = 0 . Second, a three-phase voltage source with E n 0 . In both cases, the zero component E 0 of the three-phase voltage source was different from zero. This choice was made to put into evidence the effects of source unbalancing and the source zero component. The circuit parameters were selected as follows: f = 50   Hz ,   E p = 100   V ,   E n = 0   first   case ,   E n = 30 e j π / 6   V   second   case ,   E 0 = 50 e j π / 4   V ,     R 1 = R 2 = R 3 = 1   Ω ,   R 01 = R 02 = 0.1   Ω ,   L p h = 50   mH ,   L m = 0 . Notice that the numerical values of the circuit parameters were not related to a specific application, but they were selected with the only objective of putting clearly into evidence the properties of the circuit under analysis.

3.1.1. Three-Phase Voltage Source with E n = 0

Figure 5 shows the trajectory of the Cartesian space vector i a b c . The black line corresponds to the previous steady state (i.e., the steady-state condition for t < 0 ). At t = 0 , the curve (in red color) leaves the black trajectory and shows a transient trajectory towards the new steady state. Notice that, for t < 0 , the trajectory lies on the α β plane because the zero component i 0 = 0 (in fact, all the switches were open for t < 0 ). On the contrary, since the new steady-state condition has i 0 0 (due to the connection of R 02 ), the corresponding steady-state trajectory lies on the g 1 g 2 plane defined by Equation (14) when written for i a b c .
Figure 6a shows the trajectory of the space vector i ¯ t on the α β plane. The steady-state trajectory at the end of the transient is circular since, in this case, the negative-sequence component of the current is zero (as a consequence of the assumption E n = 0 ). The corresponding time-domain behavior of the real and the imaginary parts of the space vector i ¯ t , i.e., the components i α t and i β t , is shown in Figure 6b.
Figure 7 shows the behavior of the zero-component current i 0 t , which is responsible for driving the Cartesian space vector i a b c in Figure 5 from the α β plane to the g 1 g 2 plane.
Figure 8a shows the trajectory of the Park space vector i ¯ d q t on the dq plane. The trajectory converges towards the positive-sequence phasor (blue line) corresponding to the new steady-state solution. Figure 8b shows the corresponding time-domain behavior of the d and q components of the Park space vector i ¯ d q t , i.e., the currents i d t and i q t , converging towards the real part and the imaginary part of the positive-sequence phasor in Figure 8a (blue line).

3.1.2. Three-Phase Voltage Source with E n 0

Figure 9 shows the trajectory of the Cartesian space vector i a b c . Also in this case, as in Figure 5, for t < 0 , the trajectory lies on the α β plane because the zero component i 0 = 0 , and since the new steady-state condition has i 0 0 , the corresponding steady-state trajectory lies on the g 1 g 2 plane defined by Equation (14) when written for i a b c .
Figure 10a shows the trajectory of the space vector i ¯ t on the α β plane. In this case, the steady-state trajectory at the end of the transient is elliptical since the negative-sequence component of the current is different from zero (as a consequence of the assumption E n 0 ). The corresponding time-domain behavior of the real and the imaginary parts of the space vector i ¯ t , i.e., the components i α t and i β t , is shown in Figure 10b.
Notice that the zero-component current i 0 t is not affected by the negative-sequence component in the voltage source; therefore, its time-domain behavior is still the behavior represented in Figure 7.
Figure 11a shows the trajectory of the Park space vector i ¯ d q t on the dq plane. The straight blue line shows the positive-sequence phasor. In this case, however, the new steady-state trajectory is a circular trajectory (centered on the positive-sequence phasor) corresponding to the negative-sequence phasor rotating at double (negative) angular frequency 2 ω . This is apparent in Figure 11b where, according to the theory presented in Section 2, both a transient and a steady-state component 2 ω can be distinguished.

3.2. Second Case Study: Second-Order Circuit in the Space-Vector Domain

Figure 12a shows a three-phase network including a three-phase mutual inductor, a three-phase capacitor, and resistive elements. The three-phase capacitor C and the ground resistor R 0 are inserted at t = 0 by closing the switch. The equivalent circuits in the space-vector domain and in the zero-component domain are shown in Figure 12b,c, respectively. In particular, notice that, for the space vectors, the inductance L = L p h L m , whereas for the zero components, L 0 = L p h + 2 L m . Moreover, the single-phase ground network R 0 has no effect on the space-vector equivalent circuit since each three-phase wye connection is equivalent to a short circuit for space vectors. On the contrary, in the zero-component equivalent circuit, the ground resistance R 0 must be multiplied by the factor 3, whereas the wye connection of the mutual inductor corresponds to an open circuit.
The space-vector equivalent circuit in Figure 12b is a second-order circuit whose conventional state-space representation, according to Equation (15), is given by [21]:
d d t i ¯ v ¯ = A i ¯ v ¯ + B e ¯ = R 2 / L 1 / L 1 / C 1 / R 1 C i ¯ v ¯ + 0 1 / R 1 C e ¯
The analytical solution for the inductor current i ¯ t and the capacitor voltage v ¯ t can be written in the usual form:
i ¯ t = C 1 , 1 e λ 1 t + C 1 , 2 e λ 2 t + i s ¯ t
v ¯ t = C 2 , 1 e λ 1 t + C 2 , 2 e λ 2 t + v s ¯ t
where λ 1 ,   λ 2 are the eigenvalues of the state matrix A (here, distinct eigenvalues are assumed), i s ¯ t and v s ¯ t are the steady-state solutions, and C i , k   i = 1 , 2   k = 1 , 2 are the constant coefficients depending on the initial state:
C = C 1 , 1 C 1 , 2 C 2 , 1 C 2 , 2 = x ¯ 0 A x ¯ 0 1 λ 1 1 λ 2 1
where x ¯ 0 = i ¯ 0 i s ¯ 0 v ¯ 0 v s ¯ 0 T .
The steady-state solutions can be readily obtained from Figure 12b:
i s ¯ t = e p ¯ t Z e q · Z C Z C + R 2 + Z L + e n ¯ t Z e q · Z C Z C + R 2 + Z L
v ¯ s t = e p ¯ t · Z L C Z e q + e n ¯ t · Z L C Z e q
where Z L = j ω L , Z C = 1 / j ω C , Z L C = Z C // R 2 + Z L , and Z e q = R 1 + Z L C .
The initial condition i ¯ 0 in (30) can be readily calculated by evaluating the steady-state solution at t = 0 before the operation of the switch (i.e., by considering the switch as an open circuit):
i ¯ 0 = E p R 1 + R 2 + Z L + E n R 1 + R 2 + Z L
whereas the initial condition v ¯ 0 is a given data. For the sake of simplicity, hereafter, we assume v ¯ 0 = 0 .
The zero-component equivalent circuit in Figure 12c is a conventional first-order RC circuit (in fact i 0 t = 0 ) whose analytical solution can be written as:
v 0 t = v 0 0 v 0 s 0 e t / τ 0 + v 0 s t
where v 0 0 is the initial condition ( v 0 0 = 0 in this specific case), τ 0 = R 1 + 3 R 0 C , and v 0 s t is the steady-state solution corresponding to the phasor solution:
V 0 = E 0 · Z C R 1 + 3 R 0 + Z C
Once the space vectors i ¯ t , v ¯ t and the zero components i 0 t , v 0 t were calculated, the phase variables abc can be recovered through Equation (7) when written in terms of currents and voltages. Therefore, the Cartesian space vectors i a b c , v a b c , the α and β space-vector components, and the Park space vectors i ¯ d q t , v ¯ d q t with the corresponding d and q components can be readily calculated from the relevant relationships reported in Section 2.
Similar to Section 3.1, two different source conditions were considered for the three-phase network in Figure 12. First, a three-phase voltage source with a negative-sequence component E n = 0 . Second, a three-phase voltage source with E n 0 . In both cases, the zero component E 0 of the three-phase voltage source was different from zero. This choice was made to put into evidence the effects of source unbalancing and the source zero component.
The circuit parameters were selected as follows: f = 50   Hz ,   E p = 100   V ,   E n = 0   first   case ,   E n = 30 e j π / 6   V   second   case ,   E 0 = 100 e j π / 4   V ,     R 1 = 1   Ω ,   R 2 = R 0 = 0.1   Ω ,   L p h = 1   mH ,   L m = 0 , C = 10   mF . Notice that the numerical values of the circuit parameters were not related to a specific application, but they were selected with the only objective of putting clearly into evidence the properties of the circuit under analysis. In particular, the selected parameters lead to complex conjugate eigenvalues λ 1 , 2 = 100 ± j 316 . Thus, since the imaginary part 316 is very close to the source angular frequency ω (where ω = 2 π f 314   s 1 ), according to Equation (2), the Park space vectors contain a damped component with a frequency close to zero and a damped component with nearly double frequency 2 ω . This point is clarified by the results reported in the following sections.

3.2.1. Three-Phase Voltage Source with E n = 0

Figure 13 shows the trajectory of the Cartesian space vector i a b c . The black line corresponds to the previous steady state (i.e., the steady-state condition for t < 0 ). At t = 0 , the curve (in red color) leaves the black trajectory and shows a transient trajectory towards the new steady state. Notice that the transient trajectory also lies on the α β plane because the zero component i 0 = 0 (in fact, the wye connection of the mutual inductor is not connected to the single-phase network).
Figure 14a shows the trajectory of the space vector i ¯ t on the α β plane. The steady-state trajectory at the end of the transient is circular since, in this case, the negative-sequence component of the current is zero (as a consequence of the assumption E n = 0 ). The corresponding time-domain behavior of the real and the imaginary parts of the space vector i ¯ t , i.e., the components i α t and i β t , is shown in Figure 14b.
Figure 15a shows the trajectory of the Park space vector i ¯ d q t on the dq plane. The trajectory converges towards the positive-sequence phasor (blue line) corresponding to the new steady-state solution. Figure 15b shows the corresponding time-domain behavior of the d and q components of the Park space vector i ¯ d q t , i.e., the currents i d t and i q t , converging towards the real part and the imaginary part of the positive-sequence phasor in Figure 15a (blue line). As mentioned before, due to the specific numerical values of the imaginary part of the eigenvalues, we expect transient components at 2 ω . This is apparent in Figure 15b, where an oscillating behavior with a period T = 10   ms (i.e., half period with respect to the source frequency f = 50   Hz ) is put into evidence.
Figure 16 shows the trajectory of the Cartesian space vector v a b c . Since the assumed initial condition is v a b c 0 = 0 , the transient trajectory starts at the axes’ origin. Moreover, due to the zero component in the voltage source, the capacitor voltages have a transient behavior including the zero component v 0 t . This component is responsible for a transient trajectory converging to a g 1 g 2 plane different from the α β plane.
Figure 17a shows the trajectory of the space vector v ¯ t on the α β plane. The steady-state trajectory at the end of the transient is circular since, in this case, the negative-sequence component of the voltage is zero (as a consequence of the assumption E n = 0 ). The corresponding time-domain behavior of the real and the imaginary parts of the space vector v ¯ t , i.e., the components v α t and v β t , is shown in Figure 17b.
Figure 18 shows the behavior of the zero-component voltage v 0 t , which is responsible for driving the Cartesian space vector v a b c in Figure 16 from the axes’ origin to the g 1 g 2 plane instead of the α β plane.
Figure 19a shows the trajectory of the Park space vector v ¯ d q t on the dq plane. The trajectory converges towards the positive-sequence phasor (blue line) corresponding to the new steady-state solution. Figure 19b shows the corresponding time-domain behavior of the d and q components of the Park space vector v ¯ d q t , i.e., the voltages v d t and v q t , converging towards the real part and the imaginary part of the positive-sequence phasor in Figure 19a (blue line). Also in this case, due to the specific numerical values of the imaginary part of the eigenvalues, we expect transient components at 2 ω . This is apparent in Figure 19b, where an oscillating behavior with a period T = 10   ms (i.e., half period with respect to the source frequency f = 50   Hz ) is put into evidence.

3.2.2. Three-Phase Voltage Source with E n 0

Figure 20 shows the trajectory of the Cartesian space vector i a b c . The black line corresponds to the previous steady state (i.e., the steady-state condition for t < 0 ). At t = 0 , the curve (in red color) leaves the black trajectory and shows a transient trajectory towards the new steady state. Notice that the transient trajectory also lies on the α β plane because the zero component i 0 = 0 (in fact, the wye connection of the mutual inductor is not connected to the single-phase network).
Figure 21a shows the trajectory of the space vector i ¯ t on the α β plane. The steady-state trajectory at the end of the transient is elliptical since, in this case, the negative-sequence component of the current is not zero (as a consequence of the assumption E n 0 ). The corresponding time-domain behavior of the real and the imaginary parts of the space vector i ¯ t , i.e., the components i α t and i β t , is shown in Figure 21b.
Figure 22a shows the trajectory of the Park space vector i ¯ d q t on the dq plane. The trajectory converges towards a circle centered on the positive-sequence phasor (blue line). The circular steady-state trajectory is due to the negative-sequence component of the solution. Figure 22b shows the corresponding time-domain behavior of the d and q components of the Park space vector i ¯ d q t , i.e., the currents i d t and i q t , converging towards the real part and the imaginary part of the positive-sequence phasor in Figure 22a (blue line), with a superimposed double-frequency solution due to the negative-sequence component. As mentioned before, due to the specific numerical values of the imaginary part of the eigenvalues, we expect transient components at 2 ω . This is apparent in Figure 22b, where such a transient oscillating component with a period T = 10   ms (i.e., half period with respect to the source frequency f = 50   Hz ) is superimposed on the double-frequency steady-state solution due to the negative-sequence component.
Figure 23 shows the trajectory of the Cartesian space vector v a b c . Since the assumed initial condition is v a b c 0 = 0 , the transient trajectory starts at the axes’ origin. Moreover, due to the zero component in the voltage source, the capacitor voltages have a transient behavior including the zero component v 0 t . This component is responsible for a transient trajectory converging to a g 1 g 2 plane different from the α β plane.
Figure 24a shows the trajectory of the space vector v ¯ t on the α β plane. The steady-state trajectory at the end of the transient is elliptical since, in this case, the negative-sequence component of the voltage is not zero (as a consequence of the assumption E n 0 ). The corresponding time-domain behavior of the real and the imaginary parts of the space vector v ¯ t , i.e., the components v α t and v β t , is shown in Figure 24b.
The zero-component voltage v 0 t , which is responsible for driving the Cartesian space vector v a b c in Figure 23 from the axes’ origin to the g 1 g 2 plane instead of the α β plane, is not affected by the negative-sequence component in the voltage source. Therefore, its time behavior is the same as in Figure 18.
Figure 25a shows the trajectory of the Park space vector v ¯ d q t on the dq plane. The trajectory converges towards a circular trajectory, due to the negative-sequence component, centered on the positive-sequence phasor (blue line). Figure 25b shows the corresponding time-domain behavior of the d and q components of the Park space vector v ¯ d q t , i.e., the voltages v d t and v q t , converging towards the real part and the imaginary part of the positive-sequence phasor in Figure 25a (blue line), with a superimposed double-frequency steady-state solution due to the negative-sequence component. Also in this case, due to the specific numerical values of the imaginary part of the eigenvalues, we expect transient components at 2 ω . This is apparent in Figure 25b, where both the transient oscillating behavior and the steady-state oscillating behavior with a period T = 10   ms (i.e., half period with respect to the source frequency f = 50   Hz ) can be identified.

3.3. Third Case Study: Three-Bus Power System

Figure 26 shows a three-bus medium-voltage (MV) 60 Hz power system [23]. The analytical methods outlined in Section 2 can be used to study the transient behavior of the power system once it is energized by the generator at t = 0 . Also in this case, the analytical results have been checked against numerical results obtained by means of Matlab/Simulink simulations.
Generator and loads data are reported in Figure 26. The positive-sequence and the zero-sequence impedances of the lines are Z = 0.1901 + j 0.3937   Ω / km and Z 0 = 0.368 + j 1.84   Ω / km , respectively. Notice that, in order to use the analytical results derived in Section 2, the loads must be converted into equivalent impedances. This can be easily done by recalling that the complex power S L and the impedance Z L of a given load are related by Z L = V 2 / S L , where V is the reference phase-to-phase voltage (i.e., V = 12.47   kV in this case), and the asterisk denotes a complex conjugate. Figure 27 shows the topology of the space-vector equivalent circuit. Notice that the same topology holds for the zero components in case the ground is involved in the current path.
The space-vector equivalent circuit corresponding to Figure 27 can be solved by resorting to the conventional state-space approach, leading to the dynamical model in Equation (15). Notice that the effective dynamical order of the circuit in Figure 27 is four (instead of seven) because the seven inductive currents have three independent constraints given by the Kirchhoff Laws at the three buses. By assuming that the network is energized by the generator at t = 0 , and by setting the generator with a 5% negative-sequence voltage component, the trajectory of the Cartesian space vector of the generator currents is represented in Figure 28. Notice that, since the generator currents have a null zero component, the Cartesian space vector lies on the α β plane. The space vector i ¯ g t of the generator currents is represented in Figure 29a, whereas the corresponding α and β components are represented in Figure 29b as functions of time. Notice that the steady-state trajectory in Figure 29a is elliptical because of the 5% negative-sequence component in the generator voltage. Figure 30a shows the trajectory of the Park space vector of the generator currents, whereas Figure 30b shows the time-domain behavior of the corresponding dq components. Notice that, in Figure 30a, the negative-sequence component in the steady-state currents results in the small circle around the phasor (blue line) corresponding to the positive-sequence component solution. The negative-sequence component can also be observed in Figure 30b, where oscillations at frequency 2 f = 120   Hz can be detected. Notice that, in this case, the oscillations have only the component at frequency 2 f , whereas the phenomenon of frequency shifting cannot be observed because the eigenvalues of the network in Figure 27 are real eigenvalues (complex conjugate eigenvalues would be possible only for an RLC network).
Finally, the effect of the zero-component current can be investigated by considering the ground connections of the load and of the secondary of the transformer (normally wye-connected to avoid the circulation of zero-component currents). In this case, the zero-component current has a transient behavior (see Figure 31), leading to a Cartesian space vector lying on the g 1 g 2 plane (see Figure 32) instead of the α β plane, as in Figure 28.

4. Conclusions

The analytical tools used in the paper allow a comprehensive investigation of the properties of three-phase transients with unbalanced sources and involving the fourth wire. In particular, a complete geometrical representation is given by the Cartesian space vector of a three-phase variable, whereas the space vectors based on the Clarke transformation allow us to study the original three-phase transient as a superposition of two independent transients. The dq0 representation provided by the Park transformation was also discussed, and the effects of the frequency shift with respect to the Clarke transformation were shown through specific case studies.
The methodology and the properties derived and discussed in the paper provide a general theoretical framework for three-phase transient analysis. Thus, the paper is intended to provide a theoretical contribution, whereas the selection of numerical methods for the effective solution of time-domain equations is the topic of the cited relevant literature.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Bellan, D.; Pignari, S.A.; Superti-Furga, G. Consistent circuit technique for zero-sequence currents evaluation in interconnected single/three-phase power networks. J. Electr. Syst. 2016, 12, 230–238. [Google Scholar]
  2. Fortescue, C.L. Method of symmetrical coordinates applied to the solution of polyphase networks. Trans. AIEE 1918, 37, 1027–1140. [Google Scholar]
  3. Bellan, D. Circuit modeling and statistical analysis of differential-to-common-mode noise conversion due to filter unbalancing in three-phase motor drive systems. Electronics 2020, 9, 1612. [Google Scholar] [CrossRef]
  4. Das, J.C. Understanding Symmetrical Components for Power System Modeling; Wiley: Hoboken, NJ, USA, 2017. [Google Scholar]
  5. Chicco, G.; Mazza, A. 100 Years of Symmetrical Components. Energies 2019, 12, 450. [Google Scholar] [CrossRef] [Green Version]
  6. Bollen, M.H.J.; Yu-Hua Gu, I. On the analysis of voltage and current transients in three-phase power systems. IEEE Trans. Power Deliv. 2007, 22, 1194–1201. [Google Scholar] [CrossRef]
  7. Mahseredjian, J.; Dinavahi, V.; Martinez, J.A. Simulation tools for electromagnetic transients in power systems: Overview and challenges. IEEE Trans. Power Deliv. 2009, 24, 1657–1669. [Google Scholar] [CrossRef]
  8. Van der Sluis, L. Transients in Power Systems; John Wiley & Sons Ltd: Chichester, UK, 2001. [Google Scholar]
  9. Greenwood, A. Electrical Transients in Power Systems, 2nd ed.; Wiley India (P.) Ltd.: New Delhi, India, 1991. [Google Scholar]
  10. Bellan, D. Clarke transformation solution of asymmetrical transients in three-phase circuits. Energies 2020, 13, 5231. [Google Scholar] [CrossRef]
  11. Ilic, M.; Zaborszky, J. Dynamics and Control of Large Electric Power Systems; Wiley-IEEE Press: New York, NY, USA, 2000. [Google Scholar]
  12. Baimel, D.; Belikov, J.; Guerrero, J.M.; Levron, Y. Dynamic Modeling of Networks, Microgrids, and Renewable Sources in the dq0 Reference Frame: A Survey. IEEE Access 2017, 5, 21323–21335. [Google Scholar] [CrossRef]
  13. Vega-Herrera, J.; Rahmann, C.; Valencia, F.; Strunz, K. Analysis and Application of Quasi-Static and Dynamic Phasor Calculus for Stability Assessment of Integrated Power Electric and Electronic Systems. IEEE Trans. Power Syst. 2021, 36, 1750–1760. [Google Scholar] [CrossRef]
  14. Stankovic, A.M.; Lesieutre, B.C.; Aydin, T. Modeling and analysis of single-phase induction machines with dynamic phasors. IEEE Trans. Power Syst. 1999, 14, 9–14. [Google Scholar] [CrossRef] [Green Version]
  15. Miao, Z.; Piyasinghe, L.; Khazaei, J.; Fan, L. Dynamic Phasor-Based Modeling of Unbalanced Radial Distribution Systems. IEEE Trans. Power Syst. 2015, 30, 3102–3109. [Google Scholar] [CrossRef]
  16. Levron, Y.; Belikov, J. Modeling power networks using dynamic phasors in the dq0 reference frame. Electr. Power Syst. Res. 2017, 144, 233–242. [Google Scholar] [CrossRef]
  17. Hassani, R.; Mahseradjian, J.; Tshibungu, T.; Karaagac, U. Evaluation of time-domain and phasor-domain methods for power system transients. Electr. Power Syst. Res. 2022, 212, 108335. [Google Scholar] [CrossRef]
  18. O’Rourke, C.J.; Qasim, M.M.; Overlin, M.R.; Kirtley, J.L., Jr. A geometric interpretation of reference frames and transformations: dq0, Clarke and Park. IEEE Trans. Energy Convers. 2019, 34, 2070–2083. [Google Scholar] [CrossRef] [Green Version]
  19. Mohanty, R.; Sahu, N.K.; Pradhan, A.K. Time-Domain Techniques for Line Protection Using Three-Dimensional Cartesian Coordinates. IEEE Trans. Power Deliv. 2022, 37, 3740–3751. [Google Scholar] [CrossRef]
  20. Casado-Machado, F.; Martinez-Ramos, J.L.; Barragán-Villarejo, M.; Maza-Ortega, J.M.; Rosendo-Macías, J.A. Reduced Reference Frame Transform: Deconstructing Three-Phase Four-Wire Systems. IEEE Access 2020, 8, 143021–143032. [Google Scholar] [CrossRef]
  21. Bellan, D.; Superti-Furga, G. Space-vector state-equation analysis of three-phase transients. J. Electr. Syst. 2018, 14, 188–198. [Google Scholar]
  22. Bellan, D. Probability Density Function of Three-Phase Ellipse Parameters for the Characterization of Noisy Voltage Sags. IEEE Access 2020, 8, 185503–185513. [Google Scholar] [CrossRef]
  23. Paranavithana, P.; Perera, S.; Koch, R.; Emin, Z. Global Voltage Unbalance in MV Networks Due to Line Asymmetries. IEEE Trans. Power Deliv. 2009, 24, 2353–2360. [Google Scholar] [CrossRef]
Figure 1. Example of trajectory of the Cartesian space vector v a b c in the orthogonal abc space. If no specific assumptions are made, the trajectory is not within a plane.
Figure 1. Example of trajectory of the Cartesian space vector v a b c in the orthogonal abc space. If no specific assumptions are made, the trajectory is not within a plane.
Energies 15 09122 g001
Figure 2. Example of trajectory of the Cartesian space vector v a b c in the orthogonal abc space in case of sinusoidal v a ,   v b ,   v c and null zero component (i.e., v 0 = 0 ). In this case, the trajectory is lying on the α β plane with equation a + b + c = 0 and is defined by the unit vectors e ^ α ,   e ^ β .
Figure 2. Example of trajectory of the Cartesian space vector v a b c in the orthogonal abc space in case of sinusoidal v a ,   v b ,   v c and null zero component (i.e., v 0 = 0 ). In this case, the trajectory is lying on the α β plane with equation a + b + c = 0 and is defined by the unit vectors e ^ α ,   e ^ β .
Energies 15 09122 g002
Figure 3. Trajectory of the Cartesian space vector v a b c in the orthogonal abc space in the case of sinusoidal waves v a ,   v b ,   v c . The black line ellipse refers to the case of the null zero component (see Figure 2), whereas the blue line ellipse refers to the case V 0 = V p / 2 . The zero component V 0 results in an elliptical trajectory lying on the plane defined by the vectors g 1 , g 2 instead of the α β plane.
Figure 3. Trajectory of the Cartesian space vector v a b c in the orthogonal abc space in the case of sinusoidal waves v a ,   v b ,   v c . The black line ellipse refers to the case of the null zero component (see Figure 2), whereas the blue line ellipse refers to the case V 0 = V p / 2 . The zero component V 0 results in an elliptical trajectory lying on the plane defined by the vectors g 1 , g 2 instead of the α β plane.
Energies 15 09122 g003
Figure 4. Three-phase network (a) leading to a first-order circuit in both the space-vector domain (b) and the zero-component domain (c).
Figure 4. Three-phase network (a) leading to a first-order circuit in both the space-vector domain (b) and the zero-component domain (c).
Energies 15 09122 g004
Figure 5. Trajectory of the Cartesian space vector i a b c in the case E n = 0 . The black line shows the steady-state trajectory for t < 0 , whereas the red line shows the transient behavior, starting at t = 0 , converging towards the new steady state. The black steady-state curve lies on the α β plane because i 0 = 0 for t < 0 , whereas the new steady-state trajectory lies on the new plane g 1 g 2 .
Figure 5. Trajectory of the Cartesian space vector i a b c in the case E n = 0 . The black line shows the steady-state trajectory for t < 0 , whereas the red line shows the transient behavior, starting at t = 0 , converging towards the new steady state. The black steady-state curve lies on the α β plane because i 0 = 0 for t < 0 , whereas the new steady-state trajectory lies on the new plane g 1 g 2 .
Energies 15 09122 g005
Figure 6. Trajectory of the space vector i ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n = 0 . In this case, the space-vector trajectory of the new steady state is circular.
Figure 6. Trajectory of the space vector i ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n = 0 . In this case, the space-vector trajectory of the new steady state is circular.
Energies 15 09122 g006
Figure 7. Transient behavior of the zero-component current i 0 t responsible for driving the Cartesian space vector i a b c in Figure 5 from the α β plane to the g 1 g 2 plane.
Figure 7. Transient behavior of the zero-component current i 0 t responsible for driving the Cartesian space vector i a b c in Figure 5 from the α β plane to the g 1 g 2 plane.
Energies 15 09122 g007
Figure 8. Trajectory of the Park space vector i ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n = 0 . The straight blue line in (a) shows the steady-state phasor at the end of the transient, i.e., the positive-sequence phasor.
Figure 8. Trajectory of the Park space vector i ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n = 0 . The straight blue line in (a) shows the steady-state phasor at the end of the transient, i.e., the positive-sequence phasor.
Energies 15 09122 g008
Figure 9. Trajectory of the Cartesian space vector i a b c in the case E n 0 . The black line shows the steady-state trajectory for t < 0 , whereas the red line shows the transient behavior, starting at t = 0 , converging towards the new steady state. The black steady-state curve lies on the α β plane because i 0 = 0 for t < 0 , whereas the new steady-state trajectory lies on the new plane g 1 g 2 .
Figure 9. Trajectory of the Cartesian space vector i a b c in the case E n 0 . The black line shows the steady-state trajectory for t < 0 , whereas the red line shows the transient behavior, starting at t = 0 , converging towards the new steady state. The black steady-state curve lies on the α β plane because i 0 = 0 for t < 0 , whereas the new steady-state trajectory lies on the new plane g 1 g 2 .
Energies 15 09122 g009
Figure 10. Trajectory of the space vector i ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n 0 . In this case, the space-vector trajectory of the new steady state is elliptical.
Figure 10. Trajectory of the space vector i ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n 0 . In this case, the space-vector trajectory of the new steady state is elliptical.
Energies 15 09122 g010
Figure 11. Trajectory of the Park space vector i ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n 0 . The straight blue line in (a) shows the positive-sequence phasor of the steady-state solution. In this case, the new steady-state trajectory is a circular trajectory (centered on the positive-sequence phasor) corresponding to the negative-sequence phasor rotating at double (negative) angular frequency 2 ω . This is apparent in (b) where, according to the theory, both a transient and a steady-state component at 2 ω can be distinguished.
Figure 11. Trajectory of the Park space vector i ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n 0 . The straight blue line in (a) shows the positive-sequence phasor of the steady-state solution. In this case, the new steady-state trajectory is a circular trajectory (centered on the positive-sequence phasor) corresponding to the negative-sequence phasor rotating at double (negative) angular frequency 2 ω . This is apparent in (b) where, according to the theory, both a transient and a steady-state component at 2 ω can be distinguished.
Energies 15 09122 g011
Figure 12. Three-phase network (a) leading to a second-order circuit in the space-vector domain (b) and a first-order circuit in the zero-component domain (c).
Figure 12. Three-phase network (a) leading to a second-order circuit in the space-vector domain (b) and a first-order circuit in the zero-component domain (c).
Energies 15 09122 g012
Figure 13. Trajectory of the Cartesian space vector i a b c in the case E n = 0 . The black line shows the steady-state trajectory for t < 0 , whereas the red line shows the transient behavior, starting at t = 0 , converging towards the new steady state. Both the black steady-state trajectory and the transient trajectory lie on the α β plane because i 0 = 0 .
Figure 13. Trajectory of the Cartesian space vector i a b c in the case E n = 0 . The black line shows the steady-state trajectory for t < 0 , whereas the red line shows the transient behavior, starting at t = 0 , converging towards the new steady state. Both the black steady-state trajectory and the transient trajectory lie on the α β plane because i 0 = 0 .
Energies 15 09122 g013
Figure 14. Trajectory of the space vector i ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n = 0 . In this case, the space-vector trajectory of the new steady state is circular.
Figure 14. Trajectory of the space vector i ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n = 0 . In this case, the space-vector trajectory of the new steady state is circular.
Energies 15 09122 g014
Figure 15. Trajectory of the Park space vector i ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n = 0 . The straight blue line in (a) shows the steady-state phasor at the end of the transient, i.e., the positive-sequence phasor.
Figure 15. Trajectory of the Park space vector i ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n = 0 . The straight blue line in (a) shows the steady-state phasor at the end of the transient, i.e., the positive-sequence phasor.
Energies 15 09122 g015
Figure 16. Trajectory of the Cartesian space vector v a b c in the case E n = 0 . The red line shows the transient behavior, starting at t = 0 , converging towards the new steady state on the g 1 g 2 plane.
Figure 16. Trajectory of the Cartesian space vector v a b c in the case E n = 0 . The red line shows the transient behavior, starting at t = 0 , converging towards the new steady state on the g 1 g 2 plane.
Energies 15 09122 g016
Figure 17. Trajectory of the space vector v ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n = 0 . In this case, the space-vector trajectory of the new steady state is circular.
Figure 17. Trajectory of the space vector v ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n = 0 . In this case, the space-vector trajectory of the new steady state is circular.
Energies 15 09122 g017
Figure 18. Transient behavior of the zero-component voltage v 0 t responsible for driving the Cartesian space vector v a b c in Figure 16 from the axes’ origin to the g 1 g 2 plane instead of the α β plane.
Figure 18. Transient behavior of the zero-component voltage v 0 t responsible for driving the Cartesian space vector v a b c in Figure 16 from the axes’ origin to the g 1 g 2 plane instead of the α β plane.
Energies 15 09122 g018
Figure 19. Trajectory of the Park space vector v ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n = 0 . The straight blue line in (a) shows the steady-state phasor at the end of the transient, i.e., the positive-sequence phasor.
Figure 19. Trajectory of the Park space vector v ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n = 0 . The straight blue line in (a) shows the steady-state phasor at the end of the transient, i.e., the positive-sequence phasor.
Energies 15 09122 g019
Figure 20. Trajectory of the Cartesian space vector i a b c in the case E n 0 . The black line shows the steady-state trajectory for t < 0 , whereas the red line shows the transient behavior, starting at t = 0 , converging towards the new steady state. Both the black steady-state trajectory and the transient trajectory lie on the α β plane because i 0 = 0 .
Figure 20. Trajectory of the Cartesian space vector i a b c in the case E n 0 . The black line shows the steady-state trajectory for t < 0 , whereas the red line shows the transient behavior, starting at t = 0 , converging towards the new steady state. Both the black steady-state trajectory and the transient trajectory lie on the α β plane because i 0 = 0 .
Energies 15 09122 g020
Figure 21. Trajectory of the space vector i ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n 0 . In this case, the space-vector trajectory of the new steady state is elliptical.
Figure 21. Trajectory of the space vector i ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n 0 . In this case, the space-vector trajectory of the new steady state is elliptical.
Energies 15 09122 g021
Figure 22. Trajectory of the Park space vector i ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n 0 . The straight blue line in (a) shows the positive-sequence phasor where the steady-state circular trajectory is centered.
Figure 22. Trajectory of the Park space vector i ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n 0 . The straight blue line in (a) shows the positive-sequence phasor where the steady-state circular trajectory is centered.
Energies 15 09122 g022
Figure 23. Trajectory of the Cartesian space vector v a b c in the case E n 0 . The red line shows the transient behavior, starting at t = 0 , converging towards the new steady state on the g 1 g 2 plane.
Figure 23. Trajectory of the Cartesian space vector v a b c in the case E n 0 . The red line shows the transient behavior, starting at t = 0 , converging towards the new steady state on the g 1 g 2 plane.
Energies 15 09122 g023
Figure 24. Trajectory of the space vector v ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n 0 . In this case, the space-vector trajectory of the new steady state is elliptical.
Figure 24. Trajectory of the space vector v ¯ t on the α β plane (a) and time-domain behavior of the corresponding α and β components (b) in the case E n 0 . In this case, the space-vector trajectory of the new steady state is elliptical.
Energies 15 09122 g024
Figure 25. Trajectory of the Park space vector v ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n 0 . The straight blue line in (a) shows the positive-sequence phasor where the steady-state circular trajectory is centered.
Figure 25. Trajectory of the Park space vector v ¯ d q t on the d q plane (a) and time-domain behavior of the corresponding d and q components (b) in the case E n 0 . The straight blue line in (a) shows the positive-sequence phasor where the steady-state circular trajectory is centered.
Energies 15 09122 g025
Figure 26. Three-bus MV power system. Loads are given in terms of apparent power and lagging power factor (pf).
Figure 26. Three-bus MV power system. Loads are given in terms of apparent power and lagging power factor (pf).
Energies 15 09122 g026
Figure 27. Topology of the space-vector equivalent circuit of the three-bus power system reported in Figure 26.
Figure 27. Topology of the space-vector equivalent circuit of the three-bus power system reported in Figure 26.
Energies 15 09122 g027
Figure 28. Trajectory of the Cartesian space vector of the generator currents. The trajectory lies on the α β plane because the zero component of the currents is null.
Figure 28. Trajectory of the Cartesian space vector of the generator currents. The trajectory lies on the α β plane because the zero component of the currents is null.
Energies 15 09122 g028
Figure 29. Space vector of the generator currents on the α β plane (a) and time-domain behavior of the corresponding α and β components (b).
Figure 29. Space vector of the generator currents on the α β plane (a) and time-domain behavior of the corresponding α and β components (b).
Energies 15 09122 g029
Figure 30. Park space vector of the generator currents on the d q plane (a) and time-domain behavior of the corresponding d and q components (b). The circle in (a) corresponds to the negative-sequence component in the steady-state solution.
Figure 30. Park space vector of the generator currents on the d q plane (a) and time-domain behavior of the corresponding d and q components (b). The circle in (a) corresponds to the negative-sequence component in the steady-state solution.
Energies 15 09122 g030
Figure 31. Time-domain behavior of the zero-component current.
Figure 31. Time-domain behavior of the zero-component current.
Energies 15 09122 g031
Figure 32. Trajectory of the Cartesian space vector of the generator currents. The trajectory lies on the g 1 g 2 plane instead of the α β plane because of the zero component of the generator currents.
Figure 32. Trajectory of the Cartesian space vector of the generator currents. The trajectory lies on the g 1 g 2 plane instead of the α β plane because of the zero component of the generator currents.
Energies 15 09122 g032
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bellan, D. Analytical Investigation of the Properties of Transients in Unbalanced Three-Phase Four-Wire Networks. Energies 2022, 15, 9122. https://doi.org/10.3390/en15239122

AMA Style

Bellan D. Analytical Investigation of the Properties of Transients in Unbalanced Three-Phase Four-Wire Networks. Energies. 2022; 15(23):9122. https://doi.org/10.3390/en15239122

Chicago/Turabian Style

Bellan, Diego. 2022. "Analytical Investigation of the Properties of Transients in Unbalanced Three-Phase Four-Wire Networks" Energies 15, no. 23: 9122. https://doi.org/10.3390/en15239122

APA Style

Bellan, D. (2022). Analytical Investigation of the Properties of Transients in Unbalanced Three-Phase Four-Wire Networks. Energies, 15(23), 9122. https://doi.org/10.3390/en15239122

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