# 6.3: Parallel Plate Waveguide- TE Case, Electric Field

- Page ID
- 24813

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

In Section 6.2, the parallel plate waveguide was introduced. At the end of that section, we described the decomposition of the problem into its TE and TM components. In this section, we find the electric field component of the TE field in the waveguide.

Figure \(\PageIndex{1}\) shows the problem addressed in this section. (Additional details and assumptions are addressed in Section 6.2.)

Figure \(\PageIndex{1}\): TE component of the electric field in a parallel plate waveguide. ( CC BY-SA 4.0; C. Wang)Since \(\widetilde{E}_x=\widetilde{E}_z=0\) for the TE component of the electric field, Equations 6.2.11 and 6.2.13 are irrelevant, leaving only:

\[\frac{\partial^2}{\partial x^2}\widetilde{E}_y + \frac{\partial^2}{\partial z^2}\widetilde{E}_y = - \beta^2 \widetilde{E}_y \label{m0174_eDE} \]

The general solution to this partial differential equation is:

\[\begin{align} \widetilde{E}_y =&~~~~~e^{-jk_z z} \left[ A e^{-jk_x x} + B e^{+jk_x x} \right] \nonumber \\ &+e^{+jk_z z} \left[ C e^{-jk_x x} + D e^{+jk_x x} \right] \label{m0174_eGS}\end{align} \]

where \(A\), \(B\), \(C\), and \(D\) are complex-valued constants and \(k_x\) and \(k_z\) are real-valued constants. We have assigned variable names to these constants with advance knowledge of their physical interpretation; however, at this moment they remain simply unknown constants whose values must be determined by enforcement of boundary conditions.

Note that Equation \ref{m0174_eGS} consists of two terms. The first term includes the factor \(e^{-jk_z z}\), indicating a wave propagating in the \(+\hat{\bf z}\) direction, and the second term includes the factor \(e^{+jk_z z}\), indicating a wave propagating in the \(-\hat{\bf z}\) direction. If we impose the restriction that sources exist only on the left (\(z<0\)) side of Figure \(\PageIndex{1}\), and that there be no structure capable of wave scattering (in particular, reflection) on the right (\(z>0\)) side of Figure \(\PageIndex{1}\), then there can be no wave components propagating in the \(-\hat{\bf z}\) direction. In this case, \(C=D=0\) and Equation \ref{m0174_eGS} simplifies to:

\[\widetilde{E}_y = e^{-jk_z z} \left[ A e^{-jk_x x} + B e^{+jk_x x} \right] \label{m0174_eGS2} \]

Before proceeding, let’s make sure that Equation \ref{m0174_eGS2} is actually a solution to Equation \ref{m0174_eDE}. As we shall see in a moment, performing this check will reveal some additional useful information. First, note:

\[\frac{\partial \widetilde{E}_y}{\partial x} = e^{-jk_z z} \left[-A e^{-jk_x x} + B e^{+jk_x x} \right]\left(jk_x\right) \nonumber \]

So:

\[\frac{\partial^2 \widetilde{E}_y}{\partial x^2} = e^{-jk_z z} \left[ A e^{-jk_x x} + B e^{+jk_x x} \right]\left(-k_x^2\right) \nonumber \]

Comparing this to Equation \ref{m0174_eGS2}, we observe the remarkable fact that

\[\frac{\partial^2 \widetilde{E}_y}{\partial x^2} = -k_x^2 \widetilde{E}_y \nonumber \]

Similarly, note:

\[\frac{\partial \widetilde{E}_y}{\partial z} = e^{-jk_z z} \left[ A e^{-jk_x x} + B e^{+jk_x x} \right]\left(-jk_z\right) \nonumber \]

So:

\[\begin{align} \frac{\partial^2 \widetilde{E}_y}{\partial z^2} &= e^{-jk_z z} \left[ A e^{-jk_x x} + B e^{+jk_x x} \right]\left(-k_z^2\right) \nonumber \\ &= -k_z^2 \widetilde{E}_y\end{align} \nonumber \]

