2.4 Resonance, Frequency Response & Harmonic Distortion

Resonance, Frequency Response & Harmonic Distortion

In AC power engineering, resonance and harmonic distortion are critical phenomena that directly impact equipment reliability, system losses, and power quality. While resonance principles are utilized beneficially in tuned communication filters and induction heating circuits, unintentional resonance in power distribution networks can cause catastrophic dielectric breakdown, capacitor bank explosions, and transformer overheating.

Modern power distribution systems are dominated by non-linear solid-state power electronic loads—variable frequency drives (VFDs), uninterruptible power supplies (UPS), battery chargers, arc furnaces, and LED drivers—that inject harmonic currents into the power grid. Understanding resonance conditions, the IEEE 519 standard, harmonic sequence properties, and mitigation technologies is vital for the NCEES PE Power examination.

---

1. Series RLC Resonance Mechanics

Consider a series RLC circuit excited by a sinusoidal voltage source $\mathbf{V}_s$ of variable angular frequency $\omega$:

                 +----[ R ]----[ L ]----[ C ]----+ 
                 |   v_R(t)   v_L(t)   v_C(t)    |
               ( + )                             |
              V_s(t)                             |
               ( - )                             |
                 |                               |
                 +-------------------------------+ 

The total series impedance is:

$\mathbf{Z}(j\omega) = R + j\left(\omega L - \frac{1}{\omega C}\right)$

Resonance occurs at the specific angular frequency $\omega_0$ where the imaginary component of impedance vanishes ($X_L = X_C$):

$\omega_0 L = \frac{1}{\omega_0 C} \implies \omega_0^2 = \frac{1}{LC} \implies \omega_0 = \frac{1}{\sqrt{LC}}\text{ rad/s}$

$f_0 = \frac{1}{2\pi\sqrt{LC}}\text{ Hz}$

+-----------------------------------------------------------------------------------------+
|                        SERIES RLC RESONANCE CHARACTERISTICS                             |
|                                                                                         |
|   1. IMPEDANCE IS MINIMIZED:                                                            |
|      - Z(omega_0) = R + j(0) = R (Purely resistive, minimum possible impedance).        |
|                                                                                         |
|   2. CURRENT IS MAXIMIZED:                                                              |
|      - I_0 = V_s / R (In phase with source voltage, Power Factor = 1.0 Unity).          |
|                                                                                         |
|   3. REACTIVE VOLTAGE MAGNIFICATION:                                                    |
|      - V_L(omega_0) = I_0 * (omega_0 * L) = (V_s / R) * (omega_0 * L) = Q_s * V_s       |
|      - V_C(omega_0) = I_0 * (1 / (omega_0 * C)) = Q_s * V_s                             |
|      - Inductor and capacitor voltages are 180 degrees out of phase, canceling at the   |
|        terminals, but individual component voltages can be 10x-50x the applied V_s!     |
+-----------------------------------------------------------------------------------------+

Series Quality Factor ($Q_s$) and Half-Power Bandwidth ($BW$)

The Quality Factor $Q$ measures the sharpness of resonance and the ratio of stored energy to dissipated energy per cycle:

$Q_s = \frac{\omega_0 L}{R} = \frac{1}{\omega_0 R C} = \frac{1}{R}\sqrt{\frac{L}{C}}$

The Half-Power Frequencies (cutoff frequencies $\omega_1, \omega_2$ where current drops to $I_{max}/\sqrt{2} = 0.707 I_{max}$ and power drops to $P_{max}/2$) define the Bandwidth ($BW$):

$BW = \Delta \omega = \omega_2 - \omega_1 = \frac{R}{L} = \frac{\omega_0}{Q_s}\text{ rad/s}$

$BW = \Delta f = f_2 - f_1 = \frac{f_0}{Q_s}\text{ Hz}$

For high-Q circuits ($Q_s \ge 10$), the response is symmetric around $f_0$:

$f_1 \approx f_0 - \frac{BW}{2}, \qquad f_2 \approx f_0 + \frac{BW}{2}$

---

2. Parallel RLC Resonance Mechanics

In a parallel RLC circuit excited by an AC current source $\mathbf{I}_s$, the total input admittance is:

$\mathbf{Y}(j\omega) = \frac{1}{R} + j\left(\omega C - \frac{1}{\omega L}\right)$

                             +-------+-------+-------+
                             |       |       |       |
                           ( | )   [ R ]   [ L ]   [ C ]
                          I_s(t)     |       |       |
                             |       |       |       |
                             +-------+-------+-------+

Resonance occurs when the susceptance equals zero ($B_C = B_L$):

