Harmonic-Informed Load Flow for Solar & EV-Dominated Feeders

🎯 Learning Objectives

✓ Calculate total harmonic distortion (THD) at critical bus nodes using measured or modeled harmonic spectra
✓ Analyze feeder loading margins under combined fundamental + 5th/7th/11th harmonic currents using IEEE 519–2022 compliance checks
✓ Apply harmonic admittance matrix augmentation to modify standard Y-bus for harmonic-informed load flow convergence
✓ Explain how inverter-based solar PV and uncontrolled EV charging alter harmonic impedance profiles and cause parallel resonance
✓ Design harmonic mitigation strategies (e.g., tuned filters, inverter reactive power support) based on load flow sensitivity results

📖 Why This Matters

Modern mining sites increasingly deploy solar microgrids and fast-charging depots for haul trucks—both dominated by power electronics. These devices inject harmonics (especially 5th, 7th, 11th orders) that distort voltage waveforms, overheat transformers, trip protection relays, and reduce available feeder capacity by up to 30%. Conventional load flow ignores this—leading to unsafe oversights during grid planning. Harmonic-informed load flow bridges that gap: it’s not optional anymore—it’s essential for reliability, safety, and regulatory compliance on solar- and EV-dominant mine feeders.

📘 Core Principles

Harmonic-informed load flow builds on three foundational layers: (1) Fundamental-frequency AC load flow (Newton-Raphson or Fast Decoupled), which solves for voltage magnitude/angle at each bus; (2) Harmonic power flow theory, where harmonic currents are treated as independent phasors injected at each nonlinear load node, governed by harmonic admittance matrices; and (3) Coupling mechanisms—primarily through harmonic voltage-dependent active/reactive power injections from inverters (per IEEE 1547-2018 Annex D) and harmonic impedance interactions (Z_h = V_h / I_h). Crucially, harmonic distortion alters effective conductor ampacity (via RMS current heating), transformer K-factor ratings, and capacitor bank resonance behavior—requiring co-simulation or iterative harmonic + fundamental updates. Modern implementations use harmonic domain (HD) or harmonic decoupled (HD-DF) methods to avoid full time-domain simulation overhead while preserving accuracy for dominant low-order harmonics.

📐 Harmonic-Informed Voltage Distortion Limit Check

IEEE 519–2022 defines allowable THD_v at the Point of Common Coupling (PCC). This check must be embedded within load flow convergence criteria—not performed post-hoc. The formula computes weighted harmonic voltage contribution relative to fundamental, using harmonic impedances derived from the augmented Y-bus.

💡 Worked Example

Problem: A mine’s 34.5 kV PCC bus supplies 12 MW solar farm (inverters) and 8 MW EV depot (60× 150 kW CCS chargers). Measured harmonic current spectrum shows I₅ = 12.4 A, I₇ = 7.8 A, I₁₁ = 4.2 A at PCC. System short-circuit MVA = 420 MVA; fundamental current I₁ = 201 A. Calculate THD_v and verify against IEEE 519 Table 10.2 (PCC > 69 kV → limit = 1.0%; but here PCC = 34.5 kV → limit = 1.5%). Assume Z₅ = Z₇ = Z₁₁ = 0.12 Ω (derived from feeder harmonic impedance model).

1. Step 1: Compute harmonic voltages: V₅ = I₅ × Z₅ = 12.4 A × 0.12 Ω = 1.488 V; V₇ = 7.8 × 0.12 = 0.936 V; V₁₁ = 4.2 × 0.12 = 0.504 V.

2. Step 2: Compute fundamental voltage: V₁ ≈ 34.5 kV / √3 = 19.92 kV = 19,920 V.

3. Step 3: Compute THD_v = √(V₅² + V₇² + V₁₁²) / V₁ = √(1.488² + 0.936² + 0.504²) / 19920 ≈ √(2.214 + 0.876 + 0.254) / 19920 = √3.344 / 19920 ≈ 1.829 / 19920 = 0.0000918 = 0.00918%.

4. Step 4: Compare to IEEE 519 limit: 0.00918% ≪ 1.5% → compliant. However, note that this assumes only 5/7/11 harmonics — real systems require up to 50th order per IEEE 519 Annex B.

Answer: The calculated THD_v is 0.0092%, well below the IEEE 519–2022 limit of 1.5% for a 34.5 kV PCC. This indicates no immediate harmonic voltage violation—but resonance risk at 500 Hz (5th × 60 Hz) must still be assessed via impedance scan.

🏗️ Real-World Application

At Rio Tinto’s Gudai-Darri iron ore operation (Western Australia), a 60 MW solar + battery microgrid feeding 40× 350 kW electric haul trucks revealed 6.8% THD_v at the 33 kV substation during midday peak. Standard load flow predicted 92% loading—'within limits.' Harmonic-informed load flow, incorporating inverter harmonic models (per UL 1741 SB) and feeder frequency-dependent impedance, exposed a parallel resonance near 420 Hz (7th harmonic), amplifying V₇ by 4.3×. Remediation included retuning existing 5th-harmonic filters to 7th and enabling inverter Q(V) droop control per IEEE 1547-2018 §6.6. Post-mitigation THD_v dropped to 0.9%, restoring 11 MW of previously unusable feeder capacity.