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⊥=0∇⋅E=0∇×E=−∂B∂t∇⋅B=0∇×B=1c2∂E∂t \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

[∂2∂x2+∂2∂y2+(ωc)2−k2]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

[∂2∂x2+∂2∂y2+(ωc)2−k2]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=0→∂Ex∂x+∂Ey∂y=0 E_z = 0 \rightarrow \pdv{E_x}{x} + \pdv{E_y}{y} = 0

and

Bz=0→∂Ey∂x−∂Ex∂y=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 }

[∂2∂x2+∂2∂y2+(ωc)2−k2]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:

Yd 2Xd x2+Xd 2Yd y2+[(ω/x)2−k2]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

1Xd 2Xd x2=−kx2 \frac{1}{X} \dv{^2X}{x^2} = -k_x ^2

1Yd 2Yd y2=−ky2 \frac{1}{Y} \dv{^2Y}{y^2} = -k_y ^2

where

−kx2−ky2+(ω/c)2−k2=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)=Bcos⁡kx(0)=B=0→X(a)=Asin⁡kxa=0→kxa=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)=Asin⁡kyy 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 k→0 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=ωk→0→cπ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)2→k=ω2−ωmn2c→k=ω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 ωd k=1d kd ω=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)

1xd 2Xd x2=−kx2→X(x)=Asin⁡kxx+Bcos⁡kxx \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

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

→Akxcos⁡kxx−Bkxsin⁡kxx=0→A=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 ,

Bkxsin⁡kxa=0→kxa=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 a≥b 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π2cm→fco=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 m≥1 m \geq 1 and n≥1 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=c2→fco=7.5GHz2≈5.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}