I'm guessing it is infeasible. Did you start with the feasible problem I 
mentioned in my response from Nov 28th?

-- 
Ray Zimmerman
Senior Research Associate
419A Warren Hall, Cornell University, Ithaca, NY 14853
phone: (607) 255-9645




On Dec 23, 2012, at 3:49 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
>>> 
>>>  
>> 
>> 
>> 
> 
> 
> 
> 
> 
> 
> 

Reply via email to