Received 24 February 2016; accepted 20 March 2016; published 23 March 2016
Harmonic in power system caused by highly non-linear devices degrades its performance. Forced commutated VSCs are the main building block for low and medium power application. Due to recent development in the semiconductor technology and availability of high power switches e.g. Insulated Gate Bipolar Transistor (IGBT) and Gate Turn Transistor (GTO) have widespread acceptance in for high power VSC’s, which are used for FACTS controllers. The power system harmonic analysis is the process of calculating the magnitude and phase of fundamental and higher order harmonic of system signals. The generation of harmonics in modern power system is due to large size of power converter. To reduce the harmonics in the system filter, modern switching patterns are used. The increasing prevalence of flexible AC transmission system (FACTS) devices makes these accurate model devices essential. One attracting method for modeling the steady state performance of these devices is frequency domain analysis  . Harmonic phasor contains both positive and negative frequency terms for phase dependence. FACTS devices are characterized by their switching nature. One of the important FACTS devices is SSSC. SSSC is a series connected device that is able to provide active or reactive power support at the network location  . To understand the harmonics interaction between the SSSC based on VSC and utility system, there is a need of appropriate model.
VSC based SSSC acts as source of harmonic current injection into the system and also interacts with harmonic distortions present within the system. Various modulation techniques and topologies are used for SSSC to effectively reduce generated harmonics including PWM techniques, multi-module PWM techniques, selective harmonic elimination, and multi-level topology.
Previous work conducted in dynamic harmonic domain has been primarily focus on modeling PWM, multi- module, selective harmonic elimination based SSSC  -  . The power quality index can be assessed directly from the DHD modeling. Active, reactive, apparent powers as well as power factor are some important power quality indices. For the linear circuit the indices are defined in term of fundamental frequency whereas for nonlinear circuit, when nonlinear elements are present in a circuit such as electronics devices, given on basic of Fourier coefficient or given in terms of Total Harmonic Distortion (THD). For assessment of accurate power quality indices precise calculation of harmonic component is needed during transient period. Other authors have already made their contribution regarding modeling of FACTS devices and power system element using DHD method   .
There are two methods of DHD modeling. The first one is that direct mathematical mapping of all system equations and input is in frequency domain; and the other method that does not mapped all equations and input is in frequency domain gives more accurate result during transients. In second method, transient occurs due to change in circuit parameters. Second method is more suitable for calculation of instantaneous power quality indices, protection analysis and real time application all through it is violating causality and spurious dynamics result provided enough harmonic consideration  .
To extend the results obtained by other authors in the modeling of SSSC FACTS device, the proposed SVPWM VSC based SSSC model is developed in order to be able to obtained the evolution in the of harmonic components of the SSSC signal. Dynamic harmonic analysis of PWM based SSSC is shown in  as well as multi-pulse SSSC which is shown in  . This paper presents the dynamic analysis of SVPWM based SSSC during steady state and dynamic condition. This gives the information about harmonic indices during transient condition also. This can be used for control system design for space vector based SSSC as space vector modulation scheme is more suitable for digital implementation. The paper is organized as follows. The second section provides fundamentals of DHD method. The third section is fundamental of space vector based pulse width modulation based voltage source converter; this arrangement is used in this paper. The fourth section deals with DHD modeling of SVPWM VSC used for SSSC.
2. Dynamic Harmonic Domain Basic
The system given by Ordinary Differential Equations (ODE) can be transform to an alternative arrangement called the Dynamic Harmonic Domain, based on the approximation of the system by Fourier series over a period of the fundamental frequency  .
Linear periodic system (LTP) is converted into linear time domain (LTI) system using dynamic harmonic domain (DHD).
Consider a LTP system is given by
where A(t) is given as
where h is the highest harmonics of interest and ω0 is the fundamental frequency of the system.
Equation (1) can be represented in the Fourier series, the ordinary differential equation to an alternative DHD representation. Therefore
Variable given in the Equation (1) can transform into the vectors in the Equation (3) with their coefficient related to the harmonic components of their instantaneous signals as,
S is the matrix of differential and given by
The matrices A, B, C and D has Toeplitz structure and their time domain counterpart such that
The steady state solution is obtained directly from Equation (3), considering () gives,
The matrices A, B, C and D are constant and input U is also constant. The solution of X and Y is obtained from Equation (4). The solution of Equation (4) can be used to initialize the DHD simulation in time domain.
3. Space Vector Modulation
Principles of SVPWM
SVPWM is based in such a way that there are only two independent variables in a three-phase voltage system. We use orthogonal coordinates to represent the 3-phase voltage in the phasor diagram. A three-phase voltage vector represented by complex space vector as in Equation (1) neglecting zero sequence components  
SVM use the combinations of switching states to approximate the locus of Vref. In α-β plane a hexagon centred at origin of αβ plane, identifies the space vectors shown in Figure 1. Space vector plane of Figure 1 is divide into six sectors. Each sector covers the space corresponding to 60˚.
The distinct possible switching states of the 2-level VSC are represented as eight voltage vectors, out of which six are active states (V1 - V6) and two are null states (V0, V7). The active states contribute output line voltage as +Vdc or −Vdc, where as null states does not contribute any output voltage for VSC. The eight voltage vectors are tabulated in Table 1 denotes ON state of the switch and 0 denotes OFF state of the switch.
The reference vector is synthesized by the three adjacent switching vectors. For example, when Vref falls into sector I as shown in Figure 2, it can be synthesized by V1, V2 and V0. The volt second balance equation is
Figure 1. Switching vector of 2-level converter in αβ plane.
Figure 2. Representation of reference vector.
Table 1. 2-level inverter voltage vectors voltage vectors.
where Ts is the period of the switching cycle, T1 and T2 are the switching times of the vectors V1 and V2. T1 and T2 are calculated as
Similar calculation is applied to sector II to VI Vector V8 can be used in place of V7. The choice is based on the requirement to minimize average number of switching per cycle.
The maximum value of Vref is obtain when θ = 30˚ and Vref is given by
This is the maximum value of line to line voltage injected by the converter. The maximum magnitude of Vref is also the radius of circle inscribed in th hexagon shown in Figure 1. The Square wave converter generates a space vector of magnitude Vdc the maximum value of the modulation index as
4. Dynamic Harmonic Domain Modelling of SVPWM Based VSC
SSSC is important series FACT device which inject voltage into the transmission line. It is combination of voltage source converter, capacitor and coupling transformer. Primary of transformer is connected to the VSC and secondary is connected to the transmission line. If the voltage of VSC is in quadrature of line current then some reactive power exchange is take place and if voltage is not in quadrature of line current active power exchange is take place. DHD model of SSSC is presented in [Bharat] considering selective harmonic elimination method. In preceding discussion we extent this model considering space vector modulation techniques. Figure 3 shows the equivalent circuit of SSSC connected between the transmission line. AC side of SSSC shows the sinusoidal three phase voltage source connected in series. Re + jXe shows the resistance and impedance of coupling transformer. the three-phase voltages and currents on the AC side of the SSSC are VABCl(t) and iABC(t), respectively, and can be expressed in terms of the DC side voltage vdc(t), DC side current i1(t) and the switching functions as
vABC1(t) and iABC(t) are three phase voltage and current vectors given by:
In Equation (12), ps(t) and qs(t) are transformation vectors, which are given by:
The state equation can be written for SSSC are:
Figure 3. Static series synchronus compensator (SSSC).
The linear time periodic Equations (13) & (15) in matrix form will be written as:
The Equation (16) can be transfer to linear time invariant equation considering theory of DHD analysis as:
The initial condition are obtained by considering derivatives of state variable as zero gives:
Switching Vectors in Harmonic Domain
The General switching function is obtained in time domain for space vector modulation. The harmonic content in switching function is given by Fourier series.
where Sa, Sb, Sc are switching function obtained by using SVPWM algorithm. The line switching vector is defined as
The switching vector for harmonic domain is defined as
5. Simulation Result
In order to access dynamic harmonics response including power quality indices, the per-phase inductive reactance and resistance of the coupling transformer and the capacitance of the dc capacitor are R = 0.04 Ω, L = 0.2 mH and C = 5000 μF, respectively.
Under steady state conditions the bus per phase voltages VR and VS in volts at 50 Hz are
Assume disturbances in the voltages starting at 0.04 seconds and lasting for 0.005 seconds. During disturbances voltages on the VS bus is half of the original value. The simulation was started at t0 = 0 seconds with final time tf = 0.1 seconds and an integration time step is 0.001 s. 50 harmonics are considered. System is simulated using MATLAB software.
Figure 4 exhibits the dynamic behaviour of only the 1st, 3rd, 5th and 7th harmonic components with time for the phase-a current of the SSSC and DC capacitor voltage magnitudes. It clearly shows that behaviour of these variables from the onset of the disturbance until its ending, and their post-disturbance variations. The effectiveness of SSSC in bringing the system back to original steady state is also shown in these plots.
SSSC’s dynamic power quantities of all the three phases are shown in Figure 5. Dynamic power quantities of the SSSC such as RMS current Irms (t), apparent power S(t), active power P(t), reactive power Q(t), and distortion power D(t) in all the three phases are of interest for power quality assessments and control. Their changes with time during and after the disturbance are shown in the plots.
The total harmonic distortion measurement is important information for assessment of the device under different conditions. The total harmonic distortions (THD) of the voltage and current of the SSSC in all the three phases are shown in Figure 6. These plots show that, during the disturbance, the THD changes considerably due to the dynamic response of the harmonic components. The voltage and the current THD of the SSSC in all three phases are shown in Figure 6.
This paper presents the space vector based switching strategy for a SSSC that utilizes the voltage source converter to minimize the harmonic at the point of common contact. The linear time periodic equation is converted into linear time invariant system which is done using dynamic harmonic domain for calculation of harmonic interference in the system during transient and evaluated based dynamic harmonic domain algorithms using MATLAB code.
The proposed model is used to calculate harmonic interference produced by space vector based SSSC during
Figure 4. SSSC voltage and current. (a) Current at SSSC terminal phase A; (b) Voltage on DC capacitor.
Figure 5. SSSC terminal electric quantities. (a) RMS current; (b) Apparent power; (c) Distorted power; (d) RMS voltage on capacitor; (e) Active power; (f) Reactive power.
Figure 6. THD in voltage and current. (a) THD in voltage; (b) THD in current.
steady state and during disturbances. This gives accurate result and information about harmonics indices during transient operation of system as compared to time domain simulation which is important for designing control system for the system.
 Collins, C.D., Bathurst, G.N., Watson, N.R. and Wood, A.R. (2005) Harmonic Domain Approach of STATCOM Modelling. IEE Proceeding Generation Transmission, Distribution, 152, 194-200.
 Yahyaie, F. and Lehn, P.W. (2015) On Dynamic Evaluation of Harmonics Using Generalized Averaging Techniques. IEEE Transactions on Power Systems, 30, 2216-2224.
 Lian, K.L. and Lehn, P.W. (2006) Steady-State Solution of a Voltage-Source Converter with Full Closed-Loop Control. IEEE Transactions on Power Delivery, 21, 2071-2080.
 Jesus Rico, J., Madrigal, M. and Acha, E. (2003) Dynamic Harmonic Evolution Using the Extended Harmonic Domain. IEEE Transactions on Power Delivery, 18, 587-594.
 Vyakaranam, B., Madrigal, M., Eugenio Villaseca, F. and Rarick, R. (2010) Dynamic Harmonic Evolution in FACTS via the Extended Harmonic Domain Method. Power and Energy Conference at Illinois (PECI), 29-38.
 Tripura, P., Kishore Babu, Y.S. and Tagore, Y.R. (2011) Space Vector Pulse Width Modulation Schemes for Two-Level Voltage Source Inverter. ACEEE International Journal on Control System and Instrumentation, 2.
 Nanda, B. (2014) Total Harmonic Distortion of Dodecagonal Space Vector Modulation. International Journal of Power Electronics and Drive System (IJPEDS), 4, 308-313.