12.3 Power Quality, Harmonics, THD & IEEE 519 Compliance

12.3 Power Quality, Harmonics, THD & IEEE 519 Compliance

Executive Overview: Modern power systems are heavily populated by non-linear electronic loads—including Variable Frequency Drives (VFDs), uninterruptible power supplies (UPS), solar inverters, and arc furnaces—that inject harmonic currents and distort system voltages. On the PE Power exam, power quality questions focus on distinguishing disturbance classes (sags, swells, transients, flicker, unbalance), evaluating harmonic sequence behavior (especially triplen harmonics in neutrals), computing $THD$, $TDD$, and True Power Factor, sizing K-Factor transformers, and applying IEEE 519-2022 Point of Common Coupling (PCC) limits.

---

1. Power Quality Disturbance Classifications (IEEE 1159)

Voltage Unbalance Formulation (NEMA vs. True Symmetrical)

Per NEMA MG-1, the Phase Voltage Unbalance Rate ($PVUR\%$) is defined as:

$$PVUR\% = \frac{\text{Maximum Deviation from Average Phase Voltage}}{\text{Average Phase Voltage}} \times 100\%$$

---

2. Harmonic Fundamentals & Sequence Properties

Fourier Series & Non-Linear Loads

Any periodic distorted waveform $f(t)$ with fundamental frequency $\omega_1 = 2\pi f_0$ ($60\text{ Hz}$) is expressed via Fourier expansion:

$$f(t) = a_0 + \sum_{h=1}^\infty \left[ a_h \cos(h \omega_1 t) + b_h \sin(h \omega_1 t) \right] = I_0 + \sum_{h=1}^\infty \sqrt{2} I_h \sin(h \omega_1 t + \theta_h)$$

Characteristic Harmonics of Static Converters

For a line-commutated $p$-pulse converter (where $p=6$ for standard 6-pulse bridge, $p=12$ for dual-bridge with $30^\circ$ phase-shifting transformer), the characteristic ac harmonic orders generated are:

$$h = k p \pm 1 \quad \text{where } k = 1, 2, 3, \dots$$

• 6-Pulse Converter: $h = 5, 7, 11, 13, 17, 19, 23, 25, \dots$
• 12-Pulse Converter: $h = 11, 13, 23, 25, 35, 37, \dots$ (5th and 7th harmonics are canceled by the $30^\circ$ delta-wye phase shift).
• Theoretical Current Magnitude: In ideal square-wave switching, $I_h = \frac{I_1}{h}$.

Sequence Breakdown of Harmonics in Balanced Three-Phase Systems:

 Harmonic Order (h)  | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15
---------------------+---+---+---+---+---+---+---+---+---+----+----+----+----+----+----
 Sequence Component | + | - | 0 | + | - | 0 | + | - | 0 | +  | -  | 0  | +  | -  | 0 

1. Positive Sequence ($h = 3k + 1 = 1, 4, 7, 10, 13...$): Rotate in the same forward direction as fundamental flux ($a-b-c$).
2. Negative Sequence ($h = 3k + 2 = 2, 5, 8, 11, 14...$): Rotate in reverse direction ($a-c-b$), inducing double-slip rotor heating in rotating machines.
3. Zero Sequence / Triplen Harmonics ($h = 3k = 3, 9, 15, 21...$): All three phase currents are exactly in phase with each other ($\theta_a = \theta_b = \theta_c$).

$$\mathbf{I}_{N} = \mathbf{I}_a + \mathbf{I}_b + \mathbf{I}_c = 3 I_3 + 3 I_9 + 3 I_{15} + \dots$$

---

3. Distortion Metrics: THD, TDD, True RMS & Power Factor

Total Harmonic Distortion (THD)

Quantifies the total harmonic content relative to fundamental magnitude:

$$THD_V = \frac{\sqrt{\sum_{h=2}^\infty V_h^2}}{V_1} \times 100\%, \qquad THD_I = \frac{\sqrt{\sum_{h=2}^\infty I_h^2}}{I_1} \times 100\%$$

Total Demand Distortion (TDD)

Under light-load conditions, fundamental current $I_1$ is small, causing $THD_I$ to appear alarmingly high (e.g., $80\%$) even when harmonic amps are harmlessly low. IEEE 519 resolves this by defining Total Demand Distortion ($TDD$), which normalizes harmonic currents to the facility's Maximum Demand Load Current ($I_L$) (15-to-30-minute peak demand):

$$TDD = \frac{\sqrt{\sum_{h=2}^\infty I_h^2}}{I_L} \times 100\% \quad [\%]$$

