Electronics Guide

Rotor Dynamics

The rotor of a synchronous generator obeys Newton's second law for rotation. With J the combined moment of inertia of the generator and its turbine, ω_m the mechanical angular velocity, T_m the mechanical torque of the prime mover, and T_e the electromagnetic torque of the generator,

J dω_m/dt = T_m − T_e

In steady state the torques balance and the rotor turns at synchronous speed. A fault or a switching event changes T_e almost instantly, while T_m changes only as fast as the turbine controls allow. The difference accelerates or decelerates the rotor, and its angle relative to the rest of the system changes.

The Inertia Constant H

Stability studies describe inertia by the inertia constant H, the kinetic energy stored in the rotating mass at synchronous speed divided by the machine's rated apparent power:

H = ½ J ω_sm² / S_rated

With J in kg·m², ω_sm the synchronous mechanical speed in rad/s, and S_rated in volt-amperes, H is in seconds, usually quoted as MW·s/MVA on the machine rating. A two-pole 60 Hz turbine-generator turns at 3,600 r/min, or 377 rad/s; with an illustrative J of 50,000 kg·m², it stores 3,553 MJ, an H of 5.92 s on a 600 MVA rating. ERCOT, the operator of the Texas Interconnection, gives the physical reading: a 200 MVA unit with an inertia constant of 3 s would use up all of its kinetic energy in 3 s if producing at full output. Twice H is the mechanical starting time, the time rated torque would take to accelerate the rotor from standstill to rated speed.

The Per-Unit Form and Its Conventions

Let δ be the electrical angle of the rotor measured against a reference that rotates at synchronous speed, so that δ is constant in steady state. For a machine with p poles, electrical angles are p/2 times mechanical angles, and ω_s = (p/2)ω_sm = 2πf₀ is the synchronous speed in electrical radians per second, where f₀ is the nominal frequency. Substituting H and dividing by the rating gives the swing equation in its usual form:

(2H / ω_s) d²δ/dt² = P_m − P_e

• P_m is the mechanical input power and P_e the electrical output power, both in per unit on the machine rating, the same base as H.
• H is in seconds, δ in electrical radians, and t in seconds; ω_s is 377 rad/s on a 60 Hz system and 314 rad/s on a 50 Hz system.
• The derivation takes per-unit torque as equal to per-unit power, which is accurate while speed stays close to synchronous, as it does during electromechanical swings.
• Damping is neglected. A term DΔω subtracted from the right side, with Δω = (1/ω_s) dδ/dt the per-unit speed deviation and D in per-unit power per per-unit speed, represents damper windings and the frequency sensitivity of load.

With δ in electrical degrees the equation becomes (H / 180f₀) d²δ/dt² = P_m − P_e. A positive right side, more mechanical input than electrical output, advances the rotor ahead of the synchronous reference.

Several Machines on a Common Base

A multimachine study writes one swing equation per machine on a common system base, converting H by the power ratio alone: H_sys = H_mach S_mach / S_sys. Machines that swing together, such as the units of one plant, can be combined into one equivalent machine whose H·S is the sum of theirs. The network couples the equations, so the complete model joins differential equations for the machines and controls with algebraic equations for the network.

The Classical Model

The simplest useful model for rotor angle studies, called the classical model, rests on four assumptions:

• Each generator is a constant voltage magnitude E′ behind its direct-axis transient reactance X′_d, the reactance described in Fault Analysis and Symmetrical Components, and the angle of E′ stands in for the rotor angle δ.
• Mechanical power stays constant, because turbine governors act over seconds.
• Damping and the actions of excitation systems are neglected.
• The network is lossless and in sinusoidal steady state at nominal frequency, and loads are constant impedances.

The model suits the first swing after a fault, which in the worked example below peaks within about 0.7 s, but not later swings that controls can make grow. Production programs add detailed models of field and damper windings, exciters, and governors from synchronous machine theory.

Power Transfer and Equilibrium

