DIgSILENT PowerFactory

Technical Reference Documentation

Series Reactor
ElmSind
 DIgSILENT GmbH

Heinrich-Hertz-Str. 9
72810 - Gomaringen
Germany

T: +49 7072 9168 00

F: +49 7072 9168 88

http://www.digsilent.de
info@digsilent.de
r1284

Copyright ©2011, DIgSILENT GmbH. Copyright of this document belongs to DIgSILENT GmbH.
No part of this document may be reproduced, copied, or transmitted in any form, by any means
electronic or mechanical, without the prior written permission of DIgSILENT GmbH.

Series Reactor (ElmSind) 1

Contents

Contents

1 General Description 4

1.1 Input Parameter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

1.1.1 Input mode of the rating parameter . . . . . . . . . . . . . . . . . . . . . . 4

1.1.2 Short-circuit voltage uk <=> Impedance Zrea relation . . . . . . . . . . . 5

1.1.3 Copper Losses PCU <=> Short-circuit voltage ukr, real part relation . . . 6

1.1.4 Reactance Xrea <=> Inductance Lrea relation . . . . . . . . . . . . . . . 6

1.1.5 Resistance Rrea <=> Copper Losses PCU , Short-circuit voltage ukr, real
part relation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

1.1.6 Reactance Xrea <=> short-circuit voltage, uk and ukr, Impedance Zrea . 7

1.2 Load flow model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

1.2.1 AC-Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

1.2.2 DC-Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

1.3 Short-circuit model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

1.4 RMS Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

1.4.1 AC-Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

1.4.2 DC-Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

1.5 EMT Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8

1.5.1 AC-Model without Saturation . . . . . . . . . . . . . . . . . . . . . . . . . 8

1.5.2 AC-Model with Saturation . . . . . . . . . . . . . . . . . . . . . . . . . . . 9

1.5.3 DC-Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

2 Dynamic Simulation 15

2.1 RSM and EMT Simulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15

2.1.1 AC-Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15

2.1.2 DC-Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15

A Parameter Definitions 16

B Signal Definitions 16

List of Figures 17

Series Reactor (ElmSind) 2

Contents

List of Tables 18

Series Reactor (ElmSind) 3

1 General Description

1 General Description

Terminal i Terminal j
Terminal i Terminal j
Rrea Xrea
A Rrea Xrea
A A
A

Rrea Xrea
Rrea Xrea
B B
B B

Rrea Xrea
Rrea Xrea
C C
C C

Figure 1.1: Three Phase Model

Terminal i Terminal j
Terminal i Terminal j
I:bus1 Rrea Xrea I:bus2
I:bus1 Rrea Xrea I:bus2

U:bus2 U:bus2
U:bus2 U:bus2

Figure 1.2: Single Phase Model

The series reactor can be used to limit the short-circuit current.

A three phase and single phase model is available.

1.1 Input Parameter

The resistance and the reactance of the reactor can be entered at different input options:

• Short-circuit voltage uk in % and copper losses in kW or real part of the short-circuit

voltage in %
• Impedance (magnitude) in Ohm and copper losses in kW or real part of the short-circuit
voltage in %
• Reactance X and Resistance R in Ohm
• Inductance L in mH and Resistance in Ohm

1.1.1 Input mode of the rating parameter

For the definition of the rating parameter they are two different input modes available:

Series Reactor (ElmSind) 4

1 General Description

a) Rated Voltage in kV and Rated Power in MVA

b) Rated Voltage in kV and Rated Current in kA

For a three phase reactor is the rated current calculated as:

Sn
In = √ in kA
3 · Un

for a single phase reactor:

√
3 · Sn
In = in kA
Un

and for a DC single phase reactor:

Sn
In = in kA
Un

where,

• In : Rated Current in kA
• Un : Rated Voltage (Line-Line) of the reactor in kV
• Sn : Rated Power of the reactor in MVA

1.1.2 Short-circuit voltage uk <=> Impedance Zrea relation

For a three phase reactor:

uk U2
Zrea = · n in Ohm
100% Sn

