In that case, simply subtract the amount of the reactive injection from the
reactive power load at the bus … and in that case you don’t have to do anything
about changing the bus type. At non-generator buses it should change the power
flow solution. At generator buses, it should just re-dispatch the corresponding
generator’s reactive power to offset your injection.
bus_idx = <anything from 1 to 14>;
mpc.bus(bus_idx, QD) = mpc.bus(bus_idx, QD) - VAr_injection;
Ray
> On Aug 17, 2016, at 1:11 PM, siddique tomal <[email protected]> wrote:
>
> Actually i wanted to run a program that will place a certain amount of
> reactive power in every Bus and will calculate the corresponding power loss.
> But you are telling that for case14 as it has just 5 generators so i cant
> place reactive power in all of them, isnt it??
> if so then really bad luck for me.
> Hope you have then some way to do this. Waiting to hear from you.
> Thanks
>
> On Wed, Aug 17, 2016 at 9:14 PM, Ray Zimmerman <[email protected]
> <mailto:[email protected]>> wrote:
> I think what you meant to do was something like the following …
>
> define_constants;
> mpc = loadcase('case14');
> gen_idx = 8;
> bus_idx = find(mpc.gen(gen_idx, GEN_BUS) == mpc.bus(:, BUS_I));
> mpc.bus(bus_idx, BUS_TYPE) = PQ;
> mpc.gen(gen_idx, QG) = 50;
> results = runpf(mpc);
>
> … except that case14 has only 5 generators, so you’ll have to set get_idx to
> something between 1 and 5.
>
> Ray
>
>
>
>
> > On Aug 17, 2016, at 10:08 AM, siddique tomal <[email protected]
> > <mailto:[email protected]>> wrote:
> >
> > Thanks but i got just new problem..
> >
> >
> > define_constants;
> > mpc=loadcase(case14);
> > mpc.gen(8,PV)=PQ;
> > mpc.gen(8,QG)=50;
> > rslt=runpf(mpc);
> >
> >
> >
> > and they are showing the following error:
> >
> > Subscript indices must either be real positive integers or logicals.
> >
> > Error in ext2int (line 175)
> > gs = ( mpc.gen(:, GEN_STATUS) > 0 & ... TCL MERGE ERROR (
> > 08/17/2016 11:14:46 ): "extra characters after close-quote"
> OutmailID: 120695916, List: 'matpower-l', MemberID: 76173427
> SCRIPT: "gen status > > Error in runpf (line 131) > mpc = ext2int(mpc); > >
> > > what is the wrong here going on ?? > > > On Wed, Aug 17, 2016 at 7:54
> PM, Ray Zimmerman <[email protected] <mailto:[email protected]>
> <mailto:[email protected] <mailto:[email protected]>>> wrote: >
> define_constants; > mpc.bus(BUS_NO, BUS_TYPE) = PQ; > > > >> On Aug 17,
> 2016, at 9:47 AM, siddique tomal <[email protected]
> <mailto:[email protected]> <mailto:[email protected]
> <mailto:[email protected]>>> wrote: >> >> Thank you Sir but still
> having problem. >> For setting a bus as a PQ bus i command the following: >>
> >> mpc.gen(BUS_NO,GEN_STATUS)= 1; >> >> is that command is correct for
> assigning a bus as a PQ? >> if so then i also have tried with that command
> like following: >> >>> define_constants; >>> mpc=loadcase(case14); >>>
> mpc.gen(BUS_NO,GEN_STATUS)=1; >>> mpc.gen(BUS_NO,QG)=45; >>>
> result=runpf(mpc); >> >> but still no change.. >> >> On Wed, Aug 17, 2016
> at 7:35 PM, Ray Zimmerman <[email protected] <mailto:[email protected]>
> <mailto:[email protected] <mailto:[email protected]>>> wrote: >> When running a
> power flow, the reactive injection of a generator at a PV bus is an output,
> solved by the power flow, not an input. If you want to set the reactive power
> injection as an input you have to make the bus a PQ bus and allow the voltage
> to deviate from the setpoint. >> >> Ray >> >> >>> On Aug 16, 2016, at
> 8:48 PM, siddique tomal <[email protected]
> <mailto:[email protected]> <mailto:[email protected]
> <mailto:[email protected]>>> wrote: >>> >>> Dear Sir, >>>
> define_constants; >>> mpc=loadcase(case14); >>> mpc.gen(BUS_NO,QG)=45; >>>
> result=runpf(mpc); >>> >>> After running this whenever i run the command
> losses=sum(real(get_losses(result))) it just shows the same result of power
> losses as before injecting the reactive power. >>> >>> What's the problem
> here i'm missing? Please need your help. >> >> > >
> --Apple-Mail=_C902470B-BB6B-4FAC-95F7-EE3FE09A5F85 Content-Transfer-Encoding:
> quoted-printable Content-Type: text/html; charset="utf-8"
> <html><head><meta http-equiv="Content-Type" content="text/html
> charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode:
> space; -webkit-line-break: after-white-space;" class="">I think what you
> meant to do was something like the following …<div class=""><br
> class=""></div><div class=""><div class=""><font face="Courier"
> class="">define_constants;</font></div><div class=""><font face="Courier"
> class="">mpc = loadcase('case14');</font></div><div class=""><font
> face="Courier" class="">gen_idx = 8;</font></div><div class=""><font
> face="Courier" class="">bus_idx = find(mpc.gen(gen_idx, GEN_BUS) ==
> mpc.bus(:, BUS_I));</font></div><div class=""><font face="Courier"
> class="">mpc.bus(bus_idx, BUS_TYPE) = PQ;</font></div><div class=""><font
> face="Courier" class="">mpc.gen(gen_idx, QG) = 50;</font></div><div
> class=""><font face="Courier" class="">results = runpf(mpc);</font></div><div
> class=""><br class=""></div><div class="">… except that case14 has only 5
> generators, so you’ll have to set <font face="Courier"
> class="">get_idx</font> to something between 1 and 5.</div><div class=""><br
> class=""></div><div class=""> Ray</div><div class=""><br
> class=""></div><div class=""><br class=""></div><div class=""><br
> class=""></div><div class=""><br class=""><div><blockquote type="cite"
> class=""><div class="">On Aug 17, 2016, at 10:08 AM, siddique tomal <<a
> href="mailto:[email protected] <mailto:[email protected]>"
> class="">[email protected] <mailto:[email protected]></a>>
> wrote:</div><br class="Apple-interchange-newline"><div class=""><div
> dir="ltr" class="">Thanks but i got just new problem..<div class=""><br
> class=""></div><div class=""><br class=""><b class=""><i
> class="">define_constants;<br class=""></i></b><div class=""><b class=""><i
> class="">mpc=loadcase(case14);</i></b></div><div class=""><b class=""><i
> class="">mpc.gen(8,PV)=PQ;</i></b></div><div class=""><b class=""><i
> class="">mpc.gen(8,QG)=50;</i></b></div><div class=""><b class=""><i
> class="">rslt=runpf(mpc);</i></b></div></div><div class=""><b class=""><i
> class=""><br class=""></i></b></div><div class=""><b class=""><i class=""><br
> class=""></i></b></div><div class=""><b class=""><i class=""><br
> class=""></i></b></div><div class="">and they are showing the following
> error:</div><div class=""><br class=""></div><div class=""><div
> class=""><font color="#ff0000" class="">Subscript indices must either be real
> positive integers or logicals.</font></div><div class=""><font
> color="#ff0000" class=""><br class=""></font></div><div class=""><font
> color="#ff0000" class="">Error in ext2int (line 175)</font></div><div
> class=""><font color="#ff0000" class="">
> gs = ( mpc.gen(:, GEN_STATUS) > 0 & ..."
> gen status</font></div><div class=""><font color="#ff0000" class=""><br
> class=""></font></div><div class=""><font color="#ff0000" class="">Error in
> runpf (line 131)</font></div><div class=""><font color="#ff0000" class="">mpc
> = ext2int(mpc);</font></div></div><div class=""><b class=""><i class=""><br
> class=""></i></b></div><div class=""><b class=""><i class=""><br
> class=""></i></b></div><div class=""><b class=""><i class=""><br
> class=""></i></b></div><div class="">what is the wrong here going on ??<br
> class=""><br class=""></div></div><div class="gmail_extra"><br class=""><div
> class="gmail_quote">On Wed, Aug 17, 2016 at 7:54 PM, Ray Zimmerman <span
> dir="ltr" class=""><<a href="mailto:[email protected]
> <mailto:[email protected]>" target="_blank" class="">[email protected]
> <mailto:[email protected]></a>></span> wrote:<br class=""><blockquote
> class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc
> solid;padding-left:1ex"><div style="word-wrap:break-word" class=""><font
> face="Courier" class="">define_constants;</font><div class=""><font
> face="Courier" class="">mpc.bus(BUS_NO, BUS_TYPE) = PQ;</font></div><div
> class=""><div class="h5"><div class=""><br class=""></div><div class=""><br
> class=""></div><div class=""><br class=""><div class=""><blockquote
> type="cite" class=""><div class="">On Aug 17, 2016, at 9:47 AM, siddique
> tomal <<a href="mailto:[email protected]
> <mailto:[email protected]>" target="_blank"
> class="">[email protected] <mailto:[email protected]></a>>
> wrote:</div><br class=""><div class=""><div dir="ltr" class="">Thank you Sir
> but still having problem.<div class="">For setting a bus as a PQ bus i
> command the following:</div><div class=""><br class=""></div><div class=""><b
> class=""><i class="">mpc.gen(BUS_NO,GEN_STATUS)= 1;</i></b></div><div
> class=""><br class=""></div><div class="">is that command is correct for
> assigning a bus as a PQ?<br class="">if so then i also have tried with that
> command like following:<br class=""><br class=""><blockquote type="cite"
> style="color:rgb(80,0,80);font-size:12.8px" class=""><div dir="ltr"
> class=""><div class=""><b class=""><i
> class="">define_constants;</i></b></div><div class=""><b class=""><i
> class="">mpc=loadcase(case14);</i></b></div><div class=""><b class=""><i
> class="">mpc.gen(BUS_NO,GEN_STATUS)=1;</i></b></div><div class=""><b
> class=""><i class="">mpc.gen(BUS_NO,QG)=45;</i></b></div><div class=""><b
> class=""><i
> class="">result=runpf(mpc);</i></b></div></div></blockquote></div><div
> class=""><div dir="ltr" class=""><div class="">but still no
> change..</div></div></div></div><div class="gmail_extra"><br class=""><div
> class="gmail_quote">On Wed, Aug 17, 2016 at 7:35 PM, Ray Zimmerman <span
> dir="ltr" class=""><<a href="mailto:[email protected]
> <mailto:[email protected]>" target="_blank" class="">[email protected]
> <mailto:[email protected]></a>></span> wrote:<br class=""><blockquote
> class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc
> solid;padding-left:1ex"><div style="word-wrap:break-word" class=""><div
> class="">When running a power flow, the reactive injection of a generator at
> a PV bus is an <b class="">output</b>, solved by the power flow, not an
> input. If you want to set the reactive power injection as an input you have
> to make the bus a PQ bus and allow the voltage to deviate from the
> setpoint.</div><span class=""><font color="#888888" class=""><div
> class=""><br class=""></div><div class="">
> Ray</div></font></span><div class=""><div class=""><div class=""><br
> class=""></div><br class=""><div class=""><blockquote type="cite"
> class=""><div class="">On Aug 16, 2016, at 8:48 PM, siddique tomal <<a
> href="mailto:[email protected] <mailto:[email protected]>"
> target="_blank" class="">[email protected]
> <mailto:[email protected]></a>> wrote:</div><br class=""><div
> class=""><div dir="ltr" class=""><div class="">Dear Sir,</div><div
> class=""><b class=""><i class="">define_constants;</i></b></div><div
> class=""><b class=""><i class="">mpc=loadcase(case14);</i></b></div><div
> class=""><b class=""><i class="">mpc.gen(BUS_NO,QG)=45;</i></b></div><div
> class=""><b class=""><i class="">result=runpf(mpc);</i></b></div><div
> class=""><b class=""><i class=""><br class=""></i></b></div><div
> class="">After running this whenever i run the command <i
> style="font-weight:bold" class="">losses=sum(real(get_losses(res<wbr
> class="">ult)))</i> it just shows the same result of power losses
> as before injecting the reactive power.<br class=""><br class="">What's the
> problem here i'm missing? Please need your help.</div></div>
> </div></blockquote></div><br
> class=""></div></div></div></blockquote></div><br class=""></div>
> </div></blockquote></div><br
> class=""></div></div></div></div></blockquote></div><br class=""></div>
> </div></blockquote></div><br class=""></div></div></body></html>
>
