Dear all, I'm getting hard times to solve a power flow on a real Italian medium voltage network. MATPOWER does not converge with any algorithm. I'm wondering if the branch data must have a certain order or my problem is simply the R/X ratio. The casefile is the following:
function mpc = LFMT mpc.version = 2; mpc.baseMVA = 10; mpc.bus=[ 1 3 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 2 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9; 3 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9; 4 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9; 5 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9; 6 1 0.05400 0.02615 0 0 1 1.0 0 20 1 1.1 0.9; 7 1 0.08640 0.04185 0 0 1 1.0 0 20 1 1.1 0.9; 8 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9; 9 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9; 10 1 0.02700 0.01308 0 0 1 1.0 0 20 1 1.1 0.9; 11 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 12 2 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 13 2 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 14 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 15 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 16 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 17 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 18 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 19 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 20 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 21 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 22 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 23 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 24 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 25 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 26 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 27 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; 28 1 0.00000 0.00000 0 0 1 1.0 0 20 1 1.1 0.9; ]; mpc.gen=[ 12 0.99440 0.00000 100 -100 1.0 100.000 1 1000.0 0.0 0 0 0 0 0 0 0 0 0 0 0; 13 0.99440 0.00000 100 -100 1.0 100.000 1 1000.0 0.0 0 0 0 0 0 0 0 0 0 0 0; 1 0.00000 0.00000 100 -100 1.0 100.000 1 1000.0 0.0 0 0 0 0 0 0 0 0 0 0 0; ]; mpc.branch=[ 19 27 0.01430 0.01060 0.01240 50.0 50.0 50.0 0 0 1 360 -360; 27 10 0.00220 0.00120 0.00120 50.0 50.0 50.0 0 0 1 360 -360; 27 15 0.00320 0.00240 0.00280 50.0 50.0 50.0 0 0 1 360 -360; 15 24 0.00840 0.00630 0.00740 50.0 50.0 50.0 0 0 1 360 -360; 24 16 0.00520 0.00390 0.00450 50.0 50.0 50.0 0 0 1 360 -360; 24 4 0.00040 0.00010 0.00820 50.0 50.0 50.0 0 0 1 360 -360; 16 25 0.00310 0.00230 0.00270 50.0 50.0 50.0 0 0 1 360 -360; 25 28 0.00170 0.00130 0.00150 50.0 50.0 50.0 0 0 1 360 -360; 25 5 0.00160 0.00090 0.00090 50.0 50.0 50.0 0 0 1 360 -360; 28 13 0.00400 0.00290 0.04580 50.0 50.0 50.0 0 0 1 360 -360; 28 14 0.00070 0.00020 0.01370 50.0 50.0 50.0 0 0 1 360 -360; 13 6 0.01320 0.00980 0.12280 50.0 50.0 50.0 0 0 1 360 -360; 6 26 0.00840 0.00630 0.00740 50.0 50.0 50.0 0 0 1 360 -360; 26 17 0.00130 0.00100 0.00110 50.0 50.0 50.0 0 0 1 360 -360; 26 8 0.00740 0.00420 0.00410 50.0 50.0 50.0 0 0 1 360 -360; 17 11 0.00900 0.00640 0.09160 50.0 50.0 50.0 0 0 1 360 -360; 11 18 0.00950 0.00710 0.05630 50.0 50.0 50.0 0 0 1 360 -360; 18 23 0.00060 0.00050 0.00060 50.0 50.0 50.0 0 0 1 360 -360; 23 3 0.02320 0.01320 0.01300 50.0 50.0 50.0 0 0 1 360 -360; 23 20 0.00520 0.00390 0.00450 50.0 50.0 50.0 0 0 1 360 -360; 20 21 0.01950 0.01450 0.01700 50.0 50.0 50.0 0 0 1 360 -360; 21 22 0.01490 0.01110 0.01300 50.0 50.0 50.0 0 0 1 360 -360; 22 7 0.00700 0.00520 0.06060 50.0 50.0 50.0 0 0 1 360 -360; 22 9 0.02480 0.01470 0.01490 50.0 50.0 50.0 0 0 1 360 -360; 7 2 0.01130 0.00650 0.00630 50.0 50.0 50.0 0 0 1 360 -360; 1 19 0.01970 0.01460 0.02750 50.0 50.0 50.0 0 0 1 360 -360; Thank you.
