Swing equation: rotor dynamics and power system stability

Overview

The swing equation is the fundamental mathematical model used to describe the dynamic behavior of synchronous machines within a power system. A power system consists of a number of synchronous machines operating synchronously under all operating conditions. Under normal operating conditions, the relative position of the rotor axis and the resultant magnetic field axis is fixed. The angle between the two is known as the power angle, torque angle, or rotor angle. This angular relationship is critical for maintaining synchronism across the grid. The swing equation is defined as a non-linear second order differential equation that describes the swing of the rotor of a synchronous machine. It captures the relative motion that occurs when the system is subjected to a disturbance. During any disturbance, the rotor decelerates or accelerates with respect to the synchronously rotating air gap magnetomotive force, creating relative motion. The equation describing this relative motion is known as the swing equation. It links the mechanical and electrical dynamics of the rotor, providing the basis for transient stability analysis. The power exchange between the mechanical rotor and the electrical grid due to the rotor swing is called Inertial response. This inertial response is a key characteristic of synchronous generators, allowing them to absorb or release kinetic energy to buffer sudden changes in power balance. The non-linear nature of the equation arises from the sinusoidal relationship between the power angle and the electrical power output. As the rotor swings, the torque angle changes, affecting the electromagnetic torque acting on the rotor. The differential equation governs how these torques interact with the mechanical torque and the rotor's moment of inertia. Understanding the swing equation is essential for engineers and researchers analyzing the stability of power systems. It explains how synchronous machines maintain their speed and phase relationship with the grid during transient events. The model is applied to single machines and multi-machine systems to predict whether the generators will remain in synchronism or lose stability. The equation does not account for all electrical transients but focuses on the mechanical swing of the rotor. This simplification allows for efficient analysis of the first few seconds following a disturbance. The inertial response described by the equation is a primary defense against frequency deviations in the grid. Without this response, small imbalances between generation and load could lead to rapid frequency changes. The swing equation thus provides the theoretical foundation for understanding how kinetic energy stored in rotating masses stabilizes the electrical network. It remains a central concept in power system dynamics and control theory. The model is used to design control strategies and assess the impact of new generation technologies on system inertia.

Derivation of the rotor motion equation

The derivation of the rotor motion equation begins with Newton’s second law for rotational motion. For a synchronous machine, the net torque acting on the rotor determines its angular acceleration. The fundamental relationship states that the accelerating torque is the product of the rotor’s moment of inertia and its angular acceleration. This forms the basis for understanding how mechanical and electrical forces interact to maintain synchronism or cause deviation during disturbances.

Defining Torques and Inertia

Let J represent the moment of inertia of the rotor, measured in kg⋅m2. The mechanical torque, denoted as Tm, is the torque supplied by the prime mover (such as a steam or hydro turbine). Mathematically, this is expressed as:

Ta=Tm−Te

According to Newton’s second law, this accelerating torque equals the product of the moment of inertia and the angular acceleration of the rotor. If ωr is the angular velocity of the rotor in radians per second, the angular acceleration is the time derivative of ωr. Thus, the equation of motion is:

Jdt2d2θm=Tm−Te

where θm is the angular position of the rotor in mechanical radians. This second-order differential equation describes the dynamic behavior of the rotor. The term Jdt2d2θm represents the inertial response of the machine, linking the mechanical energy stored in the rotating mass to the electrical power exchange with the grid. Any imbalance between Tm and Te causes the rotor to accelerate or decelerate, altering the power angle relative to the synchronously rotating magnetic field.

How is the power angle defined and calculated?

The power angle, also referred to as the torque angle or rotor angle, is a fundamental parameter in the analysis of synchronous machine dynamics within power systems. It is defined as the angular displacement between the rotor axis of the synchronous machine and the axis of the resultant magnetic field in the air gap. Under normal, steady-state operating conditions, the synchronous machines in a power system operate synchronously, meaning they rotate at the same electrical speed. In this state, the relative position between the rotor axis and the resultant magnetic field axis remains fixed, establishing a constant power angle that corresponds to the steady power output of the machine.

