Hi all, I want to study the effect of multiple line contingency by recursively solving pf. The procedure goes as follows: start with some initial line outage, repeatedly solve dcpf, disable the lines that exceed limit until no more overloaded lines. I have some questions when implementing this: --- I disable the overloaded lines by setting their status to 0. If they disconnect the network, then island emerges. How to determine slack bus in each island? Can I arbitrary choose one bus to be slack and will this choice affect pf? --- If power is not balanced in an island, meaning that the generation is not equal to demand, I'd like to apply some control to retain balance by loading shedding or AGC. Is it fine if I directly modify PG and PD value in the case file of that island? --- What runpf does with the case when power is not balance? As we know, pf is solved by P=B\theta where B is Laplacian. And the sum of all the entries of P needs to be zero. While if not, runpf will still return some result. How to interpret that result?
Many thanks for your help! Best, Yueyun
