Update to PSY 0.22

Project.toml

@@ -1,6 +1,6 @@
 name = "LITS"
 uuid = "86b0dc02-7903-11e9-325f-f195ca7e6c1a"
-version = "0.3.5"
+version = "0.4.0"
 
 [deps]
 DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e"
@@ -15,7 +15,7 @@ SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 [compat]
 DiffEqBase = "6"
 ForwardDiff = "~v0.10"
-InfrastructureSystems = "~0.6.3"
+InfrastructureSystems = "~0.11"
 NLsolve = "4"
-PowerSystems = "~0.17"
+PowerSystems = "~0.21"
 julia = "^1.2"

docs/src/Examples/example_OMIB.md

@@ -41,7 +41,7 @@ As mentioned earlier, we require to create a `Vector` of `PSY.Bus` to define the
 
 ```julia
 #Define the vector of buses
-nodes_OMIB = [
+buses_OMIB = [
     PSY.Bus(
         1, #number
         "Bus 1", #Name
@@ -55,7 +55,7 @@ nodes_OMIB = [
 ]
 ```
 
-Note that two buses are defined in the vector `nodes_case1`. It is important that the bus numbers are ordered from ``1`` to ``n``, since that structure will be used to construct the vector of variables. Future versions of `LITS` will allow to relax this assumption. Similarly, to define the branches (that also has some parameters that are currently not used, such as the rate and angle limits):
+Note that two buses are defined in the vector `buses_case1`. It is important that the bus numbers are ordered from ``1`` to ``n``, since that structure will be used to construct the vector of variables. Future versions of `LITS` will allow to relax this assumption. Similarly, to define the branches (that also has some parameters that are currently not used, such as the rate and angle limits):
 
 ```julia
 #Define the vector of branches
@@ -64,7 +64,7 @@ branch_OMIB = [PSY.Line(
     true, #available
     0.0, #active power flow initial condition (from-to)
     0.0, #reactive power flow initial condition (from-to)
-    Arc(from = nodes_OMIB[1], to = nodes_OMIB[2]), #Connection between buses
+    Arc(from = buses_OMIB[1], to = buses_OMIB[2]), #Connection between buses
     0.01, #resistance in pu
     0.05, #reactance in pu
     (from = 0.0, to = 0.0), #susceptance in pu
@@ -82,7 +82,7 @@ branch_OMIB_fault = [PSY.Line(
     true, #available
     0.0, #active power flow initial condition (from-to)
     0.0, #reactive power flow initial condition (from-to)
-    Arc(from = nodes_OMIB[1], to = nodes_OMIB[2]), #Connection between buses
+    Arc(from = buses_OMIB[1], to = buses_OMIB[2]), #Connection between buses
     0.02, #resistance in pu
     0.1, #reactance in pu
     (from = 0.0, to = 0.0), #susceptance in pu
@@ -101,7 +101,7 @@ Secondly, we will define devices that can inject/withdraw electric current direc
 loads_OMIB = [PSY.PowerLoad(
     "LBus1", #name
     true, #availability
-    nodes_OMIB[2], #bus
+    buses_OMIB[2], #bus
     PSY.LoadModels.ConstantPower, #type
     0.3, #P
     0.01, #Q
@@ -112,7 +112,7 @@ loads_OMIB = [PSY.PowerLoad(
 inf_gen_OMIB = [PSY.Source(
     "InfBus", #name
     true, #availability
-    nodes_OMIB[1], #bus
+    buses_OMIB[1], #bus
     1.05, #VR
     0.0, #VI
     0.000005, #Xth
@@ -157,7 +157,7 @@ pss_OMIB = PSSFixed(0.0) #No PSS without AVR
 gen_OMIB = PSY.DynamicGenerator(
     1, #Number
     "OMIB_Gen", #name
-    nodes_OMIB[2], #bus
+    buses_OMIB[2], #bus
     1.0, #ω_ref
     1.0, #V_ref
     0.5, #P_ref
@@ -184,7 +184,7 @@ sys = PSY.System(
 )
 
 #Add the buses to the system
-for bus in nodes_OMIB
+for bus in buses_OMIB
     PSY.add_component!(sys, bus)
 end
 
@@ -218,7 +218,7 @@ In here, he will use a Three Phase Fault, that is modeled by modifying the admit
 #Obtain the Ybus of the faulted system
 Ybus_fault = PSY.Ybus(
     branch_OMIB_fault, #fault set of lines
-    nodes_OMIB, #set of buses
+    buses_OMIB, #set of buses
 )[:,:]
 
 #Construct the perturbation

docs/src/Examples/example_lines.md

@@ -25,7 +25,7 @@ To start we will define the data structures for the network.
 ### Buses and Branches
 
 ```julia