Consider a generator connected to an infinite bus, a bus whose voltage magnitude V and frequency stay constant whatever the generator does, through a total transfer reactance X that includes X′_d. Its electrical power is

P_e = (E′V / X) sin δ = P_max sin δ

The curve rises to P_max = E′V/X at δ = 90° and falls beyond it. A horizontal line at P_m crosses it twice, at δ_s = arcsin(P_m/P_max) and at δ_u = 180° − δ_s. The first crossing is stable: if the rotor advances past δ_s, electrical output exceeds mechanical input and the rotor decelerates back. The second is unstable: a rotor carried past δ_u finds electrical output below mechanical input, keeps accelerating, and slips poles. The same sine law governs the power carried by a line between two buses, and Transmission and Distribution Lines uses it for line loadability, including the St. Clair curve and its 30 percent steady-state stability margin.

Synchronizing Power

The slope of the curve at the operating point is the synchronizing power coefficient, in per unit of power per electrical radian:

K_s = dP_e/dδ = P_max cos δ_s

A positive K_s pulls the rotor back after a small displacement. Heavier loading moves δ_s toward 90° and reduces K_s, and a weaker network, with a larger X and a lower P_max, does the same. At δ_s = 90°, K_s reaches zero and no steady-state margin remains.

A solved power flow, found as Power Flow Analysis describes, supplies the starting point for any stability study. From a generator's terminal voltage V_t and current I, the internal voltage of the classical model is E′ = V_t + jX′_dI, as the worked example shows.

The Criterion

Multiply both sides of the swing equation by dδ/dt and integrate. Starting at synchronous speed at δ₀, the rotor satisfies

(H / ω_s) (dδ/dt)² = ∫ (P_m − P_e) dδ, integrated from δ₀ to δ

The left side is proportional to the kinetic energy of the rotor's motion relative to the synchronous reference, and the right side is a net area on the power–angle diagram. The rotor turns back where its speed deviation, and with it the net area, returns to zero. The system therefore survives the first swing if the decelerating area available above the P_m line on the postfault curve, up to δ_u, at least equals the accelerating area gained below that line before clearing. This equal-area criterion applies exactly to one machine against an infinite bus, or to two machines, under the assumptions of the classical model.

Critical Clearing Angle and Time

If the fault clears at the angle δ_c, the accelerating area is A₁ = ∫ (P_m − P₂ sin δ) dδ from δ₀ to δ_c, and the largest decelerating area is A₂ = ∫ (P₃ sin δ − P_m) dδ from δ_c to δ_u, where δ_u = 180° − arcsin(P_m/P₃). The critical clearing angle δ_cr makes the two equal. Evaluating the integrals, with angles in radians, gives

cos δ_cr = [P_m(δ_u − δ₀) + P₃ cos δ_u − P₂ cos δ₀] / (P₃ − P₂)

The result applies only when P₃ > P_m, so that the postfault curve has an equilibrium. The textbook special case is an idealized temporary fault at the generator's high-voltage bus that removes all transfer while it lasts (P₂ = 0) and clears without any change to the network (P₃ = P₁). Then δ_u = π − δ₀ and P_m = P₁ sin δ₀, and the formula reduces to

δ_cr = arccos[(π − 2δ₀) sin δ₀ − cos δ₀]

Protection engineers need a time rather than an angle. While P_e = 0, the rotor has the constant acceleration ω_sP_m/2H, so δ = δ₀ + ω_sP_mt²/4H, and the critical clearing time is

t_cr = √[4H(δ_cr − δ₀) / (ω_sP_m)]

with angles in radians. This expression holds for any fault with P₂ = 0, whether or not the postfault network matches the prefault one. When some power still flows during the fault (P₂ > 0), the area formula still gives δ_cr, but finding the time at which the rotor reaches it requires integrating the swing equation numerically. Either way, H only rescales time in the classical model: δ_cr does not depend on H, and every critical clearing time grows with the square root of H.

Assumptions and Limits

