Why would PU bus voltages blow up to infinity in a FLAT dynamic simulation?
Hello,
I just started working extensively with the dynamic simulator in PSS/E 33, and I am currently running a flat (No Contingencies), 1 second, dynamic simulation of a simple 9 bus (3 gens, 3 loads) case. I have loaded all the .raw data and .dyr data into the case via python and the PSSE API (see applicable source code below).
Does anyone have any idea why the PU bus voltages blow up to infinity? I must be missing something somewhere - This exact case simulates perfectly in PowerWorld.
Here is some of the applicable source code:
p.case(static_case_file) #Loads a succesfully solved (with FNSL) case. (Verified against PowerWorld output)
##convert to dynamic case
p.conl(-1,1,1) #setup
p.conl(-1,1,2,[0,0],[100,0,0,100]) #convert loads
p.conl(-1,1,3) #cleanup
p.cong() # convert Gens
p.ordr() # preserve sparsity of network matrices
p.fact() # factorize the network admittance matrix
p.tysl() # switching study network solutions
## Save dynamic case binary file
p.save(dynamic_case_file)
## load dynamic model definitions
p.dyre_new([1,1,1,1], dyr_infile) # loads dynamic parameters in the .dyr file overwrites existing dynamic data
## Save binary file for dynamic information
p.snap([-1,-1,-1,-1,-1], dynamic_snap_file)
##LOAD Saved Files into working memory
p.case(dynamic_case_file)
p.rstr(dynamic_snap_file)
##Setup Dynamic Sim
p.chsb(0,1,[-1,-1,-1,1,13,0])
## strt(Option, Outfile) execute before 'run'
p.strt(0, dynamic_channel_outfile)
## run(OPTION, TPause, numTimeStepsBTWPrint, numTimeStepsOUT, NumtimeStepPLT)
p.run(tpause = .1) # initial run to stabalize (suggested in PSSe)
p.run(tpause = 1) # proceding runs...
Following are the generator models that I am using in the dyr_infile:
1, 'GENSAL', 1, 8.96, .0333, .0333, 23.64, 23.64, .146, .0969, .0608, .05, 0, 0, 0 /
2, 'GENSAL', 1, 6, .0333, .0333, 6.4, 6.4, .8958, .8645, .1198, .05, 0, 0, 0 /
3, 'GENSAL', 1, 5.89, .0333, .0333, 3.01, 3.01, 1.3125, 1.2578, .1813, .05, 0, 0, 0 /
Finally here is the channel_outfile graphed (and zoomed in...)
A few other notes:
I have already verified that complex machine impedance (Zx of ZSOURCE) matches the Xd" and Xq" value in the GENSAL model
For simplicity sake I am modeling 3 transformer branches as simple Non-Transformer Branches.
I am brand new to PSSE, and would really appreciate any thoughts on what might be causing this. I assume there is some parameter causing the differential equations in the dynamic simulation to blow up - or I am missing a command/step in allowing PSSE to completely understand the topology. Any thoughts on how to troubleshoot this beast would also be helpful. In the meantime, I will continue reading through the loads of documentation...
thanks
does adding a `p.fnsl()` after you load the case help just in case it's saved in a funny state?
It looks like that TYSL can not solve the case. Check the message after you run p.strt to see on which bus the power flow solution blows up.
@JervisW -- Thanks for the comment. Adding the "p.fnsl()" did not change anything about the output of the simulation. After some further testing and reading I was able to figure a few things out. Please see my answer below.
@yfwing -- 1/2 -- Thanks for the thought. That makes complete sense, and after looking at the PSS/E output a few times (still getting used to reading it) I did notice that bus 1 and 2 (buses with generators) had some "suspect intial conditions." Unfortunately I was not sure how to use that to
@yfwing --2/2-- troubleshoot the dynamic simulation, as it doesn't really say much about what variable in the dynamic simulation is causing the blow up. Perhaps my answer below will supply some clarification?