Neutral Shift and High-Altitude TOV in Mining: Thermal Analysis and K-Factor per NFPA 70B

Electrophysical Contextualization of High-Altitude Mining Operations

The operation of electrical power systems in mining facilities located at altitudes exceeding 3,000 meters above sea level (masl) imposes severe constraints on the dielectric strength of the insulating medium and the thermal dissipation capabilities of electrical equipment. At high elevations, the relative air density ( $\delta$ ) decreases drastically in accordance with the fundamental barometric law:

\delta = \frac{P}{P_0} \cdot \frac{T_0}{T} = e^{-\left(\frac{g \cdot M \cdot h}{R \cdot T_0}\right)}

Where $P$ and $T$ represent the absolute pressure and temperature at altitude $h$ (masl), $P_0$ and $T_0$ are the standard sea-level conditions (101.325 kPa and 288.15 K), $g$ is the acceleration due to gravity, $M$ is the molar mass of dry air, and $R$ is the universal gas constant. This reduction in density alters the mean free path of free electrons, increasing the probability of collision ionization under the application of an electric field, according to Paschen's Law:

V_b = \frac{B \cdot (p \cdot d)}{\ln(A \cdot p \cdot d) - \ln\left[\ln\left(1 + \frac{1}{\gamma_{se}}\right)\right]}

Where $V_b$ is the breakdown voltage, $p$ is the atmospheric pressure, $d$ is the gap distance between electrodes, $A$ and $B$ are composite gas constants, and $\gamma_{se}$ is Townsend's second ionization coefficient. As a direct consequence, air insulation distances (phase-to-ground and phase-to-phase) as well as creepage distances experience non-linear degradation.

In accordance with IEC 60071-2 and IEEE C37.100.1, the altitude correction factor for external insulation ( $K_a$ ) is defined according to the mathematical formulation:

K_a = e^{m \cdot \left(\frac{h - 1000}{8150}\right)}

Where $m$ is an empirical parameter dependent on the type of applied voltage (power frequency, switching impulse, or lightning impulse) and the magnitude of the withstand voltage. For mining installations situated at 4,500 masl, the factor $K_a$ typically ranges between 1.45 and 1.60, requiring the derating or oversizing of the rated power-frequency withstand voltage ( $V_w$ ) of circuit breakers, insulators, and medium-voltage switchgear assemblies according to:

V_w \ge V_{req, sea\_level} \times K_a

Simultaneously, thermal dissipation by natural convection is drastically attenuated due to the reduced volumetric mass of the cooling fluid (air). The thermal derating factor ( $kalt$ ) applied to the nominal continuous current rating of equipment such as transformers, busbars, and circuit breakers is approximated via IEEE C57.96:

kalt = 1 - \alpha \cdot \left( \frac{h - 1000}{100} \right)

Where $\alpha \approx 0.005$ to $0.008$ depending on the geometry of the radiant element and the ratio between radiation and convection losses. The simultaneous interaction of dielectric degradation and thermal derating establishes a critical scenario for the inception of temporary overvoltages (TOV) and neutral displacement phenomena in ungrounded or high-resistance grounded (HRG) networks.

Physics of Neutral Displacement and Temporary Overvoltages (TOV)

Neutral displacement (or neutral shift) is an electrodynamic phenomenon occurring in three-phase systems that are not effectively grounded when phase-to-ground admittance symmetry is disrupted or when a single line-to-ground (SLG) fault takes place. Consider a medium-voltage mining distribution system supplying non-linear loads where the system neutral $N$ is connected to ground through a neutral grounding admittance $Y_n = 1/RNGR + j\omega C_n$ .

Applying Kirchhoff's Current Law at the neutral node $N$ with respect to the reference ground potential $G$ , the neutral displacement voltage ( $VNG$ ) under an imbalance of per-phase phasor admittances ( $Y_A, Y_B, Y_C$ ) and generated phase voltages ( $E_A, E_B, E_C$ ) is rigorously expressed as:

VNG = \frac{E_A Y_A + E_B Y_B + E_C Y_C}{Y_A + Y_B + Y_C + Y_n}

Where the admittance of each phase $i \in \{A, B, C\}$ incorporates the distributed capacitance of the shielded/underground cable system ( $C_{0,i}$ ), insulation leakage admittances, and any passive filter components connected to that phase:

Y_i = G_i + j\omega C_{0,i} + \sumh \frac{1}{R_{h,i} + j\left(h\omega L_{h,i} - \frac{1}{h\omega C_{h,i}}\right)}