For a single phase reactor:

uk Un2
Zrea = · in Ohm
100% 3 · Sn

• Zrea : Impedance magnitude of the reactor in Ohm

• uk : Short-circuit voltage in %
• Un : Rated Voltage (Line-Line) of the reactor in kV
• Sn : Rated Power of the reactor in MVA

Series Reactor (ElmSind) 5

1 General Description

1.1.3 Copper Losses PCU <=> Short-circuit voltage ukr, real part relation

PCU
ukr = · 100 in %
Sn · 1000

• ukr : Short-circuit voltage, real part in %

• PCU : Copper Losses in kW

• Sn : Rated Power of the reactor in MVA

1.1.4 Reactance Xrea <=> Inductance Lrea relation

Xrea
Lrea = · 1000 in mH
2 · π · fnom

• Lrea : Inductance of the reactor in mH

• Xrea : Reactance of the reactor in Ohm
• fnom : Nominal frequency of the grid in Hz

1.1.5 Resistance Rrea <=> Copper Losses PCU , Short-circuit voltage ukr, real part
relation

For a three phase reactor:

2 2
ukr Unom PCU Unom
Rrea = · = · 2 in Ohm
100% Snom 1000 Snom

For a single phase reactor:

2 2
ukr Unom PCU Unom
Rrea = · = · 2
in Ohm
100% 3 · Snom 1000 3 · Snom

• Rrea : Resistance of the reactor in Ohm

• PCU : Copper Losses in kW

• Unom : Rated Voltage (Line-Line) of the reactor in kV

• Snom : Rated Power of the reactor in MVA
• ukr : Short-circuit voltage, real part in %

Series Reactor (ElmSind) 6

1 General Description

1.1.6 Reactance Xrea <=> short-circuit voltage, uk and ukr, Impedance Zrea

For a three phase reactor:

√ 2
uk 2 − ukr2 Unom p
2 − R2
Xrea = · = Zrea rea in Ohm
100% Snom

For a single phase reactor:

√ 2
uk 2 − ukr2 Unom p
2 − R2
Xrea = · = Zrea rea in Ohm
100% 3 · Snom

• Xrea : Reactance of the reactor in Ohm

• Unom : Rated Voltage (Line-Line) of the reactor in kV
• Snom : Rated Power of the reactor in MVA
• uk : Short-circuit voltage, magnitude in %
• ukr : Short-circuit voltage, real part in %
• Zrea : Impedance magnitude of the reactor in Ohm

1.2 Load flow model

The series reactor is completely considered in the load flow calculation. A warning message
will be printed on the output window if the rated voltage differs more than 10% as the nominal
voltage of the bus bars. An error message will be printed on the output window if the rated
voltages differs more than 50%.

1.2.1 AC-Model

For the load flow calculation is the reactor modelled by it’s nominal frequency behaviour:

Ubus1 − Ubus2 = Ibus1 · (Rrea + jω · Lrea)

Ibus1 + Ibus2 = 0

The loading of the reactor is calculated as follow:

max(Ibus1 , Ibus2 )
loading = · 100 in %
In

• loading : Loading in %
• In : Rated current of the reactor in kA

Series Reactor (ElmSind) 7

1 General Description

• Ibus1 : Magnitude of the current at terminal i

• Ibus2 : Magnitude of the current at terminal j

For an unbalanced load flow calculation is the highest current of all phases used.

1.2.2 DC-Model

Under DC-conditions, a reactor is a closed circuit, hence only the resistance is considered:

Ubus1 − Ubus2 = Ibus1 · Rrea

Ibus1 + Ibus1 = 0

1.3 Short-circuit model

For the short-circuit calculation is the reactor modelled by it’s nominal frequency behaviour:

Ubus1 − Ubus2 = Ibus1 · (Rrea + jω · Lrea)

Ibus1 + Ibus2 = 0

1.4 RMS Model

1.4.1 AC-Model

The RMS-simulation model of the series capacitance is fully compatible to the load flow model.

1.4.2 DC-Model

