12.1 Transmission Line Modeling, Voltage Drop & Conductor Sizing

12.1 Transmission Line Modeling, Voltage Drop & Conductor Sizing

Executive Overview: Transmission and distribution lines constitute the vital backbone of electric power systems. On the NCEES PE Power examination, line modeling questions demand total fluency across three interrelated competencies: computing fundamental line parameters ($R, L, C$) from physical conductor geometries and bundle configurations; formulating two-port $ABCD$ transmission matrices across Short, Medium (Nominal $\pi$), and Long line models; and evaluating voltage drop, voltage regulation ($VR%$), transmission efficiency ($\eta%$), and thermal ampacity limits.

---

1. Transmission Line Parameters & Physical Geometry

Every overhead transmission line exhibits four distributed parameters per unit length: series resistance ($r$), series inductance ($l$), shunt capacitance ($c$), and shunt conductance ($g$). In power frequency steady-state analysis, shunt conductance $g$ (representing leakage across insulators) is negligible ($g \approx 0$).

Transmission Line Distributed Parameter Per-Unit-Length Representation:

         r*dx        j*omega*l*dx
---o---/\/\/\-------UUUUUUUUUU-------------------o---
   |                                             |
  ---                                           ---
  --- (j*omega*c*dx)/2                         --- (j*omega*c*dx)/2
   |                                             |
---o---------------------------------------------o---

Series Resistance and Skin Effect

The ac resistance of a conductor exceeds its dc resistance due to skin effect, where alternating magnetic flux forces current density toward the outer conductor surface:

$R_{ac} = k_{skin} \cdot R_{dc} \quad (k_{skin} \approx 1.02 \text{ to } 1.05 \text{ at } 60\text{ Hz for standard ACSR})$

Temperature correction for conductor resistance follows the linear relation:

$R_{T2} = R_{T1} \left[ 1 + \alpha_{T1} (T_2 - T_1) \right] = R_{T1} \left( \frac{T_2 + T_0}{T_1 + T_0} \right)$

where for hard-drawn aluminum (EC grade), the inferred zero-resistance temperature is $T_0 = 228.1^\circ\text{C}$ and $\alpha_{20} = 0.0039 / ^\circ\text{C}$.

Series Inductance, GMD, and GMR

Inductance per phase per unit length for a transposed three-phase line arises from both internal and external magnetic flux linkages:

$L = 2 \times 10^{-7} \ln\left( \frac{GMD}{GMR_L} \right) \quad [\text{H/m}] \quad \Longleftrightarrow \quad X_L = 2\pi f L \quad [\Omega/\text{m}]$

1. Geometric Mean Distance ($GMD$): Represents the mutual geometric spacing between phases $a, b, c$: $GMD = \sqrt[3]{D_{ab} \cdot D_{bc} \cdot D_{ca}}$
2. Geometric Mean Radius ($GMR_L$ or $D_s$): Accounts for internal flux by replacing a solid conductor of radius $r$ with a fictitious tubular conductor of radius $r' = r e^{-1/4} = 0.7788 r$. For stranded conductors (e.g., ACSR), $GMR_L$ values are tabulated directly in the NCEES PE Power Reference Handbook.

Conductor Bundling Dynamics

High-voltage transmission lines ($\ge 230\text{ kV}$) utilize bundled conductors (multiple sub-conductors per phase held by mechanical spacers at distance $d$). Bundling increases the effective $GMR$, which:

• Drastically reduces series inductive reactance $X_L$.
• Increases shunt capacitance $C$.
• Suppresses surface electric field gradients, eliminating corona discharge, audible noise, and radio interference.

$\begin{array}{|c|c|c|} \hline \textbf{Bundle Configuration} & \textbf{Number of Sub-conductors } (n) & \textbf{Modified Bundle } GMR_{bundle} \\ \hline \text{2-Conductor Bundle} & 2 & GMR_{b} = \sqrt{D_s \cdot d} \\ \hline \text{3-Conductor Bundle (Equilateral)} & 3 & GMR_{b} = \sqrt[3]{D_s \cdot d^2} \\ \hline \text{4-Conductor Bundle (Square)} & 4 & GMR_{b} = \sqrt[4]{D_s \cdot d^3 \cdot \sqrt{2}} = 1.091 \sqrt[4]{D_s \cdot d^3} \\ \hline \end{array}$

