0% found this document useful (0 votes)
11 views5 pages

Fast Decouple Method

Fast decoupled

Uploaded by

Hema
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
11 views5 pages

Fast Decouple Method

Fast decoupled

Uploaded by

Hema
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

Exp - 09

Load Flow Solution Using Fast De-Coupled

Aim:

To determine the load flow solution of the given problem with Fast De-coupled method using
MATLAB coding.

Apparatus Required:

SI.No Apparatus Specification


Dual core, RAM 512 MB 1.2 GHz speed,
1 PC
80 GB
2 MATLAB 7.5

Theory:

Due to the weak coupling between PV and q-δ half of the elements of jacobian matrix are
neglected further the assumptions made are

Cos δij = 1

Sin δij = 0

Qi << Bij |V|2

The simplified FDLF equations are

(∆P / |V|) = [B’][∆δ]

(∆Q / |V|) = [B’’][∆|V|]

One iteration implies one solution for [∆δ] to update [δ] and one solution for [∆P/|V|] to
update [|V|] and is termed as 1-δ and 1-V iteration. The convergence for the real and reactive power is
achieved when max[∆P] < ∑P; max[∆Q]<=∑Q.

The main advantage of the decoupled load flow as compared to Newton Raphson method is its
reduced memory is storing Jacobian.
Algorithm:

Step1:Formulate Y bus matrix, then compute the bus susceptance matrices B’&B’’

Step2: Assume flat start for status V solution

 10  0 for i=1,2,…..N (all buses except slack buses)

Vi  1 0 for i = M+1…….N (for all PQ buses)

Vi  Vi spec
for all PV buses & slack bus

Step3: For load buses calculate Pi cal & Qi cal using

N
Pi cal   Vi Yij V j cos(Qij   j   i )
j 1

N
Qi cal   Vi Yij V j sin (Qij   j   i )
i 1

Step4: For PV buses check for Q limit violation if Qimin < Q i< Qi max calculate Pi cal. Qi cal > Qi max,,
Qi spec=Qi max then PV bus will act as PQ bus

Step5: Calculate mismatch vector using

ΔPi=Pi(spec)-Pi cal

ΔQi=Qi(spec)-Qi cal

Step6: compute

Pi (max)  max Pi , i  1,2,.......N & slack bus

Qi (max)  max Qi , M  1.......N

Step7: calculate

Pi
 i    B1
Vi

Vi    B " 1 Qi


Vi

Step8: This procedure is continued until

Pi    Qi   otherwise go to step 3)


Program:

clear all;
clc;
n=input('enter the number of buses:');
for i=1:n
for j=1:n
y(i,j)=input('Enter the admittance value:');
end
end
yb(n,n)=0;
for i=1:n
for j=1:n
if i==j
for k=1:n
yb(i,j)=yb(i,j)+y(i,k);
end
else
yb(i,j)=-y(i,j);
end
end
end
th(n)=0;
for i=1:n
Bus=i
mag(i)=input('Enter the voltage value:');
acp(i)=input('Enter real power value:');
acq(i)=input('Enter reactive power');
end
j=1;
for i=1:n
if(mag(i)==1)
k(j)=i;
j=j+1;
end
end
my=abs(yb);an=angle(yb);g=real(yb);b=imag(yb);
yb
mag
th
acp
acq
Chmag(n)=0;Chth(n)=0;Pp(n)=0;Qq(n)=0;
for i=1:n
for j=1:n
Pp(i)=Pp(i)+mag(i)*my(i,j)*mag(j)*cos(an(i,j)-th(i)+th(j));
Qq(i)=Qq(i)-mag(i)*my(i,j)*mag(j)*sin(an(i,j)-th(i)+th(j));
end
Pp
Qq
end
for i=2:n
for j=2:n
if i~=j
j1(i,j)=mag(i)*mag(j)*(g(i,j)*sin(th(i)-th(j))-b(i,j)*cos(th(i)-th(j)));
j3(i,j)=0;
j2(i,j)=0;
j4(i,j)=j1(i,j);
else
j1(i,j)=-Qq(i)-b(i,j)*(mag(i)^2);
j2(i,j)=0;
j3(i,j)=0;
j4(i,j)=Qq(i)-b(i,j)*(mag(i)^2);
end
end
end
ja1(1:n-1,1:n-1)=j1(2:n,2:n);
ja2(1:n-1,1:n-1)=j2(2:n,2:n);
ja3(1:n-1,1:n-1)=j3(2:n,2:n);
ja4(1:n-1,1:n-1)=j4(2:n,2:n);
jacob=[ja1 ja2;ja3 ja4]
bb(1:n-1, 1:n-1)=b(2:n,2:n)
for i=2:n
j=1;
if(i==k(j))
bbb(j,:)=b(i,i);
j=j+1;
end
end
bbb
delp(2:n)=acp(2:n)-Pp(2:n)
delq(2:n)=acq(2:n)-Qq(2:n)
for i=2:n
x(i)=delp(i)/abs(mag(i));
y(i)=delq(i)/abs(mag(i));
end
Chth(2:n)=-inv(bb)*x(2:n)';
for i=2:n
j=1;
if(i==k(j))
Chmag(i)=-inv(bbb)*y(i)';
end
end
mag=mag+Chmag;
th=th+Chth;
fprintf('the voltage value for buses');
mag
fprintf('the angle values for buses');
th
Result:

You might also like