Subjects power systems

Newton Raphson Loadflow

Step-by-step solutions with LaTeX - clean, fast, and student-friendly.

Search Solutions

Newton Raphson Loadflow


1. **Problem Statement:** We have a 3-bus power system with bus 1 as slack bus and bus 2 as PQ bus with given power demands $P_2 = -0.8$ pu and $Q_2 = -0.3$ pu. The goal is to perform one iteration of the Newton-Raphson load flow method starting from flat start ($|V_2|=1.0$ pu, $\delta_2=0$) to find updated voltage magnitude and angle at bus 2. 2. **Y-bus Matrix:** Given line impedance $j0.1$ pu between bus 1 and 2, admittance is $$Y_{12} = \frac{1}{j0.1} = -j10 = j10 \quad \text{(since admittance is imaginary)}$$ Thus, $$Y = \begin{bmatrix} -j10 & j10 \\ j10 & -j10 \end{bmatrix}$$ 3. **Power Flow Equations:** Power injection at bus 2 is $$S_2 = P_2 + jQ_2 = V_2 I_2^* = V_2 \left( \sum_{k=1}^2 Y_{2k} V_k \right)^*$$ Current at bus 2: $$I_2 = Y_{21} V_1 + Y_{22} V_2 = j10 \times 1.05 - j10 V_2$$ 4. **Mismatch Equations:** Define mismatches as $$\Delta P = P_2^{spec} - P_2^{calc}$$ $$\Delta Q = Q_2^{spec} - Q_2^{calc}$$ Using power flow formulas, $$P_2 = |V_2||V_1| B_{21} \sin(\delta_2 - \delta_1) = 10.5 |V_2| \sin \delta_2$$ $$Q_2 = |V_2||V_1| B_{21} \cos(\delta_2 - \delta_1) - |V_2|^2 B_{22} = 10.5 |V_2| \cos \delta_2 - 10.5 |V_2|^2$$ At flat start $|V_2|=1$, $\delta_2=0$, calculate: $$P_2 = 10.5 \times 1 \times \sin 0 = 0$$ $$Q_2 = 10.5 \times 1 \times \cos 0 - 10.5 \times 1^2 = 10.5 - 10.5 = 0$$ Mismatch: $$\Delta P = -0.8 - 0 = -0.8$$ $$\Delta Q = -0.3 - 0 = -0.3$$ 5. **Jacobian Matrix:** Partial derivatives at flat start: $$\frac{\partial P_2}{\partial \delta_2} = 10.5 |V_2| \cos \delta_2 = 10.5$$ $$\frac{\partial P_2}{\partial |V_2|} = 10.5 \sin \delta_2 = 0$$ $$\frac{\partial Q_2}{\partial \delta_2} = 10.5 |V_2| \sin \delta_2 = 0$$ $$\frac{\partial Q_2}{\partial |V_2|} = -10.5 \cos \delta_2 - 2 \times 10.5 |V_2| = -10.5 - 21 = -31.5$$ Jacobian: $$J = \begin{bmatrix} 10.5 & 0 \\ 0 & -31.5 \end{bmatrix}$$ 6. **Calculate Voltage Corrections:** $$\begin{bmatrix} \Delta \delta_2 \\ \Delta |V_2| \end{bmatrix} = J^{-1} \begin{bmatrix} \Delta P \\ \Delta Q \end{bmatrix} = \begin{bmatrix} \frac{1}{10.5} & 0 \\ 0 & \frac{1}{-31.5} \end{bmatrix} \begin{bmatrix} -0.8 \\ -0.3 \end{bmatrix} = \begin{bmatrix} -0.0762 \\ 0.0095 \end{bmatrix}$$ 7. **Update Voltage Magnitude and Angle:** $$|V_2|^{new} = 1.0 + 0.0095 = 1.0095$$ $$\delta_2^{new} = 0 - 0.0762 = -0.0762 \text{ radians}$$ **Interpretation:** The voltage magnitude at bus 2 slightly increases from 1.0 to approximately 1.0095 pu, and the voltage angle shifts slightly negative, indicating a small phase lag relative to the slack bus. This is consistent with the load demands and system admittance.