Re: Power Flow in Matpower

2018-03-22 Thread Ray Zimmerman
   293 293 
>> 1   0   1   -360360;
>> 30   29  0.0034410.0226660.0194  293 293 293 
>> 1   0   1   -360360;
>> 37   39  0.0040440.0266330.0228  293 293 293 
>> 1   0   1   -360360;
>> 39   40  0.0007130.0015620.1176  293 293 293 
>> 1   0   1   -360360;
>> 18   19  0.0009460.0018470.1644  375 375 375 
>> 1   0   1   -360360;
>> 74   71  0.0025810.0170000.0146  293 293 293 
>> 1   0   1   -360360;
>> 75   74  0.0027900.0183760.0157  293 293 293 
>> 1   0   1   -360360;
>> 69   70  0.0034410.0226660.0194  293 293 293 
>> 1   0   1   -360360;
>> 57   55  0.0019200.0122910.0236  293 293 293 
>> 1   0   1   -360360;
>> 20   19  0.0009090.0017750.1580  375 375 375 
>> 1   0   1   -360360;
>> 21   20  0.0010320.0020150.1793  375 375 375 
>> 1   0   1   -360360;
>> 60   61  0.270.0061330.  375037503750
>> 1.160   1   -360360;
>> 98   0.750.0110370.  200020002000
>> 1.020   1   -360360;
>> 59   58  0.0001500.0214790.  100010001000
>> 1.060   1   -360360;
>> 58   57  0.000.0232800.  100010001000
>> 1.010   1   -360360;
>> 623  0.0001500.0232630.  100010001000
>> 0.980   1   -360360;
>> 88   87  0.0001500.0232630.  100010001000
>> 1.010   1   -360360;
>> 43   44  0.0001500.0222610.  100010001000
>> 0.9 0   1   -360360;
>> 232  0.0001500.0213820.  100010001000
>> 0.950   1   -360360;
>> 428  0.0001500.0213820.  100010001000
>> 0.9 0   1   -360360;
>> 12   13  0.0001500.0238550.  100010001000
>> 1   0   1   -360360;
>> 45   46  0.0001500.0222610.  100010001000
>> 0.990   1   -360360;
>> 62   63  0.0001500.0232800.  100010001000
>> 1.020   1   -360360;
>> 79   78  0.0001500.0232800.  100010001000
>> 1.1 0   1   -360360;
>> 871  0.0001500.0232800.  100010001000
>> 1   0   1   -360360;
>> 41   37  0.0002520.0221700.  630 630 630 
>> 1   0   1   -360360;
>> ];
>> 
>> ++++++++
>> 
>> and i got some buses relatively far from the real values
>> 
>> the real values of the voltage are attached 
>> 
>> On 21 March 2018 at 17:33, Carlos E Murillo-Sanchez 
>> <ce.murillosanc...@gmail.com <mailto:ce.murillosanc...@gmail.com>> wrote:
>> The Ybus matrix computed from the data in your file has NaN's and Inf's 
>> because branch # 69 from bus 16 to bus 17 has zero series impedance.  You 
>> must collapse buses 16 and 17 into a single bus before applying any 
>> algorithm to the system because the electrical "distance" between these two 
>> buses is zero.
>> 
>> Carlos.
>> 
>> Mohammed Alhajri wrote:
>>> i got
>>> 
>>> 
>>> 
>>> 19980   NaN
>>> 19981   NaN
>>> 19982   NaN
>>> 19983   NaN
>>> 19984   NaN
>>> 19985   NaN
>>> 19986   NaN
>>> 19987   NaN
>>> 19988   NaN
>>> 19989   NaN
>>> 19990   NaN
>>> 19991   NaN
>>> 19992   NaN
>>> 19993   NaN
>>> 19994   NaN
>>> 19995   NaN
>>> 19996   NaN
>>> 19997   NaN
>>> 19998   NaN
>>

Re: Power Flow in Matpower

2018-03-21 Thread Mohammed Alhajri
3 293 293 1 0 1 -360 360;
> 23 24 0.002962 0.019509 0.0167 293 293 293 1 0 1 -360 360;
> 37 40 0.004535 0.029871 0.0256 293 293 293 1 0 1 -360 360;
> 75 77 0.001217 0.008014 0.0069 293 293 293 1 0 1 -360 360;
> 46 47 0.006760 0.044523 0.0381 293 293 293 1 0 1 -360 360;
> 71 73 0.005973 0.039342 0.0337 293 293 293 1 0 1 -360 360;
> 49 50 0.001475 0.009714 0.0083 293 293 293 1 0 1 -360 360;
> 81 78 0.000549 0.002420 0.0137 146 146 146 1 0 1 -360 360;
> 78 80 0.000725 0.003010 0.0681 293 293 293 1 0 1 -360 360;
> 63 64 0.004302 0.028333 0.0243 293 293 293 1 0 1 -360 360;
> 78 82 0.000935 0.004305 0.0019 146 146 146 1 0 1 -360 360;
> 77 78 0.000701 0.004614 0.0040 293 293 293 1 0 1 -360 360;
> 16 17 0.09 0.17 0.0060 375 375 375 1 0 1 -360 360;
> 23 25 0.001413 0.009309 0.0080 293 293 293 1 0 1 -360 360;
> 25 26 0.004265 0.028090 0.0241 293 293 293 1 0 1 -360 360;
> 47 48 0.006366 0.041933 0.0359 293 293 293 1 0 1 -360 360;
> 13 18 0.001106 0.002159 0.1921 375 375 375 1 0 1 -360 360;
> 13 15 0.000492 0.000960 0.0854 375 375 375 1 0 1 -360 360;
> 28 30 0.003196 0.021047 0.0180 293 293 293 1 0 1 -360 360;
> 30 29 0.003441 0.022666 0.0194 293 293 293 1 0 1 -360 360;
> 37 39 0.004044 0.026633 0.0228 293 293 293 1 0 1 -360 360;
> 39 40 0.000713 0.001562 0.1176 293 293 293 1 0 1 -360 360;
> 18 19 0.000946 0.001847 0.1644 375 375 375 1 0 1 -360 360;
> 74 71 0.002581 0.017000 0.0146 293 293 293 1 0 1 -360 360;
> 75 74 0.002790 0.018376 0.0157 293 293 293 1 0 1 -360 360;
> 69 70 0.003441 0.022666 0.0194 293 293 293 1 0 1 -360 360;
> 57 55 0.001920 0.012291 0.0236 293 293 293 1 0 1 -360 360;
> 20 19 0.000909 0.001775 0.1580 375 375 375 1 0 1 -360 360;
> 21 20 0.001032 0.002015 0.1793 375 375 375 1 0 1 -360 360;
> 60 61 0.27 0.006133 0. 3750 3750 3750 1.16 0 1 -360 360;
> 9 8 0.75 0.011037 0. 2000 2000 2000 1.02 0 1 -360 360;
> 59 58 0.000150 0.021479 0. 1000 1000 1000 1.06 0 1 -360 360;
> 58 57 0.00 0.023280 0. 1000 1000 1000 1.01 0 1 -360 360;
> 6 23 0.000150 0.023263 0. 1000 1000 1000 0.98 0 1 -360 360;
> 88 87 0.000150 0.023263 0. 1000 1000 1000 1.01 0 1 -360 360;
> 43 44 0.000150 0.022261 0. 1000 1000 1000 0.9 0 1 -360 360;
> 2 32 0.000150 0.021382 0. 1000 1000 1000 0.95 0 1 -360 360;
> 4 28 0.000150 0.021382 0. 1000 1000 1000 0.9 0 1 -360 360;
> 12 13 0.000150 0.023855 0. 1000 1000 1000 1 0 1 -360 360;
> 45 46 0.000150 0.022261 0. 1000 1000 1000 0.99 0 1 -360 360;
> 62 63 0.000150 0.023280 0. 1000 1000 1000 1.02 0 1 -360 360;
> 79 78 0.000150 0.023280 0. 1000 1000 1000 1.1 0 1 -360 360;
> 8 71 0.000150 0.023280 0. 1000 1000 1000 1 0 1 -360 360;
> 41 37 0.000252 0.022170 0. 630 630 630 1 0 1 -360 360;
> ];
>
> 
>
> and i got some buses relatively far from the real values
>
> the real values of the voltage are attached
>
> On 21 March 2018 at 17:33, Carlos E Murillo-Sanchez <
> ce.murillosanc...@gmail.com> wrote:
>
>> The Ybus matrix computed from the data in your file has NaN's and Inf's
>> because branch # 69 from bus 16 to bus 17 has zero series impedance.  You
>> must collapse buses 16 and 17 into a single bus before applying any
>> algorithm to the system because the electrical "distance" between these two
>> buses is zero.
>>
>> Carlos.
>>
>> Mohammed Alhajri wrote:
>>
>> i got
>>
>> 
>>
>> 19980   NaN
>> 19981   NaN
>> 19982   NaN
>> 19983   NaN
>> 19984   NaN
>> 19985   NaN
>> 19986   NaN
>> 19987   NaN
>> 19988   NaN
>> 19989   NaN
>> 19990   NaN
>> 19991   NaN
>> 19992   NaN
>> 19993   NaN
>> 19994   NaN
>> 19995   NaN
>> 19996   NaN
>> 19997   NaN
>> 19998   NaN
>> 1   NaN
>> 2   NaN
>> Gauss-Seidel power flow did not converge in 2 iterations.
>>
>> >>>>>  Did NOT converge (47.99 seconds)  <<<<<
>>
>>
>> On 20 March 2018 at 19:44, Ray Zimmerman <r...@cornell.edu> wrote:
>>
>>> Unfortunately, I do not have time to work on this myself. I was just
>>> giving a suggestion for another direction to try if you want to understand
>>> the issue that MATPOWER is having with your case. Could you post the output
>>> (using verbose set to 2) of runpf() when using a MATPOIWER case file
>>> that corresponds to the solved ca

