NivaarExam PrepOfficial exam papers ↗

22-Elec-B7 Power Systems Engineering · May 2018

Question 2 of 7: The Reactive Capability Curve and a Round-Rotor Operating Table

Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)

Notes on this paper

Paper format. National Exams May 2018, 16-Elec-B7 Power Systems Engineering — open book, 3 hours, non-communicating calculator permitted. Seven problems of equal value; a complete paper (100 per cent) corresponds to 100 points scored out of a total possible 140 points, so every problem carries 20 points. Candidates are told to attempt all parts of all questions and to state any interpretive assumptions with the answer script.

Reference texts. M. E. El-Hawary, Electrical Power Systems: Design and Analysis, revised printing, IEEE Press / Wiley — the text whose notation (the cantilever transformer model, the sequence-network reduction, the two-curve equal-area construction) this paper follows throughout. J. D. Glover, T. J. Overbye and M. S. Sarma, Power System Analysis and Design, 6th ed. — line parameters, bundling, power flow and symmetrical components. J. J. Grainger and W. D. Stevenson, Power System Analysis — fault analysis and network reduction. P. Kundur, Power System Stability and Control — the generator capability diagram and transient stability. S. J. Chapman, Electric Machinery Fundamentals, 5th ed. — transformer and synchronous-machine behaviour. For Canadian protection and insulation-coordination practice, the IEEE C37 relay standards and CSA C22.3 overhead-systems series apply.

Question 2: The Reactive Capability Curve and a Round-Rotor Operating Table (20 points)

Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.

Part (a) — the reactive capability curve and what sets each of its boundaries.

-1.2-0.8-0.40.00.40.81.20.000.250.500.751.001.25reactive power Q at the machine terminals (pu)active power P (pu)under-excited (absorbing Q)over-excited (supplying Q)operatingregionarmature (stator) heatingfield (rotor) heatingprime-mover MW limitend-region heating / stability
Figure 2.1 — reactive capability (P-Q loading) diagram of a cylindrical-rotor synchronous machine at rated terminal voltage. The permitted operating region is the area enclosed by all four limits.

The capability curve is drawn in the plane of active power against reactive power at the machine terminals, with lagging (supplied) reactive power conventionally to the right. Every boundary on it is a thermal or a stability limit, and each has a different physical origin.

The armature-current limit is the arc of a circle centred on the origin. With the terminal voltage held at its rated value the apparent power is $S=V_tI_a$, so a fixed maximum stator current is a fixed maximum $\sqrt{P^{2}+Q^{2}}$. What it protects is stator copper loss and therefore the temperature of the armature winding insulation. It is the binding constraint over the middle of the diagram, near unity power factor.

The field-current limit is the arc of a second circle, centred on the negative reactive axis at $-V_t^{2}/X_s$ and of radius $E_{max}V_t/X_s$. It comes from the same two power expressions used in part (b): eliminating $\delta$ between $P=(E V_t/X_s)\sin\delta$ and $Q=(EV_t\cos\delta-V_t^{2})/X_s$ gives $$P^{2}+\left(Q+\frac{V_t^{2}}{X_s}\right)^{2}=\left(\frac{EV_t}{X_s}\right)^{2}$$ so a ceiling on the excitation is a ceiling on $E$ and hence a circle. It protects the rotor winding from over-heating and binds in the over-excited region, where the machine is exporting reactive power.

In the under-excited region two further limits appear. The end-region heating limit arises because a weakly excited machine has a low rotor field, so the armature end-winding leakage flux is no longer opposed and penetrates the stator core end laminations in the axial direction, inducing eddy currents in the end packets and the clamping structure and producing local hot spots. Close beside it lies the steady-state stability limit: as $\delta$ approaches ninety degrees the synchronising torque coefficient $dP/d\delta$ falls to zero and the machine can no longer hold synchronism against a disturbance, so a practical margin of the order of ten per cent in power is left. Finally the horizontal prime-mover limit caps the active power at whatever the turbine can deliver, which is a mechanical rather than an electrical constraint, and a lower boundary at minimum stable turbine load often applies as well. Reducing the terminal voltage shrinks the armature-limit circle in proportion to $V_t$ and moves the field-limit circle inward and to the left, so the reactive capability of a generator is always poorest exactly when the system most needs it, during a low-voltage event.

Given. A cylindrical-rotor machine on an infinite bus, with the three partly specified operating conditions of the printed table.

Given data and the quantities to be found
Condition$P$ (pu)$Q_2$ (pu)$E$ (pu)$\delta$
A?$-0.1$1.25?
B1.350.0??
C??1.3$40^{\circ}$
Common data$V=1.00$ pu$X_s=0.375$ pu$R_a=0$saturation neglected

Find. The two missing entries in each row of the table, so that all three operating points are fully specified.