In DC circuits, network transients are considered as well. The series reactor is hence modelled
by the differential equation of a resistance in series with a inductance:

d(Ibus1 (t))
Ubus1 (t) − Ubus2 (t) = Rrea · Ibus1 (t) + Lrea ·
dt
Ibus1 + Ibus2 = 0

1.5 EMT Model

1.5.1 AC-Model without Saturation

The series reactor is modelled by the differential equation of a resistance in series with an
inductance:

Series Reactor (ElmSind) 8

1 General Description

d(Ibus1:A (t))
Ubus1:A (t) − Ubus2:A (t) = Rrea:r · Ibus1:A (t) + Lrea:r ·
dt
d(Ibus1:B (t))
Ubus1:B (t) − Ubus2:B (t) = Rrea:s · Ibus1:B (t) + Lrea:s ·
dt
d(Ibus1:C (t))
Ubus1:C (t) − Ubus2:C (t) = Rrea:t · Ibus1:C (t) + Lrea:t ·
dt
Ibus1:A (t) + Ibus2:A (t) = 0
Ibus1:B (t) + Ibus2:B (t) = 0
Ibus1:C (t) + Ibus2:C (t) = 0

1.5.2 AC-Model with Saturation

The saturation of the reactance can be simulated in the EMT model.

The model supports the following options:

• Linear: no saturation considered

• Two slope: the saturation curve is approximated by a two linear slopes

• Polynomial: the saturation curve is approximated by a polynom of user-defined order. The
polynom fits asymptotically into the piecewise linear definition.
• Current/Flux values: the user inputs current-flux values as a sequence of points and se-
lects among a piecewise-linear or spline interpolation.

Linear
The input parameters are listed in Table 1.1.

Table 1.1: Basic data of the linear characteristic

Parameter Description Unit Symbol

xreapu Magnetizing reactance for unsaturated p.u. lunsat
conditions.
In p.u. values, the linear reactance is
equal to the reciprocal of the
magnetizing current (reactive part of
the exciting current).

No saturation is considered. In this case, the current is proportional to the flux:

ψM
iM =
lunsat

Where:

• iM is the magnetizing current in p.u.

Series Reactor (ElmSind) 9

1 General Description

• ψM is the magnetizing flux in p.u.

The p.u. values used for the definition of the saturation characteristic of the positive sequence
model are referred to the same base quantities mentioned in 1.5.2.

Two slope
Figure 1.3 shows the magnetizing current-flux curves for the two slope and polynomial char-
acteristics. The input parameters of the two slope saturation characteristic are listed in Table
1.2.

Figure 1.3: Two slope and polynomial saturation curves

Table 1.2: Basic data of the two-slope saturation characteristic

Parameter Description Unit Symbol

psi0 Knee-point Flux of asymptotic p.u. ψknee
piece-wise linear characteristic.
Typical value around 1.1 to 1.2 times
the rated flux.
xreapu Magnetizing reactance for unsaturated p.u. lunsat
conditions.
In p.u. values, the linear reactance is
equal to the reciprocal of the
magnetizing current (reactive part of
the exciting current).
xmair Magnetizing reactance for saturated p.u. lsat
conditions.

In the case of the two slope curve, the series reactor saturation is approximated by two linear
slopes, which are separated by the value of the knee flux. The first slope is given by the linear
reactance, while the second is given by the saturated reactance. The saturated reactance value
is fairly low compared with the unsaturated reactance.

Series Reactor (ElmSind) 10

1 General Description

ψM
iM = for ψM ≤ ψknee
lunsat
1
iM = (ψM − ψknee ) + iknee for ψM > ψknee
lsat
ψknee
iknee =
lunsat

Where

• iM is the magnetizing current in p.u.

• ψM is the magnetizing flux in p.u.

The p.u. values used for the definition of the saturation characteristic of the positive sequence
model are referred to the same base quantities mentioned in 1.5.2.

Polynomial
Figure 1.3 shows the magnetizing current-flux curves for the polynomial characteristic. The
input parameters are listed in Table 1.3.

