]>

Spherical Waves and Applications

5
Spherical Waves and Applications

So far, we have discussed plane waves and cylindrical waves. In both of these cases, we could develop with relative ease the solutions based on TMz and TEz modes. We illustrated these solutions by solving some waveguides and cavity problems. The superscript z is added here to indicate that the transverse plane is perpendicular to the z-coordinates. The structures under consideration had a uniform cross section for all values of z. In spherical coordinates, we do not have any such transverse planes, and thus spherical wave problems are mathematically more involved. It is still possible to have mode classification as TMr and TEr.

The solution in spherical coordinates involves special functions: spherical Bessel functions and associated Legendre functions. Spherical Bessel functions are related to half-integral Bessel functions.

5.1Half-Integral Bessel Functions

In connection with the solution for a scalar Helmholtz equation in spherical coordinates, we come across the differential equation

(5.1)
udduudfdu+u2n+122f=0,

which, when expanded, gives

(5.2)
d2fdu2+1udfdu+1n+1/22u2f=0.

By comparing with Equation 2.120, we immediately realize that f can be written as a linear combination of half-integral Bessel functions:

(5.3)
f=Jn+1/2uYn+1/2uHn+1/21uHn+1/22uIn+1/2uKn+1/2u.

Half-integral Bessel functions can be shown to reduce to simpler forms:

(5.4)
Jn+1/2x=2πxn+1/21xndndxnsinxx,
(5.5)
Yn+1/2x=1n+12πxn+1/21xndndxncosxx.

In particular, we can show from Equation 5.4 that

(5.6)
J1/2x=2πxsinx,
(5.7)
J3/2x=2πxsinxxcosx.

From Equations 5.6 and 5.7, we can obtain the zeros of half-integral Bessel functions and their derivatives. These are listed in Tables 5.1 and 5.2.

TABLE 5.1
Zeros of Half-Integral Bessel Functions χn+1/2,I

J1/2 J3/2 J5/2
3.1416 4.4934 5.7635
6.2832 7.7253 9.0950

TABLE 5.2
Zeros of Derivatives of Half-Integral Bessel Functions χn+1/2,I

J1/2 J3/2 J5/2
1.1656 2.4605 3.6328
4.6042 6.0293 7.3670

5.2Solutions of Scalar Helmholtz Equation

The scalar Helmholtz equations

(5.8)
2F+k2F=0

and

1r2rr2Fr+1r2sinθθsinθFθ+1r2sin2θ2Fϕ2+k2F=0

can be solved in spherical coordinates by the usual technique of separation of variables. Let

(5.9)
F=r,θ,ϕ=f1rf2θf3ϕ.

Substituting Equations 5.9 into 5.8 and manipulating the equation, we obtain the ordinary differential equations for f1, f2, and f3 involving two separation constants m and n:

(5.10)
ddrr2df1dr+k2r2nn+1f1=0,
(5.11)
1sinθddθsinθdf2dθ+nn+1m2sin2θf2=0,
(5.12)
d2f3dϕ2+m2f3=0.

By substituting

(5.13)
r=ukandf1=fr=kuf

in Equation 5.10, we obtain Equation 5.2, with f thus given by a linear combination of half-integral Bessel functions given by Equation 5.3. The solution for f1 is given by

(5.14)
f1r=fur=fkrr.

Defining the spherical Bessel function bn in terms of the half-integral Bessel functions Bn+1/2, we obtain

(5.15)
bnkr=π2krBn+1/2kr,

where

  • bn stands for any of the spherical Bessel functions jn, yn, hn1, and hn2

  • Bn+1/2 stands for the half–integral Bessel functions, Jn + 1/2, Yn+1/2, Hn+1/21, and Hn+1/22

Equation 5.11 is called the associate Legendre equation whose solutions are linear combinations of associate Legendre polynomials of the first kind Pnm and the second kind Qnm of degree n and order m. These polynomials are given by

(5.16)
Pnmx=x21m/2dmdxmPnx,

(5.17)
Pnx=12nn!dndxnx21n,
(5.18)
Qnmx=x21m/2dmdxmQnx,
(5.19)
Qnx=12Pnxln1+x1xl=1n1lPl1xPn1x.

Associated Legendre polynomials of zero order (m = 0) are called Legendre polynomials. In Equation 5.16, Equation 5.17, Equation 5.18 and Equation 5.19, x = cos θ. Note that the order m is an integer between 0 and n. Plots and expressions for a few Legendre polynomials are given in Figure 5.1 and Table 5.3. Note that on the polar axis, that is, θ = 0 or π, Qn → ∞. So, for the problem that includes a positive or negative polar axis in the field region, the function Pnm alone is the suitable solution. The solution for Equation 5.12 is a linear combination of cos mϕ and sin mϕ. Putting together the solutions for f1, f2, and f3, the solution for Equation 5.8 may be written as

(5.20)
Fr,θ,ϕ=jnkrynkrhn1krhn1krPnmcosθQnmcosθcosmϕsinmϕ.

image

FIGURE 5.1
Polynomials Pn(cos θ). The degree n is indicated on the graph.