$\omega_0 = \frac{1}{\sqrt{LC}}\text{ rad/s}, \qquad f_0 = \frac{1}{2\pi\sqrt{LC}}\text{ Hz}$

+-----------------------------------------------------------------------------------------+
|                       PARALLEL RLC RESONANCE CHARACTERISTICS                            |
|                                                                                         |
|   1. ADMITTANCE IS MINIMIZED / IMPEDANCE IS MAXIMIZED:                                  |
|      - Y(omega_0) = 1/R -> Z(omega_0) = R (Purely resistive, maximum possible impedance)|
|                                                                                         |
|   2. LINE CURRENT IS MINIMIZED:                                                         |
|      - Total current drawn from the line reaches a minimum: I_line = V / R              |
|                                                                                         |
|   3. CIRCULATING TANK CURRENT MAGNIFICATION:                                            |
|      - Large circulating currents flow between L and C: I_L = I_C = Q_p * I_line        |
|                                                                                         |
|   4. PARALLEL QUALITY FACTOR (Q_p):                                                     |
|      - Q_p = R / (omega_0 * L) = omega_0 * R * C = R * sqrt(C / L)                      |
|      - Bandwidth: BW = omega_0 / Q_p = 1 / (R * C) rad/s                                |
+-----------------------------------------------------------------------------------------+

---

3. Power System Harmonics & Non-Linear Loads

A harmonic is a sinusoidal component of a periodic wave having a frequency that is an integral integer multiple of the fundamental power frequency:

$f_h = h \times f_1$

Where $h$ is the integer harmonic order (e.g., for $f_1 = 60\text{ Hz}$: $h=3 \rightarrow 180\text{ Hz}$, $h=5 \rightarrow 300\text{ Hz}$, $h=7 \rightarrow 420\text{ Hz}$).

+-----------------------------------------------------------------------------------------+
|                        HARMONIC CLASSIFICATION & SEQUENCE PROPERTIES                    |
|                                                                                         |
|   HARMONIC ORDER (h)        SEQUENCE COMPONENT       IMPACT ON THREE-PHASE POWER SYSTEM |
|   h = 1, 4, 7, 10, 13...   Positive Sequence (+)    Rotates forward (motor torque)      |
|   h = 2, 5, 8, 11, 14...   Negative Sequence (-)    Rotates backward (braking/heating)  |
|   h = 3, 9, 15, 21...      Zero Sequence (Triplen)  IN PHASE in all three conductors!   |
+-----------------------------------------------------------------------------------------+

The Triplen Harmonic Neutral Hazard

Triplen harmonics ($h = 3, 9, 15, \dots$) are odd multiples of the third harmonic. In a balanced 3-phase, 4-wire wye system:

• Fundamental currents ($I_a, I_b, I_c$) are displaced by $120^\circ$ and sum to zero in the neutral ($I_N = 0$).
• Triplen harmonic currents are displaced by $3 \times 120^\circ = 360^\circ = 0^\circ$. They are in phase with each other in all three phases!
• The triplen currents add arithmetically in the neutral conductor:

$I_{N(3rd)} = I_{a,3} + I_{b,3} + I_{c,3} = 3 I_{3}$

In facilities with heavy concentrations of single-phase switched-mode power supplies (PCs, servers, electronic ballasts), the neutral current can exceed the phase conductor current ($I_N > 173% I_{phase}$), causing severe neutral conductor overheating, transformer saturation, and building fires.

---

4. Total Harmonic Distortion (THD) & Power Factor with Harmonics

Total Harmonic Distortion (THD) quantifies the distortion of a voltage or current waveform relative to the fundamental component.

+-----------------------------------------------------------------------------------------+
|                        TOTAL HARMONIC DISTORTION (THD) FORMULAS                         |
|                                                                                         |
|   VOLTAGE THD (THD_V):                                                                  |
|   THD_V = sqrt( sum_{h=2}^infinity V_h^2 ) / V_1 * 100% = sqrt( V_rms^2 - V_1^2 ) / V_1 |
|                                                                                         |
|   CURRENT THD (THD_I):                                                                  |
|   THD_I = sqrt( sum_{h=2}^infinity I_h^2 ) / I_1 * 100% = sqrt( I_rms^2 - I_1^2 ) / I_1 |
|                                                                                         |
|   TRUE RMS QUANTITY:                                                                    |
|   I_rms = sqrt( I_1^2 + I_2^2 + I_3^2 + ... + I_h^2 ) = I_1 * sqrt( 1 + THD_I^2 )       |
+-----------------------------------------------------------------------------------------+

