Analysis of Comparison Challenge Problem

1
Analysis of Comparison Challenge Problem

• by
• Levent U. Gökdere
• University of South Carolina

2
Objective

• Modeling and simulation of a challenge problem
  using various languages and integration methods

Simulation languages
Integration methods
SABER ACSL VTB
Gears Stiff Runge-Kutta 2nd - 4th
Euler Predictor-Corrector Trapezoidal
Programming languages
MS C Borland C TI C
3
Challenge problem
Load1
BCB1
Transmission line
FCB
Fault
S
Turbine
Generator
Load2
BCB2
Time in seconds
0
25
15
FCB Feeder circuit breaker
Simulation starts
S
Simulation ends
BCB1 Branch circuit breaker
BCB2 Branch circuit breaker
S Fictitious switch
Fault occurs
4
Turbine and generator
Transmission line
Generator
Turbine
Single-phase synchronous generator
Current dynamics
Speed dynamics
5
Transmission line and circuit breakers
BCB1
Load1
Transmission line
FCB
Fault
BCB2
Load2
FCB
100 sign( )

-
BCB1
BCB2
6
Loads and fault
Load1
Fault
S
Load2
S
S
7
Simulation with SABER
by Power Electronics Group of USC
MAST templates
Torque source
Transmission line
Generator
FCB
Generator
BCB1
BCB2
Load1
Native models
BCB1
Torque source
Fault
FCB
Transmission line
BCB2
Load2
Load1
Load2
Fault
Schematic view of challenge problem in SABER
8
Simulation with ACSL
by ACSL User Group of USC
Load1
BCB1
Transmission line
FCB
Fault
S
R1
R
Load2
BCB2
FCB
R2
FCB
R R1 R2 100 ?
9
Simulation with VTB
by VTB Network Solver Group of USC
Torque source
Native models
Torque source
Generator
Transmission line
Circuit breakers
Load1
Load2
Fault
Schematic view of challenge problem in VTB
10
Simulation with MS C
by VTB Translator Group of USC
Load1
BCB1
Transmission line
FCB
Fault
Load2
1.0E9 ?
BCB2
Fault
1.0 ?
Predictor-Corrector Runge-Kutta 2nd
Runge-Kutta 4th Euler
Alpha 800 MHz
0.533 s 0.492 s
0.586 s 0.475 s
Time step
Predictor-Corrector Runge-Kutta 2nd
Runge-Kutta 4th Euler
Pentium III 500 MHz
0.766 s 0.776 s
0.9635 s 0.731 s
11
Simulation with Borland C and TI C
by Taganrog State University of Radio Engineering
System 1
System 2
System 3
System 4
Seconds
0
25
15
BCB1
Simulation starts
Simulation ends
Fault occurs
System 1 System 2 System 3 System 4
Time step
Integration method Processor
Execution time
Borland C
Runge-Kutta 3rd Pentium II (266
MHz) 45.81 s
Runge-Kutta 3rd TMS320C44
2 min 40 s
TI C
12
Simulation output data
Observed variables
Current through BCB1
Voltage across BCB1
Velocity of generator
Time interval for observation
15.0
14.9
15.1
0.0
Simulation time (s)
Fault occurs
Time step
Communication interval
13
Results
Maximum time step
Backward-Euler-based Timing Analysis using
Gauss-Jacobi scheme of parts' decomposition with
the modified coupling approach.
Fixed time step
Incomplete
14
Current through upper branch circuit breaker
Current through upper branch circuit breaker
15
Voltage across upper branch circuit breaker
Voltage across upper branch circuit breaker
16
Velocity of generator
Velocity of generator