Dear Dr. Ray Thanks very much for your support. I made use of the case file that you has sent me and I have changed the values of BR_R and BR_X to per unit values dividing them by Zbase and then the program converges after 5 iterations and I think the results are reasonable. I will try to run optimal power flow. Thanks
Amal Elramly ________________________________ From: Ray Zimmerman <[email protected]> To: MATPOWER discussion forum <[email protected]> Sent: Wednesday, November 28, 2012 6:16 PM Subject: Re: question Ok, I spent some time trying to get to something that converges … - there is no reference bus (set bus type to 3 for bus 1) - gencost format is incorrect, it should be polynomial (set first column to 2) - the bus matrix is missing the VM column If you fix these issues, then it is a valid, but infeasible problem. I believe the infeasibility is related to voltage and reactive power. In case it is useful, dividing BR_R by 1000, BR_X by 100 and setting RATE_A to 0 in the branch matrix and dividing QD by 10 in the bus matrix, yields a feasible problem. This may give you a useful starting point for this topology from which to begin investigating where the unfeasibility comes from. -- Ray Zimmerman Senior Research Associate 419A Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645 On Nov 28, 2012, at 3:47 AM, Amal Elramly <[email protected]> wrote: Dr. Ray > >yes all lines have non-zero reactance >The case file of the test bus system is attached below: > system MVA basempc.baseMVA = 1; >%% bus data >% bus_i type Pd Qd Gs Bs Vm Va baseKV zone Vmax >Vminmpc.bus = [ >1 2 0 0 0 0 1 0 12.6 1 1.05 0.95; >2 1 0.1 0.06 0 0 1 0 12.6 1 1.1 0.95; >3 1 0.09 0.04 0 0 1 0 12.6 1 1.05 0.95; >4 1 0.12 0.08 0 0 1 0 12.6 1 1.05 0.95; >5 1 0.06 0.03 0 0 1 0 12.6 1 1.05 0.95; >6 1 0.06 0.02 0 0 1 0 12.6 1 1.05 0.95; >7 1 0.2 0.1 0 0 1 0 12.6 1 1.05 0.95; >8 1 0.2 0.1 0 0 1 0 12.6 1 1.05 0.95; >9 1 0.06 0.02 0 0 1 0 12.6 1 1.05 0.95; >10 1 0.06 0.02 0 0 1 0 12.6 1 1.05 0.95; >11 1 0.045 0.03 0 0 1 0 12.6 1 1.05 0.95; >12 1 0.06 0.035 0 0 1 0 12.6 1 1.05 0.95; >13 1 0.06 0.035 0 0 1 0 12.6 1 1.1 0.95; >14 1 0.12 0.08 0 0 1 0 12.6 1 1.05 0.95; >15 1 0.06 0.01 0 0 1 0 12.6 1 1.05 0.95; >16 1 0.06 0.02 0 0 1 0 12.6 1 1.05 0.95; >17 1 0.06 0.02 0 0 1 0 12.6 1 1.05 0.95; >18 1 0.09 0.04 0 0 1 0 12.6 1 1.05 0.95; >19 1 0.09 0.04 0 0 1 0 12.6 1 1.05 0.95; >20 1 0.09 0.04 0 0 1 0 12.6 1 1.05 0.95; >21 1 0.09 0.04 0 0 1 0 12.6 1 1.05 0.95; >22 1 0.09 0.04 0 0 1 0 12.6 1 1.1 0.95; >23 1 0.42 0.05 0 0 1 0 12.6 1 1.1 0.95; >24 1 0.42 0.2 0 0 1 0 12.6 1 1.05 0.95; >25 1 0.06 0.2 0 0 1 0 12.6 1 1.05 0.95; >26 1 0.06 0.025 0 0 1 0 12.6 1 1.05 0.95; >27 1 0.06 0.025 0 0 1 0 12.6 1 1.1 0.95; >28 1 0.12 0.02 0 0 1 0 12.6 1 1.05 0.95; >29 1 0.2 0.07 0 0 1 0 12.6 1 1.05 0.95; >30 1 0.15 0.6 0 0 1 0 12.6 1 1.05 0.95; >31 1 0.21 0.07 0 0 1 0 12.6 1 1.05 0.95; >32 1 0.06 0.1 0 0 1 0 12.6 1 1.05 0.95; >33 1 0.06 0.04 0 0 1 0 12.6 1 1.05 0.95; >]; >%% generator data >% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 >Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apfmpc.gen >= [ >1 0.5 0 1 -0.5 1 100 1 6 0 0 0 0 0 0 0 0 0 0 > 0 0; >];%% branch data >% fbus tbus r x b rateA rateB rateC ratio angle status > angmin angmaxmpc.branch = [ >1 2 0.0922 0.0477 0 0.2 0.2 0.2 0 0 1 -360 360; >2 3 0.493 0.2511 0 0.1 0.1 0.1 0 0 1 -360 360; >3 4 0.3660 0.1864 0 0.2 0.2 0.2 0 0 1 -360 360; >4 5 0.3811 0.1941 0 0.1 0.1 0.1 0 0 1 -360 360; >5 6 0.8190 0.7070 0 0.1 0.1 0.1 0 0 1 -360 360; >6 7 0.1872 0.6188 0 0.3 0.3 0.3 0 0 1 -360 360; >7 8 1.7114 1.2351 0 0.3 0.3 0.3 0 0 1 -360 360; >8 9 1.0300 0.7400 0 0.1 0.1 0.1 0 0 1 -360 360; >9 10 1.0040 0.7400 0 0.1 0.1 0.1 0 0 1 -360 360; >10 11 0.1966 0.0650 0 0.1 0.1 0.1 0 0 1 -360 360; >11 12 0.3744 0.1238 0 0.1 0.1 0.1 0 0 1 -360 360; >12 13 1.4680 1.1550 0 0.1 0.1 0.1 0 0 1 -360 360; >13 14 0.5416 0.7129 0 0.2 0.2 0.2 0 0 1 -360 360; >14 15 0.5910 0.5260 0 0.1 0.1 0.1 0 0 1 -360 360; >15 16 0.7463 0.5450 0 0.1 0.1 0.1 0 0 1 -360 360; >16 17 1.2890 1.7210 0 0.1 0.1 0.1 0 0 1 -360 360; >17 18 0.7320 0.5740 0 0.1 0.1 0.1 0 0 1 -360 360; >2 19 0.1640 0.1565 0 0.1 0.1 0.1 0 0 1 -360 360; >19 20 1.5042 1.3554 0 0.1 0.1 0.1 0 0 1 -360 360; >20 21 0.4095 0.4784 0 0.1 0.1 0.1 0 0 1 -360 360; >21 22 0.7089 0.9373 0 0.1 0.1 0.1 0 0 1 -360 360; >3 23 0.4512 0.3083 0 0.5 0.5 0.5 0 0 1 -360 360; >23 24 0.8980 0.7091 0 0.5 0.5 0.5 0 0 1 -360 360; >24 25 0.8960 0.7011 0 0.3 0.3 0.3 0 0 1 -360 360; >6 26 0.2030 0.1034 0 0.1 0.1 0.1 0 0 1 -360 360; >26 27 0.2842 0.1447 0 0.1 0.1 0.1 0 0 1 -360 360; >27 28 1.0590 0.9337 0 0.2 0.2 0.2 0 0 1 -360 360; >28 29 0.8042 0.7006 0 0.3 0.3 0.3 0 0 1 -360 360; >29 30 0.5075 0.2585 0 0.7 0.7 0.7 0 0 1 -360 360; >30 31 0.9744 0.9630 0 0.3 0.3 0.3 0 0 1 -360 360; >31 32 0.3105 0.3619 0 0.2 0.2 0.2 0 0 1 -360 360; >32 33 0.3410 0.5302 0 0.1 0.1 0.1 0 0 1 -360 360; >]; >%%generator cost datampc.gencost = [ >1 0 0 3 0.02 2 0; >]; > >Amal Elramly > > > >________________________________ > From: Ray Zimmerman <[email protected]> >To: MATPOWER discussion forum <[email protected]> >Sent: Tuesday, November 27, 2012 4:52 PM >Subject: Re: question > > >Do all lines have non-zero reactance? Do you have a reference bus? > > >-- >Ray Zimmerman >Senior Research Associate >419A Warren Hall, Cornell University, Ithaca, NY 14853 >phone: (607) 255-9645 > > > > >On Nov 27, 2012, at 3:20 AM, Amal Elramly <[email protected]> wrote: > > >> >>I have run power flow for 33 radial test bus system, but the it doesn't >>converge and I do not know what is the problem >>could you help me ? >>Attached below the results of the program: >> >> >>MATPOWER Version 4.1, 14-Dec-2011 -- AC Power Flow (Newton) >>Warning: Matrix is singular to working precision. >>> In newtonpf at 109 >> In runpf at 220 >>Warning: Matrix is singular to working precision. >>> In newtonpf at 109 >> In runpf at 220 >>Warning: Matrix is singular to working precision. >>> In newtonpf at 109 >> In runpf at 220 >>Warning: Matrix is singular to working precision. >>> In newtonpf at 109 >> In runpf at 220 >>Warning: Matrix is singular to working precision. >>> In newtonpf at 109 >> In runpf at 220 >>Warning: Matrix is singular to working precision. >>> In newtonpf at 109 >> In runpf at 220 >>Warning: Matrix is singular to working precision. >> >>Newton's method power did not converge in 10 iterations. >>Did not converge (0.09 seconds) >>================================================================================ >>| System Summary | >>================================================================================ >>How many? How much? P (MW) Q (MVAr) >>--------------------- ------------------- ------------- ----------------- >>Buses 33 Total Gen Capacity 6.0 -0.5 to 1.0 >>Generators 1 On-line Capacity 6.0 -0.5 to 1.0 >>Committed Gens 1 Generation (actual) NaN NaN >>Loads 32 Load 3.7 2.3 >> Fixed 32 Fixed 3.7 2.3 >> Dispatchable 0 Dispatchable -0.0 of -0.0 -0.0 >>Shunts 0 Shunt (inj) -0.0 0.0 >>Branches 32 Losses (I^2 * Z) NaN NaN >>Transformers 0 Branch Charging (inj) - NaN >>Inter-ties 0 Total Inter-tie Flow 0.0 0.0 >>Areas 1 >> Minimum Maximum >> ------------------------- -------------------------------- >>Voltage Magnitude NaN p.u. @ bus 1 NaN p.u. @ bus 1 >>Voltage Angle NaN deg @ bus 1 NaN deg @ bus 1 >>P Losses (I^2*R) - NaN MW @ line 1-2 >>Q Losses (I^2*X) - NaN MVAr @ line 1-2 >>================================================================================ >>| Bus Data | >>================================================================================ >> Bus Voltage Generation Load >> # Mag(pu) Ang(deg) P (MW) Q (MVAr) P (MW) Q (MVAr) >>----- ------- -------- -------- -------- -------- -------- >> 1 NaN NaN NaN NaN - - >> 2 NaN NaN - - 0.10 0.06 >> 3 NaN NaN - - 0.09 0.04 >> 4 NaN NaN - - 0.12 0.08 >> 5 NaN NaN - - 0.06 0.03 >> 6 NaN NaN - - 0.06 0.02 >> 7 NaN NaN - - 0.20 0.10 >> 8 NaN NaN - - 0.20 0.10 >> 9 NaN NaN - - 0.06 0.02 >> 10 NaN NaN - - 0.06 0.02 >> 11 NaN NaN - - 0.04 0.03 >> 12 NaN NaN - - 0.06 0.04 >> 13 NaN NaN - - 0.06 0.04 >> 14 NaN NaN - - 0.12 0.08 >> 15 NaN NaN - - 0.06 0.01 >> 16 NaN NaN - - 0.06 0.02 >> 17 NaN NaN - - 0.06 0.02 >> 18 NaN NaN - - 0.09 0.04 >> 19 NaN NaN - - 0.09 0.04 >> 20 NaN NaN - - 0.09 0.04 >> 21 NaN NaN - - 0.09 0.04 >> 22 NaN NaN - - 0.09 0.04 >> 23 NaN NaN - - 0.42 0.05 >> 24 NaN NaN - - 0.42 0.20 >> 25 NaN NaN - - 0.06 0.20 >> 26 NaN NaN - - 0.06 0.03 >> 27 NaN NaN - - 0.06 0.03 >> 28 NaN NaN - - 0.12 0.02 >> 29 NaN NaN - - 0.20 0.07 >> 30 NaN NaN - - 0.15 0.60 >> 31 NaN NaN - - 0.21 0.07 >> 32 NaN NaN - - 0.06 0.10 >> 33 NaN NaN - - 0.06 0.04 >> -------- -------- -------- -------- >> Total: NaN NaN 3.69 2.30 >>================================================================================ >>| Branch >>Data | >>================================================================================ >>Brnch From To From Bus Injection To Bus Injection Loss (I^2 * >>Z) >> # Bus Bus P (MW) Q (MVAr) P (MW) Q (MVAr) P (MW) Q (MVAr) >>----- ----- ----- -------- -------- -------- -------- -------- >>-------- >> 1 1 2 NaN NaN NaN NaN NaN NaN >> 2 2 3 NaN NaN NaN NaN NaN NaN >> 3 3 4 NaN NaN NaN NaN NaN NaN >> 4 4 5 NaN NaN NaN NaN NaN NaN >> 5 5 6 NaN NaN NaN NaN NaN NaN >> 6 6 7 NaN NaN NaN NaN NaN NaN >> 7 7 8 NaN NaN NaN NaN NaN NaN >> 8 8 9 NaN NaN NaN NaN NaN NaN >> 9 9 10 NaN NaN NaN NaN NaN NaN >> 10 10 11 NaN NaN NaN NaN NaN NaN >> 11 11 12 NaN NaN NaN NaN NaN NaN >> 12 12 13 NaN NaN NaN NaN NaN NaN >> 13 13 14 NaN NaN NaN NaN NaN NaN >> 14 14 15 NaN NaN NaN NaN NaN NaN >> 15 15 16 NaN NaN NaN NaN NaN NaN >> 16 16 17 NaN NaN NaN NaN NaN NaN >> 17 17 18 NaN NaN NaN NaN NaN NaN >> 18 2 19 NaN NaN NaN NaN NaN NaN >> 19 19 20 NaN NaN NaN NaN NaN NaN >> 20 20 21 NaN NaN NaN NaN NaN NaN >> 21 21 22 NaN NaN NaN NaN NaN NaN >> 22 3 23 NaN NaN NaN NaN NaN NaN >> 23 23 24 NaN NaN NaN NaN NaN NaN >> 24 24 25 NaN NaN NaN NaN NaN NaN >> 25 6 26 NaN NaN NaN NaN NaN NaN >> 26 26 27 NaN NaN NaN NaN NaN NaN >> 27 27 28 NaN NaN NaN NaN NaN NaN >> 28 28 29 NaN NaN NaN NaN NaN NaN >> 29 29 30 NaN NaN NaN NaN NaN NaN >> 30 30 31 NaN NaN NaN NaN NaN NaN >> 31 31 32 NaN NaN NaN NaN NaN NaN >> 32 32 33 NaN NaN NaN NaN NaN NaN >> -------- -------- >> Total: NaN NaN >> >> >>Amal Elramly >> >> >> > > >