True Power Factor vs. Displacement Power Factor

When harmonic distortion is present, the standard displacement power factor is insufficient:

$\text{Displacement Power Factor: } \text{PF}_{disp} = \cos\theta_1$

$\text{Distortion Power Factor: } \text{PF}_{dist} = \frac{I_1}{I_{rms}} = \frac{1}{\sqrt{1 + \text{THD}_I^2}}$

$\text{True Power Factor: } \text{PF}_{true} = \frac{P}{S} = \text{PF}_{disp} \times \text{PF}_{dist} = \cos\theta_1 \times \frac{1}{\sqrt{1 + \text{THD}_I^2}}$

[!CAUTION] Capacitor Overheating from Harmonics: Capacitive reactance decreases with frequency ($X_C(h) = \frac{1}{h \omega C}$). As a result, power factor correction capacitors act as low-impedance sinks for high-frequency harmonic currents, leading to severe thermal overloads and potential dielectric rupture.

---

5. IEEE 519 Standard & Harmonic Mitigation Technologies

IEEE Standard 519 (Standard for Harmonic Control in Electric Power Systems) establishes recommended distortion limits at the Point of Common Coupling (PCC)—the interface between the utility and the industrial customer.

+-----------------------------------------------------------------------------------------+
|                        IEEE 519 VOLTAGE DISTORTION LIMITS AT PCC                        |
|                                                                                         |
|   Bus Voltage at PCC (V)            Individual Harmonic (%)    Total Harmonic Dist (THD)|
|   V <= 1.0 kV                       5.0%                       8.0%                     |
|   1.0 kV < V <= 69 kV               3.0%                       5.0%                     |
|   69 kV < V <= 161 kV               1.5%                       2.5%                     |
|   V > 161 kV                        1.0%                       1.5%                     |
+-----------------------------------------------------------------------------------------+

Harmonic Mitigation Strategies

1. Detuned Capacitor Banks (Anti-Resonance Reactors): Series inductors added to power factor capacitor banks, tuned to a sub-harmonic frequency (typically $4.2\text{nd}$ or $4.7\text{th}$ harmonic, below the lowest 5th harmonic). This prevents harmonic resonance while still supplying 60 Hz reactive power.
2. Passive Tuned LC Notch Filters: Series LC branches tuned precisely to a dominant offending harmonic ($h=5$ at $300\text{ Hz}$ or $h=7$ at $420\text{ Hz}$) to trap and divert harmonic current to ground.
3. Active Harmonic Filters (AHF): High-speed IGBT-based inverters that measure load harmonic currents in real time and inject equal-and-opposite anti-phase harmonic currents, neutralizing harmonic distortion at the source.
4. Multi-Pulse Rectifiers: Utilizing 12-pulse or 18-pulse transformer configurations with phase-shifting windings (e.g., $30^\circ$ delta-wye shifts) to cancel 5th, 7th, 11th, and 13th harmonics via magnetic cancellation ($h = kp \pm 1$).

---

6. Comprehensive Worked Mathematical Examples

Example 1: Series RLC Resonant Circuit Calculations

Problem: A series RLC circuit connected across a $120\text{ V}_{rms}$ AC source has $R = 5.0\ \Omega$, $L = 20.0\text{ mH}$, and $C = 3.166\ \mu\text{F}$.

1. Find the resonant frequency $f_0$, quality factor $Q_s$, and half-power bandwidth $BW$.
2. Calculate the circuit current $I_0$ and the voltages across the inductor ($V_L$) and capacitor ($V_C$) at resonance.

Step 1: Calculate Resonant Frequency ($f_0$)

$\omega_0 = \frac{1}{\sqrt{LC}} = \frac{1}{\sqrt{(0.020\text{ H})(3.166 \times 10^{-6}\text{ F})}} = \frac{1}{\sqrt{6.332 \times 10^{-8}}} = \frac{1}{2.5163 \times 10^{-4}} = 3974.0\text{ rad/s}$

$f_0 = \frac{\omega_0}{2\pi} = \frac{3974.0}{2\pi} = 632.48\text{ Hz}$

Step 2: Calculate Quality Factor ($Q_s$) and Bandwidth ($BW$)

$Q_s = \frac{\omega_0 L}{R} = \frac{3974.0 \times 0.020}{5.0} = \frac{79.48}{5.0} = 15.90$

$BW = \frac{f_0}{Q_s} = \frac{632.48\text{ Hz}}{15.90} = 39.78\text{ Hz}$

Step 3: Calculate Current and Component Voltages at Resonance At resonance, $\mathbf{Z} = R = 5.0\ \Omega$:

$I_0 = \frac{V_s}{R} = \frac{120.0\text{ V}}{5.0\ \Omega} = 24.0\text{ A}_{rms}$

$X_L = \omega_0 L = 3974.0 \times 0.020 = 79.48\ \Omega$

$V_L = I_0 X_L = 24.0\text{ A} \times 79.48\ \Omega = 1907.5\text{ V}_{rms} = Q_s \times V_s = 15.90 \times 120.0 = 1908.0\text{ V}_{rms}$

$V_C = I_0 X_C = Q_s \times V_s = 1908.0\text{ V}_{rms}$

Note: Although the applied line voltage is only $120\text{ V}$, the internal reactive component voltages exceed $1.9\text{ kV}$ due to resonant magnification!

---

Example 2: Non-Linear Load Harmonics, True RMS & Neutral Current

Problem: A $480\text{ V}$ 3-phase 4-wire balanced non-linear computer center load draws phase current with the following spectral distribution:

• Fundamental ($60\text{ Hz}$): $I_1 = 200.0\text{ A}{rms}$, $\text{PF}{disp} = 0.90\text{ lagging}$
• 3rd Harmonic ($180\text{ Hz}$): $I_3 = 80.0\text{ A}_{rms}$
• 5th Harmonic ($300\text{ Hz}$): $I_5 = 50.0\text{ A}_{rms}$
• 7th Harmonic ($420\text{ Hz}$): $I_7 = 30.0\text{ A}{rms}$ Calculate the True RMS phase current $I{rms}$, current total harmonic distortion $\text{THD}I$, true power factor $\text{PF}{true}$, and the total neutral conductor current $I_N$.

Step 1: Calculate True RMS Phase Current ($I_{rms}$)

$I_{rms} = \sqrt{I_1^2 + I_3^2 + I_5^2 + I_7^2} = \sqrt{200^2 + 80^2 + 50^2 + 30^2} = \sqrt{40000 + 6400 + 2500 + 900} = \sqrt{49800} = 223.16\text{ A}_{rms}$

Step 2: Calculate Current THD ($\text{THD}_I$)

$\text{THD}_I = \frac{\sqrt{I_3^2 + I_5^2 + I_7^2}}{I_1} \times 100\% = \frac{\sqrt{6400 + 2500 + 900}}{200.0} \times 100\% = \frac{\sqrt{9800}}{200.0} \times 100\% = \frac{98.995}{200.0} \times 100\% = 49.50\%$

Step 3: Calculate True Power Factor ($\text{PF}_{true}$)

$\text{PF}_{dist} = \frac{I_1}{I_{rms}} = \frac{200.0}{223.16} = 0.8962$

$\text{PF}_{true} = \text{PF}_{disp} \times \text{PF}_{dist} = 0.90 \times 0.8962 = 0.8066 \approx 0.807\text{ Lagging}$

Step 4: Calculate Neutral Conductor Current ($I_N$) In a balanced system, fundamental ($h=1$), 5th ($h=5$, negative sequence), and 7th ($h=7$, positive sequence) harmonic currents sum to zero at the neutral point. Only zero-sequence triplen harmonics ($h=3$) add arithmetically:

$I_N = 3 \times I_3 = 3 \times 80.0\text{ A} = 240.0\text{ A}_{rms}$

The neutral current ($240.0\text{ A}$) exceeds the true RMS phase current ($223.16\text{ A}$), illustrating the severe danger of triplen harmonics in 4-wire distribution systems.

---

7. Common PE Exam Traps & Pitfalls

1. Inverting Series and Parallel Quality Factor ($Q$) Formulas: For series RLC, $Q_s = \frac{\omega_0 L}{R} = \frac{1}{R}\sqrt{L/C}$ (smaller $R$ increases $Q$). For parallel RLC, $Q_p = \frac{R}{\omega_0 L} = R\sqrt{C/L}$ (larger $R$ increases $Q$).
2. Confusing Displacement and True Power Factor: Displacement PF is $\cos\theta_1$. If harmonics are present, true PF is always lower: $\text{PF}{true} = \text{PF}{disp} / \sqrt{1 + \text{THD}_I^2}$. Standard analog wattmeters measure true power, but simple V-I-PF meters without true-RMS capability will give erroneous readings.
3. Assuming Balanced 3-Phase Neutrals Carry Zero Current: While $I_N = 0$ for balanced linear sinusoidal loads, non-linear triplen harmonics add constructively in the neutral ($I_N \approx 3 I_3$), requiring double-sized ($200%$) neutral conductors per NEC standards.