# Transmission Line Losses on a Loaded Lossy Line¶

When dealing with RF power, for instance in radio, industry or scientific applications, a recurrent problem is to handle the inevitable RF losses correctly to avoid overheating of cables and components.

In this example, we will use scikit-rf to evaluate the losses in a 50 Ohm, 20 meters long, RG-8 cable (VF=0.84) charged with a resistive load $$R_L=50\Omega$$ at 13.56 MHz. The cable losses are estimated to 1.483 dB/100 meters and the source power to 400W.

First, the usual Python imports:

[1]:

%matplotlib inline

[2]:

import matplotlib.pyplot as plt
import numpy as np
import skrf as rf

[3]:

rf.stylely()


Let’s define the problem constants:

[4]:

Pin = 400  # W
z0 = 50 # Ohm
freq = rf.Frequency(13.56, npoints=1, unit='MHz')
VF = 0.84
RL = 50  # Ohm
L = 20  # m


The propagation constant of the transmission line $$\gamma=\alpha+j\beta$$ is:

[5]:

alpha = rf.db_2_np(1.483/100)  # Np/m
beta = freq.w/rf.c/VF
gamma = alpha + 1j*beta


The matched line loss (or power attenuation), $$a=e^{2\alpha L}$$, is:

[6]:

a = np.exp(2*alpha*L)  # also simply 2.84/100*20
print('Matched line loss: a=', rf.mag_2_db10(a), 'dB')

Matched line loss: a= 0.29660000000000003 dB


If the line is charged with a matched load, ie. $$R_L=50\Omega$$, then the total line loss is $$a$$. The power lost in the cable will thus be:

[7]:

print('(Forward) Power delivered to the  load:', Pin/a, 'W')
print('Power lost in the cable:', Pin *( 1 - 1/a), 'W')

(Forward) Power delivered to the  load: 373.59408479883416 W
Power lost in the cable: 26.405915201165886 W


Which can also be checked by the scikit-rf transmission line function zl_2_total_loss:

[8]:

a_skrf = rf.zl_2_total_loss(z0, zl=RL, theta=gamma*L)
print('Power lost in the cable:',
Pin * (1 - 1/a_skrf), 'W')

Power lost in the cable: [26.4059152] W


Another way to evaluate the total power dissipated in the circuit is to evaluate the power expression:

$P_{delivered} = \frac{1}{2} \Re \left[ V I^* \right]$

where $$V$$ and $$I$$ are the (peak) total voltage and current. They can be evaluate using the transmission line function voltage_current_propagation:

[9]:

# reflection coefficient and input impedance
Gamma_in = rf.zl_2_Gamma_in(z0, RL, theta=gamma*L)
Z_in = rf.zl_2_zin(z0, RL, theta=gamma*L)

# voltage and current at the line input as a function of source power
V_in = np.sqrt(2*z0*Pin)*(1 + Gamma_in)
I_in = V_in/Z_in

# voltage and current at z=L
V,I = rf.voltage_current_propagation(V_in, I_in, z0, gamma*L)
P_delivered = 1/2 * np.real(V * np.conj(I))
print('Power delivered to the load: ', P_delivered, 'W')
print('Power dissipated in the cable: ',Pin - P_delivered, 'W')

Power delivered to the load:  [373.5940848] W
Power dissipated in the cable:  [26.4059152] W


However, if the load is not perfectly matched to the line characteristic impedance $$z_0$$, for example with $$R_L=200 + 30j\Omega$$, additional losses are induced by the reflected wave. The reflection coefficient $$\Gamma_{load}$$ induced by this load is:

[10]:

z0 = 50
ZL = 200 - 30j

|Gamma_load|= [0.6075238]


while the reflection coefficient seen at the input of the transmission line $$\Gamma_{in}$$ is:

[11]:

Gamma_in = rf.zl_2_Gamma_in(z0, ZL, theta=gamma*L)
SWR = rf.Gamma0_2_swr(rf.zl_2_Gamma_in(z0, ZL, theta=gamma*L))
print('|Gamma_in|=', np.abs(Gamma_in), '(SWR=', SWR,')')

|Gamma_in|= [0.56741825] (SWR= [3.62340355] )


The total loss in dB due to SWR is often stated as:

$a_{[dB]} + 10 \log_{10} \frac{1 - |\Gamma_{in}|^2}{1 - |\Gamma_{load}|^2}$
[12]:

10*np.log10(a) + 10*np.log10((1 - np.abs(Gamma_in)**2)/(1 - np.abs(Gamma_load)**2))

[12]:

array([0.60942359])


