# Travelling Wave Line Model

## Introduction Figure 1. Travelling wave along a line (image courtesy of Wikipedia)

A travelling wave on a transmission line is a transient disturbance that that moves along the line at a constant speed yet maintains its shape (see Figure 1). Examples include lightning surges, switching transients, faults, etc.

In this article, we will derive various time-domain models for travelling waves on transmission lines, starting from the simplest case of the lossless single-phase line, then building to a generic lossy single-phase line. The derivation of the model is also looked at starting from the frequency-domain, showing the symmetry of the frequency and time-domain models. Finally, the considerations required to extend the model to a multi-conductor line is discussed.

This article provides the theoretical background for the travelling wave model. Actual implementations of the model suitable for use in computer simulations (i.e. discrete-time) can be found in the following articles:

## Lossless Single-Phase Lines Figure 2. Travelling wave line model - Small Segment from a lossless single-phase line

Consider the small segment of length $\Delta x$ metres from a lossless single-phase transmission line shown in Figure 2.

Notice that the model is similar to that of the classical distributed parameter line model, except that in this model, the voltages and currents are functions of both distance $x$ and time $t$ and are real, not phasor quantities. In the classical distributed parameter line model, the voltages and currents are evaluated at steady-state and are only functions of distance $x$ .

### Derivation of Voltage and Current Equations

Analysing this circuit using Kirchhoff's voltage law and Kirchhoff's current law, we can get the following pair of equations:

$V(x+\Delta x,t)=V(x,t)-L\Delta x{\frac {\partial I(x+\Delta x,t)}{\partial t}}\,$ $I(x,t)=I(x+\Delta x,t)+C\Delta x{\frac {\partial V(x,t)}{\partial t}}\,$ Notes:

1. The voltage across an inductor is $v=L{\frac {di}{dt}}$ and the current through a capacitor is $i=C{\frac {dv}{dt}}$ .
2. Partial derivatives are used because the voltages and currents are functions of both time and distance.

Re-arranging the above equations, we get:

${\frac {V(x+\Delta x,t)-V(x,t)}{\Delta x}}=-L{\frac {\partial I(x+\Delta x,t)}{\partial t}}\,$ ${\frac {I(x+\Delta x,t)-I(x,t)}{\Delta x}}=-C{\frac {\partial V(x,t)}{\partial t}}\,$ Taking the limit as $\Delta x\to 0\,$ yields the well known Telegrapher's equations:

${\frac {\partial V(x,t)}{\partial x}}=-L{\frac {\partial I(x,t)}{\partial t}}\,$ ... Equ. (1)
${\frac {\partial I(x,t)}{\partial x}}=-C{\frac {\partial V(x,t)}{\partial t}}\,$ ... Equ. (2)

Differentiating the voltage equation with respect to x and the current equation with respect to t:

${\frac {\partial ^{2}V(x,t)}{\partial x^{2}}}=-L{\frac {\partial ^{2}I(x,t)}{\partial x\partial t}}\,$ ... Equ. (3)
${\frac {\partial ^{2}I(x,t)}{\partial x\partial t}}=-C{\frac {\partial ^{2}V(x,t)}{\partial t^{2}}}\,$ ... Equ. (4)

Now substituting Equ. (4) into Equ. (3), we get:

${\frac {\partial ^{2}V(x,t)}{\partial x^{2}}}=LC{\frac {\partial ^{2}V(x,t)}{\partial t^{2}}}\,$ ... Equ. (5)

Similarly, we could have differentiated the current equation with respect to x and the voltage equation with respect to t and then solved for the current equation:

${\frac {\partial ^{2}I(x,t)}{\partial x^{2}}}=LC{\frac {\partial ^{2}I(x,t)}{\partial t^{2}}}\,$ ... Equ. (6)

The pair of equations (5) and (6) above are referred to as the transmission line wave equations.

The general solution to either of these equations can be found using d'Alembert's formula:

$I(x,t)=i^{+}(x-vt)+i^{-}(x+vt)\,$ $V(x,t)=v^{+}(x-vt)+v^{-}(x+vt)\,$ where $v={\frac {1}{\sqrt {LC}}}\,$ is the velocity of propagation (m/s).

$i^{+}(t)\,$ and $i^{-}(t)\,$ are arbitrary current functions of time.
$v^{+}(t)\,$ and $v^{-}(t)\,$ are arbitrary voltage functions of time.

Using the current equation above, we can solve for the voltage equation as follows (see here for the full derivation):

$v(x,t)=Z_{c}\left[i^{+}(x-vt)-i^{-}(x+vt)\right]\,$ where $Z_{c}={\frac {1}{vC}}={\frac {\sqrt {LC}}{C}}={\sqrt {\frac {L}{C}}}$ is the characteristic impedance (Ohms)

