Xem mẫu
- Spheroidal Wave Functions in Electromagnetic Theory
Le-Wei Li, Xiao-Kang Kang, Mook-Seng Leong
Copyright 2002 John Wiley & Sons, Inc.
ISBNs: 0-471-03170-4 (Hardback); 0-471-22157-0 (Electronic)
5
EM Scattering by a
Coated Dielectric
Spheroid
5.1 GEOMETRY OF THE PROBLEM
In this chapter we consider the scattering of a linearly polarized plane mono-
chromatic wave by a homogeneous lossy/lossless dielectric spheroid with a
confocal lossy/lossless dielectric coating immersed in a homogeneous isotropic
medium. It is assumed that the surrounding medium is nonconducting and
nonmagnetic. Results are presented only for the prolate spheroids, as the re-
sults for the oblate cases can be obtained by the transformations < --+ i
- 116 EM SCATTERING BY A COATED DIELECTRIC SPHEROID
Incident Wave $2
Wave
Fig. 5.1 Scattering geometry of a coated spheroid.
- INCIDENT, TRANSMITTED AND SCATTERED FIELDS 117
62
c2 = -co* (5.lc)
lr CO
5.2 INCIDENT, TRANSMITTED AND SCATTERED FIELDS
Unlike the case of the perfectly conducting spheroid, the existence of fields
inside the dielectric spheroid makes the problem more complicated, and it is
necessary to use the magnetic fields in the boundary conditions. The magnetic
field of a propagating wave is related to the electric field via
H - &VxE,
- (5 .2)
where k is the wave number and 2 is the characteristic impedance of the
medium. The incident electric fields for TE and TM polarizations are the
same as those of the conductor case, and the corresponding magnetic fields
are
00
HTE =
jETE0
MO coseo n=m x m=O
(5.3a)
HTM = (5.3b)
n=m m=O
The scattered fields are also derived in the same way as in the conductor
case, but a slight change in the notation of the scattering coefficients is now
made:
(5.4a)
(5*4b)
I
for the TE case; and
aM M+t4)
&TM = mn *,m+l,ntCO) + PmM+l,n+l M$Z+l,n+l
00
M- M-c4)
+>: Qcln oln CcO), (5.5a)
n=l
aM
mn
M- N-c4) (5.5b)
%n oln (‘0)
- 118 EM SCATTERUVG BY A COATED DIELECTRIC SPHEROID
for the TM case.
The transmitted fields inside the core (t < &) can be expressed as
(5.6a)
for the TE case, and
E2t,TM = d!Zn”~~~l,n(CZ) + dF+l 9n+l”Z(ltl+l 9 9n+l
00
n=mm=O L
+c M-M-(‘)
n=l
%t oln (‘2)7 (5.7a)
*
oM
mn N~Zil,n(c2) + 7mM+l,n+l Nfifz+l,n+l 4
1
+% OK- N,-,{Z) (c2) (5.7b)
z2 n=l
for the TM case. Note that the spheroidal radial functions of the first cind are
now used instead because they are finite near the origin, and the transmitted
wave has to be finite at the center of the spheroid.
For the region & < c < &, the transmitted waves are represented as
El t,TE = YEmn M + 6~nM~~!,(cl)
(5.8a)
(5.8b)
for the TE case; and
El t,TM = + 6mM+1 9n+l MZ’Z+l 9 9n+l(‘l))
n=mm=O L
+ KYrz M~E~l,n(cl) + P,“,l 7n+l MZ(Z-i+l9 7n+l
( )I‘1
M-M-(~)
In oln (‘1) + Aln
M-
l”oln
--w
(cl)
1 9 (5.9a)
- RELATIONSHIP BETWEEN INCIDENT AND SCATTERED FIELDS 119
+ XmMnN~?j-l,n(C1) + P,M+l 9n+lNz(z+l 7 9n+l
( )Ic1
YlM,- Niln (l)(cl) + AK- N,i’,2’(cl) (5.9b)
I
for the TM cMX3* 7mn9 6m+l,n+l7 Am727 and Pm+r,n+r are the unknown trans-
mission coefficients. Note that in the coating region both first and second
kinds of spheroidal functions are used, because down here, the wave consists
of two components, one propagating inward and one outward. This is due to
reflections at both boundary surfaces.
5.3 RELATIONSHIP BETWEEN INCIDENT AND SCATTERED
FIELDS
5.3.1 Boundary Conditions
The unknown scattering and transmission coefficients can be determined by
applying the boundary conditions (i.e., continuity of the tangential compo-
nents of the electric and magnetic fields across each of the spheroidal surfaces
at < = 51 and 5 = &). Thus, we have
where the suffixes q and 4 denote, respectively, the r)- and +-components of
the fields. These equations must hold for all allowed values of 0 - 4 - 27r
< <
and -1 - q - 1.
< <
5.3.2 TE Polarization for Nonaxial Incidence
The same method is applied as for the conductor case in Chapter 4. This
means that making use of the orthogonality of the trigonometric functions
and spheroidal angular functions, and integrating over q, we will get
-1 E E E
(COS~O) (Qm>Am = (R,)S,, m = 0,172, l l l 7 (5.11)
and
-1 E E E
(coseo) (Q+)A+ = (R+)S+, (5.12)
- 120 EM SCATTERING BY A COATED DIELECTRIC SPHEROID
where
’ -X$(Q)) -V$(co) ’
-x:,, -v:‘, (co)
(co)
-xC4) -vZ& (co)
(co) m2
-- --- --- --
(R a >=
- y%(co) -w$(co)
- Y 21 (co) -3 (co)
-yap, -wt4) (co)
m2
.
I
e2w I v:&(4
(1)
xi!(cl) I vml(cl)
x:$1, I v:&(Q)
. .
. .
. I .
-- -- --_ --_-
CR) =
b
e-&l) I w:~ccl)
ezlk4 I w:;(Q)
Y?~2h~ I w:k4
. .
. .
. I .
- RELAT/ONSH/P BETWEEN /NC/DENT AND SCATTERED FIELDS 121
CR ) =
d ------
,
Bg&l) I Bg$(Cl) ’
Bg&l) I Bg)&)
CR e > =
e!oh) I $2&I)
&!&l) I &2&l)
.
.
. .
.
.
\ I
-
-
Wf)
Bgo(cl) I l$i?&l)
~~,h) I J&?&l)
.
.. .
.
I .
- 122 EM SCATTERING BY A COATED DIELECTRIC SPHEROID
f -q&2, 1 -V$(c2)
-x(‘)ml (c2) 1 -v(l) (c2)
ml
. .
. .
. I .
(R9 >=
f 43~9~2) 1 -By,
-
G&2) I -g$(cz)
.
. .
.
. I .
CR >=
h
w(2)(c1 Cl) 0
WfyCl’< 7 1) 0
.
. .
.
. I . I
---------- -- -
(4)
-B,, (co,&) I B$)$a a I
(1)
-BiP,)(co&) I B,, &l) ‘) I
.
.
.
.
I .. I
WE>
+= _------- ---- - ------ ----me
0 (1) (2)
I w,, (Q&> WlO (Cl 3 -wg(c2 C2)
0 (1) -W(‘)(cz’ 11
I ‘w,, (cl7c2> Wl”,‘(Cl, C2) ’ C2)
. .
. .
.
. .
I I
---e--B- ---- - ----m- ---w--
0 I B$)@l ?t2> @(Cl 9 (2) 49~2 SO 3 t2)
(1) (1)
0 I q1 hb> B(2)(~~
Sl 9 t2) -B,, cc27
.
.
.
.
.
.
I
- RElATIOfUSHlP BETWEEN lNClDENT AND SCATTERED FIELDS 123
(Qa) = km _ _ i - _ c1 =
(Q
y~~o(co) -Y(l)cm0 (co)
Y i% (co) -y::, (co)
.
.
.
I
’ Bg,(co) x~~oko)
B (1) w (1)
Xfml (‘0)
am1
.
. .
.
. .
(Qb) = k, ----- 7 (Q d > =
-----
(1)
B cmow x($) (co)
(1)
B cm1 9
hd $21 (co)
.
. .
.
. .
Sirn
E’
cY,,,+l
.
.
y(l) a10 (co 7 Sl) .
(1) -----
ya,,bh)
E’
l
.
P m+l,m+l
. E’
-- ---
P m+l,m+2
.
.
.
X~$yco,h~
-----
(QE)
+= x;; (co, Cl)
7 SEml -
-
. YEImm
. E’
.
Yon,m+l
-----
.
.
0 .
----- -- ---
0 bEf
m+l,m+l
liEI
m+l,m+2
.
.
.
- 124 EM SCATTERING BY A COATED DIELECTRIC SPHEROID
AE’ \ E’
F” P 00
XE E’
m,m+l P 01
.
.
.
----- -----
PiX’+l.,m+l liEI
00
Pit’+l,m+2 liE’01
.
.
.
SE -- ----- 7 SE - -----
m2 +- 7
oE’
mm PoEd
E’
cTm,m+l POE;
.
.
.
----- -----
E’ E’
7,+l,m+l 700
E’ E’
7,+l,m+2 701
.
.
. .
k m = 2 for nz = 0 and km = 1 for nz > 0, and
The column matrices Am and A+ are the same as those defined for the
conductor case in Chapter 4. 0 is the zero matrix. Besides the row matrices
used in Chapter 4, the other row matrices are defined as follows:
i i
B 0 (c> =
amN
0
(amNO (’ > .. .
>7 (5.13)
where
- 77-a 1
+ (i)
a@) (c)
mNn IAmNn + ICmNn - cc,” _ q1/2 Rm,m+n(cO)IBmNn
- m(m + ‘> R(i)
m,m+n(CO)IDmNn;
(s,” - 1y2
B@) (c)
hN
= (b@) (c)
mN0
b@) (c)
mN1
b@) (c)
mN2
. . 0)9 (5.14)
where
bci) (c)
mNn IEmNn -IGmNn + (m + ‘j2” Rf$+l,m+n+l({ 0 )I HmNn ;
c; - 1 I
0i (5.15)
B cmN ( c) = (P
mN0 (c) cCi)
mN1 (c) P)
mN2 (c) l l 0)
9
- RELATIONSHIP BETWEEN INCIDENT AND SCATTERED FIELDS 125
where
. 1
0
cmNn
2 (c) = R:;m+n (tO)IImNn
(g - 1y2
.
0
+ (r,” _ml),/2 RG,m+n (b)IKmNn
IJmNn
t=t0
IJmNn;
t= = (pmN0 (c) pmN1 (c) pmN2 (c) 3 (5.18)
where
1
&~,(c) = - RE$2,rn+n+2 (
- 126 EM SCATTERING BY A COATED DIELECTRIC SPHEROID
. .
0 h@) hb)
x( h2N ) = (h lt0 Nl N2 l l l )7 (5.19)
where
.
h0
AZ = IAmNn +kmNn ;
m=l m=l
and finally,
0i (,(i) & SE; . . .) 7
B sN = NO (5.20)
where
& -
- 2(I es N n- &~sNn) .
Here CO= s m,m+N d% (5.21d)
I -1 -
drl
ICmNn = (5.21e)
IC2mNn = ~m+2,m+n+2~m,m+N dq, (5.21f)
IDmNn = w1g)
ID2mNn = (5.21h)
IEmNn = (1 - 172)1’2Em+l,m+n+1Sm,m+N dq, (5.21i)
IGmNn = Fm+l,m+n+lSm,m+N dq, (5.21j)
IHmNn = sm+l,m+n+lSm,m+N d% (5.21k)
IImNn = Gm,m+nSm,m+N dq, (5.211)
- RELATIONSHlP BETWEEN INCIDENT AND SCATTERED FIELDS 127
+1
(5.2lm)
b2mNn =
J-1 Gm+2,m+n+2Sm,m+N d%
+1 1
(Sm,m+n>sm,m+N d% (5.21n)
IJmNn =
J-1 (1 - q2y2
+1 1
(S m+2,m+n+2>Sm,m+N d% (5.210)
IJ2mNn =
J-1 (1 - q2y2
+l d
IKmNn =
J-1
& ( (1 - rl
q-2)1/2s
mym+n s m,m+N d% (5*2lP)
+l d
IK2mNn =
J (
-1
&j (1 _
rl
72)1/2
s m+2,m-i-n+2 Sm,m+N d%
ww
+’ d
(5.21r)
IOmNn =
J -
-1 drl
(S m+l,m+n+l )Sm,m+N d%
+1 rl
(5.21s)
IPmNn =
J-1 l-q2
P m+l,m+n+l 1s m,m+N dr]~
I esNn (5.21t)
=
J-1
+‘(I - ~2)1’2&ns1,1+N dq~
+1
(5.21~)
IGmNn =
J-1
(1 - ~2)1’2&haS1,1+N d’&
Besides the spheroidal angular functions (Sm,m+n), the other functions
used in the integrals are
B m,m+n (5.22a)
C m,m+n
(5.22b)
d
E m+l,m+n+l = r7-Sm+l,m+n+l
drl
(5.22~)
Fm+l,m+n+l = 9
(5.22d)
G m,m+n = f ((I- q2)1’2$Sm,m+n) 3 (5.22e)
Eon = r)- ivSOn$ (fiRC) 7 (5.22f)
- 128 EM SCATTERING BY A COATED DIELECTRIC SPHEROID
The functions Bm+2 m+n+2, Cm+2 m+n+2, and Gm+2 m+n+2 are obtained
from their counterparts (with subscripts m and nz + A) by substituting m
bym+Zandm+nbym+n+Z.
Equation (5.11) can be written as
SE
m = (G$)IE, m=O,l,Z ,...
where (GE) = (RE)-l(QE) represents the transformation matrix to trans-
form the incident vector 1: = Am(cos 00)-l into the column vector of the
unknown scattering and transmitting coefficients. We also note that (GE)
is dependent only on the scattering body, which is the same as that of the
conducting body.
5.3.3 TM Polarization for Nonaxial Incidence
Using the same principles as before, we will get for the TM case,
(Q,M)Am = (Rf)S,M, m = O,l,Z,. .., (5.23)
(5.24)
(Q,M) =
The matrices Qa, Qb, QC, Qd, (RE), Am, and A+ are those defined for
the TE case of the coated spheroid. The elements of Sf are the same as
those for Sz except that the superscripts are now A-4instead of E. The other
- RELATIONSHIP BETWEEN INCIDENT AND SCATTERED FIELDS 129
matrices which are required are
I
I
I
-
I
I
.
I
Pm - - ----A
1 X&Cl(1)
, b> I -x(l)610 (c2,
- 130 EM SCATTERlNG BY A COATED DIELECTRIC SPHEROlD
-- f (({’ - 1)‘/2$Rj:),,, IJllVn
c=to
R(i)
1,1+n IJllVn,
)I t=t0
(0 = 51, and/or t=& ,m=O
A00
A01
(Q&) - (Qb >
- NUMERICAL COMPUTATION AND MATHEMATICA SOURCE CODE 131
Unlike the conducting case, the truncation number required to obtain a
given accuracy in the computed cross sections depends on a number of factors:
for example, the electrical size, the properties of the coating and the core
material, and the thickness of the coating. It is very difficult to come up
with a rule that can take all these factors into account. As such, for the sizes
and permittivities of the spheroids and coatings considered here, it was found
accurate enough to consider only using m = 0, 1, 2 and 72 = m, m + 1, . . . ..
nz + 5. The reason is that due to the dimensions of the matrices involved, it
takes much longer to compute the system for a given m than in the conducting
case. Moreover, the convergence has been verified by computing until m = 4
for most cases. In all instances, the results matched at least to four significant
figures and in most cases to five significant figures. This ensures that the
scheme of truncation above is proper.
Three Mathematics packages were written for the coated dielectric spher-
oid. The first is Coatedback.nb, which contains the user module Coat-
edback [a2-, ratio2-, thickness-, epsilonrl,, epsilonr2J for computing
and plotting the normalized backscattering cross section. The second pack-
age, Coatedbistatic.nb, contains the module Coatedbistatic[a2-, ratio2,,
thickness-, epsilonr l-, epsilonr2,, thetaO-, phi-] for the computing and
plotting the normalized bistatic cross section at nonaxial incidence. The last
package, Coatedbizero.nb, with the module Coatedbistatzero[a2-, ra-
tiofl,, thickness,, epsilonrl-, epsilonr2-1, is developed for the bistatic
cross section at axial incidence. Note that the formulas for the various cross
sections are the same as those given for the conductor case.
To compute the cross sections, the user needs to load the relevant packages
and call any of the three modules. The arguments a2 and ratio2 denote
the semimajor axis length (in terms of X0) and the axial ratio of the core,
respectively. The argument thickness denotes the thickness of the dielec-
tric coating (also in terms of &), while epsilonrl and epsilonr2 are the
relative permittivities of the two regions (see Fig. 5.1). In addition, Coat-
edbistatic.m computes the bistatic cross section for a user-specified angle of
incidence (thetao-) at a user-specified azimuthal angle of observation (phi-).
For all three cases, the output is formatted in the same way as the conducting
case.
The packages mentioned above are very similar in their structures, except
for the last part, when they compute the different scattering cross sections.
So they shall be described as one entity. The program flow in the packages
can be summarized as follows:
1. Based on the user’s input, the parameters 50,
- 132 EM SCATTERING BY A COATED DlELECTRlC SPHEROID
3. The spheroidal radial functions of the first, second and fourth kinds and
their first and second order derivatives are evaluated next. All of these
are done with the following supporting modules:
l PSpheroidRl
l PSpheroidR2
l PSpheroidR2forsmallc
l PSpheroidR2forlargec
4. The various functions and integrals required for formulation of the ma-
trices are then defined. For the integrals, they are computed in a manner
similar to that for the conducting case.
5. The matrices for the TE and TM cases are then computed using the
“ Table” command.
6. Once the scattering coefficients are obtained, they are substituted into
the relevant expressions for the scattering cross sections, and a plot is
generated with the corresponding title and label.
5.5 RESULTS AND DISCUSSION
One advantage of the source code for coated dielectric spheroids is that it
allows us to model the scattering by homogeneous dielectric spheroids. This
is done by letting the relative permittivity of the coated layer be unity. Fig-
ure 5.2 shows the variation of the bistatic cross sections (at an axial incidence
of 60 = 0) of two dielectric spheroids of the same axial ratios and permittivities
but different semimajor axis lengths. The solid circles represent the results
from Asano and Yamamoto [24]. The close agreement once again verifies the
accuracy of the Mathematics source code used.
In Fig. 5.2, for the smaller spheroid (top), the scattering intensity is al-
most a constant in the plane perpendicular to the plane of polarization of
the incident wave (H-plane), whereas the scattering in the plane parallel to
the plane of polarization (E-plane) is approximately proportional to cos20.
This behavior is similar to that of the Rayleigh limit for the scattering by
very small spheres, whose radii are much smaller than for the incident wave-
length. In such cases(Rayleigh scattering), the intensity is a constant in the
plane perpendicular to the plane of polarization and varies as cos20 in the
parallel plane [132]. With the increase in size (bottom), the magnitude of the
scattered intensity increases, and the forward scattering amplitude becomes
much greater than the backscattering amplitude. Again, this is similar to the
Mie effect observed in scattering by spheres [133]. In addition, the patterns
in power intensity distributions become more and more complicated, showing
oscillating fluctuations with 0.
- RESULTS AND DISCUSSION 133
-20
-40 f
\ f
\ f
-60
0 20 40 60 80 100 120 140 160 180
8 (“1
(a) a2 = - l x0
4
0
-10
-20
-30
-40
-50 I I I ) I I 1 I 1 I 1 i I
0 20 40 60 80 100 120 140 160 180
$9 co)
(b) a2 = -X0
d
Fig. 5.2 Bistatic cross sections of dielectric spheroids with a& = 2 and a relative
permittivity of 1.769, but different sizes at an axial incidence of 80 = O”.
- 134 EM SCATTERING BY A COATED DIELECTRIC SPHEROID
Figures 5.3 and 5.4 show the variation of the bistatic cross sections with
different angles of incidence for a dielectric spheroid. For these cases, the
plane of observation is taken to be the X, z-plane. The cross sections for the
TE mode do not change dramatically as the incident angle varies, but the
minimum in the TM scattering continues to move toward the other end of
the spheroid as the incident angle increasesand very strong backscattering is
observed. The maximum in the TM scattering also increases, as the surface
area available for scattering reaches a maximum at t90= 90’. This is similar
to what was observed for a conducting spheroid in Chapter 4.
Figure 5.5 illustrates the difference between a lossy dielectric and lossless
dielectric on the backscattering cross sections. Notice that the shape of the
cross sections depends mainly on the dimensions of the spheroid. Since both
spheroids are of the sameaxial ratio and semimajor axis, they have roughly the
same shape. A smaller amplitude is observed from the lossy one as expected.
Figure 5.6 depicts the Mie scattering effect and the effect of dimensions on
the bistatic cross sections of dielectric spheres. The results were compared
with those (in bullets) by Van de Hulst [134]. In both cases, the forward-
scattering cross section amplitude is larger than the backscattering amplitude.
For a bigger sphere, the scattered intensity increases, and the pattern exhibits
more oscillations [1331.
Figure 5.7 shows the bistatic cross sections of a single dielectric spheroid
with a confocal coating for two values of the axial ratio . The results were
compared with those (in bullets) by Cooray and Ciric [33]. Here we observe
that the variation of both the E- and H-plane patterns is almost the same for
the two cases. There are signs of Rayleigh scattering, but the Mie scattering
effect is also coming into play. The difference between the two casesis that the
magnitudes of scattering cross sections are smaller for the “thinner” spheroid
(axial ratio of 5) than for the “fatter” spheroid. This is because the area
available for scattering is lesswhen the axial ratio is 5 than that when it is 2.
Figure 5.8 illustrates two different backscattering cross sections of two
coated spheroids of the same type of material but of totally different dimen-
sions. The one at the bottom has a much thicker coating relative to the inner
core, and it is also fatter (axial ratio of 2) than the one at the top (axial
ratio of 5). The one with the thinner coating exhibits a TM backscatter-
ing that is always higher than the TE backscattering. For the one with the
thicker coating, the TM backscattering starts off initially larger than the TE
backscattering, but the TE backscattering increases steadily after 00 = 45”
until it is largest at broadside incidence (00 = 900).
Here, the bistatic cross sections for dielectric-coated spheroids of oblique
angles of incidence were also computed and discussed. Figures 5.9 and 5.10
show the results for a coated dielectric spheroid, which is basically the product
of adding a confocal coating of thickness 0.02& to the dielectric spheroid of
Figs. 5.3 and 5.4. As such, the crosssections for both casesare almost identical
nguon tai.lieu . vn