Convergence of the Algorithm
As can be seen from Table 4.3 that a total number of 4 real and 3 reactive powers are known to us. We must then calculate each of these from (4.6) and (4.7) using the values of the voltage magnitudes and their angle obtained after each iteration. The power mismatches are then calculated from (4.9) and (4.10). The process is assumed to have converged when each of ΔP2 , ΔP3, ΔP4 , ΔP5 , ΔQ2 , ΔQ3 and ΔQ4 is below a small pre-specified value. At this point the process is terminated.
Sometimes to accelerate computation in the P-Q buses the voltages obtained from (4.12) is multiplied by a constant. The voltage update of bus- i is then given by
 |
(4.21) |
where λ is a constant that is known as the acceleration factor . The value of λ has to be below 2.0 for the convergence to occur. Table 4.4 lists the values of the bus voltages after the 1st iteration and number of iterations required for the algorithm to converge for different values of λ. It can be seen that the algorithm converges in the least number of iterations when λ is 1.4 and the maximum number of iterations are required when λ is 2. In fact the algorithm will start to diverge if larger values of acceleration factor are chosen. The system data after the convergence of the algorithm will be discussed later.
Table 4.4 Gauss-Seidel method: bus voltages after 1 st iteration and number of iterations required for convergence for different values of l .
l |
Bus voltages (per unit) after 1st iteration |
No of iterations
for convergence |
V2 |
V3 |
V4 |
V5 |
1 |
0.9927Ð- 2.6° |
0.9883Ð- 2.83° |
0.9968Ð- 3.48° |
1.02 Ð- 0.89° |
28 |
2 |
0.9874Ð- 5.22° |
0.9766Ð- 8.04° |
0.9918Ð- 14.02° |
1.02Ð- 4.39° |
860 |
1.8 |
0.9883Ð- 4.7° |
0.9785Ð- 6.8° |
0.9903Ð- 11.12° |
1.02Ð- 3.52° |
54 |
1.6 |
0.9893Ð- 4.17° |
0.9807Ð- 5.67° |
0.9909Ð- 8.65° |
1.02Ð- 2.74° |
24 |
1.4 |
0.9903Ð- 3.64° |
0.9831Ð- 4.62° |
0.9926Ð- 6.57° |
1.02Ð- 2.05° |
14 |
1.2 |
0.9915Ð- 3.11° |
0.9857Ð- 3.68° |
0.9947Ð- 4.87° |
1.02Ð- 1.43° |
19 |
|