5.3 Synchronous Generators: Operation, Excitation, Power-Angle Equations & V-Curves

5.3 Synchronous Generators: Operation, Excitation, Power-Angle Equations & V-Curves

Synchronous generators (also called alternators) produce virtually all commercial electrical energy worldwide across thermal, nuclear, gas turbine, and hydroelectric power stations. Unlike induction machines, synchronous machines operate at a fixed, unyielding speed directly proportional to electrical grid frequency ($N = N_s = 120 f / P$) and utilize an independent Direct Current (DC) field excitation system on the rotor.

Understanding synchronous machine modeling, phasor analysis under varying power factors, active and reactive power-angle equations, saliency effects, excitation V-curves, and P-Q capability boundaries is a primary focus of the NCEES PE Power examination.

---

1. Machine Construction: Cylindrical vs. Salient-Pole Rotors

Synchronous generators are categorized mechanically into two primary rotor topologies based on prime mover operating speed.

+-----------------------------------------------------------------------------+
|                   SYNCHRONOUS GENERATOR ROTOR TOPOLOGIES                    |
|                                                                             |
|   CYLINDRICAL / ROUND ROTOR                     SALIENT-POLE ROTOR          |
|   - High-speed (1,800 or 3,600 RPM)             - Low-speed (100 - 600 RPM) |
|   - 2 or 4 poles                                - Many poles (12 to 72 poles|
|   - Steam and gas turbines (Turbo-generators)   - Hydro-turbines and diesels|
|   - Solid steel forging, long axial length,     - Large diameter, short     |
|     small rotor diameter                          axial length              |
|   - Uniform air gap (X_d ≈ X_q = X_s)           - Non-uniform air gap       |
|                                                   (X_d > X_q)               |
+-----------------------------------------------------------------------------+

Mechanical Comparison Matrix

---

2. Equivalent Circuit & Phasor Analysis

The per-phase steady-state equivalent circuit of a cylindrical-rotor synchronous generator consists of an internal generated EMF $\mathbf{E}_f$ in series with the synchronous impedance $\mathbf{Z}_s = R_a + j X_s$.

                  PER-PHASE SYNCHRONOUS GENERATOR CIRCUIT
                  
             +------------[ R_a ]-----------[ jX_s ]------------+ A
             |            Armature          Synchronous         |
             |           Resistance          Reactance          |
           ( + )                                                |
      E_f  ( ~ ) Internal Generated EMF                  I_a    | Load Terminals
     (Excitation)                                       --->    | (Grid Bus)
           ( - )                                                |
             |                                                + | 
             |                                              V_t | Terminal Voltage
             |                                                - | 
             +--------------------------------------------------+ N

Synchronous Reactance Components

The total synchronous reactance $X_s$ accounts for two distinct physical phenomena:

$$X_s = X_{al} + X_{ar}$$

1. Armature Leakage Reactance ($X_{al}$): Leakage flux around stator slot conductors that does not cross the air gap.
2. Armature Reaction Reactance ($X_{ar}$): Magnetic flux produced by stator load current $\mathbf{I}_a$ crossing the air gap and modifying the rotor field flux.

Voltage Equation

$$\mathbf{E}_f = \mathbf{V}_t + \mathbf{I}_a (R_a + j X_s) \approx \mathbf{V}_t + j \mathbf{I}_a X_s \quad (\text{since } R_a \ll X_s)$$

                          GENERATOR PHASOR DIAGRAMS

   LAGGING PF (Overexcited, Q > 0)             LEADING PF (Underexcited, Q < 0)
            E_f                                         V_t
           /  ^                                        /  ^
          /   | j*I_a*X_s                             /   | j*I_a*X_s
         /    |                                      /    |
        /  δ  |                                  E_f/  δ  |
       +------+------> V_t                         +------+------> I_a (Leads V_t)
        \   θ |                                     \ θ
         \    |                                      \
          v   I_a (Lags V_t)                          v V_t
       |E_f| > |V_t|  ===> Supplies VARs           |E_f| < |V_t|  ===> Absorbs VARs

---

3. Power-Angle Equations & Static Stability Limit

Let the terminal voltage phasor be the reference: $\mathbf{V}_t = V_t \angle 0^\circ$. The internal generated EMF is $\mathbf{E}_f = E_f \angle \delta$, where $\delta$ is the power angle (torque angle / rotor angle).

+-----------------------------------------------------------------------------+
|                 ROUND-ROTOR POWER-ANGLE FORMULATIONS (3-PHASE)              |
|                                                                             |
|   Real Electrical Power Output (P):                                         |
|             P = (3 * E_f * V_t / X_s) * sin(δ)               [Watts]        |
|                                                                             |
|   Reactive Electrical Power Output (Q):                                     |
|             Q = (3 * E_f * V_t / X_s) * cos(δ) - (3 * V_t^2 / X_s)  [VAR]   |
+-----------------------------------------------------------------------------+

(Note: For per-phase formulations, omit the factor of 3 and use line-to-neutral voltages $V_t$ and $E_f$.)

                   POWER-ANGLE CURVE & STABILITY BOUNDARIES

   Power P
     ^
     |                     P_max = (3*E_f*V_t / X_s)
     |                        /---\
     |                       /     \
     |                      /       \
     |   P_rated           /         \
     |      *-------------*           \
     |      |            /             \
     |      |           /               \
     |      |          /                 \
     +------+---------+-------------------+-----------------------------> Angle δ
     0°     δ_rated   δ = 90°             180°
            (20°-35°) (Static Limit)
     <-- STABLE REGION --><------ UNSTABLE REGION ----->

Static Steady-State Stability Limit

• Maximum Power Transfer: Differentiating $P$ with respect to $\delta$ yields the Synchronizing Power Coefficient $P_{syn}$:

$$P_{syn} = \frac{dP}{d\delta} = \frac{3 E_f V_t}{X_s} \cos\delta$$

• Stability Criterion: For the generator to remain synchronized following a small perturbation, $P_{syn} > 0$. This restricts steady-state operation to:

$$\delta < 90^\circ$$

• Maximum Theoretical Power Capability:

$$P_{max} = \frac{3 E_f V_t}{X_s} \quad \text{at } \delta = 90^\circ$$

• Practical Operating Margin: Commercial utility generators operate at full load with $\delta$ between $20^\circ$ and $35^\circ$, preserving a massive transient stability margin to survive external grid faults and line switching.

---

4. Salient-Pole Machines & Two-Reaction Theory (Blondel)

In salient-pole generators, the non-uniform air gap causes the magnetic reluctance along the direct pole axis ($d$-axis) to be far lower than the reluctance along the interpolar quadrature axis ($q$-axis).

+-----------------------------------------------------------------------------+
|                     TWO-REACTION IMPEDANCE PARAMETERS                       |
|                                                                             |
|   Direct-Axis Synchronous Reactance (X_d):                                  |
|   - Small air gap at pole center ===> High permeance ===> Large reactance   |
|                                                                             |
|   Quadrature-Axis Synchronous Reactance (X_q):                              |
|   - Large air gap between poles ===> Low permeance ===> Smaller reactance   |
|                                                                             |
|   Typical Saliency Ratio:   X_q ≈ 0.60 to 0.75 * X_d                        |
+-----------------------------------------------------------------------------+

Armature Current Decomposition

$$\mathbf{I}_a = \mathbf{I}_d + \mathbf{I}_q$$ $$\mathbf{E}_f = \mathbf{V}_t + j \mathbf{I}_d X_d + j \mathbf{I}_q X_q$$

Salient-Pole Power Equation

Resolving the real electrical power output yields two distinct terms:

$$P = \underbrace{\frac{E_f V_t}{X_d} \sin\delta}_{\text{Excitation Power Component}} + \underbrace{\frac{V_t^2}{2} \left(\frac{1}{X_q} - \frac{1}{X_d}\right) \sin(2\delta)}_{\text{Reluctance Power Component}}$$

+-----------------------------------------------------------------------------+
|                   RELUCTANCE POWER CHARACTERISTICS                          |
|                                                                             |
|   1. Excitation Independence: Reluctance power exists even if field         |
|      excitation is completely lost (E_f = 0). It is caused by rotor iron    |
|      aligning with the stator rotating magnetic field to minimize reluctance.|
|   2. Double Frequency Angle (2δ): Reluctance power peaks at δ = 45°.         |
|   3. Total Peak Shift: Due to the sin(2δ) component, maximum power P_max    |
|      in a salient-pole generator occurs at δ ≈ 70° - 75° (less than 90°).   |
+-----------------------------------------------------------------------------+

---

5. Generator V-Curves & Inverted V-Curves

V-Curves depict armature current ($I_a$) as a function of field excitation current ($I_f$) at various constant active power output levels ($P = 0\%, 50\%, 100\%\text{ MW}$).

                        SYNCHRONOUS GENERATOR V-CURVES

   Armature Current I_a
     ^
     |   \   UNDEREXCITED   /       \    OVEREXCITED   /
     |    \  (Leading PF)  /         \   (Lagging PF) /
     |     \              /           \              /
     |      \            /   100% MW   \            /
     |       \          /               \          /
     |        \        /                 \        /     UNITY POWER FACTOR
     |         \      *-------------------*------/----  LOCUS (cos θ = 1.0)
     |          \    /    50% MW           \    /
     |           \  /                       \  /
     |            *                          *          0% MW (No Load)
     +------------+--------------------------+--------------------------> Field Current I_f
            Under-Excited                Over-Excited
            (Absorbs -Q)                 (Supplies +Q)

Key V-Curve Takeaways:

1. Minimum Armature Current Points: For any constant real power level, minimum stator current occurs at exactly Unity Power Factor ($\cos\theta = 1.0$).
2. Right Side of Unity Line ($I_f > I_{f,\text{unity}}$): Generator is overexcited, operating at a lagging power factor, supplying reactive power ($+Q$) to support grid voltage.
3. Left Side of Unity Line ($I_f < I_{f,\text{unity}}$): Generator is underexcited, operating at a leading power factor, absorbing reactive power ($-Q$) from the grid.
4. Inverted V-Curves: Plot power factor ($\cos\theta$) versus field current ($I_f$), forming inverted curves peaking at $1.0$ at nominal excitation.

---

6. Generator P-Q Capability Curve

The P-Q Capability Curve defines the safe steady-state thermal and stability operating envelope of a utility synchronous generator.

                   SYNCHRONOUS GENERATOR P-Q CAPABILITY CURVE

       Real Power P (MW)
              ^
         P_max|                  [ Turbine Real Power Rating ]
              |           /---------------------------------\
              |          /                                   \  Rotor Field
              |         /   Stator Armature                   \ Winding Heating
              |        /    Current Heating                    \ Limit
              |       /     Limit Circle                        \
              |      /      (Radius = S_rated)                   \
              |     /                                             \
              |    /                                               \
              |   / Stator End-Iron Heating                         \
              |  /  & Stability Limit (MEL)                          \
   -Q (MVAR)  +--+---------------------------------------------------+----> +Q (MVAR)
   (Leading / |                                                      (Lagging /
  Underexcited)                                                     Overexcited)

Three Primary Physical Operating Limits:

1. Armature Current Heating Limit (Stator Winding Limit):
   • Symmetrical circle centered at origin $(P=0, Q=0)$ with radius equal to rated apparent power: $S_{\text{rated}} = \sqrt{P^2 + Q^2}$.
   • Prevents overheating and insulation breakdown of the three-phase stator copper windings.
2. Field Current Heating Limit (Rotor Winding Limit):
   • Circular boundary centered at coordinates $(0, -\frac{3 V_t^2}{X_s})$ with radius $\frac{3 E_{f,max} V_t}{X_s}$.
   • Dominates the overexcited (lagging) region, limiting maximum continuous DC field current to prevent rotor thermal failure.
3. Stator End-Iron Core Heating & Steady-State Stability Limit:
   • Dominates the underexcited (leading) region.
   • In the underexcited regime, high axial leakage flux enters stator core end laminations perpendicular to the sheets, inducing massive localized eddy currents and extreme heating in the clamping fingers.
   • Protected by the Minimum Excitation Limiter (MEL) and loss-of-field relays (ANSI Device 40).

---

7. Step-by-Step Worked Mathematical Example

Problem Statement:

A 3-phase, 13.8 kV (line-to-line), 50 MVA, 60 Hz, Y-connected cylindrical-rotor synchronous generator has a synchronous reactance of $X_s = 3.20\ \Omega$ per phase (armature resistance $R_a$ is negligible). The generator delivers rated MVA at a lagging power factor of $0.80$ into an infinite power grid held at $13.8\text{ kV}$.

Calculate:

1. Rated line current $I_a$ and terminal line-to-neutral voltage $V_t$.
2. Internal generated excitation EMF per phase ($E_f$), line-to-line excitation EMF ($E_{f,LL}$), and power angle $\delta$.
3. Total active power ($P$) and reactive power ($Q$) delivered to the grid.
4. Maximum theoretical steady-state power capability ($P_{max}$) at this excitation level and the static pull-out torque at $1,800\text{ RPM}$.

---

Step-by-Step Solution:

Step 1: Terminal Quantities & Rated Current Per-phase line-to-neutral terminal voltage reference:

$$V_t = \frac{13,800\text{ V}}{\sqrt{3}} = 7,967.43\text{ V} \angle 0^\circ$$ $$\text{Rated Armature Current: } I_a = \frac{S_{\text{3phase}}}{\sqrt{3} \times V_{\text{LL}}} = \frac{50 \times 10^6\text{ VA}}{\sqrt{3} \times 13,800\text{ V}} = \frac{50,000,000}{23,902.30} = 2,091.85\text{ A}$$ $$\text{Power Factor Angle: } \theta = \arccos(0.80) = 36.87^\circ \quad (\text{lagging } \implies \mathbf{I}_a = 2,091.85 \angle -36.87^\circ\text{ A})$$ $$\mathbf{I}_a = 2,091.85 \left(\cos(-36.87^\circ) + j \sin(-36.87^\circ)\right) = 1,673.48 - j 1,255.11\text{ A}$$

Step 2: Calculate Internal Generated EMF ($E_f$) & Power Angle ($\delta$)

$$\mathbf{E}_f = \mathbf{V}_t + j \mathbf{I}_a X_s = 7,967.43 + j (3.20\ \Omega)(1,673.48 - j 1,255.11)$$ $$\mathbf{E}_f = 7,967.43 + j 5,355.14 - j^2 (4,016.35) = (7,967.43 + 4,016.35) + j 5,355.14 = 11,983.78 + j 5,355.14\text{ V}$$ $$\text{Magnitude } |E_f| = \sqrt{11,983.78^2 + 5,355.14^2} = \sqrt{143,611,000 + 28,677,525} = \sqrt{172,288,525} = 13,125.87\text{ V}$$ $$\text{Line-to-Line Generated EMF: } E_{f,LL} = \sqrt{3} \times 13,125.87\text{ V} = 22,734.7\text{ V} = 22.735\text{ kV}$$ $$\text{Power Angle: } \delta = \arctan\left(\frac{5,355.14}{11,983.78}\right) = \arctan(0.44686) = 24.08^\circ$$ $$\mathbf{E}_f = 13,125.87\text{ V} \angle 24.08^\circ$$

Step 3: Verify Active and Reactive Power Flow

$$P = \frac{3 E_f V_t}{X_s} \sin\delta = \frac{3 (13,125.87)(7,967.43)}{3.20} \sin(24.08^\circ) = \frac{313,738,400}{3.20} \times 0.4080 = 98,043,250 \times 0.4080 = 40,001,600\text{ W} = 40.00\text{ MW}$$ $$\text{Direct Check: } P = S \cos\theta = 50\text{ MVA} \times 0.80 = 40.00\text{ MW} \quad (\text{Matches exact!})$$ $$Q = \frac{3 E_f V_t}{X_s} \cos\delta - \frac{3 V_t^2}{X_s} = (98.043 \times \cos 24.08^\circ) - \frac{3 (7,967.43)^2}{3.20} = (98.043 \times 0.9130) - 59.513 = 89.513 - 59.513 = 30.00\text{ MVAR}$$ $$\text{Direct Check: } Q = S \sin\theta = 50\text{ MVA} \times 0.60 = 30.00\text{ MVAR} \quad (\text{Matches exact!})$$

Step 4: Maximum Power Capability & Pull-Out Torque

$$P_{max} = \frac{3 E_f V_t}{X_s} = \frac{3 (13,125.87)(7,967.43)}{3.20} = 98.043\text{ MW}$$ $$\text{Synchronous Speed: } \omega_s = \frac{2\pi \times 1,800}{60} = 188.496\text{ rad/s}$$ $$T_{pull-out} = \frac{P_{max}}{\omega_s} = \frac{98,043,250\text{ W}}{188.496\text{ rad/s}} = 520,134\text{ N}\cdot\text{m} = 520.13\text{ kN}\cdot\text{m}$$

$\text{Steady-State Stability Margin: } \text{SSM} = \frac{P_{max} - P_{rated}}{P_{rated}} = \frac{98.043 - 40.00}{40.00} = 145.1\%$$