Now summing these results:

\[\frac{\partial^2 \widetilde{E}_y}{\partial x^2} + \frac{\partial^2 \widetilde{E}_y}{\partial z^2} = -\left( k_x^2 + k_z^2 \right) \widetilde{E}_y \label{m0174_eDE2} \]

Comparing Equation \ref{m0174_eDE2} to Equation \ref{m0174_eDE}, we conclude that Equation \ref{m0174_eGS2} is indeed a solution to Equation \ref{m0174_eDE}, but only if:

\[\beta^2 = k_x^2 + k_z^2 \label{m0174_eBeta} \]

This confirms that \(k_x\) and \(k_z\) are in fact the components of the propagation vector

\[{\bf k} \triangleq \beta\hat{\bf k} = \hat{\bf x}k_x + \hat{\bf y}k_y + \hat{\bf z}k_z \nonumber \]

where \(\hat{\bf k}\) is the unit vector pointing in the direction of propagation, and \(k_y=0\) in this particular problem.

The solution has now been reduced to finding the constants \(A\), \(B\), and either \(k_x\) or \(k_z\). This is accomplished by enforcing the relevant boundary conditions. In general, the component of \(\widetilde{\bf E}\) that is tangent to a perfectly-conducting surface is zero. Applied to the present problem, this means \(\widetilde{E}_y = 0\) at \(x=0\) and \(\widetilde{E}_y = 0\) at \(x=a\). Referring to Equation \ref{m0174_eGS2}, the boundary condition at \(x=0\) means

\[e^{-jk_z z} \left[ A \cdot 1 + B \cdot 1 \right] = 0 \nonumber \]

The factor \(e^{-jk_z z}\) cannot be zero; therefore, \(A+B=0\). Since \(B=-A\), we may rewrite Equation \ref{m0174_eGS2} as follows:

\[\widetilde{E}_y = e^{-jk_z z} B \left[ e^{+jk_x x} - e^{-jk_x x} \right] \nonumber \]

This expression is simplified using a trigonometric identity:

\[\frac{1}{2j}\left[ e^{+jk_x x} - e^{-jk_x x} \right] = \sin{k_x x} \nonumber \]

Let us now make the definition \(E_{y0}\triangleq j2B\). Then:

\[\widetilde{E}_y = E_{y0} e^{-jk_z z} \sin k_x x \label{m0174_eGS3} \]

Now applying the boundary condition at \(x=a\):

\[E_{y0} e^{-jk_z z} \sin k_x a = 0 \nonumber \]

The factor \(e^{-jk_z z}\) cannot be zero, and \(E_{y0}=0\) yields only trivial solutions; therefore:

\[\sin{k_x a} = 0 \nonumber \]

This in turn requires that

\[k_x a = m \pi \label{m0174_ekxa} \]

where \(m\) is an integer. Note that \(m=0\) is not of interest since this yields \(k_x=0\), which according to Equation \ref{m0174_eGS3} yields the trivial solution \(\widetilde{E}_y=0\). Also each integer value of \(m\) that is less than zero is excluded because the associated solution is different from the solution for the corresponding positive value of \(m\) in sign only, which can be absorbed in the arbitrary constant \(E_{y0}\).

At this point we have uncovered a family of solutions given by Equation \ref{m0174_eGS3} and Equation \ref{m0174_ekxa} with \(m=1,2,...\). Each solution associated with a particular value of \(m\) is referred to as a *mode*, which (via Equation \ref{m0174_ekxa}) has a particular value of \(k_x\). The value of \(k_z\) for mode \(m\) is obtained using Equation \ref{m0174_eBeta} as follows:

\[\begin{align} k_z &= \sqrt{\beta^2-k_x^2} \nonumber \\ & =\sqrt{\beta^2-\left(\frac{m\pi}{a}\right)^2}\end{align} \nonumber \]