• The criterion tests the first swing only. A system can survive the first swing and lose synchronism on a later one when controls or interacting machines feed energy into the oscillation.
• A constant E′ ignores both the decay of field flux and the field forcing of a fast excitation system, which raises the internal voltage during the fault and helps.
• Neglecting damping usually makes the first-swing result conservative.
• A system of many machines has no single power–angle diagram. Numerical integration of all the machine equations is the working method, and direct methods based on transient energy functions extend the area idea to many machines.

Comparing t_cr with the actual clearing time of the protection, including the longer backup clearing when a breaker or relay fails, is the core of a transient stability check. Power System Protection describes the relays and breakers that set that time. The same area picture explains the usual remedies: faster clearing shrinks the accelerating area; single-pole tripping, additional lines, series capacitors, and fast excitation enlarge the decelerating area; and fast valving of steam turbines, braking resistors, and generator tripping schemes reduce the accelerating power.

The Prefault Operating Point

1. The reactance from the generator terminals to the infinite bus is X_T + X_L/2 = 0.37. From P = |V_t|V sin θ_t / 0.37, the terminal voltage leads the infinite bus by θ_t = arcsin(0.90 × 0.37) = 19.45°, so V_t = 0.9429 + j0.3330 pu.
2. The current is I = (V_t − V) / j0.37 = 0.9000 + j0.1543 pu, so the generator also supplies 0.154 pu of reactive power at its terminals.
3. The internal voltage is E′ = V_t + j0.30I = 0.8967 + j0.6030 pu, or 1.0806 e^j33.92° pu, so δ₀ = 33.92°.

Three Networks

• Prefault. X₁ = 0.30 + 0.12 + 0.25 = 0.67, so P₁ = 1.0806/0.67 = 1.6128 pu. As a check, 1.6128 sin 33.92° = 0.900.
• Fault on. A bolted three-phase fault at the midpoint of line 2 leaves three branches at the high-voltage bus: 0.42 toward E′, 0.50 through line 1 toward the infinite bus, and 0.25 through half of line 2 to ground. Converting that wye into a delta gives the transfer reactance X₂ = (0.42 × 0.50 + 0.50 × 0.25 + 0.25 × 0.42) / 0.25 = 1.76, so P₂ = 0.6139 pu.
• Postfault. Breakers at both ends open line 2, leaving X₃ = 0.30 + 0.12 + 0.50 = 0.92 and P₃ = 1.1745 pu. This curve crosses P_m at δ_s = 50.02° and δ_u = 129.98°, or 2.2686 rad.

Clearing the Midpoint Fault

With δ₀ = 0.5920 rad, the terms of the area formula are P_m(δ_u − δ₀) = 1.5089, P₃ cos δ_u = −0.7546, and P₂ cos δ₀ = 0.5095. Then

cos δ_cr = (1.5089 − 0.7546 − 0.5095) / (1.1745 − 0.6139) = 0.2448 / 0.5606 = 0.4367

so δ_cr = 64.11°, where both areas equal 0.2328. Integrating the fault-on swing equation numerically gives t_cr = 0.2085 s, or 12.5 cycles. The rotor then runs 4.684 electrical rad/s, or 1.24 percent, above synchronous speed, and (H/ω_s)(dδ/dt)² = (4.0/377.0) × 4.684² = 0.2328, the accelerating area, as the energy relation requires. Holding the initial acceleration of 26.27 rad/s² constant would give 0.2003 s, slightly short, because the acceleration falls as P₂ sin δ grows.

Swing Curves

Cleared at 0.100 s, the rotor peaks at 79.4° at 0.47 s and swings back to 25.9° at 1.02 s, an oscillation period of 1.11 s. Cleared at 0.204 s, just inside the critical time, it peaks at 120.1°, close to δ_u, and still returns. Cleared at 0.214 s, it passes δ_u at 0.67 s and slips poles until out-of-step protection separates it. Without damping in the model, the stable swings never decay; in a real machine, damper windings, load, and stabilizers damp them.

