Skip to content

9.5 - Guided Waves

9.5.1 Guided Waves

Moving beyond plane waves with infinite extent, now we consider waves confined to the interior to some sort of pipe, or wave guide. To make things simple, the overall geometry of the pipe should be large compared with the wavelength, and we'll assume it's made of a perfect conductor so there's no loss (perfect reflection).

Figure 9.23

Monochromatic plane wave solutions will look like

E(x,y,z,t)=E0(x,y)ei(kzωt) \vec E(x, y, z, t) = \vec{E_0} (x, y) e^{i(kz - \omega t)}

B(x,y,z,t)=B0(x,y)ei(kzωt) \vec B(x, y, z, t) = \vec{B_0} (x, y) e^{i(kz - \omega t)}

The (source-free) boundary conditions within the waveguide are

E=0B=0E=0×E=BtB=0×B=1c2Et \begin{aligned} E^\parallel & = & 0 \\ B_\perp & = & 0 \\ \div \vec E & = & 0 \\ \curl \vec E & = & - \pdv{\vec{B}}{t} \\ \div \vec B & = & 0 \\ \curl \vec B & = & \frac{1}{c^2} \pdv{\vec E}{t} \end{aligned}

So it turns out that within a waveguide we aren't necessarily limited to transverse solutions only. We include the longitudinal components of the fields when plugging in the coordinates into our boundary conditions and generic solutions:

E0=Exx^+Eyy^+Ezz^ \vec E_0 = E_x \vu x + E_y \vu y + E_z \vu z

B0=Bxx^+Byy^+Bzz^ \vec B_0 = B_x \vu x + B_y \vu y + B_z \vu z

Using the Maxwell equations and boundary conditions, eventually we can get independent expressions for the fields

[2x2+2y2+(ωc)2k2]Ez=0 \left[ \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \left( \frac{\omega}{c} \right) ^2 - k^2 \right] E_z = 0

[2x2+2y2+(ωc)2k2]Bz=0 \left[ \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \left( \frac{\omega}{c} \right) ^2 - k^2 \right] B_z = 0

For trivial solutions, we can come up with separate classes of solutions: If Ez=0 E_z = 0 , we call the solutions transverse electric (TE) waves. If Bz=0 B_z = 0 they are called transverse magnetic (TM) waves. And if both are zero, we call them TEM waves, but it turns out that TEM waves can't occur in a hollow wave guide, since

Ez=0Exx+Eyy=0 E_z = 0 \rightarrow \pdv{E_x}{x} + \pdv{E_y}{y} = 0

and

Bz=0EyxExy=0 B_z = 0 \rightarrow \pdv{E_y}{x} - \pdv{E_x}{y} = 0

imply that E0 \vec E_0 can be written as the gradient of a scalar potential that satisfies Laplace's equation. But the boundary condition on E \vec E requires that hte surface is an equipotential, so the only available potential is a constant and the field is zero everywhere, so no wave exists at all.

9.5.2 Rectangular Wave Guide

Let's look at a rectangular wave guide with height a a and width b b . As we saw, there are two chief classes of solutions, TE (Ez=0 E_z = 0 ) and TM (Bz=0 B_z = 0 ). Let's specifically take a look at the TM solutions (the process is very similar for TE waves, except that the boundary conditions are flipped).

TM solutions:

Bz=0 everywhere  B_z = 0 \text{ everywhere }

[2x2+2y2+(ωc)2k2]Ez=0 \left[ \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \left( \frac{\omega}{c} \right) ^2 - k^2 \right] E_z = 0

We've got a relatively simple situation here (really just a 2D Laplacian and constants), so start with the method of separation of variables.

Ez(x,y)=X(x)Y(y) E_z(x, y) = X(x) Y(y)

So plugging in:

Yd2Xdx2+Xd2Ydy2+[(ω/x)2k2]XY=0 Y \dv{^2 X}{x^2} + X \dv{^2 Y}{y^2} + \left[(\omega / x)^2 - k^2 \right]XY = 0

Separate out the X and Y terms

