Ferroresonance in Megawatt Charging Stations: IEEE 519-2022 Limits and I²t Thermal Verification

Phenomenology of Ferroresonance in Medium-Voltage Systems and MCS Transformers

Ferroresonance is a non-linear dynamic resonance phenomenon that occurs in electrical power systems when a circuit containing a saturable non-linear inductance (typically the ferromagnetic core of a cable-connected substation transformer, or MCS transformer) and a capacitance (associated with medium-voltage cables, transmission lines, or grading/coupling capacitors) is excited by an alternating voltage source. Unlike linear resonance, where inductive reactance remains constant and resonance takes place at a single, well-defined frequency: ferroresonance is characterized by the coexistence of multiple stable operating states for the exact same set of system parameters. This occurs because the transformer's magnetizing inductance varies drastically as a function of magnetic flux, transitioning abruptly from a high-impedance linear inductance in the unsaturated state to a very low-impedance inductance under deep core saturation. The non-linear relationship between the linked magnetic flux $\Psi$ and the magnetizing current $i_m$ within the ferromagnetic core can be analytically modeled using a higher-order odd polynomial function: Where $a$ and $b$ are characteristic coefficients of the core material and the magnetic circuit design, and $n$ is a positive integer (typically $n \ge 2$). When the transformer is subjected to an overvoltage transient or an asynchronous switching operation, the magnetic flux exceeds the saturation knee-point. At this operating point, the incremental or tangent inductance: drops by several orders of magnitude. This sharp reduction in inductance allows the inductive reactance of the transformer to match the capacitive reactance of the feeding cable at power frequency or its harmonics, triggering the ferroresonant state. The dynamic behavior of this non-linear circuit is governed by the damped, forced Duffing differential equation: Where $\gamma$ represents resistive losses (damping), $\alpha$ and \beta define the non-linear stiffness of the magnetic circuit, and $F$ is the magnitude of the external driving electromotive force. Depending on the initial load conditions, the switching inception phase angle $\theta$, and the voltage magnitude, the system can converge toward distinct phase-space attractors, such as:

• Fundamental Frequency Mode: The state variables oscillate at the exact fundamental period $T = 2\pi/\omega$ of the excitation source, exhibiting severe harmonic distortion.
• Subharmonic Mode: Oscillations settle at an integer fraction of the fundamental frequency (typically $T/3$ or $T/5$), driving sustained, high-amplitude magnetizing current peaks of prolonged duration.
• Quasi-Periodic Mode: The system exhibits oscillations at frequencies that are non-integer multiples of the driving frequency, manifesting as a continuous spectrum with discrete harmonic and interharmonic peaks.
• Chaotic Mode: Non-periodic behavior exhibiting extreme sensitivity to initial conditions, where the trajectory in phase space forms a strange attractor, inducing erratic, unpredictable, and severe overvoltages and overcurrents.

Transition of the Ferromagnetic Core into Deep Saturation

During the steady-state operation of an MCS (Medium-Voltage Cable-Connected Substation) transformer, the core operates within the linear region of its hysteresis loop, where relative magnetic permeability $\mu_r$ is extremely high ($10,000$ to $100,000$). Under these conditions, the magnetizing current required to sustain the core flux is negligible (less than 1% of the rated full-load current). However, if a single-pole switching event occurs (such as single-phase opening via an expulsion fuse or a single-phase disconnector/interrupter), the return path current is forced through the phase-to-ground and inter-phase coupling capacitances of the medium-voltage cables. This induces a severe temporary overvoltage that shifts the magnetic operating point well beyond the saturation knee ($B > 1.7 \, T$ for grain-oriented silicon steel). As the core saturates, the incremental relative permeability $\mu_r$ collapses toward unity ($\mu_r \approx 1$), approaching the magnetic behavior of air. The effective magnetizing inductance drops precipitously, allowing the cable capacitance to cyclically charge and discharge directly into the transformer winding, sustaining the ferroresonant oscillation indefinitely without requiring a high-capacity external driving source.

Circuit Topology Analysis and Mathematical Modeling of Series and Parallel Ferroresonance

The circuit configuration dictates the precise nature of the ferroresonant mode, as well as the resulting dielectric and thermal stress profiles imposed on the electrical infrastructure. A rigorous analytical distinction between series and parallel topologies is essential.

Series Ferroresonance in MCS Transformers

Series ferroresonance occurs primarily when a system capacitance is placed in direct electrical series with the non-linear magnetizing inductance of the transformer. This scenario frequently develops in MCS substations fed by long runs of shielded underground medium-voltage cables during single-phase energization/de-energization procedures, or following the operation of a single-phase primary fuse. Consider a three-phase transformer with an ungrounded wye (isolated neutral) or delta primary winding, supplied via a shielded cable of length $L_c$ exhibiting a zero-sequence capacitance $C_0$ and positive-sequence capacitance $C_1$. When Phase A opens, the magnetizing inductance of that phase becomes connected in series with the equivalent capacitive network of the cable. The Thévenin equivalent capacitance seen from the terminals of the open-phase saturable inductance is given by: Where $C_d$ represents the inter-conductor coupling capacitance of the cable system. Neglecting losses for a conservative worst-case analysis, the time-domain state equations governing this coupled electro-magnetic system are: Substituting the non-linear polynomial expression for the magnetizing current yields the coupled non-linear differential system: This mathematical formulation reveals a classic saddle-node bifurcation instability. When the supply voltage exceeds a critical threshold $Vcrit$, the operating state undergoes a discontinuous jump from the normal low-voltage stable branch to the high-voltage ferroresonant branch, where terminal voltages across the transformer and the connected cable can reach dangerous levels ranging from $3.0 \, p.u.$ to $4.5 \, p.u.$

Parallel Ferroresonance in MCS Transformers

Parallel ferroresonance develops when the non-linear magnetizing inductance is connected in parallel with the system capacitance. This condition is typical in isolated-neutral (ungrounded) networks or systems grounded through high impedance (such as Petersen coils or high-resistance grounding systems) following the initiation and clearing of a single-phase-to-ground fault. It also occurs frequently in inductive voltage transformers (VTs) connected to isolated busbar sections energized via circuit breaker grading capacitors across open contacts. In this topology, the cable-to-ground capacitance $C_g$ and the magnetizing inductance $L_m$ form a parallel resonant tank excited by zero-sequence current injection or capacitive displacement current coupled across open switching chambers. The differential equation defining the nodal voltage $v(t)$ of the parallel circuit is: Where $C_p$ is the total parallel equivalent capacitance, $R_p$ models core iron losses along with any secondary-connected resistive burden, and $i_s(t)$ is the injected excitation current. Differentiating with respect to time and reformulating in terms of the core flux $\Psi$: In parallel ferroresonance, while the total line current drawn by the system may remain relatively moderate, the neutral displacement voltage (neutral-to-ground potential) escalates to severe, sustained overvoltage levels, permanently compromising the dielectric integrity of instrument transformers and cable termination accessories.

Overvoltage Dynamics and Thermal Insulation Degradation: The I²t Thermal Boundary

The primary operational hazard of ferroresonance lies in the synergistic combination of sustained high-frequency overvoltages and a massive increase in harmonic current content. This imposes extreme electro-thermal stress simultaneously on the solid insulation of the transformer and the extruded polymeric insulation (XLPE or EPR) of the cable system.

Thermal Withstand Limits of Medium-Voltage Cables (I²t)

Under short-circuit or short-duration severe transient regimes, conductor heating in medium-voltage cables is analyzed as an adiabatic thermodynamic process. This assumes that no significant heat is transferred away from the conductor into the insulation, screen, jacket, or surrounding soil due to the ultra-short duration of the event ($t < 5 \, s$). The conductor thermal withstand limit is defined by the Joule integral, or specific let-through energy: Where:

• $I$ is the RMS value of the transient or fault current ($A$).
• $t$ is the fault or transient duration ($s$).
• $S$ is the conductor cross-sectional area ($mm^2$).
• $K$ is the thermal material constant of the conductor, determined in accordance with IEC 60287 and IEC 60949:

Where $Q_c$ is the volumetric heat capacity of the conductor material ($J/K \cdot m^3$), $\beta$ is the reciprocal of the temperature coefficient of resistance at $0 \, ^\circ C$ ($234.5$ for copper, $228$ for aluminum), $\rho_{20}$ is the conductor resistivity at $20 \, ^\circ C$, $\theta_i$ is the initial operating conductor temperature (typically $90 \, ^\circ C$ for XLPE at rated continuous load), and $\theta_f$ is the maximum permissible short-circuit temperature limit ($250 \, ^\circ C$ for XLPE). Crucially, ferroresonance is a sustained, long-duration phenomenon that can persist for minutes or hours if unmitigated by dedicated protective relays. In this non-adiabatic operational regime, the differential thermal balance equation for the cable must account for the thermal resistances of the insulation $T_1$, outer jacket $T_3$, and surrounding installation medium $T_4$: Where $Cth$ is the equivalent thermal capacitance of the cable per unit length, $Rca(\theta)$ is the AC resistance of the conductor at operating temperature $\theta$ (heavily elevated by skin and proximity effects induced by the high-frequency components of the ferroresonant current), and $\sum Rth$ is the total steady-state thermal resistance from the conductor to ambient. During sustained ferroresonance, total harmonic current distortion ($THD_I$) frequently surpasses 100%, containing high-frequency spectral components that significantly increase the effective AC resistance: This drives Joule losses to levels far exceeding rated design capacity, drastically accelerating the thermal degradation rate of the cable insulation according to the classical Arrhenius thermal aging model: Continuous exposure to temperatures exceeding emergency thermal limits ($130 \, ^\circ C$ for XLPE) causes the polymer to undergo irreversible mechanical loss-of-elasticity, accelerates thermal treeing, and ultimately culminates in catastrophic dielectric breakdown.

Thermal and Dielectric Stress in Transformer Windings

In an MCS transformer, deep core saturation forces large fractions of the total magnetic flux to spill out of the laminated iron core into surrounding structural metallic pathways, such as core clamping frames, tie-plates, and tank walls. This unguided leakage flux induces immense eddy currents in these structural steel members, producing rapid, localized extreme hot-spots. Furthermore, hysteresis and eddy-current core losses increase exponentially with magnetic flux density $B$ and oscillation frequency: Where $\eta$ is the Steinmetz coefficient, $\sigma_e$ is the electrical conductivity of the silicon steel, and $d$ is the lamination thickness. Under ferroresonant conditions, the transformer hot-spot winding temperature rise ($\theta_H$) is calculated dynamically per IEEE C57.91: Where $\theta_A$ is the ambient temperature, $\Delta \theta_{\text{TO}}$ is the top-oil temperature rise over ambient, and $\Delta \theta_{\text{H}}$ is the hot-spot temperature rise over top-oil. The compounding effect of excessive thermal elevation and severe high-frequency voltage gradients rapidly degrades the oil-impregnated Kraft paper insulation (accelerating loss of Degree of Polymerization, DP), generating combustible dissolved gases and triggering inter-turn dielectric failure.

Comparative Table of Electrical Parameters, Normative Limits, and Dielectric Consequences

The following table provides an exhaustive technical comparison of critical operational parameters, international standards boundaries, and the associated physical-chemical degradation mechanisms observed in MCS substation assets during ferroresonance:

Parameter / Scenario	Standard Limits (IEEE / IEC)	Critical Inception Condition	Dielectric & Operational Consequences
Temporary Overvoltage (TOV)	IEC 60071-1: $U_m$ (highest voltage for equipment). IEEE C62.11: Surge arrester TOV withstand curve (typically $< 1.4 \, p.u.$ for $t > 10 \, s$).	Series ferroresonance triggered by single-phase clearing or single-pole switching on long cable runs ($> 1.5 \, km$).	Puncture of solid-liquid insulation barrier in transformers; thermal overload, runaway, and explosive failure of metal-oxide surge arresters (MOVs).
Cable Thermal Limit ($I^2t$)	IEC 60287 / IEC 60949: Adiabatic limit $I^2t = K^2 S^2$. $\theta_{max} = 250 \, ^\circ C$ (short-circuit) and $130 \, ^\circ C$ (emergency overload).	Sustained fundamental or subharmonic ferroresonance with continuous RMS current exceeding $2.5 \, p.u.$	Thermal melting of metallic screens, irreversible thermo-mechanical decomposition of XLPE, conductor plastic yield, and flashover across cable terminations.
Harmonic Loss Factor in Transformers	IEEE C57.110: Harmonic Derating Factor for non-sinusoidal load currents: $FHL = \frac{\sum_{h=1}^{max} h^2 I_h^2}{\sum_{h=1}^{max} I_h^2}$.	Chaotic or subharmonic ferroresonant modes rich in low- and high-order odd harmonics (3rd, 5th, 7th, 9th).	Severe localized winding overheating driven by amplified stray and eddy-current losses; thermal embrittlement of Kraft insulation, rapid drop in DP.
Neutral Displacement Voltage	IEEE C57.105: Neutral grounding considerations. Continuous permissible neutral displacement $< 0.1 \, p.u.$	Parallel ferroresonance in ungrounded or high-impedance grounded systems during/after single-line-to-ground faults.	Extreme saturation and thermal destruction of busbar inductive voltage transformers (VTs), primary fuse blowing, and phase-to-ground secondary flashovers.
Surge Arrester Energy Absorption	IEC 60099-4: Line discharge energy classification (Thermal energy rating $Wth$, Charge transfer classes 1 to 5; max $kJ/kV$ of MCOV).	Continuous repetitive conduction during quasi-periodic ferroresonance exceeding arrester continuous operating voltage.	Thermal runaway of $ZnO$ varistor discs; internal flashover, mechanical housing rupture, and violent projection of ceramic/polymeric fragments.

Engineering Failure Forensic Analysis

Forensic failure analysis of MCS substations impacted by ferroresonant dynamics requires a multidisciplinary protocol combining metallurgy, insulation chemistry, and electromagnetic transient simulation. Distinct system components exhibit well-defined, unmistakable forensic signatures.

Forensic Examination of Transformer Core and Windings

When an MCS transformer fails due to ferroresonance, the internal physical evidence is distinct from classical symmetrical fault damage:

• Physical Core Signatures: Localized thermal discoloration and burning patterns are found concentrated at the outer core corner joints (mitred step-lap joints) and core clamping structural bolts. This damage is caused by out-of-plane magnetic flux driven into the core structure, inducing high-magnitude circular eddy currents on the lamination faces. The inter-lamination insulating coating (e.g., Carlite) undergoes pyrolysis, creating inter-laminar short circuits and permanently increasing core no-load losses.
• Dissolved Gas Analysis (DGA - IEC 60599): Dielectric oil samples taken immediately following the failure exhibit critical concentrations of characteristic diagnostic gases. Ferroresonance manifests as a "High-Temperature Thermal Fault" ($T > 700 \, ^\circ C$) coupled with high-energy electrical discharges. The dominant gases are Ethylene ($C_2H_4$) and Acetylene ($C_2H_2$), the latter confirming the presence of active electrical arcing between winding turns or from the active part to the tank wall due to the dielectric breakdown of the overheated oil. Diagnostic gas ratios (Rogers, Dornenburg, or Duval Triangle 1) place the failure definitively in the severe thermal fault ($T3$) and high-energy discharge ($D2$) zones.
• Winding Mechanical and Dielectric Collapse: Severe thermal degradation and shrinkage of the Kraft paper insulation cause complete relaxation of the axial winding clamping pressure. Under the mechanical forces exerted by the ferroresonant current peaks, the windings undergo axial collapse, radial buckling, and ultimately an unrecoverable turn-to-turn dielectric breakdown within the entry section of the primary high-voltage winding.

Forensic Evidence in Medium-Voltage Cables and Terminations

The medium-voltage cables feeding the installation serve as the distributed capacitive bank of the ferroresonant circuit. Their characteristic failure signatures include:

• Screen Melting via Zero-Sequence / Circulating Currents: During series or parallel ferroresonance involving significant neutral shifts, sustained high-frequency currents circulate through the metallic screen wires or copper tape. Because metallic shields are engineered to carry earth fault currents strictly for short clearing times (e.g., $1 \, s$ to $3 \, s$), continuous exposure to multi-ampere ferroresonant currents for several minutes melts the screen wires, burns through the outer PE/PVC oversheath, and disrupts the ground return continuity of the cable run.
• Electrical Treeing Puncture: Microscopic analysis of failed XLPE insulation specimens reveals dense, carbonized electrical tree structures initiating at points of localized electrical stress enhancement (such as semiconducting screen protrusions, contaminants, or micro-voids). These tree channels propagate rapidly under the continuous high-frequency, high-voltage stress imposed by quasi-periodic ferroresonance, resulting in full wall dielectric puncture.

