Scribd
Upload
18 views4 pages

Expt2 PSA

The document describes an experiment to perform power flow analysis of a power system network using the Gauss-Siedel method. It provides the theory behind the method and shows the system model and equations. It then outlines the procedure, including writing MATLAB code to implement the Gauss-Siedel algorithm. The code is run to iteratively solve for the voltage magnitudes and angles until the power mismatch converges to below a threshold, successfully demonstrating power flow using the method.

Uploaded by

Vinit Singh
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
18 views4 pages

Expt2 PSA

The document describes an experiment to perform power flow analysis of a power system network using the Gauss-Siedel method. It provides the theory behind the method and shows the system model and equations. It then outlines the procedure, including writing MATLAB code to implement the Gauss-Siedel algorithm. The code is run to iteratively solve for the voltage magnitudes and angles until the power mismatch converges to below a threshold, successfully demonstrating power flow using the method.

Uploaded by

Vinit Singh
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/ 4

EXPERIMENT-II

Aim: To carry out power-flow / load-flow of a given power system network


using Gauss-Siedel method
Theory:

In the above network, all the transmission lines are identical. For each line, the
shunt element is a capacitor with an admittance of Yshunt =-j 0.01, while each
series element is na inductor with na impedance of Zseries =-j0.1.
Consider Case-I of the power flow problem: given V1,S2,S3……….. Sn , find
V2,S1,V3……….. Vn. We use a form of 10.3 repeated here for convenience.

(10.3a)

(i=1,2,3,….n) (10.3b)
Note: If we know V1,V2,V3……….. Vn, we can solve for S1, explicitly using
(10.3a). Since we already do know V1, it only remains to find V2,V3……….. Vn.
These n-1 unknowns may be found from n-1 equations of (10.3b). Thus , the
heart of the problem is the solution of n-1implicit equations in the unknown
V2,V3……….. Vn, where V1 and S2,S3……….. Sn are known. Equivalently, taking
complex conjugates, we can replace (10.3) by

(i=1,2,3,….n) (10.6)

We now arrange (10.6) in a form in which a solution by iteration may be


attempted. It should be noted that there are alternative ways of setting up the
problem. Dividing (10.6) by Vi* and separating out the Yk term, we can rewrite
(10.6) as:

(i=1,2,3,….n)
Or, rearranging,

(i=1,2,3,….n)

Procedure:
• First compute the bus admittance matrix. Check that:

• Type the code given below for the Gauss-Siedel algorithm in a


MATLAB mfile.
• Debug the code and run it to get the final results
Code:
clc
%Data%
s1b=1;pvb=1;pqb=1;nb=s1b+pvb+pqb;
Y=[-j*19.98 j*10 j*10;j*10 -j*19.98 j*10;j*10 j*10 -j*19.98];
y=abs(Y); phi=angle(Y);
v(1)=1.0; th(1)=0; v(2)=1.05;
Psp(2)=0.6661; Psp(3)=-2.8653; Qsp(3)=-1.2244;
%initial guess for unknown theta and v %
th(2)=0;v(3)=1.0;th(3)=0;
m=180/pi;
%Program for Gauss Siedel algorithms%
err=1.0;t=0; %'t' is the number of iterations of algorithm
while err>1e-4
for i=2:3
Px(i)=0; Qx(i)=0;
for k=1:3
Px(i)=Px(i)+v(i)*v(k)*y(i,k)*cos(th(i)-th(k)-phi(i,k));
Qx(i)=Qx(i)+v(i)*v(k)*y(i,k)*sin(th(i)-th(k)-phi(i,k));
end
end
for i=1:3
V(i)=v(i)*(cos(th(i))+j*sin(th(i)));
end
for i=2:3
vint(i)=0;
for k=1:3
if k~=i
vint(i)=vint(i)+Y(i,k)*V(k) ;
else
end
end
end
for i=2:nb
if i<(s1b+pvb)
V(i)=(((Psp(i)-j*Qx(i))/conj(V(i))-vint(i))/Y (i,i));
else
V(i)=(((Psp(i)-j*Qsp(i))/conj(V(i))-vint(i))/Y(i,i));
end
end
for i=2:nb
if i<=(s1b+pvb)
th(i)=angle(V(i));
else
v(i)=abs(V(i));
th(i)=angle(V(i));
end
delp=Psp-Px;
delq=Qsp-Qx;
delpqx=[delp'; delq'];
err=abs(max(delpqx));
t=t+1;
end
end
[v' m*th']

RESULTS:

CONCLUSION:
Through this experiment we have successfully learnt how to carry out power-
flow / load-flow of a given power system network using Gauss-Siedel method
using Matlab/Simulink by writing desired code.

You might also like