Power Flow in Matpower
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? Branch Data.xlsx Description: MS-Excel 2007 spreadsheet Bus Data.xlsx Description: MS-Excel 2007 spreadsheet
Re: Power Flow in Matpower
any answer to this question? بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" : > 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
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: on behalf of Mohammed Alhajri Reply-To: MATPOWER discussion forum Date: Friday, March 16, 2018 at 11:26 AM To: MATPOWER discussion forum Subject: Re: Power Flow in Matpower any answer to this question? بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" 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
i tried that but unfortunately not work بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." : > 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: * on behalf of > Mohammed Alhajri > *Reply-To: *MATPOWER discussion forum > *Date: *Friday, March 16, 2018 at 11:26 AM > *To: *MATPOWER discussion forum > *Subject: *Re: Power Flow in Matpower > > > > any answer to this question? > > > > بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" : > > 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
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 wrote: > > i tried that but unfortunately not work > > بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." <mailto:abhy...@anl.gov>>: > 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: <mailto:bounce-122387760-73493...@list.cornell.edu>> on behalf of Mohammed > Alhajri mailto:u106...@gmail.com>> > Reply-To: MATPOWER discussion forum <mailto:matpowe...@list.cornell.edu>> > Date: Friday, March 16, 2018 at 11:26 AM > To: MATPOWER discussion forum <mailto:matpowe...@list.cornell.edu>> > Subject: Re: Power Flow in Matpower > > > > any answer to this question? > > > > بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" <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
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" : > 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 wrote: > > i tried that but unfortunately not work > > بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." : > >> 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: * on behalf of >> Mohammed Alhajri >> *Reply-To: *MATPOWER discussion forum >> *Date: Friday, March 16, 2018 at 11:26 AM* >> *To: MATPOWER discussion forum > >* >> *Subject: Re: Power Flow in Matpower* >> >> >> >> any answer to this question? >> >> >> >> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" : >> >> 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
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 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" <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>> wrote: >> >> i tried that but unfortunately not work >> >> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." > <mailto:abhy...@anl.gov>>: >> 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: > <mailto:bounce-122387760-73493...@list.cornell.edu>> on behalf of Mohammed >> Alhajri mailto:u106...@gmail.com>> >> Reply-To: MATPOWER discussion forum > <mailto:matpowe...@list.cornell.edu>> >> Date: Friday, March 16, 2018 at 11:26 AM >> To: MATPOWER discussion forum > <mailto:matpowe...@list.cornell.edu>> >> Subject: Re: Power Flow in Matpower >> >> >> >> any answer to this question? >> >> >> >> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" > <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
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 : > 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
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
this is the *.m file of the case On 21 March 2018 at 02:01, André Silva 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 : > >> 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
- - - > - - > 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" : > >> 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 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 <%28630%29%20252-0219> >>> >>&g
Re: Power Flow in Matpower
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 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 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" : > >> 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 wrote: >> >> i tried that but unfortunately not work >> >> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." : >> >>> 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: * on behalf of >>> Mohammed Alhajri >>> *Reply-To: *MATPOWER discussion forum >>> *Date: Friday, March 16, 2018 at 11:26 AM* >>> *To: MATPOWER discussion forum >> >* >>> *Subject: Re: Power Flow in Matpower* >>> >>> >>> >>> any answer to this question? >>> >>> >>> >>> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" : >>> >>> 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
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 : > 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 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 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" : >> >>> 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 >>> wrote: >>> >>> i tried that but unfortunately not work >>> >>> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." >> >: >>> >>>> 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: * on behalf of >>>> Mohammed Alhajri >>>> *Reply-To: *MATPOWER discussion forum >>>> *Date: Friday, March 16, 2018 at 11:26 AM* >>>> *To: MATPOWER discussion forum >>> >* >>>> *Subject: Re: Power Flow in Matpower* >>>> >>>> >>>> >>>> any answer to this question? >>>> >>>> >>>> >>>> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" : >>>> >>>> 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
yes done, thank you very much بتاريخ ٢٠١٨/٠٣/٢١ ١٢:٥٨ م، كتب "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 : > >> 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 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 >>> 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" : >>> >>>> 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 >>>> wrote: >>>> >>>> i tried that but unfortunately not work >>>> >>>> بتاريخ ٢٠١٨/٠٣/١٦ ٨:٣٠ م، كتب "Abhyankar, Shrirang G." >>> >: >>>> >>>>> 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: * on behalf of >>>>> Mohammed Alhajri >>>>> *Reply-To: *MATPOWER discussion forum >>>>> *Date: Friday, March 16, 2018 at 11:26 AM* >>>>> *To: MATPOWER discussion forum >>>> >* >>>>> *Subject: Re: Power Flow in Matpower* >>>>> >>>>> >>>>> >>>>> any answer to this question? >>>>> >>>>> >>>>> >>>>> بتاريخ ٢٠١٨/٠٢/٢٥ ٧:١٨ م، كتب "Mohammed Alhajri" : >>>>> >>>>> 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
<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,
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>:
Re: Power Flow in Matpower
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 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 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" : >> >>> 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 <%28630%29%20252-0219> >>>> >>>> www.mcs.anl.gov/~abhyshr <http://www.mcs.anl.gov/%7Eabhyshr> >>>> >>>> >>>> >>>> >>>> >>>> >>>> >>>> *From: * 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 <
Re: Power Flow in Matpower
.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 > 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 > <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 >> <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" >> <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 >>>> www.mcs.anl.gov/~abhyshr <http://www.mcs.anl.gov/%7Eabhyshr> >>>> >>>> >>>> >>>> >>>> From: >>> <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
ou 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 >>> 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" : >>> >>>> 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 <%28630%29%20252-0219> >>>>> >>>>> www.mcs.anl.gov/~abhyshr <http://www.mcs.anl.gov/%7Eabhyshr> >>>>> >>>>> >>>>> >>>>> >>>>> *From: * 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
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 >> 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
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
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
Re: Help on the reactive power flow in matpower.
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
To: matpower-l
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
Re: Help on the reactive power flow in matpower.
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
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
> *To:* matpower-l
> *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.
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
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
> *To:* matpower-l
> *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.
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
> 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
> To: matpower-l
> 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.
Ray, Thank you very much.
Jerry.
On Wed, Apr 27, 2011 at 2:18 PM, Ray Zimmerman wrote:
> 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<http://www.pserc.cornell.edu/matpower/manual.pdf>
> .
>
> --
> 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 > 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
>> *To:* matpower-l
>> *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
>
>
>
>
>
>
--
Zhengrui (Jerry) Qin
Computer Science
College of William and Mary
Williamsburg, VA 23185