Since \(k_z\) is specified to be real-valued, we require:

\[\beta^2-\left(\frac{m\pi}{a}\right)^2 > 0 \nonumber \]

This constrains \(\beta\); specifically:

\[\beta > \frac{m\pi}{a} \nonumber \]

Recall that \(\beta=\omega\sqrt{\mu\epsilon}\) and \(\omega=2\pi f\) where \(f\) is frequency. Solving for \(f\), we find:

\[f > \frac{m}{2a\sqrt{\mu\epsilon}} \nonumber \]

Therefore, each mode exists only above a certain frequency, which is different for each mode. This *cutoff frequency* \(f_c\) for mode \(m\) is given by

\[\boxed{ f_c^{(m)} \triangleq \frac{m}{2a\sqrt{\mu\epsilon}} } \label{m0174_efcm} \]

At frequencies below the cutoff frequency for mode \(m\), modes \(1\) through \(m-1\) exhibit imaginary-valued \(k_z\). The propagation constant must have a real-valued component in order to propagate; therefore, these modes do not propagate and may be ignored.

Let us now summarize the solution. For the scenario depicted in Figure \(\PageIndex{1}\), the electric field component of the TE solution is given by:

\[\boxed{ \hat{\bf y}\widetilde{E}_y = \hat{\bf y}\sum_{m=1}^{\infty} \widetilde{E}_y^{(m)} } \label{m0174_eEysum} \]

where

\[\boxed{ \widetilde{E}_y^{(m)} \triangleq \begin{cases} 0, & f<f_c^{(m)} \\ E_{y0}^{(m)} e^{-jk_z^{(m)} z} \sin k_x^{(m)} x, & f\ge f_c^{(m)} \end{cases} } \nonumber \]

where \(m\) enumerates modes (\(m=1,2,...\)),

\[\boxed{ k_z^{(m)} \triangleq \sqrt{\beta^2-\left[k_x^{(m)}\right]^2} } \nonumber \]

and

\[\boxed{ k_x^{(m)} \triangleq m\pi/a } \label{m0174_ekxma} \]

Finally, \(E_{y0}^{(m)}\) is a complex-valued constant that depends on sources or boundary conditions to the left of the region of interest.

For the scenario depicted in Figure \(\PageIndex{1}\), the electric field component of the TE solution is given by Equation \ref{m0174_eEysum} with modal components determined as indicated by Equations \ref{m0174_efcm}-\ref{m0174_ekxma}. This solution presumes all sources lie to the left of the region of interest, and no scattering occurs to the right of the region of interest.

To better understand this result, let us examine the lowest-order mode, \(m=1\). For this mode \(f_c^{(1)}=1/2a\sqrt{\mu\epsilon}\), so this mode can exist if \(f>1/2a\sqrt{\mu\epsilon}\). Also \(k_x^{(1)} = \pi/a\), so

\[k_z^{(1)} = \sqrt{\beta^2-\left(\frac{\pi}{a}\right)^2} \label{m0174_ekz1} \]

Subsequently,

\[\widetilde{E}_y^{(1)} = E_{y0}^{(1)} e^{-jk_z^{(1)} z} \sin \frac{\pi x}{a} \nonumber \]

Note that this mode has the form of a plane wave. The plane wave propagates in the \(+\hat{\bf z}\) direction with phase propagation constant \(k_z^{(1)}\). Also, we observe that the apparent plane wave is non-uniform, exhibiting magnitude proportional to \(\sin \pi x/a\) within the waveguide. This is shown in Figure \(\PageIndex{2}\) (left image).

Figure \(\PageIndex{2}\): Magnitude of the two lowest-order TE modes in a parallel plate waveguide.*Left:*\(m=1\),

*Right:*\(m=2\). ( CC BY-SA 4.0; C. Wang)

In particular, we observe that the magnitude of the wave is zero at the perfectly-conducting surfaces – as is necessary to satisfy the boundary conditions – and is maximum in the center of the waveguide.