1Xd2Xdx2=kx2 \frac{1}{X} \dv{^2X}{x^2} = -k_x ^2

1Yd2Ydy2=ky2 \frac{1}{Y} \dv{^2Y}{y^2} = -k_y ^2

where

kx2ky2+(ω/c)2k2=0 -k_x ^2 -k_y ^2 + (\omega / c)^2 - k^2 = 0

We already know the general solutions

X(x)=Asin(kxx)+Bcos(kxx) X(x) = A \sin (k_x x) + B \cos (k_x x)

The boundary conditions require that Ez=0 E_z = 0 at x=0 x = 0 and at x=a x = a , so

X(0)=Bcoskx(0)=B=0X(a)=Asinkxa=0kxa=mπkx=mπa X(0) = B \cos k_x (0) = B = 0 \\ \rightarrow X(a) = A \sin k_x a = 0 \\ \rightarrow k_x a = m \pi \\ \rightarrow k_x = \frac{m \pi}{a}

Similarly,

ky=nπbY(y)=Asinkyy k_y = \frac{n \pi}{b} \\ Y(y) = A \sin k_y y

So the TM modes are

TMmn:Ez(x,y)=E0sin(mπxa)sin(nπyb) TM_{mn}: \quad E_z (x, y) = E_0 \sin \left( \frac{m \pi x}{a} \right) \sin \left( \frac{n \pi y}{b} \right)

Note that if m=0 m = 0 or n=0 n = 0 , then Ez E_z becomes zero immediately, so the "bottom-most" term is m=n=1 m = n = 1 . How do m m and n n relate to the wavenumber and frequency of the modes? We got our separable wavenumbers from

kx2+ky2+k2=(ωc)2(mπa)2+(nπb)2+k2=(ωc)2ω=c(mπa)2+(nπb)2+k2 k_x ^2 + k_y ^2 + k^2 = \left( \frac{\omega}{c} \right) ^2 \\ \rightarrow \left( \frac{m \pi}{a} \right) ^2 + \left( \frac{n \pi}{b} \right) ^2 + k^2 = \left( \frac{\omega}{c} \right) ^2 \\ \rightarrow \omega = c \sqrt{ \left( \frac{m \pi}{a} \right) ^2 + \left( \frac{n \pi}{b} \right) ^2 + k^2 }

This is a sort of weird dispersion relation. We get a minimum ω \omega as k0 k \rightarrow 0 , since k2<0 k^2 < 0 would give an attenuated wave which we don't care about, so there is a cutoff frequency

ωcutoff=ωk0cπ2(1a2+1b2) \omega_{cutoff} = \omega_{k \rightarrow 0} \rightarrow c \sqrt{ \pi^2 (\frac{1}{a^2} + \frac{1}{b^2})}

This means that TM modes can not propagate with frequencies below cπ2(1a2+1b2) c \sqrt{ \pi^2 (\frac{1}{a^2} + \frac{1}{b^2})} .

On the flip-side, what do the group and phase velocities look like?

k=(ωc)2(mπa)2(nπb)2=(ωc)2(ωmnc)2ωmn=c(mπa)2+(nπb)2k=ω2ωmn2ck=ω1ωmn2/ω2cωk=vp=c1ωmn2/ω2 k = \sqrt{ \left( \frac{\omega}{c} \right) ^2 - \left( \frac{m \pi}{a} \right) ^2 - \left( \frac{n \pi}{b} \right) ^2 } \\ = \sqrt{ \left( \frac{\omega}{c} \right) ^2 - \left( \frac{\omega_{mn}}{c} \right) ^2} \\ \omega_{mn} = c\sqrt{\left( \frac{m \pi}{a} \right) ^2 + \left( \frac{n \pi}{b} \right) ^2 } \\ \rightarrow k = \frac{\sqrt{\omega^2 - \omega_{mn}^2}}{c} \\ \rightarrow k = \frac{\omega \sqrt{1 - \omega_{mn}^2 / \omega^2}}{c} \\ \rightarrow \frac{\omega}{k} = v_p = \frac{c}{\sqrt{1 - \omega_{mn}^2 / \omega^2}}