Table 1.3: Basic data of the polynomial saturation characteristic

Parameter Description Unit Symbol

psi0 Knee-point Flux of asymptotic p.u. ψknee
piece-wise linear characteristic.
Typical value around 1.1 to 1.2 times
the rated flux.
xreapu Magnetizing reactance for unsaturated p.u. lunsat
conditions.
In p.u. values, the linear reactance is
equal to the reciprocal of the
magnetizing current (reactive part of
the exciting current).
xmair Magnetizing reactance for saturated p.u. lsat
conditions.
ksat Saturation exponent of polynomial - ksat
representation. Typical values are
9,13,15. The higher the exponent the
sharper the saturation curve.

In the polynomial type, the saturation is also represented by two linear slopes, but smoothed by
a polynomial function, in which the higher the exponent is, the sharper the saturation becomes.

Series Reactor (ElmSind) 11

1 General Description

(  ksat )
ψM ψM
iM = · 1+ for ψM ≤ ψsat
lunsat ψ0
1
iM = (ψM − ψsat ) + iknee for ψM > ψsat
lsat
ksat + 1
ψsat = · ψknee
ksat
ψknee
iknee =
lunsat

Where

• iM is the magnetizing current in p.u.

• ψM is the magnetizing flux in p.u.
• ψ0 is a parameter automatically calculated, so that the polynomial characteristic fits the
saturated reactance in full saturation and transits steadily into the piece-wise linear char-
acteristic at the knee flux point, in p.u.

This polynomial characteristic is always inside the corresponding linear representation. In full
saturation the polynomial characteristic is extended linearly. Compared to the two-slope curve, it
does not contain a singular point at the knee flux and therefore its derivate (magnetizing voltage)
is continuously defined.

The p.u. values used for the definition of the saturation characteristic of the positive sequence
model are referred to the same base quantities mentioned in 1.5.2.

Current/Flux values
The user can also define the saturation curve in terms of measured current-flux values and
select between a piecewise linear or spline interpolation.

The current-flux values in the table are peak values in p.u.

The p.u. values used for the definition of the saturation characteristic of the positive sequence
model are referred to the same base quantities mentioned in 1.5.2.

Base Quantities
The p.u. values used for the definition of the saturation characteristic of the positive sequence
model are referred to the following bases quantities:

• Ubase is the Rated Voltage of the series reactor (ucn) in kV

• Sbase is the Rated Power of the series reactor (Sn) in M V A
√ Sbase [M V A]
• Ibase [A] = 2· √ · 1000
3 · Ubase [kV ]
√
√ Ubase [kV ]/ 3
• Ψbase [V · s] = 2 · · 1000
2πf [Hz]
2
 
Ubase [kV ] 1
• Lbase [H] = ·
Sbase [M V A] 2πf [Hz]

Series Reactor (ElmSind) 12

1 General Description

Ignoring the Saturation

The saturation is ignored when any of the following conditions are met:

• The reactance xrea is zero

• It is a DC reactor
• The input signal Xin is connected

Considerations for 3-Limb Core

The ABC per unit flux values are calculated from the dq flux components (ψ0 is ignored in the
transformation):

 
  1 √0 1  
ψM :A 
1 3
 ψd
− 1 × 
 
ψM :B  = 
  
 2 2
√   ψq 
 1 3
ψM :C 0

− − 1
2 2

The zero-sequence current in p.u. is given by:

ψ0
i0 =
2πf · lunsat:r

Considerations for 5-Limb Core

The ABC per unit flux values are calculated from the dq0 flux components:

 
  1 √0 1  
ψM :A 
1 3  ψd 

− 1 × ψq 

ψM :B  = 
 
 2 2
√ 
 1 3
ψM :C ψ0

− − 1
2 2

The zero-sequence current in p.u. is given by:

i0 = 0

Equations
The magnetizing currents in p.u. are calculated according to the type of saturation selected:

iM :A (t) = f (ψM :A )
iM :B (t) = f (ψM :B )
iM :C (t) = f (ψM :C )