Fault Location and Inertia

The closer the fault lies to the generator, the less power gets through while it lasts and the shorter the critical time: a close-in fault that also costs a line leaves only 0.115 s, about 7 cycles. The last row follows from the closed-form expressions, δ_cr = arccos 0.2626 = 74.78° and t_cr = 0.183 s, a longer time than for the close-in fault with the line removed because the prefault network returns.

Inertia acts through time, not angle. With H = 2, 3, 4, and 6 s, the midpoint fault has critical clearing times of 0.147, 0.181, 0.209, and 0.255 s, and the close-in fault 0.081, 0.100, 0.115, and 0.141 s. Both sets are exactly proportional to the square root of H. The 2021 stability paper makes the system-level point: lower total system inertia can produce larger and faster rotor swings.

Linearizing the Swing Equation

For a small deviation Δδ about δ_s, the electrical power changes by K_sΔδ, and the swing equation with damping becomes linear:

(2H / ω_s) d²Δδ/dt² + (D / ω_s) dΔδ/dt + K_sΔδ = 0

Its undamped natural frequency is ω_n = √(ω_sK_s / 2H) in rad/s, and its damping ratio is ζ = D / (4Hω_n). In the worked example, the prefault operating point has K_s = 1.338 pu/rad, so ω_n = √(377.0 × 1.338 / 8.0) = 7.94 rad/s, an oscillation of 1.26 Hz. After line 2 opens, the weaker network and the higher angle reduce K_s to 0.755 and the frequency to 0.95 Hz. The large swings of the simulation are slower still, 1.11 s per cycle, or 0.90 Hz, because the power–angle curve flattens at large angles.

Too little synchronizing power lets the angle drift away without oscillating, and too little damping lets an oscillation grow; in practice, small-disturbance rotor angle problems usually appear as poorly damped oscillations.

Eigenvalues, Mode Shapes, and Participation

A large system is linearized about its operating point into the state-space form dx/dt = Ax, with states for every machine, exciter, governor, and control. Each complex pair of eigenvalues λ = σ ± jω of A is an oscillatory mode with frequency f = ω/2π and damping ratio ζ = −σ / √(σ² + ω²). A mode at λ = −0.20 ± j6.28, for example, oscillates at 1.00 Hz with a damping ratio of 0.032, or 3.2 percent. Oscillatory instability occurs when such a pair moves into the right half of the complex plane.

Eigenvectors give each mode's shape, the relative amplitude and phase of each machine's swing, and participation factors show which states, and so which machines and controls, influence a mode most. State-Space Analysis and Control develops eigenvalues, modes, and eigenvectors.

Local and Inter-Area Modes

The North American Electric Reliability Corporation (NERC) sorts natural power system oscillations into four types in its report Interconnection Oscillation Analysis (July 2019):

• Local: one plant or unit swinging against the rest of the system, generally caused by heavy loading and generator controls.
• Intra-plant: units within one plant swinging against each other.
• Inter-area: groups of machines in one area swinging against groups in another, often where the ties between areas are relatively weak.
• Torsional: subsynchronous oscillations from resonance between highly compensated lines and the mechanical modes of a steam turbine-generator.

Inter-area modes are slow because each end carries the inertia of many machines and the tie between them is weak. From recorded events in the three major North American interconnections, the report identified the dominant modes in the table, with average damping ratios from 9.26 to 16.45 percent. It treats damping above 10 percent as well damped, without calling lower damping poor, and leaves damping criteria to each reliability coordinator's system operating limit methodology.

Against those figures, the worked example's 0.95 to 1.26 Hz is a local-mode frequency. The report also describes forced oscillations, sustained oscillations driven by external inputs such as equipment failures or control interactions; near the frequency of a natural mode, they can excite it across an interconnection. On June 17, 2016, a control valve malfunction at the Grand Gulf Nuclear Station in Mississippi drove a 200 MW oscillation near 0.27 Hz that excited an Eastern Interconnection mode, producing 40 MW swings on a tie line between New York and New England about 1,400 miles away.