Re: Power Flow in Matpower

2018-03-21 Thread Ray Zimmerman
60;
> 2 32  0.0001500.0213820.  100010001000
> 0.950   1   -360360;
> 4 28  0.0001500.0213820.  100010001000
> 0.9 0   1   -360360;
> 1213  0.0001500.0238550.  100010001000
> 1   0   1   -360360;
> 4546  0.0001500.0222610.  100010001000
> 0.990   1   -360360;
> 6263  0.0001500.0232800.  100010001000
> 1.020   1   -360360;
> 7978  0.0001500.0232800.  100010001000
> 1.1 0   1   -360360;
> 8 71  0.0001500.0232800.  100010001000
> 1   0   1   -360360;
> 4137  0.0002520.0221700.  630 630 630 
> 1   0   1   -360360;
> ];
> 
> 
> 
> and i got some buses relatively far from the real values
> 
> the real values of the voltage are attached 
> 
> On 21 March 2018 at 17:33, Carlos E Murillo-Sanchez 
> <ce.murillosanc...@gmail.com <mailto:ce.murillosanc...@gmail.com>> wrote:
> The Ybus matrix computed from the data in your file has NaN's and Inf's 
> because branch # 69 from bus 16 to bus 17 has zero series impedance.  You 
> must collapse buses 16 and 17 into a single bus before applying any algorithm 
> to the system because the electrical "distance" between these two buses is 
> zero.
> 
> Carlos.
> 
> Mohammed Alhajri wrote:
>> i got
>> 
>> 
>> 
>> 19980   NaN
>> 19981   NaN
>> 19982   NaN
>> 19983   NaN
>> 19984   NaN
>> 19985   NaN
>> 19986   NaN
>> 19987   NaN
>> 19988   NaN
>> 19989   NaN
>> 19990   NaN
>> 19991   NaN
>> 19992   NaN
>> 19993   NaN
>> 19994   NaN
>> 19995   NaN
>> 19996   NaN
>> 19997   NaN
>> 19998   NaN
>> 1   NaN
>> 2   NaN
>> Gauss-Seidel power flow did not converge in 2 iterations.
>> 
>> >>>>>  Did NOT converge (47.99 seconds)  <<<<<
>> 
>> 
>> On 20 March 2018 at 19:44, Ray Zimmerman <r...@cornell.edu 
>> <mailto:r...@cornell.edu>> wrote:
>> Unfortunately, I do not have time to work on this myself. I was just giving 
>> a suggestion for another direction to try if you want to understand the 
>> issue that MATPOWER is having with your case. Could you post the output 
>> (using verbose set to 2) of runpf() when using a MATPOIWER case file that 
>> corresponds to the solved case from the Hadi Sadat code?
>> 
>> And if you have any questions about the MATPOWER case format or MATPOWER 
>> power flow options, feel free to ask.
>> 
>> Ray
>> 
>> 
>> 
>>> On Mar 20, 2018, at 11:28 AM, Mohammed Alhajri <u106...@gmail.com 
>>> <mailto:u106...@gmail.com>> wrote:
>>> 
>>> Ok, i have attached the case information in format of Hadi Saadat code, can 
>>> you please try it in MATPOWER? 
>>> 
>>> because we have spent more than three weeks checking the format, but still 
>>> dose not converge... 
>>> 
>>> Regards,,, 
>>> 
>>> بتاريخ ٢٠١٨/٠٣/٢٠ ٦:٢٦ م، كتب "Ray Zimmerman" <r...@cornell.edu 
>>> <mailto:r...@cornell.edu>>:
>>> It’s possible that the modeling is not identical or that there is some 
>>> error in your conversion to MATPOWER format. You can check by talking the 
>>> solved case from your other software, converting that solved case to 
>>> MATPOWER and then trying the MATPOWER power flow. It should converge in a 
>>> single iteration. If it does not, then you know that there is either a 
>>> mistake somewhere or a difference in modeling.
>>> 
>>>Ray
>>> 
>>> 
>>>> On Mar 16, 2018, at 12:42 PM, Mohammed Alhajri < 
>>>> <mailto:u106...@gmail.com>u106...@gmail.com <mailto:u106...@gmail.com>> 
>>>> wrote:
>>>> 
>>>> i tried that but unfortunately not work 
>>>> 
>>>> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." < 
>>>> <mailto:abhy...@anl.gov>abhy...@anl.gov <mailto:abhy...@anl.gov>>:
>>>> See FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>
>>>>  
>>>> Thanks,
>>>> 
>>>> Shri
>>>> 
>>>> Ph: (630) 252 0219 <tel:%28630%29%20252-0219>
>>>> www.mcs.anl.gov/~abhyshr <http://www.mcs.anl.gov/%7Eabhyshr>
>>>>  
>>>>  
>>>>  
>>>> 
>>>> From: <bounce-122387760-73493...@list.cornell.edu 
>>>> <mailto:bounce-122387760-73493...@list.cornell.edu>> on behalf of Mohammed 
>>>> Alhajri < <mailto:u106...@gmail.com>u106...@gmail.com 
>>>> <mailto:u106...@gmail.com>>
>>>> Reply-To: MATPOWER discussion forum < 
>>>> <mailto:matpowe...@list.cornell.edu>matpowe...@list.cornell.edu 
>>>> <mailto:matpowe...@list.cornell.edu>>
>>>> Date: Friday, March 16, 2018 at 11:26 AM
>>>> To: MATPOWER discussion forum < 
>>>> <mailto:matpowe...@list.cornell.edu>matpowe...@list.cornell.edu 
>>>> <mailto:matpowe...@list.cornell.edu>>
>>>> Subject: Re: Power Flow in Matpower
>>>> 
>>>>  
>>>> any answer to this question?
>>>> 
>>>>  
>>>> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" < 
>>>> <mailto:u106...@gmail.com>u106...@gmail.com <mailto:u106...@gmail.com>>:
>>>> 
>>>> Hello All,
>>>> 
>>>>  
>>>> I did the power flow for a 89-bus network and it converges using Hadi 
>>>> Sadat code after 17080 iterations. The accuracy was 1e-8 and the method is 
>>>> Gauss-Seidel Method.
>>>> 
>>>>  
>>>> But when I did the power flow using matpower it does not converge! I tried 
>>>> to increase the maximum iteration, I put it 10, and still did not 
>>>> cnvarge!
>>>> 
>>>>  
>>>> I have attached the data according to Hadi Sadat Code, can any one try to 
>>>> do the power flow using matpower?
>>>> 
>>>>  
>>> 
>> 
>> 
> 
> 
> 



Re: Power Flow in Matpower

2018-03-21 Thread Mohammed Alhajri
 0.00 0.023280 0. 1000 1000 1000 1.01 0 1 -360 360;
6 23 0.000150 0.023263 0. 1000 1000 1000 0.98 0 1 -360 360;
88 87 0.000150 0.023263 0. 1000 1000 1000 1.01 0 1 -360 360;
43 44 0.000150 0.022261 0. 1000 1000 1000 0.9 0 1 -360 360;
2 32 0.000150 0.021382 0. 1000 1000 1000 0.95 0 1 -360 360;
4 28 0.000150 0.021382 0. 1000 1000 1000 0.9 0 1 -360 360;
12 13 0.000150 0.023855 0. 1000 1000 1000 1 0 1 -360 360;
45 46 0.000150 0.022261 0. 1000 1000 1000 0.99 0 1 -360 360;
62 63 0.000150 0.023280 0. 1000 1000 1000 1.02 0 1 -360 360;
79 78 0.000150 0.023280 0. 1000 1000 1000 1.1 0 1 -360 360;
8 71 0.000150 0.023280 0. 1000 1000 1000 1 0 1 -360 360;
41 37 0.000252 0.022170 0. 630 630 630 1 0 1 -360 360;
];



and i got some buses relatively far from the real values

the real values of the voltage are attached

On 21 March 2018 at 17:33, Carlos E Murillo-Sanchez <
ce.murillosanc...@gmail.com> wrote:

> The Ybus matrix computed from the data in your file has NaN's and Inf's
> because branch # 69 from bus 16 to bus 17 has zero series impedance.  You
> must collapse buses 16 and 17 into a single bus before applying any
> algorithm to the system because the electrical "distance" between these two
> buses is zero.
>
> Carlos.
>
> Mohammed Alhajri wrote:
>
> i got
>
> 
>
> 19980   NaN
> 19981   NaN
> 19982   NaN
> 19983   NaN
> 19984   NaN
> 19985   NaN
> 19986   NaN
> 19987   NaN
> 19988   NaN
> 19989   NaN
> 19990   NaN
> 19991   NaN
> 19992   NaN
> 19993   NaN
> 19994   NaN
> 19995   NaN
> 19996   NaN
> 19997   NaN
> 19998   NaN
> 1   NaN
> 2   NaN
> Gauss-Seidel power flow did not converge in 2 iterations.
>
> >>>>>  Did NOT converge (47.99 seconds)  <<<<<
>
>
> On 20 March 2018 at 19:44, Ray Zimmerman <r...@cornell.edu> wrote:
>
>> Unfortunately, I do not have time to work on this myself. I was just
>> giving a suggestion for another direction to try if you want to understand
>> the issue that MATPOWER is having with your case. Could you post the output
>> (using verbose set to 2) of runpf() when using a MATPOIWER case file
>> that corresponds to the solved case from the Hadi Sadat code?
>>
>> And if you have any questions about the MATPOWER case format or MATPOWER
>> power flow options, feel free to ask.
>>
>> Ray
>>
>>
>>
>> On Mar 20, 2018, at 11:28 AM, Mohammed Alhajri <u106...@gmail.com> wrote:
>>
>> Ok, i have attached the case information in format of Hadi Saadat code,
>> can you please try it in MATPOWER?
>>
>> because we have spent more than three weeks checking the format, but
>> still dose not converge...
>>
>> Regards,,,
>>
>> بتاريخ ٢٠١٨/٠٣/٢٠ ٦:٢٦ م، كتب "Ray Zimmerman" <r...@cornell.edu>:
>>
>>> It’s possible that the modeling is not identical or that there is some
>>> error in your conversion to MATPOWER format. You can check by talking the
>>> solved case from your other software, converting that solved case to
>>> MATPOWER and then trying the MATPOWER power flow. It should converge in a
>>> single iteration. If it does not, then you know that there is either a
>>> mistake somewhere or a difference in modeling.
>>>
>>>Ray
>>>
>>>
>>> On Mar 16, 2018, at 12:42 PM, Mohammed Alhajri < <u106...@gmail.com>
>>> u106...@gmail.com> wrote:
>>>
>>> i tried that but unfortunately not work
>>>
>>> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." <
>>> <abhy...@anl.gov>abhy...@anl.gov>:
>>>
>>>> See FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>
>>>>
>>>>
>>>>
>>>> Thanks,
>>>>
>>>> Shri
>>>>
>>>> Ph: (630) 252 0219 <%28630%29%20252-0219>
>>>>
>>>> www.mcs.anl.gov/~abhyshr <http://www.mcs.anl.gov/%7Eabhyshr>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> *From: *<bounce-122387760-73493...@list.cornell.edu> on behalf of
>>>> Mohammed Alhajri < <u106...@gmail.com>u106...@gmail.com>
>>>> *Reply-To: *MATPOWER discu

Re: Power Flow in Matpower

2018-03-21 Thread Carlos E Murillo-Sanchez

  
  
The Ybus matrix computed from the data
  in your file has NaN's and Inf's because branch # 69 from bus 16
  to bus 17 has zero series impedance.  You must collapse buses 16
  and 17 into a single bus before applying any algorithm to the
  system because the electrical "distance" between these two buses
  is zero.
  
  Carlos.
  
  Mohammed Alhajri wrote:


  i got






  19980               NaN
  19981               NaN
  19982               NaN
  19983               NaN
  19984               NaN
  19985               NaN
  19986               NaN
  19987               NaN
  19988               NaN
  19989               NaN
  19990               NaN
  19991               NaN
  19992               NaN
  19993               NaN
  19994               NaN
  19995               NaN
  19996               NaN
  19997               NaN
  19998               NaN
  1               NaN
  2               NaN
  Gauss-Seidel power flow did not converge in 2
iterations.
  
  
  >>>>>  Did NOT converge (47.99 seconds) 
<<<<<



  
  
On 20 March 2018 at 19:44, Ray
  Zimmerman <r...@cornell.edu>
  wrote:
  
Unfortunately,
  I do not have time to work on this myself. I was just
  giving a suggestion for another direction to try if you
  want to understand the issue that MATPOWER is having with
  your case. Could you post the output (using verbose set to 2) of runpf() when using a MATPOIWER
  case file that corresponds to the solved case from
  the Hadi Sadat code?
  
  
  And if you have any questions about the MATPOWER case
format or MATPOWER power flow options, feel free to ask.


  
    Ray
  

  

  
