Electromagnetic transient and

phasor domain hybrid simulation

and its application to detailed
FIDVR studies
Vijay Vittal – Ira A. Fulton Chair Professor
Qiuhua Huang – Graduate Student
John Undrill – Research Professor
Arizona State University
Work done in PSERC Project S-58

PSERC Webinar
October 4, 2016 1
 Introduction
 Significant changes in power systems
• Rapid development of • Distributed generation (PV,
renewable generation Wind)
• Shutdown of coal power • More dynamic loads (Air
plants conditioners, VFD motors)

• HVDC
• FACTS

Source: U.S. DOE 2

 Two main challenges in power system dynamic simulation

1. Modeling and representation of an increasing number of

power electronic devices in dynamic simulation
• Quasi-steady-state model or performance model in TS simulators
• Fast switching and control cannot be adequately represented

3
Two main challenges in power system dynamic simulation
2. Representation of distribution systems in power system
dynamic simulation
 Various load models: ZIP, motor + ZIP, CMPLDW 1) computational
limitation; 2) availability of distribution system data
 Developments of distributed generation, EV and storage significantly
change the behavior of distribution systems
• Past: the voltage profile gradient along the feeder was small
• Now the locations matter: A/C stalling, local volt/VAr support of DGs
• Load modeling  distribution system modeling

 Integrated transmission and distribution (T&D) systems

• EPRI proposed “the integrated grid” framework—for fully realizing the values of
distributed and central generation resources
• To analyze the increased interactions between T&D and to take advantage of them

4
Two main challenges in power system dynamic simulation

2. Representation of distribution systems in power

system dynamic simulation
Detailed 3-phase
Three-sequence, phasor representation down to
models primary feeder level

Three-sequence EMT or three-phase

transient stability dynamic simulation
5
 Development of Hybrid Simulation

• Three-phase, • Sequence,
POW model phasor model
• EMT simulator • TS simulator
Region of interest Remainder of the system
< 300 buses ~15000 buses
 Detailed system  External system

Thévenin
equivalent

Detailed system modeled in an EMT Simulator External system modeled in a TS Simulator

6
 OpenHybridSim: A new EMT-TS hybrid simulation tool
 A decoupled architecture
 Three-sequence TS simulation developed based on InterPSS
 Network equivalents:
• Three-phase Thévenin equivalent of the external system in EMT simulation
• Three-sequence current source as the equivalent of the detailed system in
three-sequence TS simulation
 TCP/IP socket communication for connecting two simulators
 A generic interface framework for integrating with different EMT
simulators, e.g., PSCAD, ATP-EMTP

7
 Equivalent of the external system in EMT simulation:
Three-phase Thévenin equivalent
Three- Three-phase
Three-phase
sequence Source
120-ABC Norton Thévenin
Norton transformation
equivalent equivalent
equivalent

120
1⁄
Y
N

1 1 1
1
1
0 0
0 0 ⁄

0 0
8
Equivalent of the detailed system in TS simulation: three-
sequence current source
• Three-sequence current source
• Seamlessly integrated into the network solution step of the
three-sequence TS simulation
• Obtained from boundary current injection waveforms
using FFT and 3-phase to 3-sequence transformation

The detailed system is represented by three-

FFT component in PSCAD
sequence current sources in TS simulation 9
 Three-sequence TS simulation

(1) (1)
I VTS ( t T )
EMT ( t ) x (t )  f ( x (t ), y((1)t ) )
x (t  T )  x (t )  x (t ) T
0  g1 ( x (t  T ), y((1)t  T ) , I EMT
(1)
(t ) )

(2)
I EMT ( t ) VTS(2)(t T )
0  g 2 ( y((2)
t  T ) , I (2)
EMT ( t ) )
I (0) VTS(0)(t T )
EMT ( t )
0  g 0 ( y((0)
t  T ) , I (0)
EMT ( t ) )

Negative sequence : =0

Zero sequence : =0
10
 Interaction protocol
• Interactions between the two simulators

I120
EMT (t )
V abc
T (t )

Serial Parallel

Combined

11
 Implementation of the Two Interaction Protocols
(1)
Serial protocol
(1)
I VTS ( t T )
x(t ) EMT ( t )
x (t )  f ( x(t ), y((1)t ) )
x(t  T )  x(t )  x (t ) T x(t  T )
V 120
TS (t  T )
0  g1 ( x(t  T ), y((1)t  T ) , I EMT
(1)
(t ) )
120
I EMT (t ) V abc
I (2)
V (2) T (t  T ) VTS120(t T )
VTS120(t ) EMT ( t ) TS ( t T )
0  g2 ( y (2)
( t  T ) ,I (2)
EMT ( t ) )
I (0) VTS(0)(t T )
120
EMT ( t )
0  g0 ( y (0)
,I (0)
)
I EMT (t )
( t  T ) EMT ( t )
V abc
T (t  T )

120 120
I EMT (t ) I EMT ( t T )

t x(t  T ) t  T
120
x(t ) V TS ( t )
120
VTS120(t ) I EMT (t )
VTS120(t T )
120
Parallel
I EMT (t )
V abc protocol
T (t )
120
120 I EMT ( t T )
I EMT ( t )

t t  T 12
 The fault-induced delayed voltage recovery problem
• What is FIDVR problem?

A 230 kV bus voltage profile

during a typical FIDVR event [1]

• Root cause:
Stalling and prolonged tripping of 1-ϕ residential air conditioner (A/C)
compressor motor
• Direct impact of distribution on transmission system
[1] D. N. Kosterev, A. Meklin, J. Undrill, B. Lesieutre, et al., "Load modeling in power system studies: WECC progress
update," in 2008 IEEE Power and Energy Society General Meeting, 2008, pp. 1-8.

13
The FIDVR Problem
• FIDVR has mainly been studied using CMPLDW and positive-sequence TS
simulation programs
• Limitations of the CMPLDW model
 Limited representation of distribution systems and DGs
 Performance-based 1- ϕ A/C compressor motor model
 Point-on-wave (POW) effects cannot be considered
 Not suitable for cases involving unbalanced conditions
• The issues above can be overcome by hybrid simulation

Load Bus

System Bus Low-side M Motor A

(230, 115, 69kV) Bus Feeder
1:T Equiv.
Rfdr +j Xfdr M Motor B
jXxf

Fb Bfdr M Motor C

Bss
(1-Fb) Bfdr M Motor D

UVLS

UFLS Electronic
DG
Static

CMPLDW Model Structure The point-on-wave effects on A/C stalling

14
Application of EMT-TS hybrid simulation to FIDVR study on the
WECC system
The WECC system Buses with a large percentage of
Transmission 1-Ф A/C load
Buses lines Generators Loads
• Bus 24151
15750 13715 3074 7787
• Bus 24138
Summary of the detailed system
Total number
238
of buses
500 kV 7
Number of 230 kV 37
buses of 161 kV 3
different 115 kV 68
voltage levels 92 kV 18
<= 66kV 105
Total Load 11.9 GW
Interface
One-line diagram of the study region 8
buses
15
Set up of the detailed system

50% 1-Ф A/C motor [2] 1 3

Z Z
4 feeder 4 feeder
25% 3-Ф induction motor
25% constant impedance
Equivalent feeder model[3]
[2] Y. Liu, V. Vittal, J. Undrill, and J. H. Eto, "Transient Model of Air-Conditioner Compressor Single Phase Induction Motor," IEEE
Transactions on Power Systems, vol. 28, pp. 4528-4536, 2013.
[3] William H. Kersting, Distribution System Modeling and Analysis (Second Edition), CRC Press, 2006, p.52-54

16
Case A: fault POW at voltage zero crossing
SLG

Voltage (pu)
Voltage (pu)
Voltage (pu)

Fault point
Voltage (pu)

Measurement
location

17
Case A: responses of A/C motors phase A phase B phase C

A/C terminal voltage (pu)

A/C terminal voltage (pu)

A/C terminal voltage (pu)

1 1 1

0.5 0.5 0.5

0 0 0
0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6
Time (s) Time (s) Time (s)
A/C Motor Speed (pu)

A/C Motor Speed (pu)

A/C Motor Speed (pu)

1 1 1

0.5 0.5 0.5

0 0 0
0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6
Time (s) Time (s) Time (s)
15 15 15
Real power
Real and reactive

Real and reactive

Real and reactive

Real power Real power
power (kW/kVar)

power (kW/kVar)

power (kW/kVar)
10 Reactive power 10 Reactive power 10 Reactive power

5 5 5

0 0 0

-5 -5 -5
0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6
Time (s) Time (s) Time (s)
Electrical/Mechanical

Electrical/Mechanical

Electrical/Mechanical
Electrical torqure Electrical torqure Electrical torqure
Torqure (N*m)

Torqure (N*m)

Torqure (N*m)
50 Mechanical torque
50 Mechanical torque
50 Mechanical torque

0 0 0

-50 -50 -50

0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6
Time (s) Time (s) Time (s)
18
 Voltage (pu) Voltage (pu) Voltage (pu) Voltage (pu) Voltage (pu)
Case A: phasor voltage magnitudes

19
 The Point-On-Wave Effects
Case B: POW at the peak of phase A voltage waveform (90 degrees)

20
Case B: responses of A/C motors

A/C terminal voltage (pu)

A/C terminal voltage (pu)

A/C terminal voltage (pu)

phase A phase B phase C

1 1 1

0.5 0.5 0.5

0 0 0
0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6
Time (s) Time (s) Time (s)
A/C Motor Speed (pu)

A/C Motor Speed (pu)

A/C Motor Speed (pu)

1 1 1

0.5 0.5 0.5

0 0 0
0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6
Time (s) Time (s) Time (s)
15 15 15
Real and reactive

Real and reactive

Real and reactive

power (kW/kVar)

power (kW/kVar)

power (kW/kVar)
Real power Real power Real power
10 Reactive power 10 Reactive power 10 Reactive power

5 5 5

0 0 0

-5 -5 -5
0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6
Time (s) Time (s) Time (s)
Electrical/Mechanical

Electrical/Mechanical

Electrical/Mechanical
Electrical torqure Electrical torqure Electrical torqure
Torqure (N*m)

Torqure (N*m)

Torqure (N*m)
50 Mechanical torque 50 Mechanical torque 50 Mechanical torque

0 0 0

-50 -50 -50

0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6 0.6 0.8 1 1.2 1.4 1.6
Time (s) Time (s) Time (s)
21
Case B: phasor voltage magnitudes

22
PSLF simulation results using CMPLDW

Vstall = 0.6 pu
Tstall = 0.033 s

23
 Application of hybrid simulation to power systems interfaced
with a LCC-HVDC system

 Test case: IEEE 39 Bus system with a LCC HVDC infeed

• HVDC system
 CIGRE HVDC model
 Rated power: 1000 MW
 The inverter is connected to bus 39
• EMT simulation
 The whole system
time step : 50 µs
• Hybrid simulation
 The HVDC system is modeled in detail in
PSCAD
 EMT time step : 50 µs
 TS time step: 5 ms
• Scenario: A SLG fault is applied at
bus 39 at 3.0 s, cleared after 4 cycles
24
 Application of hybrid simulation to power systems
interfaced with a LCC-HVDC system
1.5
10 EM T-TS EMT-TS
Phase A
Current(kA)

PSCAD 1 PSCAD
5

0 0.5

V dc (pu)
0.3
2.95 3 3.05 3.1 3.15 3.2 0 0.2
10
0.1
EM T-TS
Phase B -0.5 0
Current(kA)

PSCAD
5
3.07 3.08 3.09
-1
0 2.95 3 (a) 3.05 3.1 3.15 3.2
4
2.95 3 3.05 3.1 3.15 3.2 EMT-TS

10 2 PSCAD

Idc (pu)
EM T-TS
Phase C
Current(kA)

PSCAD
5
0
0
-2
2.95 3 3.05 3.1 3.15 3.2 2.95 3 3.05 3.1 3.15 3.2
Time (s)
(b) Time(s)

Three phase current flowing into bus 39 The DC voltage (Vdc) and current
from the HVDC inverter (Idc) of the HVDC inverter
25
 Application of hybrid simulation to power systems
interfaced with a LCC-HVDC system
EM T-TS PSCAD
5 1.5
Current (kA)

PSCAD
0 EM T-TS
1

Real power (pu)

-5
2.95 3 3.05 3.1 3.15 3.2
5
Current (kA)

0.5 1.2
0 1
-5 0.8
2.95 3 3.05 3.1 3.15 3.2 0
10 0.6
Current (kA)

3 3.05 3.1 3.15

0 -0.5
3 3.2 3.4 3.6 3.8 4 4.2 4.4 4.6 4.8 5
-10 time (s)
2.95 3 3.05 3.1 3.15 3.2 (a)
5 2
Current (kA)

0 1 PSCAD

0.8 EM T-TS
Commutation

Reactive power (pu)

-5
1.5
2.95 3 3.05 3.1 3.15 3.2
0.6
5
failure
Current (kA)

1 0.4
0

-5 3.1 3.15 3.2

2.95 3 3.05 3.1 3.15 3.2 0.5
5
Current (kA)

0
0
-5 3 3.2 3.4 3.6 3.8 4 4.2 4.4 4.6 4.8 5
2.95 3 3.05 3.1 3.15 3.2 time (s)
Time (s)
(b)

The current flows in the 6-pulse Power of the generator at bus 30

bridge of the inverter 26
 Application of hybrid simulation to power systems interfaced
with a LCC-HVDC system
Simulation differences with reference to full-blown EMT simulation
Average Maximum
Monitored parameters
difference/pu difference/pu
DC Current of the inverter 0.024 0.224
DC Voltage of the inverter 0.029 0.524
Three phase current into the network
0.049 0.351
at the inverter
Three phase voltages of bus 39 0.022 0.155

