Skip to main content

G-S load flow, j0.12 lines, losses

Short answer

Using G-S with Ybus elements j8.3333 pu (off-diagonal) and −j16.6667 (diagonal), two iterations gave V2 = 0.9576∠−8.87° pu and V3 = 0.9443∠−8.78° pu. Slack bus power: P1 = 2.432 pu, Q1 = 1.005 pu. Since line resistance is zero, real power loss = 0 pu; reactive power loss = 0.418 pu (not fully converged after two iterations).

For a network given below, series impedance of each line is 0 + j0.12 p.u. Shunt admittance of the line is negligible. Calculate: (i) Voltage at bus-2 and bus-3 by using G-S method (upto 2 iteration). (ii) Slack bus real and reactive power. (iii) Network real and reactive power losses. [Figure: three buses with lines 1-2, 1-3 and 2-3; bus 1 is the slack bus with generator, V1 = 1.0∠0°; load at bus 2 = 1.5 + j0.4 p.u; load at bus 3 = 1.2 + j0.5 p.u]

Data and Ybus

Each line: z=j0.12z = j0.12, y=1/j0.12=−j8.3333y = 1/j0.12 = -j8.3333 pu. Every bus has two lines:

Ybus=[−j16.6667j8.3333j8.3333j8.3333−j16.6667j8.3333j8.3333j8.3333−j16.6667]Y_{bus} = \begin{bmatrix} -j16.6667 & j8.3333 & j8.3333 \\ j8.3333 & -j16.6667 & j8.3333 \\ j8.3333 & j8.3333 & -j16.6667 \end{bmatrix}

Bus 1: slack, V1=1.0∠0∘V_1 = 1.0\angle 0^\circ. Buses 2 and 3 are PQ buses with net injections S2=−1.5−j0.4S_2 = -1.5 - j0.4 and S3=−1.2−j0.5S_3 = -1.2 - j0.5 pu. Flat start: V2(0)=V3(0)=1.0∠0∘V_2^{(0)} = V_3^{(0)} = 1.0\angle 0^\circ.

G-S equation (latest values used immediately):

Vi(r+1)=1Yii[Pi−jQiVi∗(r)−∑k≠iYikVk]V_i^{(r+1)} = \frac{1}{Y_{ii}}\left[\frac{P_i - jQ_i}{V_i^{*(r)}} - \sum_{k \ne i}Y_{ik}V_k\right]

(i) Iteration 1

V2(1)=1−j16.6667[−1.5+j0.41.0−j8.3333(1.0)−j8.3333(1.0)]=−1.5−j16.2667−j16.6667=0.9760−j0.0900=0.9801∠−5.2685∘V3(1)=1−j16.6667[−1.2+j0.51.0−j8.3333(1.0)−j8.3333(0.9760−j0.0900)]=−1.95−j15.9667−j16.6667=0.9580−j0.1170=0.9651∠−6.9630∘\begin{aligned} V_2^{(1)} &= \frac{1}{-j16.6667}\left[\frac{-1.5 + j0.4}{1.0} - j8.3333(1.0) - j8.3333(1.0)\right] \\ &= \frac{-1.5 - j16.2667}{-j16.6667} = 0.9760 - j0.0900 = 0.9801\angle -5.2685^\circ \\ V_3^{(1)} &= \frac{1}{-j16.6667}\left[\frac{-1.2 + j0.5}{1.0} - j8.3333(1.0) - j8.3333(0.9760 - j0.0900)\right] \\ &= \frac{-1.95 - j15.9667}{-j16.6667} = 0.9580 - j0.1170 = 0.9651\angle -6.9630^\circ \end{aligned}

Iteration 2