On Mar 20, 2018, at 11:28 AM, Mohammed
  Alhajri <u106...@gmail.com>
  wrote:


  Ok, i have attached the
case information in format of Hadi
Saadat code, can you please try it in
MATPOWER? 


because we have spent
  more than three weeks checking the
  format, but still dose not
  converge... 


Regards,,, 
  
  
بتاريخ
  ٢٠١٨/٠٣/٢٠ ٦:٢٦ م، كتب "Ray Zimmerman"
  <r...@cornell.edu>:
  
It’s
  possible that the modeling is not
  identical or that there is some
  error in your conversion to
  MATPOWER format. You can check by
  talking the solved case from your
  other software, converting that
  solved case to MATPOWER and then
  trying the MATPOWER power flow. It
  should converge in a single
  iteration. If it does not, then
  you know that there is either a
  mistake somewhere or a difference
  in modeling.
  
  
     Ray
  

  
On Mar 16, 2018, at
  12:42 PM, Mohammed Alhajri
  <

Re: Power Flow in Matpower

2018-03-21 Thread Jose Luis Marín
Hello Mohammed,

The problem is in the *branch joining buses 16 and 17* (line 191 in the
file).  This branch had zero R and zero X.  Just edit the value of X to
something reasonable (I used 0.001), and the case solves just fine, in 4
iterations.

-- 
Jose L. Marin
Grupo AIA



2018-03-21 3:51 GMT+01:00 Mohammed Alhajri <u106...@gmail.com>:

> i got
>
> 
>
> 19980   NaN
> 19981   NaN
> 19982   NaN
> 19983   NaN
> 19984   NaN
> 19985   NaN
> 19986   NaN
> 19987   NaN
> 19988   NaN
> 19989   NaN
> 19990   NaN
> 19991   NaN
> 19992   NaN
> 19993   NaN
> 19994   NaN
> 19995   NaN
> 19996   NaN
> 19997   NaN
> 19998   NaN
> 1   NaN
> 2   NaN
> Gauss-Seidel power flow did not converge in 2 iterations.
>
> >>>>>  Did NOT converge (47.99 seconds)  <<<<<
>
>
> On 20 March 2018 at 19:44, Ray Zimmerman <r...@cornell.edu> wrote:
>
>> Unfortunately, I do not have time to work on this myself. I was just
>> giving a suggestion for another direction to try if you want to understand
>> the issue that MATPOWER is having with your case. Could you post the output
>> (using verbose set to 2) of runpf() when using a MATPOIWER case file
>> that corresponds to the solved case from the Hadi Sadat code?
>>
>> And if you have any questions about the MATPOWER case format or MATPOWER
>> power flow options, feel free to ask.
>>
>> Ray
>>
>>
>>
>> On Mar 20, 2018, at 11:28 AM, Mohammed Alhajri <u106...@gmail.com> wrote:
>>
>> Ok, i have attached the case information in format of Hadi Saadat code,
>> can you please try it in MATPOWER?
>>
>> because we have spent more than three weeks checking the format, but
>> still dose not converge...
>>
>> Regards,,,
>>
>> بتاريخ ٢٠١٨/٠٣/٢٠ ٦:٢٦ م، كتب "Ray Zimmerman" <r...@cornell.edu>:
>>
>>> It’s possible that the modeling is not identical or that there is some
>>> error in your conversion to MATPOWER format. You can check by talking the
>>> solved case from your other software, converting that solved case to
>>> MATPOWER and then trying the MATPOWER power flow. It should converge in a
>>> single iteration. If it does not, then you know that there is either a
>>> mistake somewhere or a difference in modeling.
>>>
>>>Ray
>>>
>>>
>>> On Mar 16, 2018, at 12:42 PM, Mohammed Alhajri <u106...@gmail.com>
>>> wrote:
>>>
>>> i tried that but unfortunately not work
>>>
>>> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." <abhy...@anl.gov
>>> >:
>>>
>>>> See FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>
>>>>
>>>>
>>>>
>>>> Thanks,
>>>>
>>>> Shri
>>>>
>>>> Ph: (630) 252 0219 <(630)%20252-0219>
>>>>
>>>> www.mcs.anl.gov/~abhyshr
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> *From: *<bounce-122387760-73493...@list.cornell.edu> on behalf of
>>>> Mohammed Alhajri <u106...@gmail.com>
>>>> *Reply-To: *MATPOWER discussion forum <matpowe...@list.cornell.edu>
>>>> *Date: Friday, March 16, 2018 at 11:26 AM*
>>>> *To: MATPOWER discussion forum <matpowe...@list.cornell.edu
>>>> <matpowe...@list.cornell.edu>>*
>>>> *Subject: Re: Power Flow in Matpower*
>>>>
>>>>
>>>>
>>>> any answer to this question?
>>>>
>>>>
>>>>
>>>> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" <u106...@gmail.com>:
>>>>
>>>> Hello All,
>>>>
>>>>
>>>>
>>>> I did the power flow for a 89-bus network and it converges using Hadi
>>>> Sadat code after 17080 iterations. The accuracy was 1e-8 and the method is
>>>> Gauss-Seidel Method.
>>>>
>>>>
>>>>
>>>> But when I did the power flow using matpower it does not converge! I
>>>> tried to increase the maximum iteration, I put it 10, and still did not
>>>> cnvarge!
>>>>
>>>>
>>>>
>>>> I have attached the data according to Hadi Sadat Code, can any one try
>>>> to do the power flow using matpower?
>>>>
>>>>
>>>>
>>>>
>>>
>>
>


Re: Power Flow in Matpower

2018-03-20 Thread Mohammed Alhajri
i got



19980   NaN
19981   NaN
19982   NaN
19983   NaN
19984   NaN
19985   NaN
19986   NaN
19987   NaN
19988   NaN
19989   NaN
19990   NaN
19991   NaN
19992   NaN
19993   NaN
19994   NaN
19995   NaN
19996   NaN
19997   NaN
19998   NaN
1   NaN
2   NaN
Gauss-Seidel power flow did not converge in 2 iterations.

>>>>>  Did NOT converge (47.99 seconds)  <<<<<


On 20 March 2018 at 19:44, Ray Zimmerman <r...@cornell.edu> wrote:

> Unfortunately, I do not have time to work on this myself. I was just
> giving a suggestion for another direction to try if you want to understand
> the issue that MATPOWER is having with your case. Could you post the output
> (using verbose set to 2) of runpf() when using a MATPOIWER case file that
> corresponds to the solved case from the Hadi Sadat code?
>
> And if you have any questions about the MATPOWER case format or MATPOWER
> power flow options, feel free to ask.
>
> Ray
>
>
>
> On Mar 20, 2018, at 11:28 AM, Mohammed Alhajri <u106...@gmail.com> wrote:
>
> Ok, i have attached the case information in format of Hadi Saadat code,
> can you please try it in MATPOWER?
>
> because we have spent more than three weeks checking the format, but still
> dose not converge...
>
> Regards,,,
>
> بتاريخ ٢٠١٨/٠٣/٢٠ ٦:٢٦ م، كتب "Ray Zimmerman" <r...@cornell.edu>:
>
>> It’s possible that the modeling is not identical or that there is some
>> error in your conversion to MATPOWER format. You can check by talking the
>> solved case from your other software, converting that solved case to
>> MATPOWER and then trying the MATPOWER power flow. It should converge in a
>> single iteration. If it does not, then you know that there is either a
>> mistake somewhere or a difference in modeling.
>>
>>Ray
>>
>>
>> On Mar 16, 2018, at 12:42 PM, Mohammed Alhajri <u106...@gmail.com> wrote:
>>
>> i tried that but unfortunately not work
>>
>> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." <abhy...@anl.gov>:
>>
>>> See FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>
>>>
>>>
>>>
>>> Thanks,
>>>
>>> Shri
>>>
>>> Ph: (630) 252 0219 <(630)%20252-0219>
>>>
>>> www.mcs.anl.gov/~abhyshr
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> *From: *<bounce-122387760-73493...@list.cornell.edu> on behalf of
>>> Mohammed Alhajri <u106...@gmail.com>
>>> *Reply-To: *MATPOWER discussion forum <matpowe...@list.cornell.edu>
>>> *Date: Friday, March 16, 2018 at 11:26 AM*
>>> *To: MATPOWER discussion forum <matpowe...@list.cornell.edu
>>> <matpowe...@list.cornell.edu>>*
>>> *Subject: Re: Power Flow in Matpower*
>>>
>>>
>>>
>>> any answer to this question?
>>>
>>>
>>>
>>> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" <u106...@gmail.com>:
>>>
>>> Hello All,
>>>
>>>
>>>
>>> I did the power flow for a 89-bus network and it converges using Hadi
>>> Sadat code after 17080 iterations. The accuracy was 1e-8 and the method is
>>> Gauss-Seidel Method.
>>>
>>>
>>>
>>> But when I did the power flow using matpower it does not converge! I
>>> tried to increase the maximum iteration, I put it 10, and still did not
>>> cnvarge!
>>>
>>>
>>>
>>> I have attached the data according to Hadi Sadat Code, can any one try
>>> to do the power flow using matpower?
>>>
>>>
>>>
>>>
>>
>


Re: Power Flow in Matpower

2018-03-20 Thread Mohammed Alhajri
   196.0   528.4
> 92.0
>   off   -   -   -
> -   -
>   fixed  6416.2  5599.8   196.0   528.4
> 92.0
>   dispatchable  -   -   -
> -   -
> on  -   -   -
> -   -
> off -   -   -
> -   -
>   reactive (MVAr)
> dispatched   2125.7  1855.664.4   173.7
> 31.9
>   fixed  2125.7  1855.664.4   173.7
> 31.9
>   dispatchable  -   -   -
> -   -
> curtailed   -   -   -
> -   -
> nominal  2125.7  1855.664.4   173.7
> 31.9
>   on 2125.7  1855.664.4   173.7
> 31.9
>   off   -   -   -
> -   -
>   fixed  2125.7  1855.664.4   173.7
> 31.9
>   dispatchable  -   -   -
> -   -
> on  -   -   -
> -   -
> off -   -   -
> -   -
>
> Generation
>   active (MW)
> dispatched   6621.8  6234.6   387.3
> -   -
> max capacity-   -   -
> -   -
>   on-   -   -
> -   -
>   off   -   -   -
> -   -
> min capacity-   -   -
> -   -
>   on-   -   -
> -   -
>   off   -   -   -
> -   -
>   reactive (MVAr)
> dispatched   2762.6  2601.0   161.7
> -   -
> max capacity-   -   -
> -   -
>   on-   -   -
> -   -
>   off   -   -   -
> -   -
> min capacity-   -   -
> -   -
>   on-   -   -
> -   -
>   off   -   -   -
> -   -
>
> Shunt Injections
> active (MW) -   -   -
> -   -
> reactive (MVAr)  1934.0  1894.0 -   -
> 40.0
>
> 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  1   1   0
> 0   0
>   ref bus numbers  27  27
>
>
> Is this what the data is supposed to represent?
>
> Carlos.
>
> Mohammed Alhajri wrote:
>
> Ok, i have attached the case information in format of Hadi Saadat code,
> can you please try it in MATPOWER?
>
> because we have spent more than three weeks checking the format, but still
> dose not converge...
>
> Regards,,,
>
> بتاريخ ٢٠١٨/٠٣/٢٠ ٦:٢٦ م، كتب "Ray Zimmerman" <r...@cornell.edu>:
>
>> It’s possible that the modeling is not identical or that there is some
>> error in your conversion to MATPOWER format. You can check by talking the
>> solved case from your other software, converting that solved case to
>> MATPOWER and then trying the MATPOWER power flow. It should converge in a
>> single iteration. If it does not, then you know that there is either a
>> mistake somewhere or a difference in modeling.
>>
>>Ray
>>
>>
>> On Mar 16, 2018, at 12:42 PM, Mohammed Alhajri <u106...@gmail.com> wrote:
>>
>> i tried that but unfortunately not work
>>
>> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." <
>> <abhy...@anl.gov>abhy...@anl.gov>:
>>
>>> See FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>
>>>
>>>
>>>
>>> Thanks,
>>>
>>> Shri
>>>
>>> Ph: (630) 2

Re: Power Flow in Matpower

2018-03-20 Thread Mohammed Alhajri
this is the *.m file of the case

On 21 March 2018 at 02:01, André Silva <caussi...@gmail.com> wrote:

> Hi Mohammed,
>
>
> It is easier for anyone here to make an objective analysis of your case,
> if you provide the network refered in Hadi Saadat alterady converted into
> the MatPower version 2 network case format. You just have to adapt your
> excel network information sheets into an *.m file
>
>
> Typically any Power Flow execution for any given network case takes less
> than a dozen or so iterations to converge. This is verifiable not only for
> small networks, but also for networks with thousands of nodes.
>
>
> Best regards,
>
> André
>
> 2018-02-25 15:18 GMT+00:00 Mohammed Alhajri <u106...@gmail.com>:
>
>> Hello All,
>>
>> I did the power flow for a 89-bus network and it converges using Hadi
>> Sadat code after 17080 iterations. The accuracy was 1e-8 and the method is
>> Gauss-Seidel Method.
>>
>>
>>
>> But when I did the power flow using matpower it does not converge! I
>> tried to increase the maximum iteration, I put it 10, and still did not
>> cnvarge!
>>
>>
>>
>> I have attached the data according to Hadi Sadat Code, can any one try to
>> do the power flow using matpower?
>>
>>
>


OETC_PF.m
Description: Binary data


Re: Power Flow in Matpower

2018-03-20 Thread Carlos E Murillo-Sanchez
   off   -   -   -  
  -   -   
    reactive (MVAr)   
      dispatched   2762.6  2601.0   161.7
  -   -   
      max capacity    -   -   -  
  -   -   
    on    -   -   -  
  -   -   
    off   -   -   -  
  -   -   
      min capacity    -   -   -  
  -   -   
    on    -   -   -  
  -   -   
    off   -   -   -  
  -   -   
  
  Shunt Injections    
      active (MW) -   -   -  
  -   -   
      reactive (MVAr)  1934.0  1894.0 -  
  -  40.0 
  
  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  1   1   0  
  0   0   
    ref bus numbers  27 
  27   
  
  Is this what the data is supposed to represent?
  
  Carlos.
  
  Mohammed Alhajri wrote:


  Ok, i have attached the case information in format
of Hadi Saadat code, can you please try it in MATPOWER? 


because we have spent more than three weeks
  checking the format, but still dose not converge... 


Regards,,, 
  
  
بتاريخ ٢٠١٨/٠٣/٢٠ ٦:٢٦ م، كتب "Ray
  Zimmerman" <r...@cornell.edu>:
  
It’s
  possible that the modeling is not identical or that there
  is some error in your conversion to MATPOWER format. You
  can check by talking the solved case from your other
  software, converting that solved case to MATPOWER and then
      trying the MATPOWER power flow. It should converge in a
  single iteration. If it does not, then you know that there
  is either a mistake somewhere or a difference in modeling.
  
  
     Ray
  

  
On Mar 16, 2018, at 12:42 PM, Mohammed Alhajri
  <u106...@gmail.com>
  wrote:


  i tried that but unfortunately not
work 
  



  بتاريخ ٢٠١٨/٠٣/١٦
٨:٣٠ م، كتب "Abhyankar, Shrirang G." <abhy...@anl.gov>:
  

  
See FAQ #5
 

  Thanks,
  Shri
  Ph: (630) 252 0219
  www.mcs.anl.gov/~abhyshr
   

 
 

  From:
  <bounce-122387760-73493617@list.cornell.edu>
  on behalf of Mohammed Alhajri <u106...@gmail.com>
  
Reply-To:
  MATPOWER discussion forum <matpowe...@list.cornell.edu>

  Date:
Friday, Marc

Re: Power Flow in Matpower

2018-03-20 Thread André Silva
Hi Mohammed,


It is easier for anyone here to make an objective analysis of your case, if
you provide the network refered in Hadi Saadat alterady converted into the
MatPower version 2 network case format. You just have to adapt your excel
network information sheets into an *.m file


Typically any Power Flow execution for any given network case takes less
than a dozen or so iterations to converge. This is verifiable not only for
small networks, but also for networks with thousands of nodes.


Best regards,

André

2018-02-25 15:18 GMT+00:00 Mohammed Alhajri <u106...@gmail.com>:

> Hello All,
>
> I did the power flow for a 89-bus network and it converges using Hadi
> Sadat code after 17080 iterations. The accuracy was 1e-8 and the method is
> Gauss-Seidel Method.
>
>
>
> But when I did the power flow using matpower it does not converge! I tried
> to increase the maximum iteration, I put it 10, and still did not
> cnvarge!
>
>
>
> I have attached the data according to Hadi Sadat Code, can any one try to
> do the power flow using matpower?
>
>


Re: Power Flow in Matpower

2018-03-20 Thread Mohammed Alhajri
Ok, i have attached the case information in format of Hadi Saadat code, can
you please try it in MATPOWER?

because we have spent more than three weeks checking the format, but still
dose not converge...

Regards,,,

بتاريخ ٢٠١٨/٠٣/٢٠ ٦:٢٦ م، كتب "Ray Zimmerman" <r...@cornell.edu>:

> It’s possible that the modeling is not identical or that there is some
> error in your conversion to MATPOWER format. You can check by talking the
> solved case from your other software, converting that solved case to
> MATPOWER and then trying the MATPOWER power flow. It should converge in a
> single iteration. If it does not, then you know that there is either a
> mistake somewhere or a difference in modeling.
>
>Ray
>
>
> On Mar 16, 2018, at 12:42 PM, Mohammed Alhajri <u106...@gmail.com> wrote:
>
> i tried that but unfortunately not work
>
> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." <abhy...@anl.gov>:
>
>> See FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>
>>
>>
>>
>> Thanks,
>>
>> Shri
>>
>> Ph: (630) 252 0219 <(630)%20252-0219>
>>
>> www.mcs.anl.gov/~abhyshr
>>
>>
>>
>>
>>
>>
>>
>> *From: *<bounce-122387760-73493...@list.cornell.edu> on behalf of
>> Mohammed Alhajri <u106...@gmail.com>
>> *Reply-To: *MATPOWER discussion forum <matpowe...@list.cornell.edu>
>> *Date: Friday, March 16, 2018 at 11:26 AM*
>> *To: MATPOWER discussion forum <matpowe...@list.cornell.edu
>> <matpowe...@list.cornell.edu>>*
>> *Subject: Re: Power Flow in Matpower*
>>
>>
>>
>> any answer to this question?
>>
>>
>>
>> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" <u106...@gmail.com>:
>>
>> Hello All,
>>
>>
>>
>> I did the power flow for a 89-bus network and it converges using Hadi
>> Sadat code after 17080 iterations. The accuracy was 1e-8 and the method is
>> Gauss-Seidel Method.
>>
>>
>>
>> But when I did the power flow using matpower it does not converge! I
>> tried to increase the maximum iteration, I put it 10, and still did not
>> cnvarge!
>>
>>
>>
>> I have attached the data according to Hadi Sadat Code, can any one try to
>> do the power flow using matpower?
>>
>>
>>
>>
>


Re: Power Flow in Matpower

2018-03-20 Thread Ray Zimmerman
It’s possible that the modeling is not identical or that there is some error in 
your conversion to MATPOWER format. You can check by talking the solved case 
from your other software, converting that solved case to MATPOWER and then 
trying the MATPOWER power flow. It should converge in a single iteration. If it 
does not, then you know that there is either a mistake somewhere or a 
difference in modeling.

   Ray


> On Mar 16, 2018, at 12:42 PM, Mohammed Alhajri <u106...@gmail.com> wrote:
> 
> i tried that but unfortunately not work 
> 
> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." <abhy...@anl.gov 
> <mailto:abhy...@anl.gov>>:
> See FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>
>  
> 
> Thanks,
> 
> Shri
> 
> Ph: (630) 252 0219 <tel:(630)%20252-0219>
> www.mcs.anl.gov/~abhyshr <http://www.mcs.anl.gov/~abhyshr>
>  
> 
>  
> 
>  
> 
> From: <bounce-122387760-73493...@list.cornell.edu 
> <mailto:bounce-122387760-73493...@list.cornell.edu>> on behalf of Mohammed 
> Alhajri <u106...@gmail.com <mailto:u106...@gmail.com>>
> Reply-To: MATPOWER discussion forum <matpowe...@list.cornell.edu 
> <mailto:matpowe...@list.cornell.edu>>
> Date: Friday, March 16, 2018 at 11:26 AM
> To: MATPOWER discussion forum <matpowe...@list.cornell.edu 
> <mailto:matpowe...@list.cornell.edu>>
> Subject: Re: Power Flow in Matpower
> 
>  
> 
> any answer to this question?
> 
>  
> 
> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" <u106...@gmail.com 
> <mailto:u106...@gmail.com>>:
> 
> Hello All,
> 
>  
> 
> I did the power flow for a 89-bus network and it converges using Hadi Sadat 
> code after 17080 iterations. The accuracy was 1e-8 and the method is 
> Gauss-Seidel Method.
> 
>  
> 
> But when I did the power flow using matpower it does not converge! I tried to 
> increase the maximum iteration, I put it 10, and still did not cnvarge!
> 
>  
> 
> I have attached the data according to Hadi Sadat Code, can any one try to do 
> the power flow using matpower?
> 
>  
> 



Re: Power Flow in Matpower

2018-03-16 Thread Mohammed Alhajri
i tried that but unfortunately not work

بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." <abhy...@anl.gov>:

> See FAQ #5 <http://www.pserc.cornell.edu/matpower/#pfconvergence>
>
>
>
> Thanks,
>
> Shri
>
> Ph: (630) 252 0219 <(630)%20252-0219>
>
> www.mcs.anl.gov/~abhyshr
>
>
>
>
>
>
>
> *From: *<bounce-122387760-73493...@list.cornell.edu> on behalf of
> Mohammed Alhajri <u106...@gmail.com>
> *Reply-To: *MATPOWER discussion forum <matpowe...@list.cornell.edu>
> *Date: *Friday, March 16, 2018 at 11:26 AM
> *To: *MATPOWER discussion forum <matpowe...@list.cornell.edu>
> *Subject: *Re: Power Flow in Matpower
>
>
>
> any answer to this question?
>
>
>
> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" <u106...@gmail.com>:
>
> Hello All,
>
>
>
> I did the power flow for a 89-bus network and it converges using Hadi
> Sadat code after 17080 iterations. The accuracy was 1e-8 and the method is
> Gauss-Seidel Method.
>
>
>
> But when I did the power flow using matpower it does not converge! I tried
> to increase the maximum iteration, I put it 10, and still did not
> cnvarge!
>
>
>
> I have attached the data according to Hadi Sadat Code, can any one try to
> do the power flow using matpower?
>
>
>
>


Re: Power Flow in Matpower

2018-03-16 Thread Abhyankar, Shrirang G.
See FAQ #5<http://www.pserc.cornell.edu/matpower/#pfconvergence>

Thanks,
Shri
Ph: (630) 252 0219
www.mcs.anl.gov/~abhyshr<http://www.mcs.anl.gov/~abhyshr>



From: <bounce-122387760-73493...@list.cornell.edu> on behalf of Mohammed 
Alhajri <u106...@gmail.com>
Reply-To: MATPOWER discussion forum <matpowe...@list.cornell.edu>
Date: Friday, March 16, 2018 at 11:26 AM
To: MATPOWER discussion forum <matpowe...@list.cornell.edu>
Subject: Re: Power Flow in Matpower

any answer to this question?

بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" 
<u106...@gmail.com<mailto:u106...@gmail.com>>:
Hello All,

I did the power flow for a 89-bus network and it converges using Hadi Sadat 
code after 17080 iterations. The accuracy was 1e-8 and the method is 
Gauss-Seidel Method.

But when I did the power flow using matpower it does not converge! I tried to 
increase the maximum iteration, I put it 10, and still did not cnvarge!

I have attached the data according to Hadi Sadat Code, can any one try to do 
the power flow using matpower?




Re: Power Flow in Matpower

2018-03-16 Thread Mohammed Alhajri
any answer to this question?

بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" <u106...@gmail.com>:

> Hello All,
>
> I did the power flow for a 89-bus network and it converges using Hadi
> Sadat code after 17080 iterations. The accuracy was 1e-8 and the method is
> Gauss-Seidel Method.
>
>
>
> But when I did the power flow using matpower it does not converge! I tried
> to increase the maximum iteration, I put it 10, and still did not
> cnvarge!
>
>
>
> I have attached the data according to Hadi Sadat Code, can any one try to
> do the power flow using matpower?
>
>


Re: Help on the reactive power flow in matpower.

2011-04-27 Thread Ray Zimmerman
The power loss (real and reactive) in a branch is given by s_f + s_t, where ...

   s_f = v_f * conj(i_f)
   s_t = v_t * conj(i_t)

... and i_f and i_t are computed via (3.1) and (3.2) in the manual.

-- 
Ray Zimmerman
Senior Research Associate
211 Warren Hall, Cornell University, Ithaca, NY 14853
phone: (607) 255-9645



On Apr 25, 2011, at 10:59 AM, z qin wrote:

 Hello,
 
 When I calculate the power flow from bus i to bus j, I set 
 Sij=Pij+1j*Qij=V_i*(I_ij^*). From the real part, I get Pij, which matches the 
 result of matpower. From the imaginary part, I get Qij, which doesn't match 
 the result of matpower. So what is wrong with my calculation? Thanks.
 
 
 Jerry.
 
 
 On Sun, Apr 24, 2011 at 11:51 PM, jigna pandya jigna_pandya2...@yahoo.com 
 wrote:
 hello,
 Reactive power flow balance equation is as follows:
 Qgeneration+Qcharging=Qload+Qloss .
 (Due to presence of capacitive effect, transmission line injects reactive 
 power in to the system.)
 so  Qfrombus+Qtobus=Qloss doesn't hold.
 
 regards
 
 From: z qin zhengr...@gmail.com
 To: matpower-l matpowe...@cornell.edu
 Sent: Mon, April 25, 2011 8:19:11 AM
 Subject: Help on the reactive power flow in matpower.
 
 Hello,
 
 I ran runopf('case9.m') and got the result below.  In some branches (1, 4 and 
 7), Qfrombus+Qtobus=Qloss. However for the rest branches, 
 Qfrombus+Qtobus=Qloss doesn't hold. Could someone explain for me? I had 
 thought that 
 Qij=|Vi|*|Vi|*Bij+|Vi|*|Vj|*(Gij*sin(theta_i-theta_j)-Bij*cos(theta_i-theta_j),
  but it seems that it is wrong. Thanks a lot.
 
 Btw, here is how I calculate Pij and Qij. Sij=Pij+1j*Qij=V_i*(I_ij^*), from 
 the real part, I get Pij, which seems correct; from the imaginary part, I got 
 Qij, which seems wrong.
 
 Best,
 
 Jerry.
 
  
 Branch Data  |
 
 Brnch   From   ToFrom Bus Injection   To Bus Injection Loss (I^2 * Z) 
  
   # BusBusP (MW)   Q (MVAr)   P (MW)   Q (MVAr)   P (MW)   Q 
 (MVAr)
 -  -  -            
 
1  1  4 89.80 12.94-89.80 -9.02 0.000  3.92
2  4  5 35.22 -3.90-35.04-13.87 0.181  0.98
3  5  6-54.96-16.13 55.97-22.18 1.010  4.40
4  3  6 94.19-22.62-94.19 27.28 0.000  4.66
5  6  7 38.22 -5.10-38.07-18.68 0.149  1.26
6  7  8-61.93-16.32 62.21  0.82 0.279  2.36
7  8  2   -134.32  9.32134.32  0.05 0.000  9.36
8  8  9 72.11-10.14-70.72-18.94 1.394  7.01
9  9  4-54.28-31.06 54.58 12.92 0.295  2.51
 
 
 
 
 
 
 -- 
 
 
 Zhengrui (Jerry) Qin
 
 Computer Science
 College of William and Mary
 Williamsburg, VA 23185
 
 



Re: Help on the reactive power flow in matpower.

2011-04-25 Thread z qin
Jigna,

Thank you very much. Do you know the formula to calculate the reactive power
flow on a branch if we know the impedance of the branch and the voltage
vector of the two buses? Thanks.

Best,

Jerry


On Sun, Apr 24, 2011 at 11:51 PM, jigna pandya
jigna_pandya2...@yahoo.comwrote:

 hello,
 Reactive power flow balance equation is as follows:
 Qgeneration+Qcharging=Qload+Qloss .
 (Due to presence of capacitive effect, transmission line injects reactive
 power in to the system.)
 so  Qfrombus+Qtobus=Qloss doesn't hold.

 regards

 --
 *From:* z qin zhengr...@gmail.com
 *To:* matpower-l matpowe...@cornell.edu
 *Sent:* Mon, April 25, 2011 8:19:11 AM
 *Subject:* Help on the reactive power flow in matpower.

 Hello,

 I ran runopf('case9.m') and got the result below.  In some branches (1, 4
 and 7), Qfrombus+Qtobus=Qloss. However for the rest branches,
 Qfrombus+Qtobus=Qloss doesn't hold. Could someone explain for me? I had
 thought that
 Qij=|Vi|*|Vi|*Bij+|Vi|*|Vj|*(Gij*sin(theta_i-theta_j)-Bij*cos(theta_i-theta_j),
 but it seems that it is wrong. Thanks a lot.

 Btw, here is how I calculate Pij and Qij. Sij=Pij+1j*Qij=V_i*(I_ij^*), from
 the real part, I get Pij, which seems correct; from the imaginary part, I
 got Qij, which seems wrong.

 Best,

 Jerry.


 Branch Data  |

 
 Brnch   From   ToFrom Bus Injection   To Bus Injection Loss (I^2 *
 Z)
   # BusBusP (MW)   Q (MVAr)   P (MW)   Q (MVAr)   P (MW)   Q
 (MVAr)
 -  -  -          
 
1  1  4 89.80 12.94-89.80 -9.02 0.000
 3.92
2  4  5 35.22 -3.90-35.04-13.87 0.181
 0.98
3  5  6-54.96-16.13 55.97-22.18 1.010
 4.40
4  3  6 94.19-22.62-94.19 27.28 0.000
 4.66
5  6  7 38.22 -5.10-38.07-18.68 0.149
 1.26
6  7  8-61.93-16.32 62.21  0.82 0.279
 2.36
7  8  2   -134.32  9.32134.32  0.05 0.000
 9.36
8  8  9 72.11-10.14-70.72-18.94 1.394
 7.01
9  9  4-54.28-31.06 54.58 12.92 0.295
 2.51






-- 


Zhengrui (Jerry) Qin

Computer Science
College of William and Mary
Williamsburg, VA 23185




Re: Help on the reactive power flow in matpower.

2011-04-25 Thread z qin
Hello,

When I calculate the power flow from bus i to bus j, I set
Sij=Pij+1j*Qij=V_i*(I_ij^*). From the real part, I get Pij, which matches
the result of matpower. From the imaginary part, I get Qij, which doesn't
match the result of matpower. So what is wrong with my calculation? Thanks.


Jerry.


On Sun, Apr 24, 2011 at 11:51 PM, jigna pandya
jigna_pandya2...@yahoo.comwrote:

 hello,
 Reactive power flow balance equation is as follows:
 Qgeneration+Qcharging=Qload+Qloss .
 (Due to presence of capacitive effect, transmission line injects reactive
 power in to the system.)
 so  Qfrombus+Qtobus=Qloss doesn't hold.

 regards

 --
 *From:* z qin zhengr...@gmail.com
 *To:* matpower-l matpowe...@cornell.edu
 *Sent:* Mon, April 25, 2011 8:19:11 AM
 *Subject:* Help on the reactive power flow in matpower.

 Hello,

 I ran runopf('case9.m') and got the result below.  In some branches (1, 4
 and 7), Qfrombus+Qtobus=Qloss. However for the rest branches,
 Qfrombus+Qtobus=Qloss doesn't hold. Could someone explain for me? I had
 thought that
 Qij=|Vi|*|Vi|*Bij+|Vi|*|Vj|*(Gij*sin(theta_i-theta_j)-Bij*cos(theta_i-theta_j),
 but it seems that it is wrong. Thanks a lot.

 Btw, here is how I calculate Pij and Qij. Sij=Pij+1j*Qij=V_i*(I_ij^*), from
 the real part, I get Pij, which seems correct; from the imaginary part, I
 got Qij, which seems wrong.

 Best,

 Jerry.


 Branch Data  |

 
 Brnch   From   ToFrom Bus Injection   To Bus Injection Loss (I^2 *
 Z)
   # BusBusP (MW)   Q (MVAr)   P (MW)   Q (MVAr)   P (MW)   Q
 (MVAr)
 -  -  -          
 
1  1  4 89.80 12.94-89.80 -9.02 0.000
 3.92
2  4  5 35.22 -3.90-35.04-13.87 0.181
 0.98
3  5  6-54.96-16.13 55.97-22.18 1.010
 4.40
4  3  6 94.19-22.62-94.19 27.28 0.000
 4.66
5  6  7 38.22 -5.10-38.07-18.68 0.149
 1.26
6  7  8-61.93-16.32 62.21  0.82 0.279
 2.36
7  8  2   -134.32  9.32134.32  0.05 0.000
 9.36
8  8  9 72.11-10.14-70.72-18.94 1.394
 7.01
9  9  4-54.28-31.06 54.58 12.92 0.295
 2.51






-- 


Zhengrui (Jerry) Qin

Computer Science
College of William and Mary
Williamsburg, VA 23185




Help on the reactive power flow in matpower.

2011-04-24 Thread z qin
Hello,

I ran runopf('case9.m') and got the result below.  In some branches (1, 4
and 7), Qfrombus+Qtobus=Qloss. However for the rest branches,
Qfrombus+Qtobus=Qloss doesn't hold. Could someone explain for me? I had
thought that
Qij=|Vi|*|Vi|*Bij+|Vi|*|Vj|*(Gij*sin(theta_i-theta_j)-Bij*cos(theta_i-theta_j),
but it seems that it is wrong. Thanks a lot.

Btw, here is how I calculate Pij and Qij. Sij=Pij+1j*Qij=V_i*(I_ij^*), from
the real part, I get Pij, which seems correct; from the imaginary part, I
got Qij, which seems wrong.

Best,

Jerry.


Branch Data  |

Brnch   From   ToFrom Bus Injection   To Bus Injection Loss (I^2 *
Z)
  # BusBusP (MW)   Q (MVAr)   P (MW)   Q (MVAr)   P (MW)   Q
(MVAr)
-  -  -          

   1  1  4 89.80 12.94-89.80 -9.02 0.000
3.92
   2  4  5 35.22 -3.90-35.04-13.87 0.181
0.98
   3  5  6-54.96-16.13 55.97-22.18 1.010
4.40
   4  3  6 94.19-22.62-94.19 27.28 0.000
4.66
   5  6  7 38.22 -5.10-38.07-18.68 0.149
1.26
   6  7  8-61.93-16.32 62.21  0.82 0.279
2.36
   7  8  2   -134.32  9.32134.32  0.05 0.000
9.36
   8  8  9 72.11-10.14-70.72-18.94 1.394
7.01
   9  9  4-54.28-31.06 54.58 12.92 0.295
2.51


Re: Help on the reactive power flow in matpower.

2011-04-24 Thread jigna pandya
hello,
Reactive power flow balance equation is as follows:
Qgeneration+Qcharging=Qload+Qloss .
(Due to presence of capacitive effect, transmission line injects reactive power 
in to the system.)
so  Qfrombus+Qtobus=Qloss doesn't hold.

regards




From: z qin zhengr...@gmail.com
To: matpower-l matpowe...@cornell.edu
Sent: Mon, April 25, 2011 8:19:11 AM
Subject: Help on the reactive power flow in matpower.

Hello,


I  ran runopf('case9.m') and got the result below.  In some branches (1, 4  and 
7), Qfrombus+Qtobus=Qloss. However for the rest branches,  
Qfrombus+Qtobus=Qloss 
doesn't hold. Could someone explain for me? I had  thought that 
Qij=|Vi|*|Vi|*Bij+|Vi|*|Vj|*(Gij*sin(theta_i-theta_j)-Bij*cos(theta_i-theta_j), 
but it seems that it is wrong. Thanks a lot.

Btw, here is how I calculate Pij and Qij.  Sij=Pij+1j*Qij=V_i*(I_ij^*), from 
the 
real part, I get Pij, which seems  correct; from the imaginary part, I got Qij, 
which seems wrong.

Best,

Jerry.

 
Branch Data  |

Brnch   From   ToFrom Bus Injection   To Bus Injection Loss (I^2 * Z)  
  # BusBusP (MW)   Q (MVAr)   P (MW)   Q (MVAr)   P (MW)   Q (MVAr)
-  -  -            
   1  1  4 89.80 12.94-89.80 -9.02 0.000  3.92
   2  4  5 35.22 -3.90-35.04-13.87 0.181  0.98
   3  5  6-54.96-16.13 55.97-22.18 1.010  4.40
   4  3  6 94.19-22.62-94.19 27.28 0.000  4.66
   5  6  7 38.22 -5.10-38.07-18.68 0.149  1.26
   6  7  8-61.93-16.32 62.21  0.82 0.279  2.36
   7  8  2   -134.32  9.32134.32  0.05 0.000  9.36
   8  8  9 72.11-10.14-70.72-18.94 1.394  7.01
   9  9  4-54.28-31.06 54.58 12.92 0.295  2.51