Relationship between $THD_I$ and $TDD$:

$$TDD = THD_I \times \left( \frac{I_1}{I_L} \right)$$

True RMS and True Power Factor

In the presence of non-sinusoidal currents and voltages:

$$I_{rms} = \sqrt{I_1^2 + \sum_{h=2}^\infty I_h^2} = I_1 \sqrt{1 + \left( \frac{THD_I}{100} \right)^2}$$ $$PF_{true} = \frac{P}{S} = \frac{V_1 I_1 \cos\theta_1}{V_{rms} I_{rms}} = \left( \frac{V_1}{V_{rms}} \right) \left( \frac{I_1}{I_{rms}} \right) \cos\theta_1 \approx PF_{dist} \times PF_{disp}$$ $$PF_{true} = \underbrace{\frac{1}{\sqrt{1 + \left(\frac{THD_I}{100}\right)^2}}}_{\text{Distortion Power Factor } (PF_{dist})} \times \underbrace{\cos\theta_1}_{\text{Displacement Power Factor } (PF_{disp})}$$

Transformer K-Factor Rating

Harmonic currents cause severe stray eddy-current losses in transformer windings, which scale with the square of harmonic frequency ($P_{ec} \propto I_h^2 h^2$). The UL/IEEE K-Factor quantifies a transformer's ability to withstand non-sinusoidal heating without exceeding thermal insulation ratings:

$$K = \frac{\sum_{h=1}^\infty I_h^2 \cdot h^2}{\sum_{h=1}^\infty I_h^2} = \sum_{h=1}^\infty \left[ \frac{I_h}{I_{rms}} \right]^2 \cdot h^2$$

Standard Commercial K-Factor Ratings: $K-1$ (standard linear load), $K-4$ (HID lighting, light VFDs), $K-13$ (data centers, multiple VFDs), $K-20$ (dense non-linear IT/mainframe loads), $K-30$ (arc furnaces).

---

4. IEEE 519-2022 Standard Compliance at the PCC

IEEE 519 establishes harmonic distortion limits at the Point of Common Coupling (PCC)—the physical interface where the utility transmission/distribution system connects to the customer facility.

Point of Common Coupling (PCC) Interface Definition:

 Utility Grid =======[ Utility Substation ]======o (PCC - Metering Point)
                                                 |
                                       [ Customer Main Switchgear ]
                                                 |
                             +-------------------+-------------------+
                             |                                       |
                      [ Linear Loads ]                        [ Non-Linear VFDs ]

IEEE 519-2022 Voltage Distortion Limits at PCC

IEEE 519-2022 Current Distortion Limits for Systems Rated $120\text{ V} - 69\text{ kV}$

Current limits are indexed by the Short Circuit Ratio ($SCR$) at the PCC: $SCR = \frac{I_{sc}}{I_L}$, where $I_{sc}$ is maximum available short-circuit current at PCC, and $I_L$ is maximum fundamental demand load current.

---

5. Comprehensive Worked Calculation: THD, TDD, K-Factor & IEEE 519

Problem Statement

An industrial manufacturing facility connects to a utility $13.8\text{ kV}$ distribution feeder at the PCC. The available short-circuit capacity at the PCC is $150\text{ MVA}$. The facility has a contract maximum 15-minute demand load current $I_L = 350\text{ A}$. During an operating cycle, power analyzer measurements indicate fundamental current $I_1 = 280\text{ A}$ at $0.88$ displacement power factor lagging, with the following harmonic current spectrum:

• $I_5 = 42.0\text{ A}$
• $I_7 = 28.0\text{ A}$
• $I_{11} = 14.0\text{ A}$
• $I_{13} = 9.8\text{ A}$
• All other higher harmonics are negligible.

Calculate:

1. Current Total Harmonic Distortion ($THD_I\%$) and True RMS current ($I_{rms}$).
2. Total Demand Distortion ($TDD\%$).
3. True Power Factor ($PF_{true}$).
4. Transformer K-Factor.
5. Short Circuit Ratio ($SCR = I_{sc}/I_L$) and verify whether the facility complies with IEEE 519-2022 current limits.

============================== STEP-BY-STEP SOLUTION ==============================