−1.5+j0.4V2∗(1)=−1.4865+j0.5469V2(2)=(−1.4865+j0.5469)−(0.9750+j16.3167)−j16.6667=0.9462−j0.1477=0.9576∠−8.8715∘−1.2+j0.5V3∗(1)=−1.1714+j0.6650V3(2)=(−1.1714+j0.6650)−(1.2307+j16.2182)−j16.6667=0.9332−j0.1441=0.9443∠−8.7797∘\begin{aligned} \frac{-1.5 + j0.4}{V_2^{*(1)}} &= -1.4865 + j0.5469 \\ V_2^{(2)} &= \frac{(-1.4865 + j0.5469) - (0.9750 + j16.3167)}{-j16.6667} \\ &= 0.9462 - j0.1477 = 0.9576\angle -8.8715^\circ \\ \frac{-1.2 + j0.5}{V_3^{*(1)}} &= -1.1714 + j0.6650 \\ V_3^{(2)} &= \frac{(-1.1714 + j0.6650) - (1.2307 + j16.2182)}{-j16.6667} \\ &= 0.9332 - j0.1441 = 0.9443\angle -8.7797^\circ \end{aligned}

Here ∑k≠2Y2kVk=j8.3333(V1+V3(1))=0.9750+j16.3167\sum_{k \ne 2}Y_{2k}V_k = j8.3333(V_1 + V_3^{(1)}) = 0.9750 + j16.3167 and ∑k≠3Y3kVk=j8.3333(V1+V2(2))=1.2307+j16.2182\sum_{k \ne 3}Y_{3k}V_k = j8.3333(V_1 + V_2^{(2)}) = 1.2307 + j16.2182.

IterationV2V_2 (pu)V3V_3 (pu)
01.0∠0∘1.0\angle 0^\circ1.0∠0∘1.0\angle 0^\circ
10.9801∠−5.27∘0.9801\angle -5.27^\circ0.9651∠−6.96∘0.9651\angle -6.96^\circ
20.9576∠−8.87∘0.9576\angle -8.87^\circ0.9443∠−8.78∘0.9443\angle -8.78^\circ

(ii) Slack bus power (with iteration-2 voltages)

I1=Y11V1+Y12V2+Y13V3=2.4318−j1.0052S1=V1I1∗=1.0(2.4318+j1.0052)\begin{aligned} I_1 &= Y_{11}V_1 + Y_{12}V_2 + Y_{13}V_3 = 2.4318 - j1.0052 \\ S_1 &= V_1I_1^* = 1.0(2.4318 + j1.0052) \end{aligned}

So P1=2.4318P_1 = 2.4318 pu and Q1=1.0052Q_1 = 1.0052 pu.

(iii) Network losses

Line flows Sik=Vi[yik(Vi−Vk)]∗S_{ik} = V_i[y_{ik}(V_i - V_k)]^*:

LineSikS_{ik} (pu)SkiS_{ki} (pu)Loss Sik+SkiS_{ik} + S_{ki}
1-21.2307+j0.44851.2307 + j0.4485−1.2307−j0.2426-1.2307 - j0.2426j0.2059j0.2059
1-31.2011+j0.55671.2011 + j0.5567−1.2011−j0.3464-1.2011 - j0.3464j0.2103j0.2103
2-3−0.0121+j0.1068-0.0121 + j0.10680.0121−j0.10530.0121 - j0.1053j0.0015j0.0015

Answer:

  • V2=0.9576∠−8.87∘V_2 = 0.9576\angle -8.87^\circ pu, V3=0.9443∠−8.78∘V_3 = 0.9443\angle -8.78^\circ pu after 2 iterations.
  • Slack bus: P1=2.432P_1 = 2.432 pu, Q1=1.005Q_1 = 1.005 pu.
  • Losses: real power loss =0= 0 (lines have no resistance); reactive power loss =0.418= 0.418 pu.

Note: after only two iterations the solution has not fully converged, so S1S_1 minus the loads (2.432−2.7=−0.272.432 - 2.7 = -0.27) does not yet equal the losses; with more iterations P1P_1 tends to 2.7 pu.

✓ Answer checked by Claude Sonnet 5.5

All Power System Analysis II chapter 2 questions →