Skip to main content

Initial mismatch and Jacobian, 2-bus system

Short answer

Initial mismatches: ΔP2 = -1.0 pu, ΔQ2 = -0.5 pu. Initial Jacobian: ∂P2/∂δ2 = 10, ∂P2/∂|V2| = 0, ∂Q2/∂δ2 = 0, ∂Q2/∂|V2| = 10.

For the system given below, find the initial power mismatch matrix and initial Jacobian matrix. Take a flat start of V2(0) = 1.0 pu and δ2(0) = 0°. [Figure: bus 1 generator with 1.0∠0°, line j0.1 to bus 2, load 1 + j0.5 at bus 2]

Bus 1 is the slack bus (1.0∠0∘1.0\angle 0^\circ) and bus 2 is a load (PQ) bus with P2sp=−1.0P_2^{sp} = -1.0 pu and Q2sp=−0.5Q_2^{sp} = -0.5 pu. Unknowns: δ2\delta_2 and ∣V2∣|V_2|. Flat start: ∣V2∣0=1.0|V_2|^{0} = 1.0, δ20=0\delta_2^{0} = 0.

Y-bus

y12=1/j0.1=−j10y_{12} = 1/j0.1 = -j10 pu, so

Ybus=[−j10.0000j10.0000j10.0000−j10.0000]=[10∠−90∘10∠90∘10∠90∘10∠−90∘]Y_{bus} = \begin{bmatrix} -j10.0000 & j10.0000 \\ j10.0000 & -j10.0000 \end{bmatrix} = \begin{bmatrix} 10\angle -90^\circ & 10\angle 90^\circ \\ 10\angle 90^\circ & 10\angle -90^\circ \end{bmatrix}

Power equations at bus 2

P2=∣V2∣∣V1∣∣Y21∣cos⁡(θ21−δ2+δ1)+∣V2∣2∣Y22∣cos⁡θ22=10∣V2∣cos⁡(90∘−δ2)=10∣V2∣sin⁡δ2Q2=−∣V2∣∣V1∣∣Y21∣sin⁡(θ21−δ2+δ1)−∣V2∣2∣Y22∣sin⁡θ22=−10∣V2∣cos⁡δ2+10∣V2∣2\begin{aligned} P_2 &= |V_2||V_1||Y_{21}|\cos(\theta_{21} - \delta_2 + \delta_1) + |V_2|^2|Y_{22}|\cos\theta_{22} \\ &= 10|V_2|\cos(90^\circ - \delta_2) = 10|V_2|\sin\delta_2 \\ Q_2 &= -|V_2||V_1||Y_{21}|\sin(\theta_{21} - \delta_2 + \delta_1) - |V_2|^2|Y_{22}|\sin\theta_{22} \\ &= -10|V_2|\cos\delta_2 + 10|V_2|^2 \end{aligned}

Initial power mismatch

At ∣V2∣=1|V_2| = 1, δ2=0\delta_2 = 0: P20=10(1)sin⁡0=0.0000P_2^{0} = 10(1)\sin 0 = 0.0000 and Q20=−10+10=0.0000Q_2^{0} = -10 + 10 = 0.0000.

ΔP2=P2sp−P20=−1.0−0=−1.0000ΔQ2=Q2sp−Q20=−0.5−0=−0.5000\begin{aligned} \Delta P_2 &= P_2^{sp} - P_2^{0} = -1.0 - 0 = -1.0000 \\ \Delta Q_2 &= Q_2^{sp} - Q_2^{0} = -0.5 - 0 = -0.5000 \end{aligned} [ΔM0]=[ΔP2ΔQ2]=[−1.0000−0.5000][\Delta M^{0}] = \begin{bmatrix} \Delta P_2 \\ \Delta Q_2 \end{bmatrix} = \begin{bmatrix} -1.0000 \\ -0.5000 \end{bmatrix}

Initial Jacobian

∂P2∂δ2=10∣V2∣cos⁡δ2=10∂P2∂∣V2∣=10sin⁡δ2=0∂Q2∂δ2=10∣V2∣sin⁡δ2=0∂Q2∂∣V2∣=−10cos⁡δ2+20∣V2∣=−10+20=10\begin{aligned} \frac{\partial P_2}{\partial \delta_2} &= 10|V_2|\cos\delta_2 = 10 \\ \frac{\partial P_2}{\partial |V_2|} &= 10\sin\delta_2 = 0 \\ \frac{\partial Q_2}{\partial \delta_2} &= 10|V_2|\sin\delta_2 = 0 \\ \frac{\partial Q_2}{\partial |V_2|} &= -10\cos\delta_2 + 20|V_2| = -10 + 20 = 10 \end{aligned} [J0]=[10.00000.00000.000010.0000][J^{0}] = \begin{bmatrix} 10.0000 & 0.0000 \\ 0.0000 & 10.0000 \end{bmatrix}

(For use in the first iteration: Δδ2=−1/10=−0.1000\Delta\delta_2 = -1/10 = -0.1000 rad =−5.730∘= -5.730^\circ and Δ∣V2∣=−0.5/10=−0.0500\Delta|V_2| = -0.5/10 = -0.0500, giving V21=0.9500∠−5.730∘V_2^{1} = 0.9500\angle -5.730^\circ pu.)

Answer: [ΔM0]=[−1.0000, −0.5000]T[\Delta M^{0}] = [-1.0000,\ -0.5000]^T and [J0]=[100010][J^{0}] = \begin{bmatrix} 10 & 0 \\ 0 & 10 \end{bmatrix}.

✓ Answer checked by Claude Sonnet 5.5

All Power System Analysis II chapter 2 questions →