Thanks a lot. It was very useful.
Alıntı Ray Daniel Zimmerman <r...@cornell.edu>:
Hi Beyza,
Let me first mention that I did not write the radial power flow
code, so I am not as familiar with the details as the original
author, Mirko Todorovski.
1. The check for radial topology in order_radial() (called by
radial_pf()) is done without regard for the status of the branches
in the branch matrix. This means that radial_pf() expects that all
out-of-service branches have been removed before it is called. And
indeed, runpf() does call ext2int() to remove off-line equipment
before calling radial_pf().
2. Looking at the code at the top of radial_pf(), it appears that
after doing the following …
mpc = order_radial(ext2int(mpc));
… mpc.loop should contain the number of loops in the network.
As a side note, I would love to see the code redesigned so
that connected_components() , which traverses the network
efficiently, returns all of the information needed by
order_radial() to avoid duplication of network traversal code.
Making connected_components() return information about loops would
be straightforward if you (or anyone else) would like to give it a
try.
Ray
On Jun 1, 2021, at 7:37 AM, Beyza Nur Dilber <dilb...@itu.edu.tr> wrote:
Hello,
I am Beyza, from Turkey. I hope everyone is well.
I have a few question about radial power flow on Matpower.
1. I would like to run 16bus radial distribution system. There are
some branches which tie buses and cause loop. But these "branches
status" are '0'. So, they can be thought as opened switches
(sectionalizing switches). Thus, system is actually radial with
these conditions. However, when I write radial_pf, it is stated
that "There are loops". I don't understand why it says that. On the
other hand, when I write the code below,
mpopt = mpoption('pf.alg','ISUM')
results=runpf(mpc,mpopt);
radial power flow can be done. (I understand it is radial power
flow instead of newton raphson because when I change the branches
and turn into not radial situation, it shows warning that 'it is
only for radial')
2. My other question is, is there any code for radiality check?
(Instead of error, I want to get some results like 0-1)
Thanks a lot. Have a nice day.
Beyza Nur Dilber