Power Flow
Power Flow
Dr Bakary Diarra
BIUST Palapye
Power flow analysis
I. Linear system resolution using Gauss elimination
1. Using Gauss-Seidel
2. Using Newton-Raphton
( ) ( ) ( ) … ( ) ( )
( )
( ) ( ) … ( ) ( )
← −
0 =
⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮
0 ( ) ( ) … ( ) ( ) ( )
← −
0 ( ) ( ) … ( ) ( )
⋮ ⋮ = ⋮
⋮ ⋮ ⋮ ⋮ ⋮ ( )
0 … ( ) ( ) ← −
0 0
6 +2 + =3
4 + 10 + 2 = 4
3 + 4 + 14 = 2
Alternate solutions?
= [ − + + ⋯+ + + ⋯+ ]
1
= − −
with = 1,2, … ,
In matrix form
+1 = + where = ( − )
+1 = + 1.9
given 0 = 0 =1
( + 1) = 3 −
= 10
with = 1,2, … ,
+1 = + 1.9
given 0 = 0 =1
( + 1) = 3 − +1
= ⋮ =
We can write
− =0
Supposing that = + ∆ as is approximated
Electrical Energy Systems Dr B. Diarra 16
Newton-Raphson method
The equation is rewritten as
− +∆ =0
Using Taylor series expansion
+∆ = +∆ +⋯
We deduce
( + 1) = ( ) + −
Extended Gauss-Seidel
In the non-linear case, D needs to be specified
= + −
( + 1) = ( ) + () −
is the Jacobian matrix equivalent to previous
Electrical Energy Systems Dr B. Diarra 18
Newton-Raphson method
The Jacobian matrix is
…
…
= = ⋮
⋮ ⋮ ⋮ … ⋮
− 1.25
−1
DC 1
2
1 DC
2
3 3
0.08
− 12.5
4
4
=∑ for ≠
= =
0
0 =−
+ = =V Y V
+ = =V Y V
Relationship of and
=V Y V ( − − )
=V Y V ( − − )
= = + +⋯+ + ⋯+
Using Gauss-Seidel
1
= − − →
1
( + 1) = ( )− ( + 1) − ( )
∗
The calculation is done with ( ) and just after with
∗
( + 1) in each iteration
The final voltage is then obtained by
1 −
( + 1) = ∗ − ( + 1) − ( )
( + 1)
( + 1) = ( ) ( ) ( ( )− ( )− )
Finally = +
= ( ) (0 − − )
0 = 0 =1
400 + 320
̅ =− = −4 − 3.2
100
300 + 270
̅ =− = −3 − 2.7
100
⋮ ⋮ ⋮
= = and = = ; = =
⋮ ⋮ ⋮
2 = = − −
3 = =− − −
4 = = − − , = 2,3, … ,
2 = = +∑ − −
3 = = ∑ − −
4 = =− +∑ − −
for , = 2,3, … ,