Transmission, Load Flow, Fault Analysis and Protection
2. Bus admittance matrix, Gauss-Seidel and Newton-Raphson load flow
The bus admittance matrix Ybus relates bus current injections to bus voltages, I = Ybus * V, and is built directly from the network's line and shunt admittances without needing to invert anything: the diagonal entry Yii is the sum of all admittances connected to bus i, and the off-diagonal Yij is the negative of the admittance directly between buses i and j (zero if no direct line exists) — a sparse matrix for any realistic network, since most bus pairs have no direct connection. A load flow (power flow) study finds the bus voltage magnitudes and angles that satisfy the specified power injections at every bus, classified as slack (voltage and angle specified, absorbs the system's unknown loss), PV (voltage magnitude and real power specified, typically a generator bus) or PQ (real and reactive power specified, typically a load bus).
| Property | Gauss-Seidel | Newton-Raphson |
|---|---|---|
| Convergence | linear, slow, more iterations | quadratic, fast, few iterations |
| Per-iteration effort | small, simple update | large — builds and inverts a Jacobian each time |
| Sensitivity to initial guess | less sensitive | more sensitive, but rarely a problem in practice |
| Scaling to large systems | poor — iteration count grows badly with bus count | the industry standard for large networks |
3. Symmetrical components and fault analysis
Fortescue's theorem says any unbalanced set of three phasors can be resolved into three balanced sets: positive sequence (same phase order as the original, balanced), negative sequence (reverse phase order, balanced) and zero sequence (all three in phase, equal magnitude). A symmetrical fault (three-phase, all phases equally involved) has only positive-sequence current, since the system stays balanced throughout — this is why a three-phase fault is analysed with an ordinary single-phase equivalent and no symmetrical-component machinery is even needed for it, despite the topic living in the same section.
| Fault type | Sequence network connection |
|---|---|
| Single line-to-ground (LG) | positive, negative and zero networks in SERIES |
| Line-to-line (LL) | positive and negative networks in PARALLEL (zero absent) |
| Double line-to-ground (LLG) | positive, negative and zero networks all in PARALLEL |
| Three-phase (LLL / LLLG) | positive network alone; negative and zero carry no current |
4. Protection, circuit breakers and stability
| Principle | How it decides to trip |
|---|---|
| Overcurrent | current magnitude exceeds a set threshold, often with an inverse time-current curve |
| Differential | compares current entering and leaving a protected zone; a mismatch means an internal fault |
| Directional | trips only for fault current flowing in one specified direction |
| Distance | measures apparent impedance (V/I) to the fault; trips if it falls within a set reach, in zones of increasing reach and time delay |
A circuit breaker's job is to interrupt fault current — the hard part is not opening its contacts but extinguishing the arc that forms between them as they separate, which different breaker types do by cooling and de-ionising the arc in oil, air-blast, SF6 gas or vacuum; its rating includes both a continuous current rating and a much larger breaking (interrupting) capacity for fault conditions. Differential protection is the classic example of UNIT protection — it defines an exact zone (a transformer, a generator, a short bus section) and is completely selective, tripping only for a fault genuinely inside that zone, whereas overcurrent and distance protection are NON-UNIT schemes that rely on graded time delays or reach settings between successive zones to achieve selectivity (coordination) instead.
Key takeaways
- The per-unit system exists to make transformer turns ratios vanish from a one-line diagram — chosen bases on the two sides in the transformer's own ratio, so its per-unit ratio becomes 1:1.
- Newton-Raphson converges quadratically (fast, few iterations, heavy per-iteration work) and is the industry standard for large networks; Gauss-Seidel converges linearly (slow but each iteration is cheap).
- A symmetrical (three-phase) fault has only positive-sequence current; an LG fault puts the three sequence networks in series, an LL fault puts positive and negative in parallel, and an LLG fault puts all three in parallel.
- Differential protection is unit (fully selective within its own zone); overcurrent and distance protection are non-unit, achieving selectivity through graded time or reach settings between zones instead.
- The equal area criterion says a system stays stable exactly when the available decelerating area is at least as large as the accelerating area already built up — faster fault clearing shrinks the accelerating area and directly buys stability margin.
Practice questions (10)
Attempt each one before opening the answer. Every explanation names the tempting wrong option as well as the right one, because that is where marks are lost.
A long transmission line is best modelled using
Show answer
Answer: C — distributed parameters, using the exact ABCD constants from the propagation constant
A long line's parameters cannot be lumped at a few points without introducing significant error, so it is modelled with distributed parameters and exact hyperbolic-function ABCD constants; short and medium lines use progressively simpler lumped approximations because the error from lumping is small enough over their length to ignore.Classical economic load dispatch, ignoring transmission loss, allocates load among generating units such that
Show answer
Answer: B — every operating unit's incremental fuel cost is equal
Equal incremental cost across all running units is the condition for minimum total cost: if one unit's next megawatt is cheaper, shifting load to it and away from a costlier unit reduces total cost, so the optimum is reached only when no such profitable shift remains — which is exactly when every unit's incremental cost is equal.A lightly loaded, long transmission line experiences the Ferranti effect, meaning its receiving-end voltage tends to
Show answer
Answer: B — rise above the sending-end voltage
A lightly loaded long line's distributed shunt capacitance dominates over its series inductance's loading effect, so the line behaves like a capacitive load on the source and its receiving end can end up at a HIGHER voltage than the sending end — exactly the condition a shunt reactor is installed to correct.In a bus admittance matrix Ybus, the off-diagonal element Yij (i not equal to j) for two buses with no direct line between them is
Show answer
Answer: B — zero
Yij is minus the admittance of the direct line between buses i and j; with no direct line, that admittance is zero, so Yij is zero — this is exactly why Ybus is sparse for any real network, where most bus pairs are not directly connected.Which of the following are correctly true when comparing Gauss-Seidel and Newton-Raphson load flow?
Show answer
Answer: A — Newton-Raphson converges quadratically, needing far fewer iterations for large systems.; B — Gauss-Seidel's per-iteration computational effort is generally lower than Newton-Raphson's.; C — Newton-Raphson requires building and solving a Jacobian matrix at every iteration.
The first three are all correctly true — quadratic convergence, cheaper Gauss-Seidel iterations, and Newton-Raphson's Jacobian requirement are the standard comparison. Newton-Raphson, not Gauss-Seidel, is the preferred method for large modern networks, since Gauss-Seidel's iteration count grows too quickly with system size to be practical there.A symmetrical three-phase fault at a bus produces fault current in
Show answer
Answer: A — the positive-sequence network only
A three-phase fault keeps the system exactly balanced throughout, so by definition there is no negative- or zero-sequence current at all — only the positive-sequence network carries current, which is exactly why this fault type can be solved with an ordinary per-phase equivalent circuit.For a line-to-line (LL) fault, the sequence networks are connected
Show answer
Answer: B — positive and negative in parallel, zero absent
A line-to-line fault does not involve ground, so no zero-sequence current path exists at all; positive- and negative-sequence networks are connected in parallel to represent the fault between the two faulted phases — this is the standard sequence-network diagram for an LL fault.Differential protection is best described as a
Show answer
Answer: B — unit scheme, fully selective within its own defined zone
Differential protection compares current at the two ends of a precisely defined zone (a transformer, a generator winding, a bus section) and trips only for a genuine internal fault — it needs no coordination with neighbouring zones and is the classic example of a unit protection scheme, unlike overcurrent (graded time) or distance (impedance-measuring, non-unit) protection.By the equal area criterion, a power system's transient stability after a fault clearance improves when
Show answer
Answer: B — the fault is cleared faster, shrinking the accelerating area
Faster clearing directly shrinks the accelerating area accumulated while the fault is on (rotor angle advances for less time before the fault is removed), leaving more of the available decelerating area as margin — the whole point of the equal area criterion's stability condition, decelerating area at least equal to accelerating area.A generator is rated 100 MVA at 11 kV. On a system base of 100 MVA, its per-unit reactance equals its reactance in ohms multiplied by 100 divided by (base kV)^2. If the base voltage is chosen equal to the generator's own rated voltage of 11 kV, a reactance quoted as 0.15 per unit on the generator's own rating is, on this same 100 MVA base, equal to ___ per unit
Numerical answer — type the value.
Show answer
Answer: 0.15
When the chosen base MVA and base kV exactly match the machine's own rated MVA and kV, the per-unit value does not change at all when converting bases — it stays 0.15 per unit, since the conversion factor (own rating / new base) is exactly 1 in both voltage and MVA.