Series Reactor (ElmSind) 13

1 General Description

The equations for voltage and current per phase in actual values are given by:

Ibus1:A (t) = (iM :A + i0 ) · Ibase

Ibus1:B (t) = (iM :B + i0 ) · Ibase
Ibus1:C (t) = (iM :C + i0 ) · Ibase
1 d(Ψbus1:A (t))
Ubus1:A (t) − Ubus2:A (t) = Rrea:r · Ibus1:A (t) + ·
2πf dt
1 d(Ψbus1:B (t))
Ubus1:B (t) − Ubus2:B (t) = Rrea:s · Ibus1:B (t) + ·
2πf dt
1 d(Ψbus1:C (t))
Ubus1:C (t) − Ubus2:C (t) = Rrea:t · Ibus1:C (t) + ·
2πf dt
Ibus1:A (t) + Ibus2:A (t) = 0
Ibus1:B (t) + Ibus2:B (t) = 0
Ibus1:C (t) + Ibus2:C (t) = 0

1.5.3 DC-Model

The series reactor is modelled by the differential equation of a resistance in series with an
inductance:

d(Ibus1 (t))
Ubus1 (t) − Ubus2 (t) = Rrea · Ibus1 (t) + Lrea ·
dt
Ibus1 (t) + Ibus2 (t) = 0

Series Reactor (ElmSind) 14

2 Dynamic Simulation

2 Dynamic Simulation

2.1 RSM and EMT Simulation

2.1.1 AC-Model

Rin
Xin AC Model

Figure 2.1: Input/OutputRin

Definition of AC Series Reactor Model (RMS-Simulation)
Rin AC Model
Xin
Lin DC Model
2.1.2 DC-Model

Rin
Lin DC Model

Figure 2.2: Input/Output Definition of DC Series Reactor Model (RMS-Simulation)

Series Reactor (ElmSind) 15

B Signal Definitions

A Parameter Definitions

Table A.1: Series Reactor Parameters

Parameter Description Unit

loc name Name

outserv Out of Service
ucn Rated voltage kV
Sn Rated power MVA
Curn Rated current kA
pRating Thermal rating
systp System type
nphases Phases
uk Short-circuit voltage %
Zd Absolute impedance Ohm
Irea Inductance mH
xrea Reactance Ohm
Pcu Copper losses kW
ukr Short-circuit voltage %
rrea Resistance Ohm
maxLoad Thermal load limit %
iAstabint A-stable integration algorithm
psi0 Knee flux p.u.
xmair Saturated reactance p.u.
ksat Saturation exponent
satcur Current (peak) p.u.
satflux Flux (peak) p.u.
iInterPol Interpolation
smoothfac Smoothing factor %
Inom Rated current kA
xreapu Linear reactance p.u.

B Signal Definitions

Table B.1: Input/Output signals

Name Description Unit Type Model

Rin Resistance Input Ohm IN AC, DC (RMS, EMT)

Xin Reactance Input Ohm IN AC (RMS, EMT)
Lin Inductance Input Ohm IN DC (RMS, EMT)

Series Reactor (ElmSind) 16

List of Figures

List of Figures

1.1 Three Phase Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

1.2 Single Phase Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

1.3 Two slope and polynomial saturation curves . . . . . . . . . . . . . . . . . . . . 10

2.1 Input/Output Definition of AC Series Reactor Model (RMS-Simulation) . . . . . . 15

2.2 Input/Output Definition of DC Series Reactor Model (RMS-Simulation) . . . . . . 15

Series Reactor (ElmSind) 17

List of Tables

List of Tables

1.1 Basic data of the linear characteristic . . . . . . . . . . . . . . . . . . . . . . . . . 9

1.2 Basic data of the two-slope saturation characteristic . . . . . . . . . . . . . . . . 10

1.3 Basic data of the polynomial saturation characteristic . . . . . . . . . . . . . . . . 11

A.1 Series Reactor Parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

B.1 Input/Output signals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

Series Reactor (ElmSind) 18