-nodes_case9 = [
+buses_case9 = [
     PSY.Bus(1, "Bus 1", "REF", 0, 1.02, (min = 0.94, max = 1.06), 138),
     PSY.Bus(2, "Bus 2", "PV", 0, 1.00, (min = 0.94, max = 1.06), 138),
     PSY.Bus(3, "Bus 3", "PQ", 0, 1.00, (min = 0.94, max = 1.06), 138),
@@ -37,7 +37,7 @@ branch_case9 = [
         true,
         0.0,
         0.0,
-        Arc(from = nodes_case9[1], to = nodes_case9[3]),
+        Arc(from = buses_case9[1], to = buses_case9[3]),
         0.01,
         0.12,
         (from = 0.1, to = 0.1),
@@ -49,7 +49,7 @@ branch_case9 = [
         true,
         0.0,
         0.0,
-        Arc(from = nodes_case9[1], to = nodes_case9[2]),
+        Arc(from = buses_case9[1], to = buses_case9[2]),
         0.01,
         0.12,
         (from = 0.1, to = 0.1),
@@ -61,7 +61,7 @@ branch_case9 = [
         true,
         0.0,
         0.0,
-        Arc(from = nodes_case9[2], to = nodes_case9[3]),
+        Arc(from = buses_case9[2], to = buses_case9[3]),
         0.02,
         0.9,
         (from = 0.5, to = 0.5),
@@ -77,7 +77,7 @@ branch_case9_fault = [
         true,
         0.0,
         0.0,
-        Arc(from = nodes_case9[1], to = nodes_case9[3]),
+        Arc(from = buses_case9[1], to = buses_case9[3]),
         0.03,
         0.36,
         (from = 0.03, to = 0.03),
@@ -89,7 +89,7 @@ branch_case9_fault = [
         true,
         0.0,
         0.0,
-        Arc(from = nodes_case9[1], to = nodes_case9[2]),
+        Arc(from = buses_case9[1], to = buses_case9[2]),
         0.01,
         0.12,
         (from = 0.1, to = 0.1),
@@ -101,7 +101,7 @@ branch_case9_fault = [
         true,
         0.0,
         0.0,
-        Arc(from = nodes_case9[2], to = nodes_case9[3]),
+        Arc(from = buses_case9[2], to = buses_case9[3]),
         0.02,
         0.9,
         (from = 0.5, to = 0.5),
@@ -115,15 +115,15 @@ branch_case9_fault = [
 
 ```julia
 loads_case9 = [
-    PowerLoad("Bus1", true, nodes_case9[1], PSY.LoadModels.ConstantPower, 0.5, 0.1, 1.5, 0.8),
-    PowerLoad("Bus2", true, nodes_case9[2], PSY.LoadModels.ConstantPower, 1.0, 0.3, 1.5, 0.8),
-    PowerLoad("Bus3", true, nodes_case9[3], PSY.LoadModels.ConstantPower, 0.3, 0.1, 0.5, 0.3),
+    PowerLoad("Bus1", true, buses_case9[1], PSY.LoadModels.ConstantPower, 0.5, 0.1, 1.5, 0.8),
+    PowerLoad("Bus2", true, buses_case9[2], PSY.LoadModels.ConstantPower, 1.0, 0.3, 1.5, 0.8),
+    PowerLoad("Bus3", true, buses_case9[3], PSY.LoadModels.ConstantPower, 0.3, 0.1, 0.5, 0.3),
 ]
 
 inf_gen_case9 = [PSY.Source(
     "InfBus", #name
     true, #availability
-    nodes_case9[1],#bus
+    buses_case9[1],#bus
     1.00, #VR
     0.0, #VI
     0.000005,  #Xth
@@ -185,7 +185,7 @@ case9_avr = PSY.AVRTypeI(
 case9_gen = PSY.DynamicGenerator(
     1, #Number
     "Case9Gen",
-    nodes_case9[2], #bus
+    buses_case9[2], #bus
     1.0, # ω_ref,
     1.0124, #V_ref
     0.6, #P_ref
@@ -252,7 +252,7 @@ vsc = PSY.CombinedVIwithVZ(
 case9_inv = PSY.DynamicInverter(
     2, #number
     "DARCO", #name
-    nodes_case9[3], #bus location
+    buses_case9[3], #bus location
     1.0, #ω_ref
     0.8, #V_ref
     0.5, #P_ref
@@ -277,7 +277,7 @@ sys = PSY.System(
 )
 
 #Add the buses to the system
-for bus in nodes_case9
+for bus in buses_case9
     PSY.add_component!(sys, bus)
 end
 
@@ -318,7 +318,7 @@ sys2 = PSY.System(
 )
 
 #Add the buses to the system
-for bus in nodes_case9
+for bus in buses_case9
     PSY.add_component!(sys2, bus)
 end
 

src/LITS.jl

@@ -17,6 +17,8 @@ export small_signal_analysis
 export get_state_series
 export get_voltagemag_series
 export print_init_states
+export print_device_states
+export get_dict_init_states
 
 ####################################### Package Imports ####################################
 import Logging
@@ -37,6 +39,8 @@ include("base/definitions.jl")
 include("base/ports.jl")
 include("base/perturbations.jl")
 include("base/small_signal_results.jl")
+include("base/simulation_initialization.jl")
+include("base/file_system.jl")
 include("base/simulation.jl")
 
 #Common Models
@@ -65,6 +69,24 @@ include("models/inverter_models/converter_models.jl")
 include("models/load_models.jl")
 include("models/source_models.jl")
 
+#Initialization Parameters
+include("initialization/init_device.jl")
+
+#Initialization Generator
+include("initialization/generator_components/init_machine.jl")
+include("initialization/generator_components/init_shaft.jl")
+include("initialization/generator_components/init_avr.jl")
+include("initialization/generator_components/init_tg.jl")
+include("initialization/generator_components/init_pss.jl")
+
+#Initialization Inverter
+include("initialization/inverter_components/init_filter.jl")
+include("initialization/inverter_components/init_DCside.jl")
+include("initialization/inverter_components/init_converter.jl")
+include("initialization/inverter_components/init_frequency_estimator.jl")
+include("initialization/inverter_components/init_inner.jl")
+include("initialization/inverter_components/init_outer.jl")
+
 #System Model
 include("models/system.jl")
 
@@ -73,5 +95,6 @@ include("utils/plot_utils.jl")
 include("utils/immutable_dicts.jl")
 include("utils/print.jl")
 include("utils/kwargs_check.jl")
+include("utils/logging.jl")
 
 end # module

src/base/definitions.jl

@@ -66,6 +66,7 @@ Vi_cnv_var :: Voltage supplied from the converter in the I-component
     VI_inv_var = 11
     Vr_cnv_var = 12
     Vi_cnv_var = 13
+    P_ES_var = 14 #Power from energy source
 end
 
 Base.to_index(ix::inverter_inner_vars) = Int(ix)
@@ -100,5 +101,17 @@ const VOLTAGE_BUSES_NO = "voltage_buses_no"
 const CURRENT_BUSES_NO = "current_buses_no"
 const TOTAL_SHUNTS = "total_shunts"
 const AUX_ARRAYS = "aux_arrays"
+const LOOKUP = "lookup"
 
-const SIMULATION_ACCEPTED_KWARGS = [:initial_guess, :initialize_simulation]
+const SIMULATION_ACCEPTED_KWARGS = [:initial_guess, :initialize_simulation, :system_to_file]
+
+PSY.get_V_ref(value::PSY.AVRFixed) = value.Vf
+PSY.set_V_ref!(value::PSY.AVRFixed, val::Float64) = value.Vf = val
+
+get_V_ref_control(value::PSY.DynamicGenerator) = PSY.get_V_ref(PSY.get_avr(value))
+get_V_ref_control(value::PSY.DynamicInverter) =
+    PSY.get_V_ref(PSY.get_reactive_power(PSY.get_outer_control(value)))
+
+get_P_ref_control(value::PSY.DynamicGenerator) = PSY.get_P_ref(PSY.get_prime_mover(value))
+get_P_ref_control(value::PSY.DynamicInverter) =
+    PSY.get_P_ref(PSY.get_active_power(PSY.get_outer_control(value)))

src/base/file_system.jl

@@ -0,0 +1,9 @@
+function check_folder(folder::String)
+    !isdir(folder) && throw(IS.ConflictingInputsError("Specified folder is not valid"))
+    try
+        mkdir(joinpath(folder, "fake"))
+        rm(joinpath(folder, "fake"))
+    catch e
+        throw(IS.ConflictingInputsError("Specified folder does not have write access [$e]"))
+    end
+end

src/base/perturbations.jl

@@ -1,23 +1,25 @@
 abstract type Perturbation end
 
-struct ThreePhaseFault <: Perturbation
+mutable struct ThreePhaseFault <: Perturbation
     time::Float64
     Ybus::SparseMatrixCSC{Complex{Float64}, Int64}
 end
 
-get_affect(pert::ThreePhaseFault) =
+get_affect(::PSY.System, pert::ThreePhaseFault) =
     (integrator) -> PSY.get_ext(integrator.p)["Ybus"] = pert.Ybus
 
-struct ControlReferenceChange <: Perturbation
+mutable struct ControlReferenceChange <: Perturbation
     time::Float64
     device::PSY.DynamicInjection
     signal_index::Int64
     ref_value::Float64
 end
 
-function get_affect(pert::ControlReferenceChange)
+function get_affect(system::PSY.System, pert::ControlReferenceChange)
+    device = PSY.get_component(typeof(pert.device), system, PSY.get_name(pert.device))
+    pert.device = device
     return (integrator) -> begin
-        control_ref = PSY.get_ext(pert.device)[CONTROL_REFS]
+        control_ref = PSY.get_ext(device)[CONTROL_REFS]
         return control_ref[pert.signal_index] = pert.ref_value
     end
 end