TABLE 5.3
Legendre Functions

image

5.3Vector Helmholtz Equation

The electric and magnetic fields in a sourceless region satisfy the vector Helmholtz equations

(5.21)
××E˜k2E˜=2E˜+k2E˜=0,
(5.22)
××H˜k2H˜=2H˜+k2H˜=0.

Each Cartesian component of E˜ and H˜ satisfies the scalar Helmholtz equation (5.8), but the spherical components do not satisfy Equation 5.8.

It is convenient once more to think of simple solutions of these equations in terms of TEr, TMr, and TEMr modes in the context of spherical geometry.

5.4TMr Modes

These modes may be constructed by considering a vector field

(5.23)
M˜=×r^rF˜e,

where F˜e satisfies Equation 5.8. It can be shown that M˜ satisfies the vector Helmholtz equation. Identifying M˜ with H˜ and noting from Equation 5.23 that

(5.24)
Hr=0,
(5.25)
H˜θ=1sinθFeϕ,
(5.26)
H˜ϕ=Feθ,

we recognize all the above as the magnetic field components of TMr modes in terms of the electric Debye potential F˜e. The corresponding electric field components can be obtained from Maxwell’s equation

(5.27)
×H˜=jωεE˜

and are given as

(5.28)
E˜r=jωε2r2rF˜e+k2rF˜e,
(5.29)
E˜θ=jωε1r2rθrF˜e,
(5.30)
E˜ϕ=jωε1rsinθ2rϕrF˜e.

5.5TEr Modes

These modes may be constructed by considering the vector function

(5.31)
N˜=×r^rF˜m,

where F˜m satisfies Equation 5.8. It can be shown that N˜ satisfies the vector Helmholtz equation. Identifying N˜ with the E˜ field and noting from Equation 5.31 that

(5.32)
E˜r=0,

(5.33)
E˜θ=1sinθFmϕ,
(5.34)
E˜ϕ=Fmθ,

we recognize the above as the electric field expressions of TEr modes. The corresponding magnetic field components can be obtained from Maxwell’s equation

(5.35)
×E˜=jωμH˜

and are given by

(5.36)
H˜r=jωμ2r2rF˜m+k2rF˜m,
(5.37)
H˜θ=jωμ1r2rθrF˜m,
(5.38)
H˜ϕ=jωμ1rsinθ2rϕrF˜m.

5.6Spherical Cavity

We illustrate the construction of solutions in spherical coordinates by determining the resonant frequencies of a spherical cavity with PEC boundary conditions at r = a (Figure 5.2).

image

FIGURE 5.2
Spherical cavity with PEC boundary.

For TM modes F˜e is given by Equation 5.20.

Since θ = 0 or π are part of the field region, we choose Pnm for θ variation. Since we need an oscillatory function that is finite at the origin, we choose jn(kr) for n variation.

The r variation of rFe is given by

(5.39)
rjnkr=rπ2krJn+1/2kr=rπ2kJn+1/2kr.

To suit such applications as above, Schelkunoff defined another set of spherical Bessel functions:

(5.40)
B^nx=πx2Bn+1/2x.

The roots of J^nζ are given in Table 5.4.

TABLE 5.4
Zeros of J^n

ζnp n = 1 n = 2 n = 3
p = 1 4.493 5.763 6.988
p = 2 7.725 9.095 10.417
p = 3 10.904 12.323 13.698

We may as well now list in Table 5.5 the roots of J^ny=0.

TABLE 5.5
Zeros of J^n

ζnp n = 1 n = 2 n = 3
p = 1 2.744 3.87 4.973
p = 2 6.117 7.443 8.722
p = 3 9.317 10.713 12.064

The PEC boundary conditions are E˜θ=0andE˜ϕ=0 at r = a, which leads to the boundary condition

(5.41)
rrF˜e|r=a=0,that is,J^nka=0,
(5.42)
knmpTMr=ζnpa,
(5.43)
frnmpTMr=12πμεζnpa,n,p=1,2,,.

Note that n cannot be zero since

(5.44)
P00=1.

And from Equations 5.29 and 5.30, E˜θ and E˜ϕ are zero not only on the boundary but also everywhere inside the cavity. Thus, n = 0 gives a trivial solution. The lowest resonant frequency of TM modes is given by

(5.45)
fr101TMr=12πμεζ11a=2.7442πμεa.

Let us next consider the TEr modes. The boundary conditions in this case, from Equations 5.33 and 5.34, are

(5.46)
Fmθ|r=a=0,
(5.47)
Fmϕ|r=a=0,

which translate to

(5.48)
J^nka=0,
(5.49)
knmpTEr=ζnp,
(5.50)
frnmpTEr=12πμεζnpa.

The lowest value of the above occurs when n = 1 and p = 1,

(5.51)
fr101TEr=12πμε4.493a.

The lowest resonant frequency of all the modes is thus given by Equation 5.45. The modes in a spherical cavity are highly degenerate, that is, for the same resonant frequency we can have many different field distributions.

..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset
18.227.161.132