When a disturbance occurs in the power system, such as a sudden change in mechanical input or electrical load, the synchronism between the rotor and the rotating magnetic field is temporarily perturbed. The rotor may accelerate or decelerate relative to the synchronously rotating air gap magnetomotive force. This relative motion causes the power angle to vary dynamically. The swing equation, a non-linear second-order differential equation, mathematically describes this relative motion and the subsequent swing of the rotor. The equation captures the interplay between the mechanical torque applied to the rotor and the electrical torque produced by the interaction with the magnetic field.

The variation in the power angle is directly linked to the inertial response of the synchronous machine. As the rotor swings, energy is exchanged between the mechanical rotor and the electrical grid. This exchange of power, driven by the change in the rotor angle, constitutes the inertial response, which is critical for maintaining frequency stability during transient events. The power angle thus serves as a key indicator of the stability and dynamic behavior of the synchronous machine within the broader power system network.

What is the role of the inertia constant H?

The inertia constant H is a fundamental parameter in the swing equation, quantifying the rotational inertia of a synchronous machine relative to its electrical rating. It is defined as the ratio of the total stored kinetic energy in the rotor at synchronous speed to the machine's apparent power rating in megavolt-amperes (MVA). This constant allows engineers to normalize the inertial response of different machine types, facilitating comparison and per-unit system calculations across the power grid.

Definition and Calculation

The inertia constant H is calculated using the following relationship:

H = Ek / Sbase

Where Ek is the stored kinetic energy in megajoules (MJ) and Sbase is the machine rating in MVA. The resulting unit for H is typically seconds, representing the time the rotor can sustain the rated output power using only its stored kinetic energy, assuming constant speed. For a typical synchronous generator, H values range from 2 to 8 seconds, depending on the rotor design and cooling method.

Role in Per-Unit System Calculations

In per-unit system analysis, H simplifies the swing equation by normalizing mechanical and electrical powers. The swing equation in per-unit form is:

(2H / ωs) · (d2δ / dt2) = Pm - Pe

Here, ωs is the synchronous angular velocity, δ is the rotor angle, Pm is the mechanical power input, and Pe is the electrical power output. The constant H directly scales the acceleration of the rotor angle in response to power imbalances. A higher H value indicates greater stored kinetic energy relative to the machine's rating, resulting in a slower change in rotor speed during disturbances. This property is critical for assessing the inertial response of the power system, which helps stabilize frequency deviations following sudden changes in load or generation.

Impact on System Stability

The inertia constant H plays a crucial role in transient stability analysis. Machines with higher H values provide more inertia, which helps dampen oscillations and maintain synchronism during faults. In modern power systems with increasing penetration of inverter-based resources, the effective H of the grid may decrease, affecting the overall inertial response. Understanding and accurately modeling H is essential for predicting system behavior and designing control strategies to maintain stability.

Worked examples

Basic Per-Unit Conversion

The swing equation is often expressed in per-unit (p.u.) values to simplify calculations across different machine ratings. The fundamental relationship between mechanical power, electrical power, and acceleration power is defined by the variable Pa, where Pa = Pm - Pe (per grounding). In this context, Pm represents the mechanical input power and Pe represents the electrical output power. The difference creates the accelerating power that drives the rotor angle change.

Numerical Example: Calculating Accelerating Power

Consider a synchronous generator operating under normal conditions where the mechanical power input Pm is 1.0 p.u. and the electrical power output Pe is 1.0 p.u. Under these steady-state conditions, the accelerating power Pa is calculated as 1.0 - 1.0 = 0 p.u. This indicates that the rotor speed is constant, and the relative position of the rotor axis and the resultant magnetic field axis is fixed, as described in the grounding. The angle between these axes is known as the power angle or rotor angle.

If a disturbance occurs, such as a sudden increase in load, the electrical power Pe might increase to 1.2 p.u. while the mechanical power Pm remains at 1.0 p.u. The new accelerating power Pa becomes 1.0 - 1.2 = -0.2 p.u. This negative value indicates that the rotor is decelerating with respect to the synchronously rotating air gap magnetomotive force. The equation describing this relative motion is the swing equation, a non-linear second order differential equation.

Inertial Response and Torque