Failure Modes in Switchgear and Surge Arresters

Surge arresters (metal-oxide varistors, $ZnO$) represent the most thermally vulnerable components exposed to ferroresonance:

• Arrester Thermal Runaway Mechanism: Zinc-oxide arresters are designed to clamp short-duration surge events (microsecond-scale lightning impulses and millisecond-scale switching surges); they are physically incapable of dissipating long-duration continuous energy. During ferroresonance, the sustained overvoltage pushes the arrester past its Maximum Continuous Operating Voltage ($MCOV$), forcing the varistor discs into continuous heavy conduction on every half-cycle. The cumulative energy exceeds the thermal absorption capacity (expressed in $kJ/kV$), leading directly to thermal runaway: as the $ZnO$ block temperature rises, its non-linear resistance drops, increasing leakage current and internal dissipation until the blocks crack, puncture, and trigger an internal line-to-ground flashover resulting in housing explosion.

Practical Mitigation Strategies, Design Rules, and Engineering Controls

Mitigating ferroresonance requires proactive engineering measures implemented during the front-end engineering design (FEED) and detailed design phases of the MCS substation, utilizing topology hardening, operating constraints, and damping systems.

Elimination of Single-Pole Switching Devices

The single most effective engineering design measure to prevent series ferroresonance is the complete prohibition of single-phase switching devices—such as drop-out expulsion fuses and single-pole disconnectors—on the medium-voltage side of the transformer. In their place, three-phase, gang-operated circuit breakers or load-break switches must be strictly specified. This ensures that all three phases are made or broken simultaneously with a maximum inter-pole opening discrepancy of less than $2 \, ms$, eliminating single-phasing conditions where one or two phases remain energized through cable series capacitance.

Transformer Neutral Grounding Architecture

The primary winding neutral connection dictates whether a zero-sequence ferroresonant return circuit can form:

• Solidly Grounded Primary Neutral: Direct grounding of the primary wye neutral provides a low-impedance return path for zero-sequence currents, preventing the cable capacitance from becoming trapped in series with the saturable magnetizing branch. While highly effective against series ferroresonance, it increases single-line-to-ground short-circuit currents.
• Wye-Wye Transformer with Delta Tertiary Winding: Implementing a closed-delta tertiary stabilizing winding provides a circulating path for zero-sequence third-harmonic currents, lowering the effective zero-sequence magnetizing impedance and substantially increasing the threshold voltage required to drive the core into ferroresonant saturation.

Passive Damping Systems (Damping Resistors)

For existing installations where switchgear replacement or grounding modifications are technically unfeasible, passive damping systems must be designed and deployed. For inductive voltage transformers (VTs) vulnerable to parallel ferroresonance, a dedicated damping resistor ($R_d$) is installed across the broken-delta (open-delta) secondary winding: Where $C_0$ is the equivalent zero-sequence capacitance of the medium-voltage bus and connected cable system as seen from the primary side, referred to the secondary winding using the turns ratio $a_t^2$: The resistor absorbs and dissipates the stored energy in the resonant tank, forcing the non-linear phase-space trajectories to collapse back to the normal steady-state fundamental attractor. The power rating of the damping resistor must be engineered to withstand continuous circulating current resulting from normal network unbalance, as well as the short-time thermal rating during sustained single-phase-to-ground faults without thermal failure:

Implementation of Tuned Harmonic Filters and Active Electronic Suppressors

In large industrial MCS substations subject to high background harmonic distortion, active protection schemes such as Electronic Ferroresonance Suppressors (EFS) can be deployed. These micro-controller-based units continuously monitor the neutral displacement voltage. Upon identifying the distinct spectral signature of ferroresonance (notably subharmonics at $20 \, Hz$ or $16.67 \, Hz$ in $50 \, Hz$ systems), the EFS immediately triggers anti-parallel power thyristors to insert a low-ohmic damping resistor across the secondary circuit for several cycles, rapidly damping the resonance without incurring continuous operational power losses.

Practical Application and System Simulation Using Vexten Suite

