Let's take this off-list. -- Ray Zimmerman Senior Research Associate 419A Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645
On Jan 8, 2013, at 3:41 AM, Amal Elramly <[email protected]> wrote: > Dr. Ray > > I'm so sorry as I did not see your reply on Jan 2nd. I recieve a message > containing "You are out of office". Sorry again for the inconvenience. > > Concerning my question, Yes I started with the feasible problem you > mentioned. It runs for power flow but it does not run for optimal power flow. > > Thanks, > > Amal Elramly > > > From: Ray Zimmerman <[email protected]> > To: MATPOWER discussion forum <[email protected]> > Cc: Amal Elramly <[email protected]> > Sent: Monday, January 7, 2013 6:45 PM > Subject: Re: question > > As I said in my reply on Jan 2nd … > > I'm guessing it is infeasible. Did you start with the feasible problem I > mentioned in my response from Nov 28th? > > Please follow through on my suggestions before continuing to post the same > question. > > -- > Ray Zimmerman > Senior Research Associate > 419A Warren Hall, Cornell University, Ithaca, NY 14853 > phone: (607) 255-9645 > > > > > On Jan 6, 2013, at 6:09 AM, Amal Elramly <[email protected]> wrote: > >> >> >> Dear Dr. Ray >> >> I have tried to run optimal power flow for the 33 radial test bus system but >> it does not converge. Knowing that the case file works for power flow >> (runpf). >> I can not determine the source of the problem. Could you help me? >> >> Attached below the program warnings and results: >> >> >> runopf('case_amal5') >> MATPOWER Version 4.1, 14-Dec-2011 -- AC Optimal Power Flow >> MATLAB Interior Point Solver -- MIPS, Version 1.0, 07-Feb-2011Warning: >> Matrix is close to singular or badly scaled. >> Results may be inaccurate. RCOND = 5.247763e-018. >> > In mips at 422 >> In mipsopf_solver at 145 >> In opf_execute at 106 >> In opf at 225 >> In runopf at 96 >> Warning: Matrix is close to singular or badly scaled. >> Results may be inaccurate. RCOND = 3.349638e-022. >> > In mips at 422 >> In mipsopf_solver at 145 >> In opf_execute at 106 >> In opf at 225 >> In runopf at 96 >> Numerically Failed >> Did not converge in 13 iterations. >> Did not converge (0.28 seconds) >> Objective Function Value = 4.40 $/hr >> ================================================================================ >> | 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) 2.2 0.5 >> 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) 0.13 0.09 >> Transformers 0 Branch Charging (inj) - 0.0 >> Inter-ties 0 Total Inter-tie Flow 0.0 0.0 >> Areas 1 >> Minimum Maximum >> ------------------------- -------------------------------- >> Voltage Magnitude 1.506 p.u. @ bus 18 1.532 p.u. @ bus 2 >> Voltage Angle -0.16 deg @ bus 18 0.09 deg @ bus 30 >> P Losses (I^2*R) - 0.04 MW @ line 1-2 >> Q Losses (I^2*X) - 0.02 MVAr @ line 1-2 >> Lambda P -6674.13 $/MWh @ bus 1 15413474426196144000.00 $/MWh @ >> bus 18 >> Lambda Q 4445.12 $/MWh @ bus 1 8466416009517438000.00 $/MWh @ >> bus 33 >> ================================================================================ >> | Bus Data >> | >> ================================================================================ >> Bus Voltage Generation Load >> Lambda($/MVA-hr) >> # Mag(pu) Ang(deg) P (MW) Q (MVAr) P (MW) Q (MVAr) P >> Q >> ----- ------- -------- -------- -------- -------- -------- ------- >> ------- >> 1 1.526 0.000* 2.15 0.48 - - >> -6674.1314445.122 >> 2 1.532 -0.001 - - 0.10 0.06 >> 482471203417005.810249602520478259.500 >> 3 1.523 0.008 - - 0.09 0.04 >> 15072572746478563000.0008185918922953839600.000 >> 4 1.522 0.017 - - 0.12 0.08 >> 15104025286329487000.0008207876480607969300.000 >> 5 1.520 0.026 - - 0.06 0.03 >> 15134463712170088000.0008229215530480105500.000 >> 6 1.516 0.015 - - 0.06 0.02 >> 15206029152644973000.0008282433372989881300.000 >> 7 1.515 -0.023 - - 0.20 0.10 >> 15225286235917232000.0008290320354694113300.000 >> 8 1.512 -0.060 - - 0.20 0.10 >> 15296769731624745000.0008314283903063420900.000 >> 9 1.511 -0.079 - - 0.06 0.02 >> 15328004445576665000.0008322680628071578600.000 >> 10 1.510 -0.096 - - 0.06 0.02 >> 15355208175207971000.0008329745091343159300.000 >> 11 1.510 -0.097 - - 0.04 0.03 >> 15359180129512268000.0008330590821724995600.000 >> 12 1.510 -0.099 - - 0.06 0.04 >> 15365913344698071000.0008331654368793930800.000 >> 13 1.515 -0.122 - - 0.06 0.04 >> 15286706390876719000.0008252011468596544500.000 >> 14 1.508 -0.136 - - 0.12 0.08 >> 15376665231265518000.0008357579409955088400.000 >> 15 1.508 -0.142 - - 0.06 0.01 >> 15386528993294068000.0008360997889585863700.000 >> 16 1.507 -0.146 - - 0.06 0.02 >> 15395689234238870000.0008364574717577316400.000 >> 17 1.506 -0.160 - - 0.06 0.02 >> 15409536653083820000.0008370346040024593400.000 >> 18 1.506 -0.162 - - 0.09 0.04 >> 15413474426196144000.0008372131740058532900.000 >> 19 1.526 -0.003 - - 0.09 0.04 >> 483409380339560.500250432571149656.120 >> 20 1.525 -0.021 - - 0.09 0.04 >> 483931362873987.250250466952173264.060 >> 21 1.525 -0.025 - - 0.09 0.04 >> 484020789441837.310250453373003567.310 >> 22 1.531 -0.029 - - 0.09 0.04 >> 481787513781547.190247382801145038.620 >> 23 1.528 0.002 - - 0.42 0.05 >> 14979643414630095000.0008118178854098206700.000 >> 24 1.521 0.004 - - 0.42 0.20 >> 15106484506421404000.0008219389687429146600.000 >> 25 1.521 0.014 - - 0.06 0.20 >> 15109715780803389000.0008229638628413411300.000 >> 26 1.516 0.022 - - 0.06 0.03 >> 15211131578805922000.0008289737720668774400.000 >> 27 1.521 0.031 - - 0.06 0.03 >> 15098406719421895000.0008238958818964934700.000 >> 28 1.512 0.047 - - 0.12 0.02 >> 15249229614526636000.0008384065566448683000.000 >> 29 1.510 0.064 - - 0.20 0.07 >> 15283185053805144000.0008427418723589714900.000 >> 30 1.509 0.086 - - 0.15 0.60 >> 15296412867248689000.0008449006203481891800.000 >> 31 1.508 0.077 - - 0.21 0.07 >> 15317078797228669000.0008462351132649569300.000 >> 32 1.507 0.077 - - 0.06 0.10 >> 15319642074435744000.0008465341544258864100.000 >> 33 1.507 0.076 - - 0.06 0.04 >> 15321213791940397000.0008466416009517438000.000 >> -------- -------- -------- -------- >> Total: 2.15 0.48 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 -11.78 -6.13 11.83 6.16 0.044 >> 0.02 >> 2 2 3 3.41 1.85 -3.39 -1.84 0.020 >> 0.01 >> 3 3 4 0.78 0.55 -0.78 -0.55 0.001 >> 0.00 >> 4 4 5 0.72 0.51 -0.72 -0.51 0.001 >> 0.00 >> 5 5 6 0.69 0.51 -0.69 -0.51 0.002 >> 0.00 >> 6 6 7 0.44 0.17 -0.44 -0.17 0.000 >> 0.00 >> 7 7 8 0.35 0.11 -0.35 -0.11 0.001 >> 0.00 >> 8 8 9 0.25 0.06 -0.25 -0.06 0.000 >> 0.00 >> 9 9 10 0.22 0.05 -0.22 -0.05 0.000 >> 0.00 >> 10 10 11 0.19 0.04 -0.19 -0.04 0.000 >> 0.00 >> 11 11 12 0.17 0.03 -0.17 -0.03 0.000 >> 0.00 >> 12 12 13 -0.45 -0.46 0.46 0.46 0.002 >> 0.00 >> 13 13 14 1.15 1.35 -1.14 -1.34 0.005 >> 0.01 >> 14 14 15 0.13 0.05 -0.13 -0.05 0.000 >> 0.00 >> 15 15 16 0.10 0.04 -0.10 -0.04 0.000 >> 0.00 >> 16 16 17 0.07 0.03 -0.07 -0.03 0.000 >> 0.00 >> 17 17 18 0.04 0.02 -0.04 -0.02 0.000 >> 0.00 >> 18 2 19 4.67 4.35 -4.65 -4.34 0.018 >> 0.02 >> 19 19 20 0.09 0.01 -0.09 -0.01 0.000 >> 0.00 >> 20 20 21 0.05 -0.01 -0.05 0.01 0.000 >> 0.00 >> 21 21 22 -0.72 -0.99 0.72 0.99 0.003 >> 0.00 >> 22 3 23 -1.85 -1.35 1.86 1.36 0.006 >> 0.00 >> 23 23 24 1.21 0.97 -1.21 -0.96 0.006 >> 0.00 >> 24 24 25 0.03 0.10 -0.03 -0.10 0.000 >> 0.00 >> 25 6 26 0.21 0.33 -0.21 -0.33 0.000 >> 0.00 >> 26 26 27 -3.76 -1.69 3.77 1.70 0.013 >> 0.01 >> 27 27 28 1.15 1.11 -1.15 -1.10 0.007 >> 0.01 >> 28 28 29 0.34 0.43 -0.34 -0.43 0.001 >> 0.00 >> 29 29 30 0.24 0.39 -0.24 -0.39 0.000 >> 0.00 >> 30 30 31 0.16 0.10 -0.16 -0.10 0.000 >> 0.00 >> 31 31 32 0.06 0.07 -0.06 -0.07 0.000 >> 0.00 >> 32 32 33 0.03 0.02 -0.03 -0.02 0.000 >> 0.00 >> -------- >> -------- >> Total: 0.130 >> 0.09 >> ================================================================================ >> | Voltage Constraints >> | >> ================================================================================ >> Bus # Vmin mu Vmin |V| Vmax Vmax mu >> ----- -------- ----- ----- ----- -------- >> 1 6246.856 0.950 1.526 1.0501272740518625929200.000 >> 2 6302.525 0.950 1.532 1.100355570291.153 >> 3 6286.023 0.950 1.523 1.050234659548248.324 >> 4 6304.242 0.950 1.522 1.05070211654568.559 >> 5 6322.540 0.950 1.520 1.05029853950572.320 >> 6 6369.737 0.950 1.516 1.0506097293404.070 >> 7 6379.032 0.950 1.515 1.0504859694314.642 >> 8 6415.832 0.950 1.512 1.0502062405473.305 >> 9 6433.196 0.950 1.511 1.0501426960223.656 >> 10 6449.101 0.950 1.510 1.050971125114.997 >> 11 6451.543 0.950 1.510 1.050907364464.729 >> 12 6455.813 0.950 1.510 1.050785716565.193 >> 13 6529.122 0.950 1.515 1.10085464084.004 >> 14 6480.352 0.950 1.508 1.0502838355.232 >> 15 6484.456 0.950 1.508 1.050272204085.475 >> 16 6488.425 0.950 1.507 1.050279697729.872 >> 17 6494.291 0.950 1.506 1.050327620045.383 >> 18 6496.046 0.950 1.506 1.050333944214.087 >> 19 6253.843 0.950 1.526 1.0503008201144580546.000 >> 20 6262.180 0.950 1.525 1.0505934815489843.966 >> 21 6263.675 0.950 1.525 1.0501316284247402.528 >> 22 6315.013 0.950 1.531 1.100325795191.182 >> 23 6345.784 0.950 1.528 1.100251631455.982 >> 24 6307.177 0.950 1.521 1.05048890166499.626 >> 25 6310.559 0.950 1.521 1.05038714021061.191 >> 26 6374.732 0.950 1.516 1.0502324115755.694 >> 27 6434.834 0.950 1.521 1.100139848954.055 >> 28 6413.029 0.950 1.512 1.0501912131557.960 >> 29 6434.480 0.950 1.510 1.0501156823231.291 >> 30 6442.778 0.950 1.509 1.050978580603.251 >> 31 6452.337 0.950 1.508 1.050818185087.477 >> 32 6453.938 0.950 1.507 1.050795031950.793 >> 33 6454.696 0.950 1.507 1.050784387575.451 >> ================================================================================ >> | Generation Constraints >> | >> ================================================================================ >> Gen Bus Active Power Limits >> # # Pmin mu Pmin Pg Pmax Pmax mu >> ---- ----- ------- -------- -------- -------- ------- >> 1 1 10275.446 0.00 2.15 6.00 3600.015 >> Gen Bus Reactive Power Limits >> # # Qmin mu Qmin Qg Qmax Qmax mu >> --- --- ------- -------- -------- -------- ------- >> 1 16364.618 -0.50 0.48 1.0010809.675 >> ================================================================================ >> | Branch Flow Constraints >> | >> ================================================================================ >> Brnch From "From" End Limit "To" End To >> # Bus |Sf| mu |Sf| |Smax| |St| |St| mu Bus >> ----- ----- ------- -------- -------- -------- ------- ----- >> 1 1 4.887 13.28 0.10 13.33 4.342 2 >> 2 2 0.136 3.88 0.10 3.85444868379062962690.000 >> 3 >> 3 38857013.775 0.96 0.20 0.968980790.433 4 >> 4 4190975.730 0.89 0.10 0.89196199.302 5 >> 5 5 95220.586 0.86 0.10 0.85100930.224 6 >> 6 6 28641.260 0.47 0.30 0.47 28714.067 7 >> 7 7 19558.624 0.37 0.30 0.37 19668.295 8 >> 8 8 5819.117 0.26 0.10 0.26 5825.381 9 >> 9 9 5534.336 0.23 0.10 0.23 5538.127 10 >> 10 10 5302.956 0.20 0.10 0.20 5303.362 11 >> 11 11 5150.910 0.17 0.10 0.17 5151.381 12 >> 12 12 2172.596 0.64 0.10 0.64 1984.649 13 >> 13 13 2250.949 1.77 0.20 1.76 2649.396 14 >> 14 14 5063.615 0.14 0.10 0.14 5064.110 15 >> 15 15 4945.511 0.11 0.10 0.11 4945.803 16 >> 16 16 4854.824 0.08 0.10 0.08 4855.038 17 >> 17 17 4795.020 0.05 0.10 0.05 4795.043 18 >> 18 2 22.331 6.38 0.10 6.36 161.524 19 >> 19 19 4882.909 0.10 0.10 0.10 4883.192 20 >> 20 20 4782.472 0.05 0.10 0.05 4782.467 21 >> 21 21 1531.158 1.22 0.10 1.23 1351.125 22 >> 22 3 3636.219 2.29 0.50 2.30 3132.169 23 >> 23 23 7914.348 1.55 0.50 1.54 9907.350 24 >> 24 24 13363.654 0.10 0.30 0.10 13364.250 25 >> 25 6 6880.566 0.39 0.10 0.39 6885.089 26 >> 26 26 483.983 4.12 0.10 4.14 333.084 27 >> 27 27 17981.684 1.60 0.20 1.59176578.511 28 >> 28 28 45384.414 0.55 0.30 0.55 45881.876 29 >> 29 29 35066.816 0.46 0.70 0.46 35143.551 30 >> 30 30 14458.864 0.19 0.30 0.19 14465.501 31 >> 31 31 9401.831 0.09 0.20 0.09 9401.992 32 >> 32 32 4778.881 0.04 0.10 0.04 4778.887 33 >> Thanks for your patience and time >> Amal Elramly >> >> From: Amal Elramly <[email protected]> >> To: MATPOWER discussion forum <[email protected]> >> Sent: Sunday, December 9, 2012 11:06 AM >> Subject: Fw: question >> >> >> >> >> 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 data >>> mpc.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 >>>> >>>> >>> >>> >>> >> >> >> >> >> >> >> >> >> > > >