vg=dωdk=1dkdω=11c2ω2ω2ωmn2=cω2ωmn2ωvg=c1ωmn2ω2 v_g = \dv{\omega}{k} \\ = \frac{1}{\dv{k}{\omega}} = \frac{1}{\frac{1}{c} \frac{2 \omega}{2 \sqrt{\omega^2 - \omega_{mn}^2}}} \\ = \frac{c \sqrt{\omega^2 - \omega_{mn}^2}}{\omega} \\ v_g = c \sqrt{1 - \frac{\omega_{mn}^2}{\omega^2}}

What about the TEmn TE_{mn} modes?

Bz(x,y)=X(x)Y(y) B_z(x, y) = X(x) Y(y)

1xd2Xdx2=kx2X(x)=Asinkxx+Bcoskxx \frac{1}{x} \dv{^2 X}{x^2} = - k_x ^2 \\ \rightarrow X(x) = A \sin k_x x + B \cos k_x x

Now, the boundary conditions say that

B=0 B^ \perp = 0

Bxx=0 at x=0 and x=a \pdv{B_x}{x} = 0 \text{ at } x = 0 \text{ and } x = a

AkxcoskxxBkxsinkxx=0A=0 \rightarrow A k_x \cos k_x x - B k_x \sin k_x x = 0 \\ \rightarrow A = 0

Plugging in x=a x = a ,

Bkxsinkxa=0kxa=mπ Bk_x \sin k_x a = 0 \\ \rightarrow k_x a = m \pi

Bz(x,y)=B0cos(mπxa)cos(nπyb) \rightarrow B_z(x, y) = B_0 \cos \left( \frac{m \pi x}{a} \right) \cos \left( \frac{n \pi y}{b} \right)

k=(ω/c)2π2[(m/a)2+(n/b)2] k = \sqrt{(\omega / c)^2 - \pi ^2[(m / a)^2 + (n / b)^2] }

If

ω<cπ(m/a)2+(n/b)2ωc,mn \omega < c \pi \sqrt{(m / a)^2 + (n / b)^2} \equiv \omega_{c, mn}

then the wave number is imaginary, and instead of a traveling wave we have exponentially attenuated fields, so this is our cutoff frequency for TE waves. Given the convention of always associating the first index with the larger dimension (assume ab a \geq b ), the lowest frequency for a given wave guide is

ω1,0=cπ/a \omega_{1,0} = c \pi / a

Example 9.1

We have a rectangular waveguide of dimensions 1cm by 2cm. (a) What is the lowest mode? Find the cutoff frequency. (b) If the waveguide were filled with lossless plastic with ϵr=2 \epsilon_r = 2, how would the cutoff frequency change?

What is the lowest allowed mode? We should check the lowest TE and TM modes to find out!

The lowest TM mode is

ωc,TM=cπ2(1a2+1b2) \omega_{c, TM} = c \sqrt{\pi^2 \left( \frac{1}{a^2} + \frac{1}{b^2} \right) }

and the lowest TE mode is

ω1,0=cπ2cmfco=c4cm=3×1010cm/s4=7.5GHz \omega_{1, 0} = c \frac{\pi}{2 \text{cm}} \\ \rightarrow f_{co} = \frac{c}{4 \text{cm}} = \frac{3 \times 10^{10} \text{cm/s}}{4} = 7.5 \text{GHz}

It's worth noting that the lowest TM mode will always be higher, since it requires both m1 m \geq 1 and n1 n \geq 1 , so the lowest frequency will be the TE mode.

(b) For vacuum,

1c2=μ0ϵ0 \frac{1}{c^2} = \mu_0 \epsilon_0

In plastic,

v=cϵ/ϵ0=c2fco=7.5GHz25.4GHz v = \frac{c}{\sqrt{\epsilon / \epsilon_0}} = \frac{c}{\sqrt{2}}\\ \rightarrow f_{co} = \frac{7.5 \text{GHz}}{\sqrt{2}} \approx 5.4\text{GHz}