However, this expression is only correct if either properties are verified: - (i) the characteristic impedance of the line is real (distorsionless line) - (ii) reflection coefficient is real (ie real $$Z_L/Z_0$$) [1].

The 1st condition is met here, however it will not be the case in the next section.

Anyway, the scikit-rf transmission line function zl_2_total_loss is correct in all conditions:

[13]:

a = rf.zl_2_total_loss(z0, zl=ZL, theta=gamma * L)
print('Total power loss: ', rf.mag_2_db10(a), 'dB' )
print('Delivered power:', Pin/a, 'W')
print('The total power loss is the cable:', Pin*(1 - 1/a), 'W')

Total power loss:  [0.60942359] dB
Delivered power: [347.63030721] W
The total power loss is the cable: [52.36969279] W

[14]:

# reflection coefficient and input impedance
Gamma_in = rf.zl_2_Gamma_in(z0, ZL, theta=gamma*L)
Z_in = rf.zl_2_zin(z0, ZL, theta=gamma*L)

# voltage and current at the line input as a function of source power
V_in = np.sqrt(2*z0*Pin)*(1 + Gamma_in)
I_in = V_in/Z_in

# voltage and current at z=L
V,I = rf.voltage_current_propagation(V_in, I_in, z0, gamma*L)
P_delivered = 1/2 * np.real(V * np.conj(I))
print('Power delivered to the load: ', P_delivered, 'W')
print('Power dissipated in the cable: ',Pin - P_delivered, 'W')

Power delivered to the load:  [235.70604719] W
Power dissipated in the cable:  [164.29395281] W

[15]:

rf.Gamma0_2_swr(Gamma_in)

[15]:

array([3.62340355])

[16]:

10*np.log10(P_delivered/Pin)

[16]:

array([-2.29689267])


This example reproduces the example presented in reference [1].

Let’s assume a coaxial line (Wireman #551, 450-Ohm) loaded with a complex load $$Z_L=R_L + jX_L$$, with the following parameters:
- line length: 100 feet - frequency: 1.83 MHz - attenuation constant: $$\alpha=$$ 0.095 dB/100 feet - coaxial relative permittivity: $$\epsilon_r=1.194$$ (Velocity Factor VF=0.915) - real part of the characteristic impedance: $$R_0 = \Re \left[Z_0\right]$$=402.75 Ohm - Load resistance: $$R_L$$ = 4.5 Ohm - Load reactance: $$X_L$$ = -1673 Ohm
[17]:

Z_L = 4.5 - 1673j
R_0 = 402.75
freq = rf.Frequency(1.83, npoints=1, unit='MHz')
VF = 0.915
L = rf.feet_2_meter(100)


First, we can derive the propagation constant $$\gamma=\alpha+j\beta$$ with $$\beta=\frac{\omega}{c }\sqrt{\epsilon_r}$$ from the problem parameters:

[18]:

alpha = rf.db_2_np(rf.db_per_100feet_2_db_per_100meter(0.095)/100)
beta = freq.w/rf.c/VF
gamma = alpha + 1j*beta
print(gamma)

[0.00035883+0.0419169j]

However, the transmission line characteristic reactance is not given in the problem parameters and must be determined. It can be approximated from a high-frequency, low-loss approximation [1]:
$Z_0 = R_0 + j X_0 \approx R_0 - j \frac{\alpha}{\beta}R_0$
ie,
$X_0 \approx - \frac{\alpha}{\beta}R_0$
[19]:

X_0 = -alpha/beta*R_0
Z_0 = R_0 + 1j*X_0
print('X_0=', X_0)

X_0= [-3.44778935]


Now that we have both the characteristic impedance and the propagation constant of the line, the reflection coefficients, input impedance and total loss can be deduced:

[20]:

print('Gamma at load:', np.abs(rf.zl_2_Gamma0(Z_0, Z_L)))
print('Gamma at input:', np.abs(rf.zl_2_Gamma_in(Z_0, Z_L, theta=gamma*L)))

print('SWR at input:', rf.Gamma0_2_swr(rf.zl_2_Gamma_in(Z_0, Z_L, theta=gamma*L)))
print('Input impedance:', rf.input_impedance_at_theta(Z_0, Z_L, theta=gamma*L ), 'Ohm')

total_loss_db = rf.mag_2_db10(np.abs(rf.zl_2_total_loss(z0=Z_0, zl=Z_L, theta=gamma*L)))
print('Total loss:', total_loss_db, 'dB')

Gamma at load: [0.99489308]
Gamma at input: [0.97336654]

[ ]: