I guess, that's what Dr. Zimmerman has already described. You can follow
that. But the goal you want to achieve will not be fulfilled just by
checking. What you need to do, you need to put some constraints which will
deter the loop and islanding to happen. There might be literature
describing them, you can find those things there.
Another option might be that you run every possible scenarios (combination
of tie switches) and only compare among the cases where power flow
converges.
On Wed, Apr 10, 2019 at 12:48 PM Nitya Kirana <nityakirana...@gmail.com>
wrote:
> Dear Jubeyer Rahman,
>
> Yes you are right, I'm using SPSO to find the tie switches with the goals
> of reducing power losses and using the Newton Raphson method to analyze the
> power flow. But, I actually don't know how to check it for every iteration,
> please mind to help me.
>
> *Best Regards,*
>
> *Putu Diah Nitya Kirana,*
>
> Assistant and Researcher at Electric Energy Conversion Laboratory B101
> Electrical Engineering Department
> Faculty of Electrical Technology
> Institut Teknologi Sepuluh Nopember
> Surabaya, Indonesia
>
> +62 8170 3178 248
>
>
>
> On Thu, 11 Apr 2019 at 00:36, Jubeyer Rahman <jube...@hawaii.edu> wrote:
>
>> If I understand it correctly, you're running some particle swarm
>> optimization to find the best combination of tie switch states (or network
>> configuration) which will minimize your line loss. So, when particle swarm
>> optimization tries to change the value of the particles (may be the state
>> of tie switches), not every random combination will give you a feasible
>> solution. Let me know if I have got this right or not.
>>
>> On Wed, Apr 10, 2019 at 12:27 PM Jubeyer Rahman <jube...@hawaii.edu>
>> wrote:
>>
>>> No, What I have said is, you need to check at every step of your
>>> iteration. Since you're opening some tie switches, that may result in
>>> islands sometimes.
>>>
>>> On Wed, Apr 10, 2019 at 12:19 PM Nitya Kirana <nityakirana...@gmail.com>
>>> wrote:
>>>
>>>> Dear Jubeyer Rahman,
>>>>
>>>> I just check the connectivity and this is the result :
>>>>
>>>>
>>>>
>>>> Checking connectivity ... single fully connected network
>>>> Elapsed time is 9.566667 seconds.
>>>>
>>>> ================================================================================
>>>> Full
>>>> System
>>>> Number of: ----------
>>>> buses 64
>>>> loads 58
>>>> on 58
>>>> off -
>>>> fixed 58
>>>> dispatchable -
>>>> on -
>>>> off -
>>>> generators 3
>>>> on 3
>>>> off -
>>>> shunt elements -
>>>> branches 66
>>>> on 66
>>>> off -
>>>> ties (off) -
>>>>
>>>> Load
>>>> active (MW)
>>>> dispatched 25.4
>>>> fixed 25.4
>>>> dispatchable -
>>>> nominal 25.4
>>>> on 25.4
>>>> off -
>>>> fixed 25.4
>>>> dispatchable -
>>>> on -
>>>> off -
>>>> reactive (MVAr)
>>>> dispatched 29.1
>>>> fixed 29.1
>>>> dispatchable -
>>>> nominal 29.1
>>>> on 29.1
>>>> off -
>>>> fixed 29.1
>>>> dispatchable -
>>>> on -
>>>> off -
>>>>
>>>> Generation
>>>> active (MW)
>>>> dispatched -
>>>> max capacity -
>>>> on -
>>>> off -
>>>> min capacity -
>>>> on -
>>>> off -
>>>> reactive (MVAr)
>>>> dispatched -
>>>> max capacity -
>>>> on -
>>>> off -
>>>> min capacity -
>>>> on -
>>>> off -
>>>>
>>>> Shunt Injections
>>>> active (MW) -
>>>> reactive (MVAr) -
>>>>
>>>> Branch Losses
>>>> active (MW) -
>>>> reactive (MVAr) -
>>>>
>>>> DC line
>>>> export (MW)
>>>> dispatch -
>>>> max capacity -
>>>> on -
>>>> off -
>>>> min capacity -
>>>> on -
>>>> off -
>>>>
>>>> Reference Buses
>>>> num of ref buses 3
>>>> ref bus numbers 1
>>>> 2
>>>> 3
>>>>
>>>> Is there any wrong with this? because I think everything is fully
>>>> connected.
>>>> *Best Regards,*
>>>>
>>>> *Putu Diah Nitya Kirana,*
>>>>
>>>> Assistant and Researcher at Electric Energy Conversion Laboratory B101
>>>> Electrical Engineering Department
>>>> Faculty of Electrical Technology
>>>> Institut Teknologi Sepuluh Nopember
>>>> Surabaya, Indonesia
>>>>
>>>> +62 8170 3178 248
>>>>
>>>>
>>>>
>>>> On Thu, 11 Apr 2019 at 00:09, Jubeyer Rahman <jube...@hawaii.edu>
>>>> wrote:
>>>>
>>>>> Hi,
>>>>> As per the suggestion of Dr. Zimmerman, check the connectivity of the
>>>>> system by using:
>>>>>
>>>>> case_info('your case')
>>>>>
>>>>> As Ray said, you have isolation in your system, this kind of system
>>>>> usually don't converge if you don't assign different swing generators for
>>>>> different islands.
>>>>>
>>>>> -Jubeyer
>>>>>
>>>>>
>>>>> On Wed, Apr 10, 2019 at 11:37 AM Nitya Kirana <
>>>>> nityakirana...@gmail.com> wrote:
>>>>>
>>>>>> Dear all,
>>>>>>
>>>>>> Hi, sorry for contacting again, I just make the P and Q loads of the
>>>>>> bus 10 times bigger because I put the data wrong at first, but it can not
>>>>>> run and I don't know why. Can anyone help me? thank you. Here I attach my
>>>>>> program again.
>>>>>> *Best Regards,*
>>>>>>
>>>>>> *Putu Diah Nitya Kirana,*
>>>>>>
>>>>>> Assistant and Researcher at Electric Energy Conversion Laboratory B101
>>>>>> Electrical Engineering Department
>>>>>> Faculty of Electrical Technology
>>>>>> Institut Teknologi Sepuluh Nopember
>>>>>> Surabaya, Indonesia
>>>>>>
>>>>>> +62 8170 3178 248
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Tue, 9 Apr 2019 at 17:49, Ilias Sarantakos <
>>>>>> il.saranta...@gmail.com> wrote:
>>>>>>
>>>>>>> Hi Nitya,
>>>>>>>
>>>>>>> I tried your network and it converged; I am using MATPOWER 6.1.
>>>>>>>
>>>>>>> Kind regards,
>>>>>>>
>>>>>>> Ilias Sarantakos
>>>>>>> PhD Student
>>>>>>> Newcastle University
>>>>>>> UK
>>>>>>>
>>>>>>> Στις Τρί, 9 Απρ 2019 στις 9:19 π.μ., ο/η Nitya Kirana <
>>>>>>> nityakirana...@gmail.com> έγραψε:
>>>>>>>
>>>>>>>> Dear Fiaz Ahmad,
>>>>>>>>
>>>>>>>> I actually have run the simulation for IEEE 33 bus, 16 bus, and 69
>>>>>>>> bus using Newton Raphson method and the result is convergen. And also
>>>>>>>> in
>>>>>>>> IEEE 33 bus it has the ratio R/X greater than mine, so I just think
>>>>>>>> that
>>>>>>>> the problem is not the Newton Raphson method. Anw, thank you for your
>>>>>>>> recommendation.
>>>>>>>>
>>>>>>>>
>>>>>>>> *Best Regards, *
>>>>>>>>
>>>>>>>> *Putu Diah Nitya Kirana,*
>>>>>>>>
>>>>>>>> Assistant and Researcher at Electric Energy Conversion Laboratory
>>>>>>>> B101
>>>>>>>> Electrical Engineering Department
>>>>>>>> Faculty of Electrical Technology
>>>>>>>> Institut Teknologi Sepuluh Nopember
>>>>>>>> Surabaya, Indonesia
>>>>>>>>
>>>>>>>> +62 8170 3178 248
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> On Tue, 9 Apr 2019 at 11:54, Fiaz Ahmad <fiazah...@sabanciuniv.edu>
>>>>>>>> wrote:
>>>>>>>>
>>>>>>>>> Dear Nitya use forward-backward sweep algorithm for the power flow
>>>>>>>>> of this system.
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Have a nice day
>>>>>>>>>
>>>>>>>>> Regards,
>>>>>>>>> DR. FIAZ AHMAD
>>>>>>>>> PhD Mechatronics Engineering (Sabanci University, Istanbul
>>>>>>>>> Turkey),
>>>>>>>>> Mob: 0092-332-9124525,
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Tue, Apr 9, 2019 at 9:51 AM Nitya Kirana <
>>>>>>>>> nityakirana...@gmail.com> wrote:
>>>>>>>>>
>>>>>>>>>> Dear Fiaz Ahmad,
>>>>>>>>>>
>>>>>>>>>> Okay thank you, so what's your recommendation for this power
>>>>>>>>>> system analysis method?
>>>>>>>>>> *Best Regards,*
>>>>>>>>>>
>>>>>>>>>> *Putu Diah Nitya Kirana,*
>>>>>>>>>>
>>>>>>>>>> Assistant and Researcher at Electric Energy Conversion Laboratory
>>>>>>>>>> B101
>>>>>>>>>> Electrical Engineering Department
>>>>>>>>>> Faculty of Electrical Technology
>>>>>>>>>> Institut Teknologi Sepuluh Nopember
>>>>>>>>>> Surabaya, Indonesia
>>>>>>>>>>
>>>>>>>>>> +62 8170 3178 248
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> On Tue, 9 Apr 2019 at 01:04, Fiaz Ahmad <
>>>>>>>>>> fiazah...@sabanciuniv.edu> wrote:
>>>>>>>>>>
>>>>>>>>>>> Normally Newton Raphson method doesn't converge for distribution
>>>>>>>>>>> systems due to high r/x ratios
>>>>>>>>>>>
>>>>>>>>>>> On Mon, Apr 8, 2019, 22:48 Carlos A. Castro <ccas...@unicamp.br>
>>>>>>>>>>> wrote:
>>>>>>>>>>>
>>>>>>>>>>>> Dear Nitya
>>>>>>>>>>>>
>>>>>>>>>>>> I would do two things:
>>>>>>>>>>>>
>>>>>>>>>>>> 1. I would replace all commas "," by points "." in the data
>>>>>>>>>>>> file.
>>>>>>>>>>>>
>>>>>>>>>>>> 2. I would try to draw the one line diagram of the network
>>>>>>>>>>>> based on the data, to check whether the network is all connected,
>>>>>>>>>>>> that is,
>>>>>>>>>>>> that there are no isolated buses or group of buses,
>>>>>>>>>>>>
>>>>>>>>>>>> Good luck,
>>>>>>>>>>>>
>>>>>>>>>>>> Carlos
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> Em seg, 8 de abr de 2019 às 14:42, Nitya Kirana <
>>>>>>>>>>>> nityakirana...@gmail.com> escreveu:
>>>>>>>>>>>>
>>>>>>>>>>>>> Dear all,
>>>>>>>>>>>>>
>>>>>>>>>>>>> I have run the power flow of 64 radial test bus system using
>>>>>>>>>>>>> Newton Raphson in Matpower 4.1. It shows the result of power
>>>>>>>>>>>>> losses but it doesn't converge, whenever I run the matlab, the
>>>>>>>>>>>>> result is
>>>>>>>>>>>>> always different, and I don't know what's the problem. Can anyone
>>>>>>>>>>>>> help me?
>>>>>>>>>>>>> because this is my final project. It shows "Warning: Matrix
>>>>>>>>>>>>> is singular to working precision."
>>>>>>>>>>>>>
>>>>>>>>>>>>> So, here is my data :
>>>>>>>>>>>>>
>>>>>>>>>>>>> bus data
>>>>>>>>>>>>> bus type Pd Qd
>>>>>>>>>>>>> Gs Bs area Vm Va baseKV zone Vmax
>>>>>>>>>>>>> Vmin
>>>>>>>>>>>>> 1 3 0 0 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 2 3 0 0 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 3 3 0 0 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 4 1 0,0136 0,016256 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 5 1 0,0204 0,024384 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 6 1 0 0 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 7 1 0,07684 0,091846 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 8 1 0,02856 0,034138 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 9 1 0,03808 0,045517 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 10 1 0,0238 0,028448 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 11 1 0,01768 0,021133 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 12 1 0,017 0,02032 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 13 1 0,022304 0,02666 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 14 1 0,0204 0,002438 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 15 1 0,05848 0,069901 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 16 1 0,26384 0,315366 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 17 1 0,03774 0,04511 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 18 1 0,0204 0,002438 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 19 1 0,05576 0,06665 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 20 1 0,0136 0,016256 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 21 1 0,0306 0,016256 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 22 1 0,0204 0,024384 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 23 1 0,0272 0,032512 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 24 1 0,11764 0,140614 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 25 1 0,0068 0,008128 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 26 1 0,0136 0,016256 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 27 1 0,0238 0,028448 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 28 1 0,112948 0,135006 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 29 1 0,02176 0,02601 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 30 1 0 0 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 31 1 0,01768 0,021133 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 32 1 0,05168 0,061773 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 33 1 0,03468 0,041453 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 34 1 0,02856 0,034138 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 35 1 0,0272 0,032512 0 0 1 1 0 20 1 1 0,9
>>>>>>>>>>>>> 36 1 0,01768 0,021133</fon
>>>>>>>>>>>>>
>>>>>>>>>>>>
