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 <[email protected]<mailto:[email protected]>> 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