Step 1: Compute Harmonic RMS Current Sum and THD_I
  Sum of squared harmonic currents:
  Sum(I_h^2) = 42.0^2 + 28.0^2 + 14.0^2 + 9.8^2
             = 1,764.0 + 784.0 + 196.0 + 96.04 = 2,840.04 A^2
  
  Total harmonic current magnitude:
  I_harm_rms = sqrt(2,840.04) = 53.292 A
  
  Current THD relative to fundamental (I_1 = 280 A):
  THD_I% = (I_harm_rms / I_1) * 100% = (53.292 / 280.0) * 100% = 19.03%
  
  True RMS Current:
  I_rms = sqrt(I_1^2 + I_harm_rms^2) = sqrt(280.0^2 + 2,840.04) 
        = sqrt(78,400 + 2,840.04) = sqrt(81,240.04) = 285.03 A

Step 2: Compute Total Demand Distortion (TDD)
  Contract Maximum Demand Current I_L = 350 A
  TDD% = (I_harm_rms / I_L) * 100% = (53.292 / 350.0) * 100% = 15.23%
  
  Check via conversion formula:
  TDD% = THD_I% * (I_1 / I_L) = 19.033% * (280 / 350) = 15.23% (Exact match)

Step 3: Compute True Power Factor
  Displacement Power Factor PF_disp = 0.88
  Distortion Power Factor:
  PF_dist = 1 / sqrt(1 + (THD_I/100)^2) = 1 / sqrt(1 + 0.19033^2) 
          = 1 / sqrt(1 + 0.03622) = 1 / 1.01795 = 0.98237
  
  True Power Factor:
  PF_true = PF_disp * PF_dist = 0.88 * 0.98237 = 0.8645 lagging

Step 4: Compute Transformer K-Factor
  Calculate Sum(I_h^2 * h^2) including fundamental (h=1):
  h=1:   280.0^2 * 1^2   = 78,400.0 * 1   = 78,400.0
  h=5:   42.0^2  * 5^2   = 1,764.0  * 25  = 44,100.0
  h=7:   28.0^2  * 7^2   = 784.0    * 49  = 38,416.0
  h=11:  14.0^2  * 11^2  = 196.0    * 121 = 23,716.0
  h=13:  9.8^2   * 13^2  = 96.04    * 169 = 16,230.76
  --------------------------------------------------
  Sum(I_h^2 * h^2) = 78,400 + 44,100 + 38,416 + 23,716 + 16,230.76 = 200,862.76
  Sum(I_h^2) = I_rms^2 = 81,240.04
  
  K-Factor = 200,862.76 / 81,240.04 = 2.472
  
  *Selection:* A standard K-4 rated transformer (rated for K <= 4.0) is required.

Step 5: Evaluate IEEE 519-2022 PCC Compliance
  Available 3-Phase Short-Circuit Current at 13.8 kV:
  I_sc = S_sc / (sqrt(3) * V_LL) = 150,000,000 / (sqrt(3) * 13,800) = 6,275.5 A
  
  Short Circuit Ratio:
  SCR = I_sc / I_L = 6,275.5 A / 350 A = 17.93
  
  For SCR < 20 at 13.8 kV PCC:
  - IEEE 519 Individual Harmonic Limit for h < 11 is 4.0%
    Calculated 5th harmonic: I_5 / I_L = 42 / 350 = 12.0% (VIOLATION! Max is 4.0%)
    Calculated 7th harmonic: I_7 / I_L = 28 / 350 = 8.0%  (VIOLATION! Max is 4.0%)
  - IEEE 519 Max Total Demand Distortion (TDD) limit is 5.0%
    Calculated TDD = 15.23% (VIOLATION! Exceeds 5.0% limit)
  
  *Conclusion:* The facility FAILS IEEE 519 compliance. A 5th/7th passive harmonic 
  trap filter or an active harmonic filter (AHF) is mandatory.
===================================================================================

---

6. Common Exam Traps & Strategic Pitfalls

• Conflating $THD_I$ with $TDD$: Using fundamental current $I_1$ in the denominator instead of maximum demand current $I_L$ when evaluating IEEE 519 compliance. IEEE 519 strictly limits $TDD$, not $THD_I$.
• The True Power Factor Omission: Assuming $PF_{true} = \cos\theta_1$. If current distortion exists, you must multiply displacement PF by distortion PF ($PF_{dist} = 1/\sqrt{1+THD_I^2}$).
• Triplen Neutral Current Cancellation Error: Assuming neutral current is zero in a balanced 3-phase system with triplens. 3rd harmonics are zero-sequence and add directly in the neutral ($I_N = 3 I_3$).
• Neglecting Fundamental ($h=1$) in K-Factor: Calculating $K = \sum_{h>1} I_h^2 h^2 / I_{rms}^2$. The numerator must include the fundamental component $I_1^2 \times 1^2$.