Plane Wave Diffraction by Two Oppositely Placed, Parallel Two-Part Planes
Plane Wave Diffraction by Two Oppositely Placed, Parallel Two-Part Planes
Abstract: The problem of plane wave diffraction by two oppositely placed, parallel two-part
planes, consisting of the junction of perfectly conducing and resistive half-planes, is solved with the
Fourier transform technique. The problem is formulated as two simultaneous Wiener–Hopf
equations that are uncoupled by the analytical properties of the functions that occur. The solution
involves an infinite number of expansion coefficients satisfying an infinite system of linear algebraic
equations. The expansion coefficients are determined numerically, and some graphical results
showing the influence of the surface resistivity and the distances between the two-part planes on the
diffracted and transmitted waves are presented.
168 IEE Proc.-Sci. Meas. Technol., Vol. 153, No. 4, July 2006
Authorized licensed use limited to: ULAKBIM UASL - GEBZE INST OF TECH LIBRARY. Downloaded on March 23, 2009 at 05:49 from IEEE Xplore. Restrictions apply.
with R being the surface resistivity, given by account the following radiation and edge conditions
i [8, 9]:
h
R¼ si ð2cÞ
pffiffiffi @u pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
2ot ðe e0 Þ þ i r iku ! 0 r ¼ x2 þ y 2 ! 1 ð5gÞ
o @r
þ
In the above relationships, ½ denotes the discontinuity
across the transmissive sheet, and ~ n is the unit normal u ¼ Oðr1=2 Þ r ! 0 ð5hÞ
vector directed to the side indicated by the + sign. The By inserting (4a–c) into (5a–f ) and then inverting the
resistivities of the half-planes S12 and S21, are denoted by
pffiffiffiffiffiffiffiffiffiffiffi resulting integral equations, we obtain
z1;2 ¼ R1;2 =Z0 , with Z0 ¼ m0 =e0 being the characteristic
impedance of the free space. AðaÞ ¼ BðaÞeiKd þ CðaÞeiKd ð6aÞ
For analysis purposes, it is convenient to express the total
field as follows: AðaÞ ¼ þ
1 ðaÞ ð6bÞ
8
< ui þ ur þ u1 y4d
uT ðx; yÞ ¼ u2 y 2 ð0; dÞ ð3aÞ iAðaÞ½KðaÞ k=z1 iKðaÞ½BðaÞeiKd CðaÞeiKd
:
u3 yo0
2ik sin f0 eikd sin f0
Here, ur denotes the field that would be reflected if the ¼ þ
1 ðaÞ ð6cÞ
2piða k cos f0 Þ
whole plane y ¼ d was perfectly conducting
ur ðx; yÞ ¼ expfik½x cos f0 ðy 2dÞ sin f0 g ð3bÞ BðaÞ þ CðaÞ ¼ DðaÞ ð6dÞ
In (2a), un ðx; yÞ ðn ¼ 1; 2; 3Þ satisfy the Helmholtz equation
Dun þ k 2 un ¼ 0 n ¼ 1; 2; 3 DðaÞ ¼
2 ðaÞ ð6eÞ
in the regions where they are defined. For the sake of and
mathematical convenience, we assume that k has a small
positive imaginary part. The results corresponding to media DðaÞ½KðaÞ k=z2 þ KðaÞ½BðaÞ CðaÞeiKd
with negligibly small conductivity can be obtained by 2ik sin f0 eikd sin f0
making Imk ! þ0 at the end of the analysis. ¼ þ
1 ðaÞ ð6f Þ
2piða k cos f0 Þ
It is appropriate to use the following Fourier integral
representations: with
Z Z 1
u1 ¼ AðaÞeiKðaÞðydÞiax da ð4aÞ 1
þ
1 ðaÞ ¼ u1 ðx; dÞeiax dx ð6gÞ
L 2p 0
Z
Z
u2 ¼ ½BðaÞeiKðaÞy þ CðaÞeiKðaÞy eiax da ð4bÞ 1 0 @ @ ik
L 1 ðaÞ ¼ u1 ðx; dÞ u2 ðx; dÞ u1 ðx; dÞ eiax dx
2p 1 @y @y z1
and
Z ð6hÞ
u3 ¼ DðaÞeiKðaÞyiax da ð4cÞ Z
L 0
1
Here, K(a) denotes the square-root function
2 ðaÞ ¼ u3 ðx; 0Þeiax dx ð6iÞ
2p 1
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
KðaÞ ¼ k 2 a2 ð4dÞ and
which is defined in the complex a-plane, cut along a ¼ k to Z 1
1 @ @ ik
a ¼ k þ i1 and a ¼ k to a ¼ k i1, such that þ
2 ðaÞ ¼ u2 ðx; 0Þ u3 ðx; 0Þ u3 ðx; 0Þ eiax dx
Kð0Þ ¼ k. The integration line L is a straight line lying in 2p 0 @y @y z2
the strip ImðkÞoImðaÞoImðk cos f0 Þ. ð6jÞ
The spectral coefficients AðaÞ, BðaÞ, CðaÞ and DðaÞ are to In the region x o 0, where the sources are absent, the field
be determined through the following boundary and u1 behaves like a cylindrical wave that emanates from the
continuity relationships: origin and propagates towards infinity with the wave
u1 ðx; dÞ ¼ u2 ðx; dÞ x 2 ð1; 1Þ ð5aÞ number k (radiation condition):
ikx
u1 ðx; dÞ ¼ 0 xo0 ð5bÞ e
u1 ðx; yÞ ¼ O pffiffiffiffiffiffi x ! 1 ð7aÞ
@ @ @ x
ui ðx; dÞ þ ur ðx; dÞ þ u1 ðx; dÞ In the region x 4 0, where the sources are present, relying
@y @y @y
upon the boundary condition (5c), it is an easy matter to
@ ik
u2 ðx; dÞ ¼ u1 ðx; dÞ x40 ð5cÞ show that the asymptotic behaviour of u1 is the same as that
@y z1 of the incident field, namely
u2 ðx; 0Þ ¼ u3 ðx; 0Þ; x 2 ð1; 1Þ ð5dÞ
u1 ðx; yÞ ¼ Oðeikx cos f0 Þ x ! 1 ð7bÞ
u3 ðx; 0Þ ¼ 0 x40 ð5eÞ If we consider (7a) and (7b), we observe that
þ
1;2 ðaÞ and 1;2 ðaÞ are regular functions of a in the
@ @ ik
u2 ðx; 0Þ u3 ðx; 0Þ ¼ u3 ðx; 0Þ xo0 ð5f Þ half-planes ImðaÞ4Imðk cos f0 Þ and ImðaÞoImðkÞ,
@y @y z2 respectively.
To obtain a unique solution to the mixed boundary value The elimination of A(a), B(a), C(a) and D(a) between
problem stated by (3) and (5a–f ), we have to take into (6a) and (6f ) yields the following Wiener–Hopf matrix
IEE Proc.-Sci. Meas. Technol., Vol. 153, No. 4, July 2006 169
Authorized licensed use limited to: ULAKBIM UASL - GEBZE INST OF TECH LIBRARY. Downloaded on March 23, 2009 at 05:49 from IEEE Xplore. Restrictions apply.
equations: with
2 3 2 3
1 1 8 9
0 " þ # 1 > k >
6 MðaÞN1 ðaÞ 7 1 ðaÞ 6 MðaÞ 7 i1 ðaÞ >
> i$ d 2
sin ð$dÞ2 >
>
6 7 6 7 >
>d z1;2 >
>
4 1 5 þ 4 1 5 i >
> 2 3>>
1
2 ðaÞ
0 2 ðaÞ >
> 2 >
>
Z 1>
>p k >
>
MðaÞ MðaÞN2 ðaÞ < 2
6 $ 2 sin ð$dÞ
2
7=
q1;2 ðaÞ ¼ PV 6 z1;2 7
k sin f0 eikd sin f0 1 0 >
> ip6
6
7>
7>
¼ ð8aÞ >
>
> 4 þ i2$ k sinð$dÞ cosð$dÞ 5> >
>
ipða k cos f0 Þ 0 >
> 2 >
>
>
> z1;2 >
>
with >
> >
>
: ;
sin½KðaÞd
MðaÞ ¼ ð8bÞ
KðaÞ Kð$Þ þ a
ln d$ ð10bÞ
and Kð$Þ a
KðaÞ As to the explicit expression for M þ ðaÞ appearing in (9b), it
N1;2 ðaÞ ¼ is given in [11] as
½k=z1;2 KðaÞ sin½KðaÞd iKðaÞ cos½KðaÞd
ð8cÞ
sinðkdÞ 1=2 iad 2p p
M þ ðaÞ ¼ exp 1 C ln þi
k p kd 2
Y1
2.1 Solution of simultaneous Wiener–Hopf a iab
1þ exp ð11aÞ
equations m¼1
a m mp
As is well known, the first and most important step in
solving a matrix functional equation of the Wiener–Hopf and
type is the factorisation of the kernel matrix. To this end, let
us multiply (8a) from the left by the following matrix: M ðaÞ ¼ M þ ðaÞ
N1;2 þ
ðaÞ ¼ N1;2 ðaÞ ð11bÞ
L1 ðaÞ 0 In (11a), C is the Euler constant given by C ¼ 0:57721 . . ..
ð9aÞ
Lþ After the Wiener–Hopf decomposition have been
0 2 ðaÞ
performed, the matrix equation given by (9b) can be rear-
to obtain ranged as
" #" #
1=Lþ
1 ðaÞ 0 þ
1 ðaÞ þ
1 ðaÞ L ðk cos f0 Þk sin f0 eikd sin f0 L ðaÞ
þ 1 ¼ i 1 ðaÞ
Lþ
2 ðaÞ=MðaÞ Lþ
2 ðaÞ þ
2 ðaÞ L1 ðaÞ ipða k cos f0 Þ MðaÞ 2
" #" #
L
1 ðaÞ L1 ðaÞ=MðaÞ i
1 ðaÞ k sin f0 eikd sin f0
þ iL
1 ðaÞ1 ðaÞ þ
1=L i ipða k cos f0 Þ
0 2 ðaÞ 2 ðaÞ
k sin f0 eikd sin f0 1 ½L
1 ðaÞ L1 ðk cos f0 Þ ð12aÞ
¼ L ðaÞ ð9bÞ
ipða k cos f0 Þ 1 0 and
Here, Lþ
1;2 ðaÞ
denotes certain functions that regular and
2 ðaÞ iLþ ðaÞ þ
free of zeros in the half-planes ImðaÞ4ImðkÞ and ¼ 2 ðaÞ þ iLþ þ
2 ðaÞ2 ðaÞ ð12bÞ
ImðaÞoImðkÞ, respectively, resulting from the Wiener–
L2 ðaÞ MðaÞ 1
Hopf factorisation of The regularity of the right-hand side of (12a) in the lower
half-plane can be violated by the simple poles occurring at
L1;2 ðaÞ ¼ MðaÞN1;2 ðaÞ ð9cÞ the zeros of MðaÞ ¼ sin½KðaÞd=KðaÞ lying in the lower
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
as half-plane, namely at a ¼ am ¼ k 1 ðmp=kdÞ2 ,
m ¼ 1; 2; . . . , and the regularity of the right-hand side of
L1;2 ðaÞ ¼ Lþ
1;2 ðaÞL1;2 ðaÞ (12b) in the upper half-plane can be violated by the simple
þ
¼ N1;2 ðaÞM þ ðaÞN1;2
ðaÞM ðaÞ ð9dÞ poles occurring at the zeros ofqMðaÞ lying in the upper half-
ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
The explicit expression of N1;2
ðaÞ can be obtained by plane, namely at a ¼ am ¼ k 1 ðmp=kdÞ2 , m ¼ 1; 2; . . ..
P P1
following the procedure outlined in [10] If the infinite series of poles 1m¼1 am =ða þ am Þ and m¼1
bm =ða am Þ are subtracted from both sides of (12a) and
" #1=2 (12b), respectively, we obtain
þ KðaÞ
N1;2 ðaÞ¼
k=z1;2 KðaÞ sin½KðaÞdiKðaÞ cos½KðaÞd þ L ikd sin f0 X1
1 ðaÞ 1 ðk cos f0 Þk sin f0 e am
1=4 þ 2
kþa KðaÞd a þ iK ðaÞ L1 ðaÞ ipða k cos f0 Þ m¼1
a þ am
exp ln
ka p k L X1
1 ðaÞ am
¼ i 2 ðaÞ þ iL
1 ðaÞ1 ðaÞ
MðaÞ a þ am
iad 2a m¼1
0:5iKðaÞd þ 0:5q1;2 ðaÞ exp ln
p k k sin f0 eikd sin f0
þ ½L ðaÞ L
1 ðk cos f0 Þ ð13aÞ
ð10aÞ ipða k cos f0 Þ 1
170 IEE Proc.-Sci. Meas. Technol., Vol. 153, No. 4, July 2006
Authorized licensed use limited to: ULAKBIM UASL - GEBZE INST OF TECH LIBRARY. Downloaded on March 23, 2009 at 05:49 from IEEE Xplore. Restrictions apply.
and a ¼ k cos t, x ¼ r cos f and y d ¼ r sin f, with
f 2 ð0; pÞ. The saddle point occurs at t ¼ f, whose
X1
iLþ contribution is
2 ðaÞ bm þ þ 2 ðaÞ þ
¼ iL ðaÞ ðaÞ 1 ðaÞ
L ðaÞ a a m
2 2
MðaÞ eip=4 2k sin f0 sin feikd sin f0 L
1 ðk cos fÞL1 ðk cos f0 Þ
2 m¼1 u1 ðr; fÞ¼ pffiffiffiffiffiffi
X
1 2p ðcos f þ cos f0 Þ
bm X1
ð13bÞ am eikr
a am þ 2pik sin fL
1 ðk cos fÞ pffiffiffiffiffi
m¼1
m¼1 m
a k cos f kr
To cancel out the unwanted poles, the expansion coeffi- ð16Þ
cients am and bm should satisfy the following equations:
Similarly, if we use the change of variables a ¼ k cos t,
x ¼ r cos f and y ¼ r sin f, with f 2 ð0; pÞ, the saddle
Lþ
1 ðam Þ2 ðam Þ
am ¼ i m ¼ 1; 2; . . . ð13cÞ point of the integrand occurs at t ¼ 2p f, with
_
Mða mÞ f 2 ð0; pÞ, and the asymptotic evaluation of the transmitted
and field (4c) gives
pffiffiffiffiffiffi
u3 ðr; fÞ ¼ 2peip=4 k sin fLþ2 ðk cos fÞ
Lþ þ
2 ðam Þ1 ðam Þ X1
bm ¼ i m ¼ 1; 2; . . . ð13dÞ bm eikr
_ mÞ
Mða pffiffiffiffiffi ð17Þ
m¼1 m
a þ k cos f kr
In (13b) and (13c), the dot ( ) stands for the derivation with
respect to a. The application of the analytical continuation
principle together with the Liouville theorem to the
equations (13a) and (13b) produces 4 Computational results
and
Lþ ðam ÞLþ
2 ðam Þ
bm ¼ i 1
_
Mðam Þ
" #
L1 ðk cos f0 Þk sin f0 eikd sin f0 X
1
an
þ
ipðam k cos f0 Þ a þ am
n¼1 n
ð15bÞ
IEE Proc.-Sci. Meas. Technol., Vol. 153, No. 4, July 2006 171
Authorized licensed use limited to: ULAKBIM UASL - GEBZE INST OF TECH LIBRARY. Downloaded on March 23, 2009 at 05:49 from IEEE Xplore. Restrictions apply.
Fig. 5 Diffracted field amplitude against observation angle for
Fig. 3 Diffracted field amplitude against z1 for different values of d different values of distance d between two-part planes
Fig. 4 Transmitted field amplitude against z2 for different Fig. 6 Transmitted field amplitude against observation angle for
values of d different values of distance d between two-part planes
172 IEE Proc.-Sci. Meas. Technol., Vol. 153, No. 4, July 2006
Authorized licensed use limited to: ULAKBIM UASL - GEBZE INST OF TECH LIBRARY. Downloaded on March 23, 2009 at 05:49 from IEEE Xplore. Restrictions apply.
In this case, the two infinite sets of linear algebraic
equations in (15a) and (15b) are reduced to
½Lþ ðam Þ X
2 1
bn
am ¼ ð23aÞ
Mðam Þ n¼1 an þ am
and
" #
½Lþ ðam Þ2 L ikd sin f0 X1
1 ðk cos f0 Þk sin f0 e an
bm ¼ þ
Mðam Þ pðam k cos f0 Þ a þ am
n¼1 n
ð23bÞ
6 References
1 Abrahams, I.D.: ‘Scattering of sound by two parallel semi-infinite
screens’, Wave Motion, 1987, 9, pp. 289–300
Fig. 7 Transmitted field amplitude against observation angle for 2 Abrahams, I.D., and Wickham, W.G.: ‘The scattering of water waves
different values of d by two semi-infinite opposed vertical walls’, Wave Motion, 1991, 14,
pp. 145–168
3 Idemen, M.: ‘A new method to obtain exact solutions of vector
Wiener–Hopf equations’, Z. Angew. Math. Mech. (ZAMM), 1979,
59, pp. 656–658
and 4 B.uy.ukaksoy, A., Topsakal, E., and Idemen, M.: ‘Plane wave
pffiffiffiffiffiffi ip=4 diffraction by a pair of parallel soft and hard overlapped half-planes’,
u3 ðr; fÞ ¼ 2pe k sin fLþ ðk cos fÞ Wave Motion, 1994, 20, pp. 273–282
5 B.uy.ukaksoy, A., and Serbest, A.H.: ‘Matrix Wiener-Hopf factoriza-
X1
bm eikr tion methods and applications to some diffraction problems’, in
pffiffiffiffiffi ð20Þ Hashimoto, Idemen, Tretyakov (Eds.): ‘Analytical and numerical
m¼1 m
a þ k cos f kr techniques in electromagnetic wave theory’ (Science House, Tokyo,
Japan, 1993, chap. 6)
respectively. The explicit expression of L ðaÞ is given 6 B.uy.ukaksoy, A., Uzg.oren, G., and Birbir, F.: ‘The scattering of a
plane wave by two parallel semi-infinite over-lapping screens with
in [11] as dielectric loading’, Wave Motion, 2001, 34, pp. 375–389
pffiffiffi Kd aþiK iad 2p p
7 Senior, T.B.A., and Volakis, J.L.: ‘Sheet simulation of a thin dielectric
Lþ ðaÞ ¼ d e½ p lnð k Þ e p 1Cþln kd þi2 layer’, Radio Sci., 1987, 22, pp. 1261–1272
2sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 3 8 Idemen, M.: ‘Confluent edge conditions for the electromagnetic wave
Y1 2 at the edge of a wedge bounded by material sheets’, Wave Motion,
4 1 kd ad ad 2000, 32, pp. 37–55
i 5ei np ð21Þ 9 Braver, I.M., Fridberg, P.SH., Garb, K.L., and Yakover, I.M.: ‘The
n¼1
np np behavior of the electromagnetic field near the edge of a resistive half-
plane’, IEEE Trans. Antennas Propagat., 1988, AP-36, pp. 1760–1768
10 Mittra, R., and Lee, S.W.: ‘Analytical techniques in the theory of
L ðaÞ ¼ Lþ ðaÞ ð22Þ guided waves’ (Macmillan, New York, 1971)
11 Noble, B.: ‘Methods based on the Wiener–Hopf technique’ (Pergamon
with C ¼ 0:577215 being the Euler constant. Press, New York, 1958)
IEE Proc.-Sci. Meas. Technol., Vol. 153, No. 4, July 2006 173
Authorized licensed use limited to: ULAKBIM UASL - GEBZE INST OF TECH LIBRARY. Downloaded on March 23, 2009 at 05:49 from IEEE Xplore. Restrictions apply.