Ferroresonance vulnerability assessment and $I^2t$ thermal boundary validation are performed using the integrated engineering architecture of Vexten Suite. The analytical engineering workflow is executed across three dedicated calculation modules. ``` +-----------------------------------------------------------------+ | VEXTEN SUITE | | | | +---------------------+ +-------------------------------+ | | | MODULE 1: | | MODULE 2: | | | | Vexten Short- | | Vexten Cable Sizing | | | | Circuit | | & Ampacity | | | | - IEC 60909 | | - IEC 60287 | | | | - S_sc, X_s/R_s | | - I^2t adiabatic limit | | | +----------+----------+ +---------------+---------------+ | | | | | | +----------------+---------------+ | | | | | v | | +-----------------------+ | | | MODULE 3: | | | | Vexten Power Factor | | | | & Resonance | | | | - Harmonic scan | | | | - De-tuned filters | | | +-----------------------+ | +-----------------------------------------------------------------+ ```

Module 1: Vexten Short-Circuit (IEC 60909 / IEEE 141)

The engineering procedure begins by determining the equivalent driving grid impedance and system strength at the Point of Common Coupling (PCC) of the MCS substation. Utilizing the Vexten Short-Circuit engine compliant with IEC 60909, the upstream utility grid is modeled using its short-circuit apparent power ($Ssc$) and $X/R$ ratio: Where $c$ is the normative voltage factor per IEC 60909. A low short-circuit capacity (weak grid) results in a high source impedance $Z_k$. This condition dramatically amplifies susceptibility to series ferroresonance because the weak grid is unable to clamp the terminal voltage during sharp transitions in transformer magnetizing current. This module computes the complex transfer impedance matrix, which is directly exported to the electromagnetic transient solver within Vexten Suite to define the initial boundary conditions for dynamic simulation.

Module 2: Vexten Cable Sizing & Ampacity (IEC 60287 / NEC 310)

To evaluate the thermal degradation of the substation's medium-voltage cables under ferroresonant overcurrents, the Vexten Cable Sizing & Ampacity module is deployed. This solver implements the steady-state and transient heat transfer equations defined in IEC 60287. The engineer specifies cable construction parameters (e.g., $15 \, kV$ three-core XLPE cable, $150 \, mm^2$ copper phase conductors, with individual $16 \, mm^2$ copper wire screens), installation topography (direct buried, conduit, or free air), and soil thermal properties (soil thermal resistivity $g = 1.0 \, K \cdot m/W$, ambient soil temperature $\theta_s = 25 \, ^\circ C$). Under ferroresonant operating regimes identified during dynamic simulation, the computed harmonic current spectrum is imported. The module calculates the Harmonic Derating Factor, adjusting the effective AC resistance for each harmonic order $h$: The software integrates the actual Joule heating over the full transient duration: and compares the result against the adiabatic withstand limit $K^2 S^2$ as well as the non-adiabatic thermal jacket deformation limit. If the calculated thermal energy exceeds safe boundaries, Vexten Suite flags a non-compliance alert and calculates the required conductor up-sizing or metallic screen re-specification needed to prevent insulation failure.

Module 3: Vexten Power Factor & Resonance Mitigation

This module is designed for the synthesis and tuning of harmonic and ferroresonant mitigation equipment. It performs an automated Harmonic Impedance Scan across a wide frequency range (fundamental to the 50th harmonic). If the impedance scan uncovers a high-impedance parallel resonance peak aligned with common ferroresonant harmonic frequencies (such as the 3rd or 5th harmonic), the module calculates the required parameters for de-tuned capacitor banks. The de-tuning factor $p$ is determined via: Typically, a de-tuning factor of $p = 7\%$ (tuning frequency $f_r = 189 \, Hz$ for $50 \, Hz$ systems) or $p = 5.67\%$ is selected. This shifts the system resonant frequency safely below the 3rd harmonic, ensuring that the harmonic magnetizing currents generated during transformer core saturation do not encounter a high-impedance parallel tank that could ignite or sustain a parallel ferroresonant mode. Furthermore, the software simulates the dynamic insertion of the engineered broken-delta damping resistor on the voltage transformer secondary, rendering the phase-plane trajectories of flux versus magnetizing current to verify the rapid, stable collapse of the ferroresonant state back to the fundamental steady-state operating point (achieving full extinction in less than $150 \, ms$).