Now let us examine the \(m=2\) mode. For this mode, \(f_c^{(2)}=1/a\sqrt{\mu\epsilon}\), so this mode can exist if \(f>1/a\sqrt{\mu\epsilon}\). This frequency is higher than \(f_c^{(1)}\), so the \(m=1\) mode can exist at any frequency at which the \(m=2\) mode exists. Also \(k_x^{(2)} = 2\pi/a\), so

\[k_z^{(2)} = \sqrt{\beta^2-\left(\frac{2\pi}{a}\right)^2} \label{m0174_ekz2} \]

Subsequently,

\[\widetilde{E}_y^{(2)} = E_{y0}^{(2)} e^{-jk_z^{(2)} z} \sin \frac{2\pi x}{a} \nonumber \]

In this case, the apparent plane wave propagates in the \(+\hat{\bf z}\) direction with phase propagation constant \(k_z^{(2)}\), which is less than \(k_z^{(1)}\). For \(m=2\), we find magnitude is proportional to \(\sin 2\pi x/a\) within the waveguide (Figure \(\PageIndex{2}\), right image). As in the \(m=1\) case, we observe that the magnitude of the wave is zero at the PEC surfaces; however, for \(m=2\), there are *two* maxima with respect to \(x\), and the magnitude in the center of the waveguide is zero.

This pattern continues for higher-order modes. In particular, each successive mode exhibits higher cutoff frequency, smaller propagation constant, and increasing integer number of sinusoidal half-periods in magnitude.

Consider an air-filled parallel plate waveguide consisting of plates separated by \(1\) cm. Determine the frequency range for which one (and only one) propagating TE mode is assured.

###### Solution

Single-mode TE propagation is assured by limiting frequency \(f\) to greater than the cutoff frequency for \(m=1\), but lower than the cutoff frequency for \(m=2\). (Any frequency higher than the cutoff frequency for \(m=2\) allows at least 2 modes to exist.) Calculating the applicable cutoff frequencies, we find:

\[\begin{aligned} f_c^{(1)} &= \frac{1}{2a\sqrt{\mu_0\epsilon_0}} \cong 15.0 \: ~\mbox{GHz} \\ f_c^{(2)} &= \frac{2}{2a\sqrt{\mu_0\epsilon_0}} \cong 30.0 \:~\mbox{GHz} \end{aligned} \nonumber \]

Therefore, \(15.0 \: \mathrm{GHz} \leq f \leq 30.0 \: \mathrm{GHz}\).

Finally, let us consider the phase velocity \(v_p\) within the waveguide. For lowest-order mode \(m=1\), this is

\[v_p = \frac{\omega}{k_z^{(1)}} = \frac{\omega}{\sqrt{\omega^2\mu\epsilon-\left(\frac{\pi}{a}\right)^2}} \label{m0174_evp} \]

Recall that the speed of an electromagnetic wave in unbounded space (i.e., not in a waveguide) is \(1/\sqrt{\mu\epsilon}\). For example, the speed of light in free space is \(1/\sqrt{\mu_0\epsilon_0}=c\). However, the phase velocity indicated by Equation \ref{m0174_evp} is *greater* than \(1/\sqrt{\mu\epsilon}\); e.g., faster than light would travel in the same material (presuming it were transparent). At first glance, this may seem to be impossible. However, recall that *information* travels at the group velocity \(v_g\), and not necessarily the phase velocity. (See Section 6.1 for a refresher.) Although we shall not demonstrate this here, the group velocity in the parallel plate waveguide is always *less* than \(1/\sqrt{\mu\epsilon}\), so no physical laws are broken, and signals travel somewhat slower than the speed of light, as they do in any other structure used to convey signals.

Also remarkable is that the speed of propagation is different for each mode. In fact, we find that the phase velocity increases and the group velocity decreases as \(m\) increases. This phenomenon is known as *dispersion*, and sometimes specifically as *mode dispersion* or *modal dispersion*.