When an Oscillation Breaks a System

The Western Interconnection breakup of August 10, 1996, summarized in NERC's 2019 report, shows how such a mode can end. During heavy transfers from Canada through the Pacific Northwest into California, erroneous relay operations removed all 13 generating units at McNary, and a poorly damped inter-area oscillation began and became negatively damped. When the swings reached about 1,000 MW and 60 kV peak to peak at the Malin 500 kV substation, voltage collapsed, and the interconnection separated into four islands.

Synchronizing and Damping Components

For small oscillations, the change in electrical power splits into a component in phase with the angle deviation and a component in phase with the speed deviation:

ΔP_e = K_sΔδ + K_DΔω

K_s is the synchronizing coefficient and K_D the damping coefficient. A voltage regulator contributes to both, and the phase lag of the exciter and field winding can make its contribution to K_D negative.

How a Stabilizer Works

A power system stabilizer (PSS) adds a supplementary signal at the voltage regulator's summing point so that the excitation produces a component of electrical power in phase with speed deviation, which is pure damping. A typical design has four parts:

1. An input that reflects rotor speed deviation: shaft speed, terminal frequency, electrical power, or a combination, such as dual-input designs that synthesize the integral of accelerating power from speed and electrical power.
2. A washout, a high-pass filter that passes the oscillations but blocks steady offsets, so that the stabilizer does not fight slow changes in speed or voltage.
3. Lead–lag phase compensation to offset the lag between the stabilizer output and the electrical power it produces, a lag introduced by the exciter, the field winding, and the network.
4. A gain and output limits, with filters that keep shaft torsional frequencies in a speed signal from reaching the exciter.

Because the lag to be compensated varies with frequency, one stabilizer must provide useful phase compensation across the modes it should damp, from inter-area modes of a few tenths of a hertz to local modes near and above 1 Hz. IEEE Std 421.5-2016, IEEE Recommended Practice for Excitation System Models for Power System Stability Studies, gives standard stabilizer models along with its exciter models. Hydroelectric Power Electronics describes stabilizers, excitation systems, and governors in hydro plants, where dual-input designs avoid spurious output when gate movements change mechanical power.

Tuning and Placement

Engineers tune a stabilizer from the measured or computed frequency response of the generator, exciter, and network and confirm the settings with field tests, and participation factors show where stabilizers will most affect a poorly damped mode. Power oscillation damping controls on FACTS devices, such as static VAR compensators and STATCOMs, and on HVDC links also damp oscillations by modulating reactive or real power. The 2021 paper sees similar potential in supplemental controls on converter-interfaced generation, and it warns that displacing synchronous generators can remove machines whose stabilizers are crucial.

The Rate of Change of Frequency

Writing each machine's swing equation in megawatts with speed expressed as frequency, and summing over all synchronized machines, gives the relation for the frequency f of the center of inertia:

df/dt = f₀ (P_m − P_e) / (2∑H_iS_i)

Here f₀ is nominal frequency in hertz, P_m − P_e is the imbalance in megawatts, and ∑H_iS_i is the kinetic energy, in megawatt-seconds, of the machines still synchronized, each H_i on its own rating S_i. A loss of generation makes the imbalance negative, so frequency falls. The relation gives the initial rate of change of frequency (RoCoF), before load relief and governors act; tripped units are left out of the sum, because their kinetic energy leaves with them. If a 60 Hz system loses 1,200 MW with 240 GW·s still synchronized, such as 60,000 MVA of machines at an average H of 4.0 s, the initial RoCoF is 60 × (−1,200) / (2 × 240,000) = −0.15 Hz/s.

Three Periods of Response

The 2021 paper illustrates three periods in the response of a system dominated by synchronous generators to a major loss of generation, with typical times taken from a 2013 IEEE Power & Energy Society report on turbine-governor models:

1. Inertial response. Kinetic energy flows out of the rotating machines at once, with no control action, and the system inertia sets the initial slope of the decline.
2. Primary frequency response. Governors open valves and gates in proportion to the frequency deviation, and frequency-sensitive load, such as motors, draws less power. Together they arrest the decline at the frequency nadir, typically 5 to 10 s after the event, and hold frequency at a settling value below nominal, typically reached 20 to 30 s after it.
3. Automatic generation control. Secondary control changes generator setpoints to return frequency to nominal and restore scheduled interchange between areas, typically within 5 to 10 min.

Droop and Deadband

Droop is the per-unit frequency change that moves a unit across its full output range: in the example of the Federal Energy Regulatory Commission (FERC), 5 percent droop means that 3 Hz, 5 percent of 60 Hz, moves the prime mover's control from fully closed to fully open. In Order No. 842, issued February 15, 2018, FERC required newly interconnecting generating facilities, synchronous and non-synchronous, to have governors or equivalent controls with a maximum 5 percent droop, based on nameplate capacity, and a maximum deadband of ±0.036 Hz, the band of small deviations that draws no response. The order also requires a deadband without a step, so response starts from zero at its edge. With 5 percent droop, a 0.1 Hz fall on a 60 Hz system then raises a 100 MW unit's output by 100 × [(0.1 − 0.036)/60] / 0.05 = 2.1 MW, provided the unit has headroom to rise.

Once governors settle, the steady-state deviation depends on the load damping D, the per-unit change in load per per-unit change in frequency, and on the governors' combined 1/R_eff, both on the system base:

Δf / f₀ = −ΔP / (D + 1/R_eff)

where ΔP is the lost generation in per unit.

An Illustrative Event

A lumped model shows how inertia and governors divide the work. Take the system above and assume, for illustration, D = 1 and governors at 5 percent droop with no deadband on 30 percent of the capacity, each responding with a first-order lag of 8 s, while the other units sit at their limits. The combined 1/R_eff is 0.30/0.05 = 6, so frequency settles at −(1,200/60,000) / (1 + 6) × 60 = −0.171 Hz from nominal, a response of 700 MW per 0.1 Hz, the unit North American practice uses. At that frequency, load relief supplies 171 MW of the loss and the governors 1,029 MW.

Halving the inertia doubles the initial RoCoF and deepens the nadir but leaves the settling frequency unchanged, because inertia stores energy and supplies none in steady state. Faster excursions, the 2021 paper notes, put more emphasis on fast-acting controllers that arrest frequency drops as soon as they are detected.

Operating Criteria and Inertia Limits

Commission Regulation (EU) 2017/1485, the EU guideline on electricity transmission system operation, sets the reference incident for continental Europe, the largest instantaneous imbalance for which frequency containment reserves are dimensioned, at 3,000 MW. For imbalances up to that size it sets a maximum instantaneous frequency deviation of 800 mHz and a maximum steady-state deviation of 200 mHz, and for deviations of 200 mHz or more it requires at least half of the full reserve within 15 s and all of it within 30 s.

ERCOT's white paper Inertia: Basic Concepts and Impacts on the ERCOT Grid also sets a floor on system inertia. Its critical inertia is the minimum at which load resources on under-frequency relays, set no lower than 59.7 Hz and responding in about 0.416 s, can act before frequency falls below 59.3 Hz, where under-frequency load shedding begins, after the simultaneous loss of its two largest units, 2,750 MW. Simulations put that level at 94 GW·s, which ERCOT raised to 100 GW·s for margin; its lowest inertia by then had been 130 GW·s, on October 27, 2017. If 100 GW·s remained synchronized after that loss, the relation above would give an initial RoCoF of 0.825 Hz/s.

Under-frequency load shedding is the last line of defense, disconnecting blocks of load in stages as frequency falls. Converter-interfaced generation provides no inherent inertial response, but the 2021 paper notes that it can deliver primary response faster and with smaller droop settings, since electronics rather than boilers limit its speed, and that battery storage can contribute decisively. Grid Synchronization and Control describes fast frequency response and synthetic inertia from converters.

The Nose of the PV Curve

The mechanism appears in the simplest case: a source of fixed voltage V_s feeding a constant-power load P + jQ through a lossless reactance X. Solving for the load voltage V gives

V² = [V_s² − 2QX ± √((V_s² − 2QX)² − 4X²(P² + Q²))] / 2

Below a maximum load there are two voltages for each load, a high one at which the system normally operates and a low one. Plotting V against P at constant power factor traces the PV curve, or nose curve, whose two branches meet at the tip where the square root vanishes. For a load of power factor cos φ, with φ positive for a lagging (inductive) load, the maximum power and the voltage at which it occurs are

P_max = V_s² cos φ / [2X(1 + sin φ)], at V = V_s / √[2(1 + sin φ)]

Reactive support, represented here by a leading power factor at the load bus, raises the maximum transfer, but it also raises the voltage at the nose: at 0.90 leading, the limit arrives at 0.94 pu, so a heavily compensated system can near collapse with voltages that still look acceptable. Near the nose, a small load increase causes a large voltage drop, and at the nose the power flow Jacobian becomes singular, as Power Flow Analysis shows along with the continuation power flow that locates the limit.

QV Curves and Reactive Reserves

A QV curve plots the reactive power that a fictitious synchronous condenser at a bus would have to supply to hold each voltage, with real loads fixed. The bus operates where the curve crosses zero, the bottom of the curve marks the voltage stability limit, and its depth below zero is the bus's reactive power margin. For the two-bus system with a load of 0.5 pu at 0.95 lagging, the curve crosses zero at 0.862 pu, matching the table, and reaches a minimum of −0.211 pu at 0.559 pu. The bus could therefore take on about 0.21 pu more reactive load, or lose that much reactive support, before no operating point remains.

In practice, reactive reserves are the capability of generators, synchronous condensers, static VAR compensators, and STATCOMs that remains after a contingency. Shunt capacitors help least when most needed, because their output falls with the square of voltage, and generators that reach their field or armature current limits further limit voltage support and the maximum transfer. Static VAR Compensators describes the fast-acting compensators.

Short-Term and Long-Term Mechanisms

Voltage stability is subdivided by the size of the disturbance and by time frame.

• Short-term voltage stability involves fast-acting loads and devices over several seconds: induction motors, electronically controlled loads, HVDC links, and inverter-based generators. In the typical case, induction motors slow during a fault and draw more current and reactive power; after clearing, a motor that has fallen below a critical speed stalls and draws starting current until protection disconnects it, holding voltage low and possibly stalling nearby motors in turn. HVDC links with line-commutated converters on weak AC systems can become voltage unstable within about a second; voltage-source converters have significantly raised those limits.
• Long-term voltage stability involves slower equipment over several minutes: tap-changing transformers, thermostatically controlled loads, and generator current limiters. It is usually set off not by a fault itself but by the outage of transmission or generation that follows it. Tap changers and thermostats restore load power after a voltage drop, and if the restored demand exceeds the maximum transfer of the weakened network, voltages decline progressively toward collapse.

Overvoltage instability has also occurred in a few cases: lines below their surge impedance loading, shunt capacitors, and HVDC filter banks can drive voltage up while underexcitation limiters keep generators and synchronous condensers from absorbing the excess reactive power.

Planners express long-term voltage stability as a margin, the load increase from the operating point to the maximum transfer, measured along a defined direction of stress: which loads grow, and which generators pick up the increase. Countermeasures include dynamic reactive reserves, fast switching of capacitors, blocking of tap changers during low-voltage emergencies, and undervoltage load shedding.

Resonance Stability

Resonance stability covers subsynchronous resonance in both of its forms. Torsional resonance exchanges energy between a series-compensated network and the torsional modes of a turbine-generator shaft, and its oscillations can be poorly damped or even growing, threatening the shaft. Electrical resonance, the induction generator effect, is a purely electrical self-excitation between a series capacitor and a machine that presents negative resistance at subsynchronous frequencies. According to the 2021 paper, it has never been observed with conventional synchronous generation, but wind turbines with doubly fed induction generators are susceptible: the first field event occurred in ERCOT in 2009, and similar events followed in the Xcel Energy network in Minnesota. Such events are called subsynchronous control interaction because the converter controls supply much of the negative damping. Fast controls on HVDC links, static VAR compensators, STATCOMs, and power system stabilizers can also interact with shaft torsional modes, although such interaction can improve torsional damping as well as degrade it.

Converter-Driven Stability

Converter-driven stability concerns oscillations that the controls of power electronic converters drive through their coupling with the network, with machines, and with each other. The 2021 paper divides it by frequency:

• Slow interaction, typically below 10 Hz, involves outer power and voltage loops and phase-locked loops interacting with slow system dynamics. It resembles voltage instability in that a weak connection limits power transfer, but converter controls rather than loads drive it. The paper reports that the phase-locked loop and inner current loop of a converter plant can become oscillatory below a short-circuit ratio typically between 1.5 and 2, depending on the vendor and the network, and that a high-gain phase-locked loop can inject current at the wrong phase angle during nearby faults on weak networks.
• Fast interaction, typically tens to hundreds of hertz and possibly into the kilohertz range, involves fast inner current loops interacting with the transmission network, filters, and other converters, a family that power electronics engineers call harmonic instability. Reported cases include oscillations between 500 Hz and 2 kHz in large wind plants connected through VSC HVDC links, and oscillations at 2.5 Hz and 97.5 Hz between a STATCOM and weak AC and DC grids in the China Southern Grid.

The paper lists remedies: lower gains in the phase-locked and current loops, supplemental controls, emerging control strategies, and equipment that strengthens the network, such as synchronous condensers. Grid-forming control, which regulates an internal voltage instead of injecting a current aligned by a phase-locked loop, is one such strategy. Grid-Forming Inverters explains its control laws and its own stability questions.

Converters also change rotor angle stability. The 2021 paper notes that they alter flows on major ties, displace large synchronous generators and their stabilizers, and can raise or lower the damping of nearby machines. It finds no general consensus on their net effect on electromechanical modes and reports that they can help or harm transient stability, depending on grid layout, location, and control.

Time-Domain Simulation

Time-domain simulation is the workhorse of stability analysis. A program initializes every dynamic model from a solved power flow, applies a scripted sequence of events, such as a fault, its clearing, and generator tripping, and integrates the machine and control equations with the network equations, over seconds for angle stability and minutes for long-term voltage and frequency studies. Engineers judge the resulting swing curves, voltages, and frequencies against criteria for synchronism, damping, and voltage recovery.

RMS and Electromagnetic-Transient Models

Root-mean-square (RMS), or phasor, simulation represents the network by fundamental-frequency phasors, which makes it efficient enough for studies of whole interconnections. Electromagnetic-transient (EMT) simulation solves instantaneous voltages and currents with much shorter time steps and captures converter controls, harmonics, and switching. Co-simulation tools combine electromechanical and electromagnetic models when fast converter phenomena matter. The Power System Analysis overview describes programs of both kinds, such as PSS/E, PowerFactory, and PSCAD.

Other Methods

Eigenvalue analysis finds modes, damping ratios, and participation factors; energy-function methods estimate transient stability margins without full simulation; PV and QV analysis, often with continuation power flow, finds voltage stability margins; and synchrophasor measurements track modes and forced oscillations in real time.

Planning Criteria

In North America, NERC reliability standard TPL-001-5.1 requires transmission planners to show that the system remains stable, without cascading or uncontrolled islanding, for its planning events. Those events range from the loss of a single element to faults with stuck breakers or failed protection, which lengthen clearing times and so test critical clearing directly. Utilities validate the generator, exciter, governor, stabilizer, and inverter models behind such studies against recorded disturbances and field tests.