Shunt Capacitance

The line-to-neutral capacitance per phase for a transposed three-phase line is:

$C_n = \frac{2\pi \varepsilon_0}{\ln\left( \frac{GMD}{GMR_C} \right)} \quad [\text{F/m}] \quad \text{where } \varepsilon_0 = 8.854 \times 10^{-12} \text{ F/m}$

[!IMPORTANT] Critical Inductance vs. Capacitance GMR Distinction: For inductance calculations, $GMR_L$ incorporates the internal magnetic flux penalty ($r' = 0.7788r$). For capacitance calculations, all electric charge resides strictly on the outer conductor surface (no internal electric field), so the actual physical outer radius is used: $GMR_C = r$. For bundled conductors, compute $GMR_{C,bundle}$ using physical radius $r$ instead of $D_s$.

---

2. Line Classification & Two-Port ABCD Parameters

Transmission lines are modeled as two-port, four-terminal linear networks relating sending-end voltage and current $(\mathbf{V}_S, \mathbf{I}_S)$ to receiving-end quantities $(\mathbf{V}_R, \mathbf{I}_R)$:

$\begin{bmatrix} \mathbf{V}_S \\ \mathbf{I}_S \end{bmatrix} = \begin{bmatrix} A & B \\ C & D \end{bmatrix} \begin{bmatrix} \mathbf{V}_R \\ \mathbf{I}_R \end{bmatrix} \quad \Longleftrightarrow \quad \begin{aligned} \mathbf{V}_S &= A\mathbf{V}_R + B\mathbf{I}_R \\ \mathbf{I}_S &= C\mathbf{V}_R + D\mathbf{I}_R \end{aligned}$

All passive, bilateral, symmetrical transmission lines satisfy two fundamental network identities:

1. Reciprocity: $AD - BC = 1$
2. Symmetry: $A = D$

For long lines, the propagation constant $\boldsymbol{\gamma}$ and characteristic (surge) impedance $Z_c$ are:

$\boldsymbol{\gamma} = \sqrt{z \cdot y} = \alpha + j\beta \quad [\text{m}^{-1} \text{ or km}^{-1}], \qquad Z_c = \sqrt{\frac{z}{y}} \quad [\Omega]$

where $\alpha$ is the attenuation constant (neper/km) and $\beta = \omega \sqrt{lc}$ is the phase shift constant (rad/km).

Medium Line Nominal Pi Circuit Model:

         I_S ------------[ Z = R + jX ]------------> I_R
          |                                           |
         --- I_C,S                                   --- I_C,R
   V_S  | Y/2                                  V_R  | Y/2      Load
         ---                                         --- 
          |                                           |
  --------o-------------------------------------------o--------

---

3. Voltage Drop Formulations & Voltage Regulation

Exact Phasor Voltage Drop

For a three-phase system evaluated on a per-phase (line-to-neutral) basis:

$\mathbf{V}_{S,LN} = \mathbf{V}_{R,LN} + \mathbf{I}_R (R + jX_L)$

where receiving-end phase voltage is taken as the reference: $\mathbf{V}{R,LN} = |V{R,LN}|\angle 0^\circ$, and load current is $\mathbf{I}_R = |I_R|\angle -\theta$ (for lagging power factor $\cos\theta$).

Approximate Voltage Drop Formula

For short lines and radial distribution feeders where $|\mathbf{I}_R \mathbf{Z}| \ll |\mathbf{V}_R|$, the phase-to-neutral and phase-to-phase voltage drops are approximated by:

$\Delta V_{LN} \approx |I_R| (R \cos\theta \pm X_L \sin\theta) \quad [\text{Volts line-to-neutral}]$

$\Delta V_{LL} \approx \sqrt{3} |I_R| (R \cos\theta \pm X_L \sin\theta) \quad [\text{Volts line-to-line}]$

Sign convention: Use $(+)$ for lagging power factor (inductive loads); use $(-)$ for leading power factor (capacitive loads).

Percent Voltage Regulation ($VR%$)

Voltage regulation quantifies the change in receiving-end voltage magnitude from full-load to no-load with sending-end voltage held constant:

$VR\% = \frac{|V_{R,NL}| - |V_{R,FL}|}{|V_{R,FL}|} \times 100\%$

From the general $ABCD$ equation, at no-load ($\mathbf{I}_R = 0$):

$\mathbf{V}_S = A \mathbf{V}_{R,NL} \implies |V_{R,NL}| = \frac{|V_S|}{|A|}$

$\therefore VR\% = \frac{\frac{|V_S|}{|A|} - |V_{R,FL}|}{|V_{R,FL}|} \times 100\%$

[!NOTE] For short lines, $A = 1$, so $|V_{R,NL}| = |V_S|$. For medium and long lines, $|A| < 1$ under normal operating conditions, making $|V_{R,NL}| > |V_S|$ due to line charging.

Transmission Line Efficiency ($\eta%$)

$\eta\% = \frac{P_{R,3\phi}}{P_{S,3\phi}} \times 100\% = \frac{3 |V_{R,LN}| |I_R| \cos\theta_R}{3 |V_{S,LN}| |I_S| \cos\theta_S} \times 100\% = \frac{P_{R,3\phi}}{P_{R,3\phi} + 3 |I_{line}|^2 R} \times 100\%$

---

4. Conductor Thermal Ampacity vs. Voltage Drop Constraints

Conductor selection for transmission and distribution lines is governed by two distinct engineering limits:

1. Thermal Ampacity (Current Carrying Limit): Determined by the conductor's maximum operating temperature (typically $75^\circ\text{C}$ for standard ACSR, up to $250^\circ\text{C}$ for high-temperature low-sag HTLS conductors). The steady-state thermal rating is computed via IEEE Standard 738 heat balance: $q_c + q_r = q_s + I^2 R(T_c)$ where $q_c$ is convective cooling rate, $q_r$ is radiative cooling rate, $q_s$ is solar heat gain, and $I^2 R$ is internal Joule heating.
2. Voltage Drop Limit (System Operability Limit): For long lines, voltage drop or steady-state stability limits are reached long before thermal ampacity is exceeded. Utility criteria generally restrict steady-state voltage variation to $\le \pm 5%$ of nominal nominal voltage under peak transfer.

---

5. Comprehensive Worked Calculation: Medium Line Analysis

Problem Statement

A 3-phase, $60\text{ Hz}$, $138\text{ kV}$ (line-to-line), $100\text{ km}$ transmission line supplies a balanced 3-phase load of $90\text{ MVA}$ at $0.85$ power factor lagging and rated voltage ($138\text{ kV}$). The line parameters per phase per km are:

• $r = 0.125; \Omega/\text{km}$
• $x = 0.400; \Omega/\text{km}$
• $y = j 3.20 \times 10^{-6}; \text{S/km}$

Calculate:

1. Total line series impedance $\mathbf{Z}$ and shunt admittance $\mathbf{Y}$.
2. Nominal $\pi$ $ABCD$ parameters.
3. Sending-end line-to-line voltage magnitude $|V_{S,LL}|$.
4. Sending-end active power $P_S$ and transmission efficiency $\eta%$.
5. Percent voltage regulation $VR%$.

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

Step 1: Compute Total Line Parameters Z and Y
  Total length l = 100 km
  Z = (0.125 + j0.400) * 100 = 12.5 + j40.0 ohms = 41.908 /_ 72.65 deg ohms
  Y = (j 3.20e-6) * 100 = j 3.20e-4 S = 3.20e-4 /_ 90.0 deg S

Step 2: Determine ABCD Parameters for Nominal Pi Model
  Z * Y = (41.908 /_ 72.65 deg) * (3.20e-4 /_ 90.0 deg) 
        = 0.01341 /_ 162.65 deg = -0.01280 + j0.00400
  
  A = D = 1 + (Z * Y)/2 = 1 + (-0.00640 + j0.00200) 
        = 0.99360 + j0.00200 = 0.99360 /_ 0.115 deg
  
  B = Z = 12.5 + j40.0 ohms = 41.908 /_ 72.65 deg ohms
  
  C = Y * (1 + Z*Y/4) 
    = (j 3.20e-4) * (1 - 0.00320 + j0.00100) 
    = (j 3.20e-4) * (0.99680 + j0.00100) 
    = -3.20e-7 + j 3.1898e-4 S = 3.1898e-4 /_ 90.06 deg S

Step 3: Establish Receiving-End Reference Phasors
  V_R,LL = 138,000 V => V_R,LN = 138,000 / sqrt(3) = 79,674.3 V /_ 0 deg
  Load S_3ph = 90 MVA at 0.85 lag => theta = arccos(0.85) = 31.788 deg
  I_R magnitude = S_3ph / (sqrt(3) * V_R,LL) = 90,000,000 / (sqrt(3) * 138,000) = 376.53 A
  I_R phasor = 376.53 /_ -31.788 deg A = 320.05 - j 198.37 A

Step 4: Compute Sending-End Voltage Phasor
  V_S,LN = A * V_R,LN + B * I_R
  Term 1: A * V_R,LN = (0.99360 + j0.00200) * 79,674.3 = 79,164.4 + j159.35 V
  Term 2: B * I_R = (41.908 /_ 72.65 deg) * (376.53 /_ -31.788 deg)
                  = 15,779.6 /_ 40.862 deg = 11,934.3 + j10,323.2 V
  V_S,LN = (79,164.4 + 11,934.3) + j(159.35 + 10,323.2)
         = 91,098.7 + j10,482.6 V = 91,699.6 /_ 6.566 deg V

  Sending-end Line-to-Line Voltage:
  |V_S,LL| = sqrt(3) * 91,699.6 V = 158,828 V = 158.83 kV

Step 5: Compute Sending-End Current, Power, and Efficiency
  I_S = C * V_R,LN + D * I_R
  Term 1: C * V_R,LN = (3.1898e-4 /_ 90.06 deg) * 79,674.3 = 25.415 /_ 90.06 deg A = -0.027 + j25.415 A
  Term 2: D * I_R = (0.99360 /_ 0.115 deg) * (376.53 /_ -31.788 deg) = 374.12 /_ -31.673 deg A
                  = 318.39 - j196.44 A
  I_S = 318.36 - j171.03 A = 361.42 /_ -28.245 deg A

  Sending-end power factor angle theta_S = 6.566 deg - (-28.245 deg) = 34.811 deg
  cos(theta_S) = cos(34.811 deg) = 0.8210 lagging
  
  P_S,3ph = sqrt(3) * |V_S,LL| * |I_S| * cos(theta_S)
          = sqrt(3) * 158,828 * 361.42 * 0.8210 = 81,595,000 W = 81.60 MW
  P_R,3ph = 90 MVA * 0.85 = 76.50 MW
  
  Transmission Efficiency:
  eta% = (P_R / P_S) * 100% = (76.50 MW / 81.60 MW) * 100% = 93.75%

Step 6: Compute Percent Voltage Regulation
  |V_R,NL| = |V_S,LN| / |A| = 91,699.6 / 0.99360 = 92,290.3 V
  |V_R,FL| = 79,674.3 V
  
  VR% = [ (92,290.3 - 79,674.3) / 79,674.3 ] * 100% 
      = (12,616.0 / 79,674.3) * 100% = 15.83%
===================================================================================

---

6. Common Exam Traps & Strategic Pitfalls

• The Per-Phase vs. Line-to-Line Voltage Confusion: Forgetting to divide line-to-line voltage by $\sqrt{3}$ when setting up $ABCD$ matrix calculations. All matrix equations operate strictly on line-to-neutral quantities.
• The No-Load Voltage Regulation Blunder: Setting $|V_{R,NL}| = |V_S|$ for a medium or long line instead of $|V_{R,NL}| = |V_S|/|A|$. Because $|A| < 1$ due to shunt capacitance, failing to divide by $|A|$ underestimates voltage regulation significantly.
• Inductive vs. Capacitive $GMR$ Swapping: Using $0.7788r$ for capacitance calculations. Capacitance involves electrostatic charge on the surface; use physical radius $r$, not $r'$.
• Approximate Voltage Drop Sign Inversion: Adding $X_L \sin\theta$ for leading power factors. For capacitive loads, the reactive drop is subtractive ($R \cos\theta - X_L \sin\theta$), which can lead to a voltage rise.