Under a high-impedance or bolted single line-to-ground fault condition on Phase $A$ ( $R_f \to 0$ ), the admittance $Y_A \to \infty$ . The limit of the neutral displacement equation yields $VNG \to -E_A$ . Consequently, the neutral-to-ground voltage phasor equals the magnitude of the pre-fault phase-to-neutral voltage but with opposite polarity. The phasor voltages of the unfaulted phases ( $VBG$ and $VCG$ ) experience a sudden voltage rise, reaching the pre-fault line-to-line voltage magnitude ( $VLL$ ):

VBG = E_B - VNG = E_B + E_A = \sqrt{3} E_B e^{-j 30^\circ}

VCG = E_C - VNG = E_C + E_A = \sqrt{3} E_C e^{+j 30^\circ}

The temporary overvoltage factor ( $kTOV$ ) for the healthy phases in an HRG system is defined by the ratio between the peak phase-to-ground voltage during the fault and the nominal peak phase-to-neutral voltage:

kTOV = \frac{|VBG|}{|E_B|} = \sqrt{3} \approx 1.732 p.u.

In the symmetrical components domain, temporary overvoltage due to neutral displacement is evaluated by analyzing the positive sequence ( $Z_1$ ), negative sequence ( $Z_2$ ), and zero sequence ( $Z_0$ ) networks. The phase voltage expression during an SLG fault is given by:

VBG = a^2 E_A - \left( \frac{Z_2 + Z_0 + 3R_f}{Z_1 + Z_2 + Z_0 + 3R_f} a^2 + \frac{Z_2 - a Z_0}{Z_1 + Z_2 + Z_0 + 3R_f} a + \frac{Z_0 (1-a^2)}{Z_1 + Z_2 + Z_0 + 3R_f} \right) E_A

When the ratios $X_0 / X_1 \gg 3$ and $R_0 / X_1 \gg 1$ (characteristic of ungrounded or high-resistance grounded systems), the system is classified as non-effectively grounded per IEEE 142. Under these conditions, the coefficient $kTOV$ is not constrained to $1.38$ p.u. (the threshold for effectively grounded networks), but reaches steady-state levels exceeding $1.732$ p.u., with high-frequency LC oscillating transients capable of reaching peak values of:

V_{peak, transient} = Vprefault \cdot \left( 1 + \eta \cdot \sin(\omega_r t) e^{-\frac{t}{\tau}} \right) \approx 2.5 to 3.0 p.u.

The presence of high distributed capacitance in underground or open-pit mining feeders (such as SHD-GC trailing cables with high per-unit capacitance $C_0 \approx 0.3 - 0.6 \, \mu F/km$ ), combined with the non-linear inductance of potential transformers or unloaded power transformers, initiates ferroresonant coupling. Parallel ferroresonance is triggered when the zero-sequence capacitive reactance ( $XC0 = 1 / (\omega C_0)$ ) equals the non-linear magnetizing reactance of the saturated transformer ( $X_m(i)$ ):

f_r = \frac{1}{2\pi \sqrt{L_m(i) \cdot C_0}}

Under ferroresonance, the neutral potential oscillates chaotically, producing extreme sustained overvoltages ( $TOV > 3.0 p.u.$ ) with subharmonic ( $1/3 f_1$ , $1/2 f_1$ ) or harmonic spectrum components that frequently exceed the altitude-corrected dielectric withstand voltage of air ( $V_b / K_a$ ).

Mining Load Harmonics and K-Factor per NFPA 70B and IEEE C57.110

Mining processing plants operate high-power non-linear loads, such as medium-voltage variable frequency drives (VSI VFDs with 3-level Neutral Point Clamped NPC topology, 6/12/18-pulse or AFE rectifiers), SAG mill drives, and ball mill drives driven by cycloconverters or LCI-fed synchronous motors. These loads inject harmonic currents into the power system with spectral orders:

h = p \cdot k \pm 1

Where $p$ is the pulse number of the converter and $k \in \{1, 2, 3, \dots\}$ . For 6-pulse rectifiers, $h \in \{5, 7, 11, 13, 17, 19, \dots\}$ . In systems with phase unbalance or 4-wire configurations, triplen harmonics emerge ( $h = 3, 9, 15, 21, \dots$ ), which possess zero-sequence ( $Z_0$ ) characteristics.

Triplen harmonic currents do not cancel at the neutral node of wye-connected transformers with accessible neutrals; on the contrary, they add up arithmetically in the neutral conductor:

I_{N, rms} = \sqrt{3 \sum_{k=1}^{\infty} \left( I_{3(2k-1)} \right)^2 + I_{N, fundamental}^2}

The flow of harmonic currents through transformer windings causes severe overheating due to the non-linear increase in winding eddy current losses ( $PEC$ ) and stray losses in structural components ( $POS$ ). The K-Factor ( $K-Factor$ ), defined in IEEE C57.110 and adopted in NFPA 70B for insulation integrity assessment in preventive and predictive maintenance, quantifies the thermal heating effect of harmonic currents relative to the fundamental current:

K = \sum_{h=1}^{hmax} h^2 \left( \frac{I_h}{I_1} \right)^2 = \frac{\sum_{h=1}^{hmax} (h \cdot I_h)^2}{\sum_{h=1}^{hmax} I_h^2}

Where $I_h$ is the RMS current of harmonic order $h$ , and $I_1$ is the fundamental RMS current. The eddy current loss factor ( $FHL$ ) used to derate transformer nominal capacity is expressed as:

FHL = \frac{\sum_{h=1}^{hmax} h^2 \left( \frac{I_h}{I} \right)^2}{\sum_{h=1}^{hmax} \left( \frac{I_h}{I} \right)^2} = \frac{K}{1 + THD_I^2}

Where $THD_I = \frac{\sqrt{\sum_{h=2}^{\infty} I_h^2}}{I_1}$ is the Total Harmonic Distortion of current. The derated loading capacity of the transformer in per-unit ( $Pmax(pu)$ ) under a combined harmonic and high-altitude environment is recalculated by uniting IEEE C57.110 guidelines with the altitude derating factor:

Pmax(pu) = kalt \cdot \sqrt{ \frac{PLL(pu)}{1 + FHL \cdot P_{EC-R}(pu)} }

Where $PLL(pu)$ represents the total full-load losses calculated in per-unit, and $P_{EC-R}(pu)$ represents the winding eddy current losses calculated at fundamental frequency under rated conditions.

According to NFPA 70B (Standard for Electrical Equipment Maintenance), the accumulation of severe harmonic distortion combined with elevated operating temperatures accelerates the degradation process of thermal insulation following the Arrhenius equation for expected insulation lifetime ( $\mathcal{L}$ ):

\mathcal{L} = \mathcal{L}_0 \cdot \exp\left( \frac{E_a}{R \cdot Thotspot} \right)

Where $E_a$ is the dielectric activation energy, $R$ is the universal gas constant, and $Thotspot$ is the winding hot-spot temperature, which rises dramatically according to:

Thotspot = T_{amb, alt} + \Delta T_{A-R} \cdot \left[ \frac{PDC + FHL \cdot P_{EC-R}}{PDC + P_{EC-R}} \right]^y + \Delta THVR

NFPA 70B mandates for high-altitude industrial facilities compulsory infrared thermography adjusted for atmospheric transmittance corrected for altitude, dissolved gas analysis (DGA) in oil with lowered action thresholds for ethylene and acetylene due to TOV-induced partial discharges, and periodic verification of the distribution system K-Factor.

Forensic Analysis of Integrated Electromechanical Failures

The combination of neutral displacement, temporary overvoltages (TOV), low-pressure dielectric attenuation, and harmonic thermal stress triggers multi-phase failure mechanisms across critical mining assets. The following section provides a detailed forensic analysis categorized by active asset class:

Power Cable Insulation (XLPE / EPR)

In medium- and high-voltage cables (e.g., 13.8 kV or 34.5 kV XLPE-insulated systems), a neutral shift sustains a phase-to-ground voltage of $\sqrt{3} VLN$ across the insulation of healthy phases for extended hours if the HRG protection system does not isolate the fault rapidly. The peak voltage applied across the internal dielectric is:

V_{peak, actual} = \sqrt{2} \cdot \left( \sqrt{3} VLN \right) \cdot \left( 1 + THD_V \right)

This overvoltage significantly lowers the Partial Discharge Inception Voltage (PDIV / CIV). Micro-voids embedded within the XLPE matrix generate a localized electric field gradient $Evoid = \epsilon_{r} \cdot E_0$ . When $Evoid$ exceeds the breakdown voltage of the gas trapped within the void (whose dielectric strength drops per Paschen's Law under pressure migration or altitude gradients in cable terminations), sustained partial discharge activity initiates. These discharges erode the polymer matrix, driving electrical treeing and accelerating catastrophic dielectric breakdown of the cable line.

Medium-Voltage Switchgear and Air Insulation

Air within medium-voltage switchgear enclosures functions as the primary insulating medium between phase busbars and the grounded structural enclosure. At 4,200 masl, the dielectric strength of air drops by a factor of approximately $\delta \approx 0.62$ . During a TOV event associated with an SLG fault in an HRG system with high 3rd and 9th harmonic content:

The clearance distance required to withstand switching surges or arc-reignition transients must be increased by $45\% - 60\%$ pursuant to IEC 60664-1 Table A.2.
If the switchgear was rated for sea-level operation ( $Vdisruption = 95 kV BIL$ ), at 4,200 masl its actual BIL withstand capability degrades to $95 \times 0.62 = 58.9 kV BIL$ .
A TOV transient peaking at 48 kV phase-to-ground overlaid with high-frequency harmonics exceeds the reduced dielectric limit of the rarefied air, initiating phase-to-ground flashover that rapidly escalates into a destructive three-phase arc flash event via massive air ionization within the cubicle.

Power Transformers and K-Factor Degradation

In power transformers, an elevated neutral voltage shifts the operating point on the core magnetization curve ( $B-H$ ). If the TOV elevates the flux density $B$ beyond the saturation threshold ( $Bsat \approx 1.7 - 1.8 Tesla$ ):

B(t) = \frac{1}{N \cdot A} \int0^{t} VTOV(\tau) d\tau

The magnetizing current transforms into sharp current pulses rich in odd and triplen harmonics ( $Imag \propto B^{n}, n > 9$ ). This exponentially escalates core iron losses ( $P_Fe \propto f \cdot B^2 + f^2 \cdot B^2$ ) and induces elevated winding eddy current losses ( $PEC$ ). Trapped thermal energy, unvented due to $kalt$ , thermally degrades the cellulosic paper insulation, releasing carbon monoxide and carbon dioxide ( $CO, CO₂$ ) and inducing localized hot-spots that trigger furanic compound generation. This reduces the Degree of Polymerization (DP) of the Kraft paper from 1,000 to critical thresholds (< 200), signaling mechanical end-of-life and imminent turn-to-turn winding insulation failure.

Electromechanical / Insulation Parameter	Standard Operating Condition (0 masl, $THD_I < 5\%$ )	Severe Mining Operating Condition (4500 masl, $K \ge 13, TOV = 1.73 \, pu$ )	Imposed Thermal and Dielectric Consequence
Dielectric Strength of Air ( $E_b$ )	$3.0 kV/mm$	$\sim 1.85 kV/mm$ ( $\delta \approx 0.61$ )	Phase-to-ground dielectric breakdown in switchgear and post insulators.
Healthy Phase Phase-to-Ground Voltage ( $VLG$ )	$1.0 p.u.$ ( $VLL / \sqrt{3}$ )	$1.732 p.u.$ sustained (or $> 2.5 p.u.$ transient)	Accelerated aging of XLPE insulation via partial discharge and electrical treeing.
K-Factor of Harmonic Loss Derating	$K = 1.0$	$K = 13.0 - 20.0$	Winding eddy losses ( $PEC$ ) increased by $1300\% - 2000\%$ .
Hot-Spot Temperature ( $Thotspot$ )	$98^\circ C - 110^\circ C$	$> 145^\circ C$ (without specific derating)	Destruction of Kraft paper polymeric chains (DP < 200).
Partial Discharge Inception Voltage (PDIV / CIV)	$100\%$ nominal value	$\sim 60\% - 65\%$ of nominal value	Sustained corona discharge generation and external insulation erosion.
Radiator Dissipation Efficiency ( $h_c$ )	$100\%$ convective efficiency	$\sim 70\%$ convective efficiency	Systemic overheating of insulating and cooling liquid medium.

Advanced Mitigation Strategies and Design Criteria

To ensure operational continuity and complete engineering integrity across high-altitude mining electrical networks under combined conditions of high altitude, harmonic distortion, and TOV risks, the following detailed design methodologies must be executed:

Neutral Grounding Resistor (NGR) Sizing for TOV Control and Altitude Derating

In high-altitude HRG systems, the current selected for the Neutral Grounding Resistor ( $INGR$ ) must be strictly greater than the total system zero-sequence capacitive charging current ( $IC0$ ):

INGR \ge IC0 = 3 \cdot \omega \cdot C_{0, total} \cdot VLN

Where $C_{0, total}$ is the aggregated charging capacitance of all cables, capacitor banks, and stray capacitances across the network. By enforcing $RNGR \le \frac{1}{3 \omega C_{0, total}}$ , capacitive energy stored during an SLG fault is safely dissipated, fully damping arcing ground overvoltage spikes and maintaining $kTOV \le 1.732 p.u.$ .

Additionally, the physical resistor element $R$ (stainless steel or nickel-chromium alloy bank) must be sized considering high-altitude thermal derating to handle continuous or short-time fault current dissipation ( $tfault = 10 s$ or continuous per NFPA 70B):

R_{NGR, derated} = R_{NGR, 20^\circ C} \left[ 1 + \alpha_{\text{T}} \cdot (\Delta Talt) \right]

P_{NGR, nominal} \ge \frac{VLN^2}{RNGR} \times \frac{1}{k_{alt, NGR}}

Surge Arrester Selection and Insulation Coordination

Selecting zinc oxide (ZnO) gapless surge arresters for high-altitude mining applications subject to TOV requires rigorous margin evaluation between Maximum Continuous Operating Voltage (MCOV or $V_c$ ) and the arrester's TOV withstand capability curve ( $VTOV(t)$ ). The MCOV must exceed the maximum operating phase-to-ground voltage during a sustained neutral shift:

VMCOV \ge V_{LL, max} \quad (For HRG networks with clearing delays)

The surge arrester must absorb transient overvoltage energy without undergoing thermal runaway. The specific energy absorption capability ( $E_K$ in kJ/kV of MCOV) is derated for altitude due to reduced heat transfer from the polymeric housing:

E_{K, min} = \frac{Etransient + ETOV}{k_{alt, arrester}}

Furthermore, the external creepage distance of the arrester housing is sized using a minimum metric of $31 mm/kV$ of maximum system voltage, multiplied by the $K_a$ factor to withstand heavy mining pollution (conductive metallic or saline dust).

Harmonic Filtering and Detuning to Prevent Parallel Resonance

To eliminate harmonic spectrum injection and prevent transformer overheating due to high K-Factor levels, detuned passive filter banks or C-type damped passive filters are deployed. The filter tuning frequency ( $fst$ ) must be placed below the 5th harmonic ( $h=5$ ), typically tuned to $h_r = 4.2 - 4.7$ , to eliminate parallel resonance risk between power factor correction capacitors and system source inductance ( $Lsys + Ltransformer$ ):

f_r = f_1 \cdot \sqrt{ \frac{X_C}{X_L + Xsys} }

When feeding large variable frequency drives, parallel-connected Active Power Filters (APF) utilizing multilevel converter topologies inject harmonic cancellation currents in opposite phase ( $IAPF(t) = -Ih(t)$ ). This reduces current Total Harmonic Distortion ( $THD_I < 5\%$ ) at the Point of Common Coupling (PCC) per IEEE 519, effectively restoring the transformer K-Factor to $K \approx 1.0$ .

Application and Calculation Methodology with Vexten Suite

The Vexten Suite electrical engineering software integrates numerical and analytical simulation modules designed to simultaneously solve complex short-circuit calculations, insulation altitude derating, harmonic power flow, and cable thermal ampacity within high-altitude mining systems. The following calculation methodology details the workflow implemented within the Vexten Suite architecture.

Vexten Short-Circuit & Grounding Module (IEC 60909 / IEEE 141)

The Vexten Short-Circuit & Grounding Analyzer calculates the zero-sequence impedance matrix ( $Z_1, Z_2, Z_0$ ) adjusted for high-altitude operating temperatures and ground resistivity. During a single line-to-ground fault in an HRG mining network, the software calculates the initial symmetrical short-circuit current ( $Ik1''$ ), neglecting conventional load impedance effects per IEC 60909-0:

Ik1'' = \frac{\sqrt{3} \cdot c \cdot V_n}{Z_1 + Z_2 + Z_0 + 3Z_f}

Because the zero-sequence impedance in an HRG system is dominated by the neutral resistor ( $Z_0 \approx 3 RNGR$ ), Vexten solves the ground fault current via:

Ik1'' \approx \frac{c \cdot V_n}{\sqrt{3} RNGR}

Simultaneously, the Vexten core calculates phase voltage vectors across healthy buses, evaluating the exact $k_{TOV, bus}$ factor and neutral displacement $VNG$ across all operating topologies.

Vexten Cable Sizing & Harmonic Derating Engine (IEC 60287 / NEC 310)

The Vexten Cable Sizing & Harmonic Derating Engine determines the effective ampacity ( $Iamp$ ) of single-core and multi-core medium/low voltage cables (e.g., SHD-GC or XLPE types) subject to combined altitude and neutral harmonic currents using a multi-dimensional derating equation:

I_{amp, corrected} = Itabulated \cdot ktemp \cdot kgrouping \cdot kalt \cdot kharm

Where the neutral harmonic current derating factor ( $kharm$ ) is computed internally by evaluating phase and neutral conductor heating per NEC 310.15(E) and IEC 60287-1-1. If triplen harmonic content ( $h=3$ ) exceeds $33\%$ , neutral current surpasses phase current; Vexten automatically increases neutral conductor sizing to $150\% - 200\%$ of the phase area, recalculating the thermal gradient across concentric layers:

\Delta \theta = (Pcond + Pdiel) \cdot T_1 + \left[ Pcond(1 + y_s + y_p) + Pdiel \right] \cdot n \cdot (T_2 + T_3 + T_4)

Where $T_1, T_2, T_3, T_4$ represent thermal resistances of cable components (insulation, sheath, bedding, surrounding medium) adjusted within Vexten based on atmospheric density at the user-specified altitude.

Vexten K-Factor & Resonance Module (IEEE C57.110 / IEEE 519)

The Vexten Harmonic & Resonance Analyzer performs system frequency sweeps ( $Z(f)$ ) from $50 Hz$ to $2500 Hz$ in $1 Hz$ steps. The algorithm identifies zero-derivative inflection points ( $\frac{d|Z|}{df} = 0$ ), identifying exact parallel ( $f_p$ ) and series ( $f_s$ ) resonance frequencies.

Based on user-defined or imported harmonic current spectra (PQDIF/COMTRADE format), Vexten calculates:

System K-Factor per IEEE C57.110.
Winding eddy current loss factor $FHL$ .
Derated transformer load capacity in kVA ( $kVA_{derated} = kVA_{nominal} \cdot Pmax(pu)$ ).
Hot-spot temperature elevation profile ( $Thotspot$ ) and loss-of-life rate based on NFPA 70B Arrhenius modeling.

Integrated Computational Algorithm in Vexten Suite

The structural flowchart below illustrates the computational logic executed by Vexten Suite to resolve the interconnected physics of altitude insulation breakdown, neutral displacement overvoltages, and K-Factor thermal losses:

\begin{matrix} \text{\textbf{Step 1: Input Processing}} & \longrightarrow & h (masl) , Tamb, Harmonic Spectrum \{h, I_h\}, Neutral Topology (HRG/Solid) \\ \Downarrow & & \\ \text{\textbf{Step 2: Dielectric Correction}} & \longrightarrow & \delta = e^{-\frac{g M h}{R T_0}} \quad \implies \quad K_a = e^{m \frac{h-1000}{8150}} \quad \implies \quad V_{BIL, req} = V_{BIL, std} \cdot K_a \\ \Downarrow & & \\ \text{\textbf{Step 3: Displacement Analysis}} & \longrightarrow & VNG = \frac{\sum E_i Y_i}{\sum Y_i + Y_n} \quad \implies \quad kTOV = \frac{|VBG|}{|E_B|} \quad \implies \quad Check Arrester MCOV \\ \Downarrow & & \\ \text{\textbf{Step 4: K-Factor Evaluation}} & \longrightarrow & K = \frac{\sum (h I_h)^2}{\sum I_h^2} \quad \implies \quad FHL = \frac{K}{1+THD_I^2} \quad \implies \quad Pmax(pu) = kalt \sqrt{\frac{PLL}{1+FHL P_{EC-R}}} \\ \Downarrow & & \\ \text{\textbf{Step 5: Design Optimization}} & \longrightarrow & NGR Sizing (INGR \ge 3 IC0) \quad \& \quad Detuned Filter Design (fst < 5th) \end{matrix}

Through this integrated computational engine, Vexten Suite delivers a high-precision platform that eliminates catastrophic failure risks driven by temporary overvoltages, dielectric arc flashover, and thermal insulation collapse across complex high-altitude mining power networks.