Quick Start Guide

To follow along, you can download this tutorial as a Julia script (.jl) or Jupyter notebook (.ipynb).

The data for these tutorials is provided in PowerSystemCaseBuilder.jl. If you want to build your own case, take a look at the tutorial Creating and Handling Data for Dynamic Simulations

For more details about loading data and adding more dynamic components check the tutorial on Adding Data for Dynamic Simulations in PowerSystems.jl.

For a detailed tutorial about this case visit One Machine against Infinite Bus (OMIB) Simulation

Loading data

Data can be loaded from a pss/e raw file and a pss/e dyr file.

using PowerSystems
using PowerSimulationsDynamics
using PowerSystemCaseBuilder
using Sundials
using Plots

omib_sys = build_system(PSIDSystems, "OMIB System")

Define the Simulation

time_span = (0.0, 30.0)
perturbation_trip = BranchTrip(1.0, Line, "BUS 1-BUS 2-i_1")
sim = Simulation!(ResidualModel, omib_sys, pwd(), time_span, perturbation_trip)

Explore initial conditions for the simulation

x0_init = read_initial_conditions(sim)

Dict{String, Any} with 5 entries:
  "generator-102-1" => Dict(:ω=>1.0, :δ=>0.168525)
  "V_R"             => Dict(102=>1.03973, 101=>1.05)
  "Vm"              => Dict(102=>1.04, 101=>1.05)
  "θ"               => Dict(102=>0.0228958, 101=>0.0)
  "V_I"             => Dict(102=>0.0238095, 101=>0.0)

show_states_initial_value(sim)

Voltage Variables
====================
BUS 1
====================
Vm 1.05
θ 0.0
====================
BUS 2
====================
Vm 1.04
θ 0.0229
====================
====================
Differential States
generator-102-1
====================
δ 0.1685
ω 1.0
====================

Obtain small signal results for initial conditions

Show eigenvalues for operating point

small_sig = small_signal_analysis(sim)
summary_eigenvalues(small_sig)

Execute the simulation

execute!(sim, IDA(); dtmax = 0.02, saveat = 0.02, enable_progress_bar = false)

SIMULATION_FINALIZED::BUILD_STATUS = 6

Make a plot of the results

results = read_results(sim)
angle = get_state_series(results, ("generator-102-1", :δ));
plot(angle; xlabel = "time", ylabel = "rotor angle [rad]", label = "gen-102-1")