### Significance of + and - Equations

So we have the pair of voltage and current equations:

$I(x,t)=i^{+}(x-vt)+i^{-}(x+vt)\,$ $V(x,t)=Z_{c}\left[i^{+}(x-vt)-i^{-}(x+vt)\right]\,$ What is the significance of $i^{+}(x-vt)\,$ and $i^{-}(x+vt)\,$ ?

Consider that time t=0s, the function $i^{+}(x)\,$ represents an arbitrary waveform / shape distributed in space along the transmission line. For example, it could look like this: Figure 3. Arbitrary waveform $i^{+}(x)$ at time t = 0s

Now consider what happens at time t=1s, the function $i^{+}(x-v)\,$ has the same shape as the function above, but is translated to the right by $v\,$ metres: Figure 4. Arbitrary waveform $i^{+}(x-v)$ at time t = 1s

What we can conclude is that the function $i^{+}(x-vt)\,$ is an arbitrary waveform that maintains its shape but moves in the direction of the receiving end with increasing time at speed $v\,$ m/s. This is called a forward travelling wave.

We can use similar logic to conclude that the function $i^{-}(x+vt)\,$ is an arbitrary waveform that moves in the opposite direction (i.e. in the direction of the sending end). This is referred to as a backward travelling wave.

## Generic Lossy Single-Phase Lines

Consider the generic single-phase distributed parameter line model of length $l\,$ metres shown in Figure 5. Analysing the circuit using Kirchhoff's laws as we did previously in the lossless case, we can get:

${\frac {\partial V(x,t)}{\partial x}}=-I(x,t)R-L{\frac {\partial I(x,t)}{\partial t}}\,$ ${\frac {\partial I(x,t)}{\partial x}}=-V(x,t)G-C{\frac {\partial V(x,t)}{\partial t}}\,$ Where $G={\frac {1}{R_{sh}}}\,$ is the shunt conductance.

Differentiating the voltage equation with respect to x and the current equation with respect to t:

${\frac {\partial ^{2}V(x,t)}{\partial x^{2}}}=-R{\frac {\partial I(x,t)}{\partial x}}-L{\frac {\partial ^{2}I(x,t)}{\partial x\partial t}}\,$ ... Equ. (7)
${\frac {\partial ^{2}I(x,t)}{\partial x\partial t}}=-G{\frac {\partial V(x,t)}{\partial t}}-C{\frac {\partial ^{2}V(x,t)}{\partial t^{2}}}\,$ ... Equ. (8)

Substituting Equ. (8) into Equ. (7), we get:

${\frac {\partial ^{2}V(x,t)}{\partial x^{2}}}=-R{\frac {\partial I(x,t)}{\partial x}}+LG{\frac {\partial V(x,t)}{\partial t}}+LC{\frac {\partial ^{2}V(x,t)}{\partial t^{2}}}\,$ ... Equ. (9)

If we were to set R = G = 0 in the equation above, then we would get the same transmission line wave equation as Equ. (5) in the lossless case.

However, while there is a general solution to the lossless line using d'Alembert's formula (as shown earlier), the generic lossy line does not have a neat, closed-form solution.

## Frequency Domain Derivation

In the preceding sections, the travelling wave model was developed purely in the time-domain. Another way of looking at the travelling wave model is to derive it first in the frequency domain and then convert it to the time-domain. Recall from the steady-state (frequency domain) distributed parameter model that the sending end voltage and current can be related to the receiving end quantities as follows:

$\left[{\begin{matrix}{\boldsymbol {V_{s}}}\\\\{\boldsymbol {I_{s}}}\end{matrix}}\right]=\left[{\begin{matrix}\cosh({\boldsymbol {\gamma }}l)&{\boldsymbol {Z}}_{c}\sinh({\boldsymbol {\gamma }}l)\\\\{\frac {1}{{\boldsymbol {Z}}_{c}}}sinh({\boldsymbol {\gamma }}l)&\cosh({\boldsymbol {\gamma }}l)\end{matrix}}\right]\left[{\begin{matrix}{\boldsymbol {V_{r}}}\\\\{\boldsymbol {I_{r}}}\end{matrix}}\right]\,$ where ${\boldsymbol {\gamma }}={\sqrt {\boldsymbol {zy}}}$ is the propagation constant ($m^{-1}$ )

${\boldsymbol {Z}}_{c}={\sqrt {\boldsymbol {\frac {z}{y}}}}$ is the characteristic impedance ($\Omega$ )
${\boldsymbol {z}}=R+j\omega l$ and ${\boldsymbol {y}}=G+j\omega C$ are the series impedance and shunt admittance respectively

With some straightforward algebraic manipulation, the steady-state equations above can be represented in the following form:

${\boldsymbol {V_{r}}}+{\boldsymbol {Z}}_{c}{\boldsymbol {I_{r}}}=({\boldsymbol {V_{s}}}+{\boldsymbol {Z}}_{c}{\boldsymbol {I_{s}}})e^{-{\boldsymbol {\gamma }}l}$ ... Equ. (10)

The time-domain model can be derived by calculating the inverse Fourier transform of Equ. (10).

In the lossless case (i.e. R = G = 0.), the propagation constant and characteristic impedance are as follows:

${\boldsymbol {\gamma }}={\sqrt {\boldsymbol {zy}}}={\sqrt {(j\omega L)(j\omega C)}}=j\omega {\sqrt {LC}}$ ${\boldsymbol {Z}}_{c}={\sqrt {\boldsymbol {\frac {z}{y}}}}={\sqrt {\frac {(j\omega L)}{(j\omega C)}}}={\sqrt {\frac {L}{C}}}$ A property of Fourier transforms is that multiplying by a complex exponential in the frequency domain results in a time shift in the time domain. Therefore, taking the inverse Fourier transform of Equ. (10) for the lossless case, we get:

$v_{r}(t)+{\sqrt {\frac {L}{C}}}i_{r}(t)=v_{s}(t-{\sqrt {LC}}l)+{\sqrt {\frac {L}{C}}}i_{s}(t-{\sqrt {LC}}l)\,$ The lossless case above shows that the receiving end quantities are essentially time-shifted (past) versions of the sending end quantities, as was derived earlier using d'Alembert's formula.

In the general lossy case, the characteristic impedance is not a scalar constant (but rather a function of frequency ${\boldsymbol {Z}}_{c}(\omega )\,$ ) and the propagation constant does not lead to a simple time shift in the time-domain. Therefore, the inverse Fourier transform of Equ. (10) leads to a series of messy convolutions:

$v_{r}(t)+z_{c}(t)*i_{r}(t)=\left[v_{s}(t)+z_{c}(t)*i_{s}(t)\right]*a(t)\,$ where $a(t)\,$ is the inverse Fourier transform of $e^{-{\boldsymbol {\gamma }}(\omega )l}\,$ $z_{c}(t)\,$ is the inverse Fourier transform of ${\boldsymbol {Z}}_{c}(\omega )\,$ As can be imagined, the convolution integrals in the equation above are not neatly solvable.

## Multi-conductor Lines

The travelling wave line model for the single-phase lossless line can be extended to a multi-conductor line of n conductors by replacing the voltage and current quantities with $n\times 1$ vectors, e.g. for a three-phase, three-conductor line:

${\boldsymbol {V}}(x,t)=\left[{\begin{matrix}{\boldsymbol {V_{a}}}(x,t)\\{\boldsymbol {V_{b}}}(x,t)\\{\boldsymbol {V_{c}}}(x,t)\end{matrix}}\right]\,$ , ${\boldsymbol {I}}(x,t)=\left[{\begin{matrix}{\boldsymbol {I_{a}}}(x,t)\\{\boldsymbol {I_{b}}}(x,t)\\{\boldsymbol {I_{c}}}(x,t)\end{matrix}}\right]\,$ And the inductance and capacitance (per unit length) with $n\times n$ matrices to capture the mutual coupling between phases, e.g. for a three-phase, three-conductor line:

$[L]=\left[{\begin{matrix}L_{aa}&L_{ab}&L_{ac}\\L_{ba}&L_{bb}&L_{bc}\\L_{ca}&L_{cb}&L_{cc}\end{matrix}}\right]\,$ $[C]=\left[{\begin{matrix}C_{aa}&C_{ab}&C_{ac}\\C_{ba}&C_{bb}&C_{bc}\\C_{ca}&C_{cb}&C_{cc}\end{matrix}}\right]\,$ The transmission line wave equations for the single-phase model (see equations (5) and (6) above) can now be re-written for the multi-conductor model as follows:

${\frac {\partial ^{2}{\boldsymbol {V}}(x,t)}{\partial x^{2}}}=[L][C]{\frac {\partial ^{2}{\boldsymbol {V}}(x,t)}{\partial t^{2}}}\,$ ... Equ. (10)
${\frac {\partial ^{2}{\boldsymbol {I}}(x,t)}{\partial x^{2}}}=[C][L]{\frac {\partial ^{2}{\boldsymbol {I}}(x,t)}{\partial t^{2}}}\,$ ... Equ. (11)

Unlike the single-phase case, there is no general solution to the second-order differential equations (10) and (11) since the matrix product [L][C] is typically full (i.e. off-diagonals are non-zero). To solve these differential equations, a modal transformation is required to decouple the phases (refer to this section of the distributed parameter line article for more details on the modal transformation).

The voltage and current vectors are transformed from the phase domain to the modal domain by the $n\times n$ transformation matrices $[T_{v}]\,$ and $[T_{i}]\,$ , such that:

${\boldsymbol {V}}(x,t)=[T_{v}]{\boldsymbol {V'}}(x,t)$ ${\boldsymbol {I}}(x,t)=[T_{i}]{\boldsymbol {I'}}(x,t)$ where ${\boldsymbol {V'}}(x,t)$ and ${\boldsymbol {I'}}(x,t)$ are the modal voltage and current vectors, i.e.

${\boldsymbol {V'}}(x,t)=\left[{\begin{matrix}{\boldsymbol {V_{0}}}(x)\\{\boldsymbol {V_{1}}}(x)\\{\boldsymbol {V_{2}}}(x)\end{matrix}}\right]\,$ , ${\boldsymbol {I'}}(x,t)=\left[{\begin{matrix}{\boldsymbol {I_{0}}}(x)\\{\boldsymbol {I_{1}}}(x)\\{\boldsymbol {I_{2}}}(x)\end{matrix}}\right]\,$ The transformation matrices $[T_{v}]\,$ and $[T_{i}]\,$ are selected such that they are the eigenvectors of [L][C] and [C][L] respectively (note that they have the same eigenvalues).

In the modal domain, equations (10) and (11) can be re-written as sets of decoupled equations, i.e.

$\left[{\begin{matrix}{\frac {d^{2}{\boldsymbol {V_{0}}}(x,t)}{dx^{2}}}\\\\{\frac {d^{2}{\boldsymbol {V_{1}}}(x,t)}{dx^{2}}}\\\\{\frac {d^{2}{\boldsymbol {V_{2}}}(x,t)}{dx^{2}}}\end{matrix}}\right]=\left[{\begin{matrix}\lambda _{0}&&\\&\lambda _{1}&\\&&\lambda _{2}\end{matrix}}\right]\left[{\begin{matrix}{\frac {d^{2}{\boldsymbol {V_{0}}}(x,t)}{dt^{2}}}\\\\{\frac {d^{2}{\boldsymbol {V_{1}}}(x,t)}{dt^{2}}}\\\\{\frac {d^{2}{\boldsymbol {V_{2}}}(x,t)}{dt^{2}}}\end{matrix}}\right]\,$ ... Equ. (12)
$\left[{\begin{matrix}{\frac {d^{2}{\boldsymbol {I_{0}}}(x,t)}{dx^{2}}}\\\\{\frac {d^{2}{\boldsymbol {I_{1}}}(x,t)}{dx^{2}}}\\\\{\frac {d^{2}{\boldsymbol {I_{2}}}(x,t)}{dx^{2}}}\end{matrix}}\right]=\left[{\begin{matrix}\lambda _{0}&&\\&\lambda _{1}&\\&&\lambda _{2}\end{matrix}}\right]\left[{\begin{matrix}{\frac {d^{2}{\boldsymbol {I_{0}}}(x,t)}{dt^{2}}}\\\\{\frac {d^{2}{\boldsymbol {I_{1}}}(x,t)}{dt^{2}}}\\\\{\frac {d^{2}{\boldsymbol {I_{2}}}(x,t)}{dt^{2}}}\end{matrix}}\right]\,$ ... Equ. (13)

These modal differential equations can then be solved individually using d'Alembert's formula. For example, the solution for mode 0 is:

$I_{0}(x,t)=i_{0}^{+}(x-v_{0}t)+i_{0}^{-}(x+v_{0}t)\,$ $V_{0}(x,t)=v_{0}^{+}(x-v_{0}t)+v_{0}^{-}(x+v_{0}t)\,$ where $v_{0}={\frac {1}{\sqrt {\lambda _{0}}}}\,$ is the velocity of propagation for mode 0 (m/s)

$i^{+}(t)\,$ and $i^{-}(t)\,$ are current functions for mode 0
$v^{+}(t)\,$ and $v^{-}(t)\,$ are voltage functions for mode 0

Once modal domain quantities are calculated, then the phase domain quantities can be found by applying the inverse transformation.

Finally, it should be noted that the general solution using the modal transformation only applies to lossless multi-conductor lines. The generic lossy line suffers from the same issues as described earlier for the lossy single-phase line, i.e. multiple cross-coupling terms that render d'Alembert's formula unusable.

## Frequency Dependency of Line Parameters

So far, we have described the per unit length inductance (L), capacitance (C), resistance (R) and conductance (G) parameters of the line as fixed and constant for the full range of frequencies. However, this is not a valid assumption, particularly at high frequencies. For example, the skin effect leads to higher resistances at higher frequencies. It is therefore more correct to characterise the line parameters as having frequency dependencies, i.e. $L(\omega )\,$ , $C(\omega )\,$ , $R(\omega )\,$ and $G(\omega )\,$ .