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 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
>>
>>
>> 
>
>
>  

Reply via email to