Approach. Both unknowns in every row follow from the two round-rotor power transfer equations written at the infinite bus; each condition supplies two of the four quantities $P$, $Q_2$, $E$, $\delta$, and the pair of equations then closes algebraically with no iteration.

  1. Part (b) — write the two governing equations at the bus. With armature resistance neglected the machine sees the infinite bus through the synchronous reactance alone, so the complex power delivered to the bus is $S_2=V\,I_a^{*}$ with $I_a=(E\,/\,\delta-V)/(jX_s)$. Separating real and imaginary parts, $$P=\frac{EV}{X_s}\sin\delta \qquad Q_2=\frac{EV\cos\delta-V^{2}}{X_s}$$ With $V=1.00$ and $X_s=0.375$ these become $P=2.66667\,E\sin\delta$ and $Q_2=(E\cos\delta-1)/0.375$. It is worth restating the two of them in the form used below: $E\sin\delta=0.375P$ and $E\cos\delta=1+0.375Q_2$, which makes the pair a simple rectangular decomposition of the phasor $E$.
  2. Condition A: the reactive equation fixes the angle first. Here $E=1.25$ and $Q_2=-0.1$ are known, so $$E\cos\delta=1+0.375(-0.1)=0.9625 \quad\Longrightarrow\quad \cos\delta=\frac{0.9625}{1.25}=0.7700$$ which gives $$\boxed{\delta_A=39.6461^{\circ}}$$ Only the positive root is admissible because the machine is generating.
  3. Condition A: the active power then follows directly. With $\sin(39.6461^{\circ})=0.638044$, $$P_A=\frac{(1.25)(1.00)}{0.375}\sin\delta_A=3.33333\times0.638044=\boxed{P_A=2.1268\ \text{pu}}$$ The pull-out power at this excitation is $EV/X_s=3.3333$ pu, so the machine is loaded to sixty-four per cent of its steady-state limit — comfortable, but the slightly negative $Q_2$ shows it is drawing a little reactive power from the system, that is, running under-excited.
  4. Condition B: zero reactive power collapses the second equation. Setting $Q_2=0$ makes $$E\cos\delta=V=1.000$$ which says that the tip of the phasor $E$ lies on the vertical line through $V$ in the chart below. The active power then supplies the other component, $$E\sin\delta=0.375P=0.375(1.35)=0.50625$$
  5. Condition B: recover the excitation and the angle from those two components. Dividing and then combining, $$\tan\delta_B=\frac{0.50625}{1.000} \quad\Longrightarrow\quad \boxed{\delta_B=26.8508^{\circ}}$$ $$E_B=\sqrt{(1.000)^{2}+(0.50625)^{2}}=\sqrt{1.256289}=\boxed{E_B=1.12084\ \text{pu}}$$ Delivering 1.35 pu at exactly unity power factor therefore needs only about twelve per cent more excitation than the bus voltage itself.
  6. Condition C: both equations are explicit. With $E=1.3$ and $\delta=40^{\circ}$ nothing has to be inverted: $$P_C=\frac{(1.3)(1.00)}{0.375}\sin40^{\circ}=3.46667\times0.642788=\boxed{P_C=2.2283\ \text{pu}}$$ $$Q_{2C}=\frac{(1.3)\cos40^{\circ}-1}{0.375}=\frac{0.995858-1.000000}{0.375}=\boxed{Q_{2C}=-0.01105\ \text{pu}}$$ The reactive output is very slightly negative, and the reason is visible in the numerator: $E\cos\delta=0.9959$ falls just short of $V=1.000$, so the machine sits a hair on the under-excited side of the unity-reactive-power line. This is not a rounding artefact — the sign is the answer, and it means the machine is absorbing about eleven kVAr per MVA of rating.
  7. Read the three points off the operating chart. Plotting $E\sin\delta$ against $E\cos\delta$ turns the two governing equations into a pair of rulings: horizontal lines are constant $P$ and vertical lines are constant $Q_2$, with the line $E\cos\delta=V$ separating over- from under-excited operation. Conditions A and C both sit marginally to the left of that line and condition B sits exactly on it, which is the geometric statement of the three answers.
  8. Closure check on the armature current. Recomputing $I_a=(E\,/\,\delta-V)/(jX_s)$ from each completed row and forming $S_2=VI_a^{*}$ reproduces the tabulated $P$ and $Q_2$ to six decimal places, with $|I_a|$ equal to 2.1292, 1.3500 and 2.2284 pu for A, B and C respectively. Condition B has the smallest current for a substantial power output, which is the expected reward for operating at unity power factor.
0.000.250.500.751.001.250.000.250.500.751.00E cos δ (pu)E sin δ (pu)E cos δ = V(Q = 0)ABChorizontal lines are constant P; vertical lines are constant Q
Figure 2.2 — the three operating points plotted as the tip of the excitation phasor. The vertical dashed line is the locus of zero reactive power at the machine terminals.
Final results — completed operating table
Condition$P$ (pu)$Q_2$ (pu)$E$ (pu)$\delta$$|I_a|$ (pu)
A2.1268$-0.1$ (given)1.25 (given)$39.6461^{\circ}$2.1292
B1.35 (given)0.0 (given)1.12084$26.8508^{\circ}$1.3500
C2.2283$-0.01105$1.3 (given)$40^{\circ}$ (given)2.2284
Check: the question says to complete the table “neglecting armature reaction”, which cannot be taken literally — the synchronous reactance $X_s=X_l+X_{ar}$ already contains the armature-reaction reactance, and removing it would leave no reactance for power transfer at all. The instruction is read here in its usual textbook sense: neglect magnetic saturation (so $E$ is proportional to field current and $X_s$ is constant) and neglect armature resistance. Every number above follows from that reading.