The power exchange between the mechanical rotor and the electrical grid due to the rotor swing is called Inertial response. This response is driven by the torque difference. The mechanical torque Tm and electrical torque Te are directly related to the powers Pm and Pe. The inertia constant, often related to the moment of inertia J, determines how quickly the rotor accelerates or decelerates for a given power imbalance. The swing equation models this dynamic behavior, capturing the non-linear nature of the synchronous machine's operation during disturbances. The rotor angle changes as the machine accelerates or decelerates, altering the power angle and thus the power exchange with the grid.

Applications in power system stability analysis

The swing equation serves as the fundamental mathematical model for evaluating transient stability in power systems. As a non-linear second order differential equation, it describes the relative motion of the synchronous machine rotor with respect to the synchronously rotating air gap magnetomotive force during disturbances. Engineers utilize this equation to determine whether the rotor will maintain synchronism or fall out of step following a sudden change in mechanical or electrical torque.

Transient Stability Analysis

In transient stability studies, the swing equation models the dynamic behavior of the rotor angle, also known as the power angle, torque angle, or rotor angle. Under normal operating conditions, the relative position of the rotor axis and the resultant magnetic field axis remains fixed. However, during a disturbance, the rotor decelerates or accelerates, creating relative motion. The equation quantifies this swing, allowing analysts to predict the rotor's trajectory. If the rotor angle diverges excessively, the synchronous machines may lose synchrony, leading to potential system-wide instability or cascading failures.

Inertial Response and Rotor Dynamics

The power exchange between the mechanical rotor and the electrical grid due to the rotor swing is defined as inertial response. This phenomenon is critical for frequency stability. When the swing equation is applied, it reveals how kinetic energy stored in the rotating masses of synchronous machines is released or absorbed to counteract imbalances between mechanical input and electrical output. The acceleration or deceleration of the rotor directly influences the system frequency. By solving the swing equation, engineers can assess the magnitude of the inertial response and its effectiveness in mitigating frequency deviations during the initial seconds following a disturbance.

Modeling Disturbances

Disturbances such as short circuits, generator tripping, or sudden load changes alter the balance of forces acting on the rotor. The swing equation captures these dynamics by relating the rotor's angular acceleration to the difference between mechanical torque and electrical torque. This relationship enables the simulation of various operating conditions, ensuring that the power system consists of a number of synchronous machines operating synchronously under all operating conditions. Accurate modeling using this equation is essential for designing control strategies and protection schemes that maintain grid reliability.

What distinguishes linear and non-linear solutions?

The swing equation is fundamentally a non-linear second-order differential equation. This non-linearity arises because the electrical power output of a synchronous machine depends on the sine of the rotor angle (power angle, torque angle, or rotor angle). The relationship is expressed as P_e = (E V / X) sin(δ), where E is the internal voltage, V is the terminal voltage, X is the reactance, and δ is the rotor angle. Consequently, the restoring torque acting on the rotor is not directly proportional to the displacement angle δ, but rather to sin(δ). This sinusoidal relationship means that for large disturbances, the rotor's motion is complex and the system's stability cannot be determined by simple linear approximations.

Linearization for Small Angles

For small disturbances, engineers often employ a linearized approximation to simplify stability analysis. In this approach, the rotor angle δ is assumed to be small enough that sin(δ) ≈ δ. This transforms the non-linear differential equation into a linear second-order differential equation, analogous to a simple harmonic oscillator. The linearized model is highly effective for analyzing steady-state stability and small-signal stability, where the rotor oscillates with small amplitude around the equilibrium point. However, this approximation loses accuracy during large transients, such as a three-phase fault or a sudden load change, where the rotor angle may swing significantly, making the difference between sin(δ) and δ substantial.

Numerical Solution Methods

Because the exact analytical solution to the non-linear swing equation is complex, numerical methods are frequently used to determine the rotor's trajectory over time. One of the most common and effective techniques is the fourth-order Runge-Kutta algorithm. This method provides a high degree of accuracy by estimating the slope of the solution curve at multiple points within each time step. By iteratively applying the Runge-Kutta formulas, engineers can simulate the rotor angle and angular velocity as functions of time, capturing the non-linear dynamics of the inertial response. This numerical approach allows for a detailed analysis of transient stability, helping to determine whether the synchronous machine will remain synchronized with the grid or lose synchronism following a significant disturbance.

References

#Power System Stability #synchronous machine #rotor dynamics #inertial response #power angle