Computational times of hybrid simulation and EMT simulation

Simulation method Total computation time*

EMT simulation using PSCAD 352 s

EMT-TS hybrid simulation 81 s
* 5-second simulation
27
 Applied to power systems interfaced with HVDC

 IEEE 39 bus + VSC-HVDC

 VSC-HVDC
• Two-level, PWM, decoupled dq
vector control
• Carrier frequency is 1980 Hz
 EMT simulation using PSCAD
• The whole system
• Time step : 5 µs
 Hybrid simulation
• The part encircled by the dashed line
is modeled in detail in PSCAD
• Boundary buses: Buses 26 and 8
• EMT time step : 5 µs
• TS time step: 5 ms

28
 Applied to power systems interfaced with HVDC

 Scenario: a single-line-to-
ground(SLG) fault is applied on
bus 29 (AC bus of the rectifier)
at 1.0 s and cleared after 0.05 s
 Response of VSC-HVDC
rectifier to the fault:
(a) Real power
(b) Reactive power flowing into
the rectifier
(c) DC voltage
(d) DC current

29
 Applied to power systems interfaced with HVDC

• Positive sequence
voltage of bus 25
(within the external
system)

Computation times for a 2-second simulation

Computation
Simulation method time

EMT using PSCAD 1152 s

Hybrid simulation 164 s

7 times

30
 Integrated T&D system modeling
• Physically:
• Distribution systems : in general, 3-phase unbalanced
• Transmission system: 3-phase reasonably balanced
• Conditions at boundary between T&D: 1) reasonably balanced under normal operating
conditions; 2) could be significantly unbalanced during and post contingency
• Modeling

Source: Jason Fuller, PNNL, 2015

Bus i Transmission Bus k

(3-sequence)

Distribution Distribution
(3-phase) (3-phase)

31
 Integrated T&D power flow (TDPF)
 Master-slave splitting method
• Iteratively solve power flow for the transmission and the distribution
systems
 Transmission system power flow
• Positive sequence: conventional power flow
• Negative- and zero- sequence: network solution (I=YV)
 Distribution system power flow
• 3-phase power flow: backward/forward sweep algorithm

Boundary data exchange

between T&D in TDPF
algorithm LoadPQk(1)
LoadPQi(1)
Iinjk(2)
Iinji(2)
Iinjk(0)
Iinji(0)

32
 T&D dynamic simulation based on the Multi-Area Thévenin
Equivalent (MATE) approach
• Partitioned dynamic simulation method: used in PSS/E and PSLF
,  Integration step
,  Network Solution step
• Main challenges lie in the network solution step
• The MATE [4] approach is employed in the network solution step

Reconcile the subsystem solution

results at the link subsystem level at
steps 2 and 3
[4] Martí, José R., Luis R. Linares, Jorge A. Hollman, and Fernando A. Moreira. "OVNI: Integrated software/hardware
solution for real-time simulation of large power systems." In Proceedings of the PSCC, vol. 2. 2002. 33
 Three-phase dynamic simulation
• Developed by extending existing three-sequence system
modeling and TS simulation
• Modeling: Inheritance and the adapter design pattern
• Simulation procedure: the same as the positive-sequence TS except for
the three-phase oriented network solution

Development of three-phase machine dynamic model based

Transformation from three-sequence modeling on the corresponding three-sequence model using the
to three-phase modeling adapter pattern 34
Flowchart of the integrated T&D dynamic simulation

Integrated T&D
power flow

Consider unbalanced
fault(s) in transmission
system

The MATE approach

based T&D network
solution

35
 Conclusions
• A normally cleared SLG fault could result in A/C motor stalling
and propagation to the non-faulted phase, depending on the
connection of step-down transformers
• A/C compressor motors could take a much longer time than 2-5
cycles (typical Tstall value) to stall when the equivalent
impedances between the fault point and A/C motors are large.
• The point-on-wave effects deserve more attention as different
POWs could lead to significantly different results in terms of
A/C motor stalling.
• OpenHybridSim, the first open-source tool for EMT and phasor
domain hybrid simulation, has been developed and is available
from: https://github.com/OpenHybridSim
• A modeling framework and power flow and dynamic
simulation algorithms for integrated T&D systems have been
developed
36
Publication
Qiuhua Huang, V. Vittal, “Application of
Electromagnetic Transient-Transient Stability Hybrid
Simulation to FIDVR Study,” IEEE Transactions on
Power Systems, Vol. 31, No. 4, pp. 2634-2646, July
2016.

37
 Questions?

Vijay Vittal
([email protected])

38