Firstly I would like to excuse for the triple posting. I loosed the reactive limits of every generator to [-2.5 MVar , 2.5 MVar] and it converges.
Thank you for your answers. Alexandros Gesios 2016-11-24 14:33 GMT+02:00 Alexandros Gesios <alexandrosges...@gmail.com>: > I uploaded a file which shows the network configuration > > 2016-11-24 13:45 GMT+02:00 Alexandros Gesios <alexandrosges...@gmail.com>: > >> The bus 1 has constant voltage and it is considered as an infinite bus. >> The negative active power means that PVs send power to the grid. >> >> Στις 24 Νοε 2016 13:33, ο χρήστης "Saranya A" <asar...@gmail.com> έγραψε: >> >>> I dont see in the case file where you have connected the Bus 1 to the >>> grid. >>> Can you clarify on this? >>> >>> On Thu, Nov 24, 2016 at 1:31 PM, Alexandros Gesios < >>> alexandrosges...@gmail.com> wrote: >>> >>>> Bus 1 represents the grid connection which has negative active power >>>> because the active power from PV is larger that the load active power. >>>> >>>> Στις 24 Νοε 2016 07:38, ο χρήστης "Saranya A" <asar...@gmail.com> >>>> έγραψε: >>>> >>>>> I ran the "casePV3" and observed that bus 1 which is slack bus has Pg >>>>> as -3MW. Can a bus which has photovoltaic unit as generator observe real >>>>> power? I am not pretty aware of the power electronics of it. I think if >>>>> you >>>>> can reason that out you may get some direction. >>>>> >>>>> Regards, >>>>> Saranya A >>>>> >>>>> On Thu, Nov 24, 2016 at 4:08 AM, Alexandros Gesios < >>>>> alexandrosges...@gmail.com> wrote: >>>>> >>>>>> Ι made some changes in the matrices and OPF completes but I do not >>>>>> think that it allocates the reactive power among the PVs in an optimal >>>>>> way >>>>>> regarding the losses. >>>>>> The last three PVs absorb Qm reactive power while the other operate >>>>>> in a unity power factor. >>>>>> I attach the case in order to be more specific. Also the network is >>>>>> radial. >>>>>> >>>>>> >>>>>> 2016-11-23 20:46 GMT+02:00 Alexandros Gesios < >>>>>> alexandrosges...@gmail.com>: >>>>>> >>>>>>> Dr. Zimmerman >>>>>>> >>>>>>> >>>>>>> I applied the change you have mentioned and the result is the following: >>>>>>> >>>>>>> >>>>>>> MATPOWER Version 6.0b1, 01-Jun-2016 -- AC Optimal Power Flow >>>>>>> MATLAB Interior Point Solver -- MIPS, Version 1.2.1, 01-Jun-2016 >>>>>>> (using built-in linear solver) >>>>>>> Numerically Failed >>>>>>> >>>>>>> Did not converge in 9 iterations. >>>>>>> >>>>>>> >>>>> Did NOT converge (0.11 seconds) <<<<< >>>>>>> >>>>>>> >>>>>>> >>>>>>> >>>>>>> 2016-11-23 19:54 GMT+02:00 Ray Zimmerman <r...@cornell.edu>: >>>>>>> >>>>>>>> I think you want to change the first column of your gencost entries >>>>>>>> to 2 (polynomial) instead of 1 (piecewise linear). Hopefully, that >>>>>>>> fixes >>>>>>>> your problem and it will solve. >>>>>>>> >>>>>>>> Ray >>>>>>>> >>>>>>>> > On Nov 23, 2016, at 10:30 AM, Alexandros Gesios < >>>>>>>> alexandrosges...@gmail.com> wrote: >>>>>>>> > >>>>>>>> > Hello to everyone >>>>>>>> > >>>>>>>> > I am new to MATPOWER and I am trying to operate an OPF analysis >>>>>>>> in a MV grid constisted of two feeders. >>>>>>>> > The first feeder has 10 inductive loads and the second feeder has >>>>>>>> 10 photovoltaic units and their are connected to the grid by the slack >>>>>>>> bus. >>>>>>>> > What I want to achieve is loss minimization by allocating the >>>>>>>> reactive power among the photovoltaic units in a optimal way. The >>>>>>>> nominal >>>>>>>> active power of the PV units is 10.6 MW and the load is 6 MW so the >>>>>>>> extra >>>>>>>> active power is being absorbed by the HV/MV Grid. >>>>>>>> > Also I want the slack bus to operate at constant 1.05 pu voltage >>>>>>>> so in the mpc.bus slack bus has voltage limits Vmax=Vmin=1.05. >>>>>>>> > Also I want the PV units to operate in nominal power so in >>>>>>>> mpc.gen their active power limits are Pmax=Pmin=Pnom. >>>>>>>> > >>>>>>>> > Finally the mpc.gencost matrix constists of 11 rows and every row >>>>>>>> has these elements 1 0 0 2 2 4; because I >>>>>>>> read >>>>>>>> in the archive that in order to achieve loss minimization you need the >>>>>>>> same >>>>>>>> linear costs in the gencost matrix. >>>>>>>> > >>>>>>>> > When I run runopf(mycase) I take the following errors >>>>>>>> > >>>>>>>> > Index exceeds matrix dimensions. >>>>>>>> > >>>>>>>> > Error in opf_setup (line 84) >>>>>>>> > x1 = mpc.gencost(pwl1, COST+2); >>>>>>>> > >>>>>>>> > Error in opf (line 198) >>>>>>>> > om = opf_setup(mpc, mpopt); >>>>>>>> > >>>>>>>> > Error in runopf (line 75) >>>>>>>> > [r, success] = opf(casedata, mpopt); >>>>>>>> > >>>>>>>> > Also, I would like to ask whether my problem can be formulated in >>>>>>>> MATPOWER. >>>>>>>> > >>>>>>>> > >>>>>>>> > Thank you in advance >>>>>>>> > >>>>>>>> > >>>>>>>> > Alexandros Gesios >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>> >>>>>> >>>>> >>> >
