Uniform asymptotic theory of diffraction by a plane screen
Citation for published version (APA):Ahluwalia, D. S., Lewis, R. M., & Boersma, J. (1968). Uniform asymptotic theory of diffraction by a plane screen. SIAM Journal on Applied Mathematics, 16(4), 783-807. https://doi.org/10.1137/0116065
DOI:
10.1137/0116065
Document status and date: Published: 01/01/1968
Document Version:
Publisher’s PDF, also known as Version of Record (includes final page, issue and volume numbers)
Please check the document version of this publication:
• A submitted manuscript is the version of the article upon submission and before peer-review. There can be important differences between the submitted version and the official published version of record. People interested in the research are advised to contact the author for the final version of the publication, or visit the DOI to the publisher's website.
• The final author version and the galley proof are versions of the publication after peer review.
• The final published version features the final layout of the paper including the volume, issue and page numbers.
Link to publication
General rights
Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain
• You may freely distribute the URL identifying the publication in the public portal.
If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, please follow below link for the End User Agreement:
www.tue.nl/taverne
Take down policy
If you believe that this document breaches copyright please contact us at:
openaccess@tue.nl
providing details and we will investigate your claim.
SIAM J. APPI. MATI-I. Vol.16,No.4,July 1968 Printed inU.S.A.
UNIFORM
ASYMPTOTIC
THEORY
OFDIFFRACTION
BY
A
PLANE
SCREEN*D. S. AHLUWALIAf, R. M.
LEWIS
AND,.
BOERSMA:[:1. Introduction. The study of diffraction phenomena requires the solu-tionofanappropriateboundaryvalueproblemforthe reducedwave equa-tion or Maxwell’s equations. With few exceptions these problems cannot be solved exactly. Often useful approximate solutions are given by
geo-metricaloptics, but these solutions failto account for diffraction, i.e., the existence of nonzero fields in the shadow regions.
It
is now known thatgeometrical optics yields the leading termof a high-frequency asymptotic expansion of thesolutionof the boundaryvalue problem, and that higher
order terms account for diffraction. Keller’s "geometrical theory of dif-fraction"
[3]
providesasystematic meansof computing such terms.Keller’s theory has not onlybeen ofgreatpracticalvaluebuthas formed
the foundation for important further developments in the asymptotic theory of diffraction.
Many
of these developments have been motivatedby the
attempt
to overcome someofthedefects of thegeometrical theory ofdiffraction.Thesedefects,
suchasthesingularitiesat causticsandshadowboundaries,
arelistedat theendof3.
In
arecent paper[4]
Lewis andBoersma
presentedamethod ofobtaining a "uniform" asymptotic solution of problemsinvolving diffractionbythin screens. That work was largely motivated by an earlier paper of Wolfe[8],
who treated specialcasesinvolvingplane and sphericalwavesincidenton aplane screen, byasomewhatdifferentmethod.
More
recentlyBoersma
and
Kersten [1]
have extended the method of[4]
to the electromagneticcase, and Wolfe
[9]
has introduced a new method for the scalar problem basedonthe representationofthesolutionas anintegralovertheaperture.
In
severalrespects
the work of Lewis andBoersma
[4]
is incomplete. Only the first two terms of the asymptotic expansion were actuallyob-tained,
and it was conjectured that all terms could be obtained by the samemethod.Itowever
thecalculationswereprohibitively complex.It
was also conjectured thatallterms would beregular at theshadow boundaries,but thiswas proved only fortheleading term.
In
this paperwe complete*Received by the editors October 3, 1967. This researchwassupported by theAir
ForceOfficeof ScientificResearch underGrantAF-AFOSR-684-64andthe Office of Naval Researchunder Contract Nonr 285(48).
) CourantInstitute of Mathematical Sciences,NewYork University,NewYork,
NewYork10012.
:
DepartmentofMathematics, Technological University, Eindhoven, The Nether-lands.784 Do
.
AIILSWALIA
R. M. LEWI AND $o BOERSMAtheworkof
[4]
forthe special caseof screenswhichare portionsof planes.We
begin withthesameAnsat
introduced in[4],
butourtreatment of theAnsat
is significantly simpler. This enables us to obtain allterms of the expansion and to prove the conjectures.Except
for one reference to are-sultobtained in
[4]
ourworkhere is essentially self-contained.In
2
we formulate the boundary valueproblem,
and in3
we briefly summarize Keller’ssolution.In
4
wereduce theboundaryvalue problemto the determination of a certain double-valued function. This device,
which was first introduced by Sommerfeld
[6],
simplifies the remaining work.In
5
we introduce ourAnsatz
and derive the consequences of in-serting it into the reduced wave equation. There we state two theoremswhichassert the existence oftheintegrals that definetheterms of the ex-pansion and the regularity of the solution. These theorems are proved in
Appendix 2.
In
6
we present alternate forms of the solution, and in7
we compare ourresultswith Keller’s theory. Therewe obtain all terms of the expansion of the "diffracted
wave".
Keller’s theory yields only theleading term and involves a "diffraction coefficient"
D. We
find that ourleadingtermagrees withKeller’s andalltheterms canbedescribed simply in terms of successive diffraction coefficients
Do
D,
D1,
D2,Ex-plicit formulas for the coefficients
D.
are given. Appendix 1 contains abrief summary of a basic method for obtaining asymptotic solutions of thereducedwaveequation.
2. Formulationof theproblem.
We
considerproblemsofdiffractionbyascreen
S
whichlies inthe planex
0.Thescreenmay have one or moreapertures
of arbitrary shape or may consist of a collection of disjoint regionsofarbitraryshape. Thecbmplications
of thegeometry of thescreen will not concernus because our considerations willbe local.We
shall con-struct the diffracted field in a certain neighborhoodN
of the edge of a typical portion of the screen and shall ignore contributions from otherportions of thescreen aswellasthosedue tointeractionsbetween portions
of the screen. Such contributions will be considered in a later paper.
We
shallrequirethat theedgecurvex
x0(v)
beregular, i.e.,have derivative ofall orders. The parameter denotesarc length along theedge.An
incident fieldu0(x)
whichisasolution of the reducedwave equationThe neighborhoodNextendsuptothefirstcaustic pointalongeach"diffracted
ray"emanatingfromthe edge (see 5).
Thisrequirement canbe weakened.Weshallconstruct ourasymptotic solution to allordersandshow thatthefunctions inevery term areregular. Itoweverit can beshown that theconstruction canbecarriedoutto any given finite order andthe terms willhave anyspecified number ofderivatives if the edge function x0(v) has sufficiently many derivatives. Infact therequired order of differentiability ofx0()
DIFFRACTION BY A ILANE SCREEN 785
(2.1)
is prescribed. The total fieldu(x)
must then satisfy the followingconditions"
(2.1)
Au
-[-]u
0;(2.2a)
u=0 onS
or
(2.2b)
Ou/Ox
0 onS;
(2.3)
u hasa finite limitat the edge;(2.4)
u u0 is outgoing fromS.
Thuswe are in fact simultaneously consideringtwo problems correspond-ing to the two boundary conditions
(2.2a)
and(2.2b).
Condition(2.4)
isa form of the "radiation condition" which is more convenient for our asymptotic method. The definition ofthe condition is given in Appendix
1. The "edge condition"
(2.3)
is an essential part of the problem.It
iswellknown thatwithout it the solution isnot unique.
We
assume that the incident field has an asymptotic expansion, of the form(2.5)
u02
(i)-z(x),
Then
(see
Appendix1)
the phasefunctions(x)
satisfies the eiconalequa-tion
(2.6)
(Vs)
1,and theamplitude functions
z,(x)
satisfy the recurivesystemofransporl equations(2.7)
2V.s. Vz,+
ZmAS
--Az,,_..,
mO,
1, 2,..., z_ : O. Thesolutionsof theseequationsare discussed inAppendix 1.3. Keller’sasymptotic solution. According to Keller’sgeometrical theory of diffraction
[3],
the asymptotic solution of our diffraction problem isgiven by
(3.1)
where
(3.2)
(3.3)
and(3.4)
uu+u.+,
786 D. S. AHLJWALIA, I. M. LEWIS AND J. BOERSMA
The factor is one in the illuminated region of the incident wave and zeroin the (complementary) shadowregion.
We
assumethatthis wave is incidentfrom the regionx
<
0. Then theilluminatedregionincludes theregion
x
<
0 and that portion of the regionx
>
0reached by incident rays. Similarly is one in the illuminated region of the reflected wave(the
region reached by the reflected rays of geometrical optics) and zeroin the corresponding shadow region. The upper sign in
(3.3)
corresponds totheboundarycondition(2.2)
andthelower sign to(2.2b).
From (3.1)
we see that, in addition to the incident and reflected waves, there is a
"diffracted
wave"
givenby(3.4). In
order to describe thisfunctionwe must first discuss thetwo-parameter
family of "diffractedrays".
Theserys emanate from the edge. The diffracted rays througha point
x0(7)
ofthe edge generate a cone ofsemiangle
fl
fl(7)
withvertex atx0(7)
andaxistangent to the edge.
Thus,
for each fixed 7, 6, adiffracted ray is given by(3.5)
xx(,
,
)
x0()
+
u(, ),
where
U
is the unitvector(3.6)
U
cosfltl-t-sin/cosCt-
sintsinCh, -r=<
_-<
r.Here
h+/-o(7)
dxo/dr
is theunit tangent vector to the edge;t(7)
is the unitvector orthogonal to the edge, inthe plane ofthe screen, pointingawy from thescreen; and
t
isaunitvectorinthe direction of the negativex-axis.
Thesevectors are illustrated in. Fig. 1. The positive direction ofalong the edge is so chosen that h t.
;
t.
In (6), /()
is the ngle between theincident ray and the tangent to theedge at the pointx0().
Thus,
since Vs is the unit vectorinthedirectionof theincident ry, cosVs.h
In
fact(3.7)
Vs cosh
sin/cos0
t
sin/sin0
ta.
This equation merely determines the angle
0(v).
(See
Fig.1.)
If ndenotes the unit nor.malto the edge, then
t
=t=n, and the upperor lower signholdswhen thescreenis locallyconcave or convex.
In
eithercasethe curvatureis givenby 0 n
I,
where -h.1
=
0 is the "signedcurvature."
Since 0nandh
-0tl,itfollows that(3.s)
t,
t,
o.
Equation
(3.5)
defines a transformation from"ray
coordinates",
7, to Cartesian coordinatesx, x2,
xa.
The Jacobian(3.9)
jO(x,
x2,x)
Ox.Ox
X
O_x
DIFFRACTION BY A PLANE SCREEN 787
FxG. 1. Anglesand vectors atanedgeofthe screen. The vectorstl t t3,A,and B
areofunit length:tl istangent to theedgeofthescreenandpoints outoftheplaneof the figure,t.lies in theplaneofthe screenandpoints awayfromthescreen, and t3points in the directionofthenegativex3-axis. Theprojectionsofincidentand diffracted rays into theplaneofthe figure are shown. r go isthe angle between these pro-jections.The incident wavepropagatestothe right,i.e.,0-<
o
_
’.can be obtained
rom
(3.5),
(3.6)
and(3.8). A
brief calculation yields(.3.10)
j= sin.o-
1-t-where
(3.11)
sinIn
order to complete the description of Keller’s solution(3.1)
we must specify thefunctionsthat appear in(3.4).
Along thediffracted ray(3.5),
78 D. S. AHLUWALIA, R. M. LEWIS AND $. BOERSMA
(x)
is givenby(3.12)
s[xo(/)]-t-
a,wheres isthe phasefunction oftheincidentwave
(2.5).
The functionsaregivenrecursivelyalongthe diffracted rays
(see
(Al.14)
of Appendix bym(/, 6)
10
y(a’)/,,-1(’
(3.13)
2re(z)
y(z)-
-y()
)dz,
m=0,1,2,
...,
wheresin-
1-4-part integral
:f
in(3.13)
is defined in Appendix 1. Keller’s The finitemethodyields
(/, 6)
only form 0, hence only the leading term 20 of(3.4). It
isgiven by(
(3.
5)o
Dzo
Ix0(,)]
1+
where
D
is Keller’s "dractioncoefficient",
ei/4
(3.16)
D
--2
sifl
[sec
(
+
0)
sec(
A-
60)].
The upperorlower sigholds fortheboundarycondition
(2.2a)
or(2.2b).
Since increasesthstance
fromtheedgealongthedfractedrays, the lastterm in(3.1)
is clearly outgoing fromS.
Thereflectedwave u,is also clearly outgoing. Then,sinceu
u0(1
)u0
is nonzero only in theshadow regionof theincidentwave, we seethat
(3.1)
satisfies the outgoing condition(2.4).
Keller’s solution has been very useful and yields excellent agreement
with experimental results.
It
also agrees perfectly with the asymptoticexpansion of the few exact solutions that are known.
However
it suffersfromthe following defects"
(a)
As
canbe seenfrom(3.2)
and(3.3), u
isdiscontinuous across theshadow boundary of the incident wave
(the
surface that separates theillunated and shadow regions). Similarly
u
is discontinuous across the shadow boundary of the reflectedwave.(b)
Thediffracted wave becomesinfinite at both shadow boundaries, where v 0 and -v+
0, because the diffraction coefficientDIFFRACTION BY A PLANE SCREEN ?89
(c)
From
(3.15)
we see that the diffracted wave becomes infinite at the edge where 0; thus the edge conditionisviolated.(d)
The higher order terms,
m 1, 2,...,
in(3.4)
cannot bedetermined.
(e)
The value(3.16)
of the diffraction coefficient does not arise as an integralpart
of Keller’smethod;
ratheritis obtainedby comparison withthe asymptotic expansion of theexact solution of a "canonical problem," theproblemof diffractionof aplanewave byahalf-plane.
(f)
Thesolutionbecomes infiniteat thecaustic -pof thediffracted wave(see
(3.15))
as well as atany caustics of theincident and reflected waves.(g)
A
rigorous proof of the asymptotic nature of the formal solutionhas not been given.
Buchal and Keller
[2]
haveovercomedefects(a)-(e)
by boundary layermethods.
However
these methods yieldseparate
expansions in variousregions and requirerelativelycomplicated computations.
In
the succeeding sections we shall obtain, by relatively simple means, a single (uniform) asymptotic expansion which is free of defects(a)-(e).
However (f)
and (g) remain.Our
expansion isthesameasthat obtained byamorecompli-cated method in
[4].
The present method enables us to prove thecon-jectures madein
[4].
4. Thedouble-valuedsolution.Thesolutionofourdiffractionproblemis
facilitated by the introduction of a double-valued solution of the reduced wave equation.
A
similar device was used by Sommerfeld[6]
for the solu-tion of the half-plane diffraction problem.We
shall attempt to constructafunction
U
ofthe ray coordifiates z,,
whichsatisfies the conditions(corresponding to
(2.1)-(2.3))
.(4.1)
AU-I--kU
0 for>
0,(4.2)
U(a,
7,+
4r)
-
U(z,
7,),
(4.3)
limU(z,
7,)
Uo(7)
existsandis finitefor all7.-0
From
thetransformation(3.5), (3.6)
we seethat the periodicity condition(4.2)
makesU
a double-valued function of x.We
now define a single-valued functionu(x)
by setting(4.4)
uU(a,
7,)
U(a,
7, 2r),
-r_-<
_-<
r;and we observe that if
(4.1)-(4.3)
are satisfied, then u satisfies the con-ditions(2.1)-(2.3)
of the diffraction problem. (Condition(2.4)
will beverified
later.) In fact,
conditions(2.1)
and(2.3)
are clearly satisfied0 D. So
AHLUWALIA
R. Mo LEWIS AND Jo BOERSMA(3.6)
weseethatonS,
i.e.,
for =i=r, 00X.v
-sint
0(4.5)
We
assume th.at0<
fl
<:
r.Hence
(2.2a)
and(2.2b)
areequivalent to(4.6a)
u(=i=r)
0,
(4.6b)
u, (=i=r)
0.Usingtheupper sign in
(4.4)
fortheboundarycondition(4.6a)
we seefrom(4.2)
that(4.7)
u(r)
V(r)-
V(r)
0,
u(-r)
U(-r)
V(3r)
0.Similarly forthe other boundarycondition
u,()
U()
U.(2r
and
u,()
u()
u,()
o,
(4.8)
u(--r)
U(-r)
U(3r)
0.Thus theboundarycondition is verified inbothcases.
5. The uniform asymptotic solution.
We
shall construct the functionU
(asymptotically) in a neighborhoodN
of the edge defined as follows:Y
{x
x0()
+
u,
o
_<_
<
},
where cx isthe smallest positive value of such thatx
xo
+
IU
isa caustic point of the incident or diffracted wave. Thus that segment of
each diffracted ray
(3.5),
beginning at the edge and terminating at the nearestcausticpoint,lies inN.
We
shallalso referlaterto the neighborhoodNo
{x
x0()
+
u,
o
<
<
},
from whichthe edgeitselfhas been deleted.
In
order to find the functionU
we introduce theAnsatz
g
(5.1)
where --1/2ei.14
(5.2)
f(x)
--ice-
edr,
c rThis form wassuggested by the study of the exact solutionofthehalf-plane diffraction problem (see [4]).
DIFFRACTION BY A PLANE SCREEN 91
and
(5.3)
2
s.The functions s and
z
are the phase and amplitude functions of the in-cident wave(2.5),
and is the phasefunction ofKeller’s diffracted wave.It
is givenby(3.12).
The functionsv
areto be determined.It
is easyto show that s>-
0inN (see [4,
Section2, Lemma 1]).
Hence
is real anddouble-valuedinN.
We
note that ifweset(5.4)
then
(see
Fig.1)
sin(’/2)
cos((
+
0)/2)
wnishes he shadow boundary,wheret
2nr,
n 0, 1, :t:2,....
Furthermore
he ineidennddiffracted rays coincideonheshadow boundary,ndhence here.
I
follows t,ha # wnishes on he shadow boundary, nd w my choose2
Then#sisfies he periodicity condition
(4.2).
In
fe he firsterm
in(5.1)
sisfies he sme condition. This follows from he fc ht ndz
re single-wlued functions ofx,
hence hve periodLter
we shll verify h he seconderm
in(.1)
lso sisfies(4.2).
First however we inser(.1)
into he reduced wve equation, using(5.6)
f’(x)
--ic-2ixf(x)
to eliminate derivatives of
f.
The calculation is simplified if we set(5.7)
gea’
f k20
hck-l/2e
ik$and
(5.s)
v
Here
we sum over all integer values of the repeated index m, and it isunderstood that
z
andv
vanish identically for m -1, -2,....
In
computing derivatives ofU
wenote that(5.9)
Og ik(
OsOx
- -
h--O
x i-
hThenitis easyto show that
(4.1)
issatisfied, provided(5.10)
(5.11)
and(5.12)
(V)
1, v0.(vs
+
v)
o
792 D. S.
AHLI:IWALIA
R. M. LEWIS AND J. BOERSMAwhere
(5.13)
qm 2V0.Vzm
-t-
zA0.
In
verifying(4.1)
we also madeuseof(2.6)
and(2.7).
Now (5.10)
isiust
the eiconal equation for,
and is clearly satisfiedby(3.12).
Furthermore(5.11)
is satisfied because, from(5.3),
20 V0 VVs,
and2 v0.
(vs
+
v)
v
vs). (v
+
v)
(5.14)
(v)
-
(vs)
1- 1o.
Thuswe areleftwith
(5.12)
whichweshallusetodeterminethefunctionsv. We
firsttransform(5.12)
byusingtheidentityd dy
[j
1112
(5.1.5)
Atelog[jl
2Y-:d--,
Y-sin
which follows from
(A1.9).
Here,
since is the phase function of the dif-i’racted wave,j is the Jacobian of the transformation defined by the dif-fracted rays.It
is given by(3.1.0).
Now,
since V.Vv, dv,/da,(5.1_2)
becomesd Y
(Ave_:
-t-
q)
(5.16)
d--
(yv,)
From
(4.3)
we see thatv
must befinite at a 0, and from(3.10)
and(5.15)
we see that y vanishes at a 0. Therefore integration of(5.16)
yields(5.17)
vm(a)
2y(a)
y(--A_:
+
q)
da’,
m 0, 1,2,...,
provided the integral exists.
In
(5.17)
the dependence onthe ray coordi-natesn
and isnot explicitlyindicated.In
Appendix 2 we shall prove the following theorems.(The
definitions of a andNo
aregivenat thebeginningof thissection.)THEOREM
1,For
everymO,
1, 2, the integral(5.17)
existsfor
0<
a: andlim
v(a)
0.-0
TEOnEM
2.U
is a regularfunction
of
x inNo
and itsatisfies
(4.3).
Since(4.1)
is satisfied by construction and Theorem 2 establishes thevalidity of
(4.3),
it remains to verify(4.2).
We
have already seen thatDIFFRACTION BY A PLANE SCREEN 793
v(b
-t-
4r)
-
vm
(b).
This can be proved by induction on m beginningwithra 1.(v_l clearlysatisfiesthe periodicity conditionsince itvashes identically.) Since
z
is a single-valued function ofx,
it is 2-perioc in;
hence it follows from(5.13)
that q(4)
q(4).
If now we make theinductionassumptionv,_( % 4)
v_()
wesee from(5.17)
that
v(
+
4)
v().
According to
(4.4)
our uform asymptotic solution of the diffractionproblem
(2.1)-(2.4)
is nowgiven by(5.18)
u(x)
where
U(a,
,
)is given by(5.1),
(5.2),
(5.3), (5.5),
(5.15), (5.17)
and
(5.13).
The presentsolution(5.18)
satisfiesthecontions(2.1)-(2.3).
It
only remains to be verified thatthe outgoing condition(2.4)
issutis-fled.
For
that purposeweshallshowthataway fromtheshadowboundariesand fromthe edge thesolution
(5.18)
reduces to(3.1).
At
thesame timewe shall verify Keller’s theory and obtain the higher order terms in the
expansions
(3.4). We
begin with the asymptotic expansion off(x),
whichcanbe obtained from
(5.2)
by integration by parts"(5.19)
f(x)
’vo(x)
cx-
(-)(ix:)-’,
x.
Here
(5.20)
()0
1,
(),
(
+
1)... (+n- 1),
n=1,2,3,
...,
and
v0(x)
istheunitstepfunction.Thusv0(x)
1for x>
0andv0(x)
0 for x 0.Except
near the shadow boundary andthe edge,where 0 0,kO
islarge, andwe mayuse(5.19)
in(5.1).
Thisyields(5.21)
U
v
cs
+
2
U0+
-1l ei$
m-o
(i)-mOm,
whereu0istheincidentfield., givenby
(2.5),
andE-(5.22)
In
theinterval<
0, which(see
Fig.1)
coincides withtheilluminated region of theincidentwave.Silarly, inthesameinterval,yo cos
2 yo --cos 2
illumi-794 D. S. AHL’UWALIA: It. M. LEWIS AND ft. BOERMA
natedregion ofthe reflectedwave. Therefore
(5.23)
Thus away from the shadow boundaries and from the edge we see
from
(5.21)
and(5.23)
that(5.18)
reducestothe (nonuniform) asymptotic solution(3.1)-(3.4),
where(5.24)
and
9m
is given by(5.22). Hence
the outgoing condition is satisfied. 6. Alternate forms of the uniform expansion.We
first obtain a usefulalternateexpression for qm which is given by
(5.13). From. (5.3)
we see that 1(v.-
vs)
0=(6.1)
V0 A--As 1 V.Vs 20 20Hence
(6.2)
qm(V
Vs).Vz
1<A
AS 1V.Vs)
0+
z
0 03But
d()
1 1--Vs.V(6.3)
d-
V0.V
20Thus,
from(5.15)
and(6.3),
(6.4)
y-1
d(y.)_zmA
Vz,.V(1--
Vs.V)
20
+
0 z, 20Now
from(6.2),
(6.4)
and(2.7)
we seethat(6.5)
q Ifwe insert(6.5)
in(5.17),
weobtain 1fo
I
Az_d,+
1Fyzl
(6.6)
vm
=
y --Av_+.20
A
-L
0A0"
By
expanding 0 and y for small(see (7.19)
and(7.20))
we find that(6.7)
limL21/2
Sill-
sia.
DIFFRACTION BY A PLANE SCREEN
Hence
lf
on"
[
AZm_ldo.!
ZmZm[X0()]
6.8 20:
2-
2ysin(’/2)sin"
For
m 0,(6.8)
becomes z0zo[x()]
(6.9)
Vo-
+
23/.ysin(/2)sin"
In
theimportant specialcaseofanincidentplane wave,z0-
1 andzm
-
0form 1, 2,
....
Then(6.8)
simplifies to1
yAv,_t
do/,
m 1,2,..-(6.o)
-7. The nonuniform expansion.
In
5
we obtained the nonuniform ex-pansion(5.21),
(5.22)
forU
valid away fromthe shadow boundariesandfrom the edge. Using the results of Appendix 1 we shall now derive a
simple recursive formula for the coefficients
Ore.
According to(Al.14),
(Al.15), 9()
can berepresented by(7.1)
()
y(a)
2y(a)
yAO,_do",
mO, 1,
2,...,
v_=---0,
where
(7.2)
h
finy(z)O(z).
Here,
yJ I//sin
is given by(3.10).
Using(5.22),
the initialvalue,
can be expressed in terms of the known coefficientsv
and zm, viz.,E
(7.3)
cfinyv--
z_,aO nO/n
Since
v
O(Theorem 1)
and y O as a 0, thefitepart
(7.3)
re-duces to(7.)
Here
the,()
arelinear operatorsdefined by c()
finO-,,-X
yz).
(7.5)
.
z.
0
For
example, from(6.7),
c y lim )oz=
o
(7.6)
cz(xo)
sin sin(/2)
sec+
?
z(xo).
27
D. S. AHLUWALIA, R. M. LEWIS AND J. BOERSMAThus 5)0is a multiplication operator.
However,
forn>
0, 0. is a differ-entialoperator,as we shallseeshortly.Ifwe nowinsert
(7.1)
into(5.24),
we seethatwhere
(7.8)
()
(2v
)
D,,z_,,.Here
thediffraction coefficients
D,,are linearoperators definedby(7.5)
and(7.9)
D--
D(b)
:t:D.(2r-
),
--r=<
-<
r.Thus from
(7.6),
(7.10)
Dozo
Dzo(xo),
where
D
isKeller’s diffraction coefficient(3.16),
and(7.11)
Dzo(xo)y
-.
We
note that(7.7)
and(7.11)
agreeexactly with(3.13)
and(3.15).
ThuswehaveverifiedKeller’s theory.
The higher order termsin the expansion of the diffracted wave cannot beobtainedby Keller’smethod.
Here
we seethat theyaregivenrecursively by(7.7),
(7.8), (7.9)
and(7.5).
In
conclusion we may state that theuniform asymptotic solution as derivedin
4
and5
is not only of great value in itself, but it is also fundamental for the completion of Keller’snonuniform asymptoticsolution. The initial value/tin
(3.13)
and(7.7),
which was unknown until now (except for m0),
is directly obtainedfromtheuniformasymptoticsolution.
To
illustrate the application of this nonuniform asymptotic solution, we complete thecorrectionterm2.
This requires theevaluation ofc
fin0
-
yz.(7.12)
)zTo
evaluate thefinitepart
we expand 0, y, andz forsmall a. First wesee from(3.5)
and(3.6)
that(7.3)
zz(x0)
+
U. vz(x0)
+
0(
)
and
DIFFRACTION BY A PLANE SCREEN 797
Here
(7._5)
(7.16)
1 b-
,I OyOyi(x0)
U
U,
U
U, U,
U)
(cos
,
si cos,
sin sin),
and theyare Cartesian coordinates corresponding to the basevectors
t,
i 1, 2,3.
From (3.7)
we seethat(7.17)
U.
Vs(x0)
cos sinfl
cos(
+
o)
cos+
sin2/
cos.
Since
s(xo)
+
,
and 1 cosi"
sin(/2),
(7.18)
s 2sin:/
sir 12sinS3 sin(./2)
+
0(a
2)
Now,
(5.3)
and(5.5)
yield(7.19)
0-
(2a)
-3/ (sinfl
sin 1--sirt2flsin2(./2)
Furthermore(3.10)
yields(7.20)
yJ
1/2sin
rl+pTO(cr
2)
We
nowform theproduct of(7.13),
(7.19)
and(7.20).
Then we delete thesingularterms (negative powers of)
and thenlet-
0. This yields fin,*o 0-
yz, and(7.1.2)
becomesI)1Z "--2--7/2c(sin
f
sin(/2))
-s(7.21)
[(1
3b)
z(x0)
+
U.Vz(xo)l"
+
4sia/Here
ro,
b is given by(7.15)
andU
isgiven by(3.6).
The last term in(7.21)
illustrates the fact that the0
are in general differ-entialoperators.We
shall not complete the evaluation of in general, because the in-tegral in(7.7)
form 1 cannot be explicitly evaluatedin general.How-ever there are two important special cases which can be evaluated.
We
consider first thecaseinwhich
o
0 (grazing incidencetoward thescreen: see Fig.1)
and the second boundary condition(2.2b)
holds.In
this case we see from(3.16)
that thediffraction coeificientDo
D
vanishes. Then0
----
0 and it isespeciallyimportanttoevaluate because it nowprovidestheleading termin
(3.4). From (7.7),
(7.8)
and(7.9)
we seethat(7.22)
1
I)ZoD
:D (ok)
+
(2
),
r<-
<.
798 D. S.
AHLUWALIA
R. M. LEWIS AND J. BOERSMAHere
ff)lz5)1()z
is givenby(7.21)
with0 0, i" r,
andsin(i’/2)
cos
(/2).
Since cos((2r
)/2)
-cos(/2)
and p(2r)
p()
(see
(3.11)),
the first term in(7.21)
contributes nothing to the sum in(7.22). Furthermore,
since the incident rays are tangent to the screen, Os/Oya 0onS
and02S(Xo)/OyOya
0,i 1, 2,3.It
follows from(7.15)
and
(7.16)
thatb(2r
)
b();
hence thesecondterm in(7.21)
alsodoesnot contribute.
Now
from(3.6)
we seethat(7.23)
It
follows thatU()
.VzU(2r
).Vz
-2sin/3 sin Vz.t3 OZ 2sin/ sinb aox--"
e/4sinOzo
(7.24)
91
z0 4%/-
sin2/3
cos(q/2)
(x)
and if we iert
(7.22)
and(7.24)
into(3.4)
we obtain, for the leadingterm of the ffracted wave,
(
(7.25)
l(/e
1+
D’
OZo
(Xo),
where
(7.26)
D’
-e-
sin(4/2)
2/-
sin cos(/2)"
Thisresultwasalsoobtainedby Keller by expandingtheexactsolution ofa special diffraction problem.
It
is easily seen that(7.25)
and(7.26)
agree exactly with(12)
of[3]. (We
must firstcorrect anerror in the last part of(12)
which has the wrong sign. Then the results agree because 0The second specialcase occurswhen 0 r (grazing incidence from the screen: see Fig.
1)
and thefirst boundary condition(2.2a)
holds.In
this case we see from(3.16)
that the diffraction coefficientDo
D
again vanishes. Again40
0and1
provides theleadingtermin(3.4). Now (7.7),
(7.8)
and(7.9)
yieldD1
z0D1
)()
)1(2r-
),
--r <q
< r.(7.27)
z
Y
Here
lz l(q)z is given by(7.21)
with0 r,-
andsin(/2)
-sin
(/2).
Sincesin((2r
)/2)
sin(/2)
andp(2)
p(),
the first termin(7.21)
contributes nothing to the sum in(7.22).
Since againb(2r
)
b(),
the second termin(7.21)
also doesnot contribute.It
DIFFRACTION BY A PLANE SCREEN ?99
then follows from
(7.23)
thate’1sin
Oz--2"
(Xo)
(7.28)
D1
z04/-
sia sin(/2)
Ox3The leadingtermofthe diffractedwaveisnowgivenby
(7.25),
4the
-
cos(/2)
(7.29)
n’
2
/-
sia sin(/2)"
This result wasalso obtained by Keller.
(If
we set 0-
/2,
n2,
and correctsomeerrorsin(19)
of[3],
itthenagrees with(29).)
Appendix 1. Asymptotic solutions of the reduced wave equation.
We
consider solutionsu(x)
of(AI.1)
Au
ku
0whichhave an. asymptoticexpansion of theform
(A1.2)
ue’kS(’_,
(ik)-mZm(X),
/C--.
m----0
By
formally,substituting(A1.2)
into(A1.1)
we findthat (A1.1)
is satis-fied if(Vs)
1,(A1.3)
and
(A1.4)
2Vs.Vzm
z,As -Azm_l, m 0, 1, 2,..., z_l-=
0. The solutions of(A1.3)
and(A1.4)
may be described conveniently by introducing two-parameterfamily of straight lines (rays)(A1.5)
xx(a,
.,
a)
whichare orthogonal to a level surface
(wave
front) s(x)
So ofs. Thelabeling
parameters
a2, aarefixedon a ray andadenotesarclength along the ray from the givenwavefront in the direction of increasings.Thenwe seefrom(A1.3)
that(A1.6)
six(a,
z,a)]
so
+
.
This provides the solution of
(A1.3).
It
is easily seen that the rays are orthogonal to every wavefront s const.An
asymptotic solution of(AI.1)
of the form(A1.2)
is saidto be out-goingfromama,nifo]dM
ifallof theraysof thefamily associated withthesolution, emanate from
M
and o each ray, in a neighborhood ofM,
the phasefunction sincreases with distance from 214 alongthe ray.800 D. S.
AHLUWALIA
R. M. LEWIS ,AND J. BOERSMFor
eachm,
(A1.4)
is an ordinary differential equation alonga raybe-cause
Vs.
Vz
dz/da.
This equation can be conveniently solved by introducing theJacobianofthe"ray
transformation" xOx
Ox
Here
wehave used theexpansionof thedeternantinter of cofactomoftheithrow, i 1, 2 or3. Since thedeternantvanishes if
two
rows areidentical,
wehave(Al.8)
Ox
Ox
j=t
cof
where isthe Kronecker symbol.
It
foows thatz
co 0X =z
0rox
(A1.9)
0 dx
Thus,
from(A1.4),
d
(
jl
J
da
+
2[2Vs.Vz+ zhs]
( 1.1o)
Zm-I
2
By
integration (longrays)
we obtain thereoursiveformulsforthez’s"
z,()
J()
"
1i
j(,){l
m= 0,1,2,....
Here
wehavenot indicated the dependence of all quantities ona and a.In
generalwe canof course takea0 0 in(AI.ll).
However
ifj(0)0,
the pointa 0 iscalledacaustic pointanditcanbeshownthat the integral
in
(AI.ll)
wouldthen divergeat the lower endpointa’
0.To
avoid this difficultywe introduce afini
partintegral defined asfollows"For
e 0 letf(e)
have an asymptotic expansion, in powers (perhaps fractional) of e as e 0. Leftf(e)
denote the singular terms (negative powers ofe)
ofthisexpansion.We
definethefinite
partoff(e)
as 0by(Al.12)
finf(e)
lim[f(e)
The present method of solution of the equations (A1.4) is different from the
method usedin [5]andelsewhere. Thelattermethod led toasolution containingthe
expansion ratioda(ao)/da(a), wheredastands forthecross-sectional area ofa tube ofrays. Thesolutionsareequivalentbecause j(ao)/j() da(ao)/da(a).
DIFFRACTION BY A PLANE SCREEN 801
Now
ifg(x)dx
is divergent or convergent atx 0, we define thefinite
partof
the integral as(Al.13)
g(x)
dx fin.g(x)
dx.-Ifa 0 isacaustic point,thesolution
(A1.11)
is meaninghflfora0 0.Let
usnowtakethefinitepartof(A1.11)
asa0-*0.Then"
1o
o"
(hl.14)
z,,,()
lj()I
3--
az.,,_,(’)
d’,
m 0, 1,2,...,
,(z, .)
fin]j(z0)
a0-0
The initialwlue my be chosen
o
mee
heboun.dryconditions ofhe problemfor(A1.1).
For 0 he integralerm
inA1.1]
is mising. If 0 ino
cusicpoint, he inegrM inA1.14
inordiaary integralndhe finite
pr
of(A1.15)
reduceso
nordinarylimit,(A1.16)
j(O)
]l[2Zm(O).
It
isthen clear that(A1.14)
reduces to(A1.11)
with a0replaced by zero. Appendix 2. Proofs of theorems.In
this Appendix we shall proveTheorems iand2whicharestatedin
5.
In
thebodyofthepaperwemade heavy useofthe"ray
coordinates",
v, defined by the transformation(A2.1)
xx0()
+
U(,, ),
where
U
is aut
vector inthe rection of the dracted ray. ThU
is givenby(3.6)
or,intermsoftheut
vectorst,
A,
B (illustrated
in Fig.1),
by(A2.2)
U
cosflhW,sincosAW sinflsin
Here
it is convenient to introducea new set of coordinates w, w,w
de-finedby
(A2.3)
w
,
v
(2a)
/sin(/2),
va
(2a)/
cos(/2).
Thus
(A2.4)
2a W,
2cosw
asin
W3,and,
from(A1.1)
and(A1.2),
02 D. S. AHLUWALIA R. M. LEWIS AND ft. ]3OERSMA
Here
x0 andthe orthogonal unitvectors h,A,
B
are functionsof nl,and
(A2.5)
defines a transformation xx(,n,
/,/.).
This transforma-tion maps the/-space onthedoubly-sheeted x-space.Two
points(m,
=t:/.,j) have thesameimage in x-space.
In
ordertocompute
thegradientand Laplacian operatorsin the newco-ordinateswe first
note
that(A2.6)
A
-cos0 t2 sin0
hence
(3.8)
yields(A2.7)
.
coso
h
+
0 B,
It
follows thatB
sin0t. coso
t
sino
h6o
A.
(A2.8)
(A2.9)
(A2.10)
where(A2.11)
(A2.12)
(A2.13)
:
X10X//0/1
(1
-
el)h
-t-
e A
A-
ea
B,
x2
Ox/O12
V2cosf
h WsinA
A-
/sinfl
B,
xa
0x/0
v.
cos/ hA-
/.sinfl
A
-t-
/sinfl
B,
I"
el -sin
3[-(
+
,)
+
,
cos q,o(sin
e2 cos
(,
,:)
+
-
cos cos0(,2
+
o
sin2a,ea
cos fl2ya cos sin0(y2+
W:)
+
o
sin(Va
The Jacobian
J
O(x,
x2,xa)/O(m,
va,72)
of the transformation(A2.5)
can be computed directly from
(A2.8)-(A2.13). However
it is simplertouse
(A2.2), (A2.3)
and(3.10),
whichyieldsin./
(1
+
-
jO(x,x2,
x)_
O(x,
x2,x)
\
o(,
,
)
o(,
,
)
(A2.14)
j
0(,)
1j.0(,
)
2The metric coefficients g of
(A2.5)
are defined by or(g)
(A2.15)
gx.x
0 0
kO]Here
the accentdenotes the transposedmatrix. Clearly,[
DIFFRACTION BY A PLANE SCREEN
The reciprocalcoefficients
g’
aredefinedby(A2.17)
(g)
(g)-i
org"
gwhere
G
;G
is th cofactor ofg. Then(see,
e.g.,[7])
for arbitrary functions,
7,(A2.18)
and(A2.19)
1 0(
..OOZe)=
j_l 0(
..00)
.0,0_
1 0,O
We
now introduce two classes of functions g(nl, 2,3). We
shall saythatgis anoddor evenfunction ifitisregularinaneighborhood of theedge
m
n
0(i.e.,
can beexpressed as apower seriesinn
and n3withco-efficients thatareregular functions of
v)
and if(A .20)
respectively.The definitionshavesome immediate andusefulconsequences"
If g isodd, then g(, 0,
0)
0. Theproductoftwooddfunctions iseven, etc.From (A2.5)
weseethatx iseven; henceifg(x) isregularinaneigh-borhood of the edge, theng[x(w,
,
w)]
is even. From(A2.14),
(A2.4)
and
(3.11)
it is easyto show that(A2.21)
zJ-
iseven.In
ordertoproveanimportantlemmaabout theregularityofthefunction defined by(5.3)
and(5.5),
we introduce that segment $ of the shadow boundary thatlies inthe neighborhoodN
whichwas defined at thebegin-ning of
5.
In
termsofthecoordinates(w,
,.,
va)
weseefrom(A2.3)
thatLEMMA
1. 0 iS a regularfunction
of
(,
,
)
in a neighborhoodM
of
$. Furthermore 0is odd.
Proof.
Let
(A2.22)
U
cos/h
-{-sinfl A
cos/h
sint
cos0h
sin/
sin0
(A2.23)
U
B
sin0t-
cosO0t,(A2.24)
U
U, X
U
sin/h-cosfA.
Then
U,
has thedirectionof theincident ray(see (3.7)),
andfrom(A2.2)
we seethat inthe
U,,
U,,
U
basis0 Do S. AHLUWALIA R. M. LEWIS AND J. BOERSMA
We
consideranarbitrarypointP
x0-
zU1
ontheshadow boundary andaneighboring pointx x0
-t-
aU.
The differenceisx-
(u- )
(A2.26)
g[sin2(cos
1),
sinfl
sin’,
cos/ sinfl(1
cos)].
ttence,
from(A2.4),
(A2.27)
h (hi,h., h3)
(-
sin2
2,
sinfl
.3, cos sinfl
.).
Now (x)
s(x0)
W
s(P);
therefore by Taylor’stheorem,
providedP
isnot a caustic point of the incidentwave,1
(A2.28)
(x)
s(x)
s(P)
s(x)
--
.
sil’"n(P)
h,
hn.
Sinceat
P, (s,
s,s)
Vs(1,
0,0),
weseefrom(A2.27)
and(A2.28)
that
(A2.29)
(x)
s(x)
(
sin)
s
.
(
sinf)
-t-
r,whereeveryterminrcontainsafactor
w.
.
Thus(A2.30)
(x)
s(x)
(,
sinwherep 1 s2v3 isregular ina neighborhoodof
v.
0andeven.Furthermorewesee from
(A2.4)
that,onthe shadowboundary where 0, v3 2 and(A2.31)
p 1We
now use the followingidentity which is givenby[4,
(18),
Appendix2]"
(A2.32)
(p
-t-
)(pa-- )
s2+
o2-t-
o3 PHere
p.,p3 are the principal radii ofcurvature
of theincident wavefrontat
x0(v). It
followsfrom(A2.31)
thatonthe shadow boundary1
+
(A2.33)
P(1
-
a/p)(1
-t-At
the edge,a 0andp 1. Since p can vanishonly at the caustic pointa -p and is continuous except at the caustic points a -p and a
-p3, weseethatp is finiteandpositive in$,hencein aneighborhood
M
of$.
From
(5.3), (5.5),
(A2.3)
and(A2.30)
we nowseethatDIFFRACTION BY A PLANE SCREEN 805
Since p isregular and positive atr/,. 0,weseethatt isaregular function
of
(r/l,
r/s,r/3)
ina neighborhoodof r/,. 0;and since p is even, tis odd.COROLLARY
1. tis aregularfunction
of
x(xl
x2x3)
inNo.
Proof.
From (5.3)
we see that 0 is a regular function of x except atcaustic
(where
sor failsto be regular) andperhapsattheshadowbound-ary where s and 72 0.
But
fromLemma
1, in a neighborhoodM
ofthe shadow boundary segment $, t is a regular function of (r/z,r/,
r/3),
hence of
x,
except where the JacobianJ
vanishes.From (A2.14)
we see thatJ
vanishes only atthe caustic -pand attheedge 0.Hence
0isaregularfunctionofx inNo.
Proof
of
Theorem 2. Thefunctionf
defined by(5.2)
is entire and thez
and areregularfunctionsofx
except
atcaustics.Hence
from Corollary 1 the first termin(5.1) is
regularinNo.
Theregularity of the second termcan beproved by induction:Ifvm_isregular in
No,
thenAv_1 isregular, andfromCorollary1 and(5.13)
wesee thatqm is regular.Thusfrom(5.17),
(3.10)
and the formulayJ I"/sin/,
v
isregularinNo.
Condition(4.3)
followsfromTheorem 1.
The proofofTheorem 1 isbased onthreemorelemmas.
LEMMA
2.(i)
Ifi 1,j 2,3orj 1,i 2, 3, thenJ-G
isodd.(ii)
If
i j, thenJ-G
is even.(iii)/f
i 2,j 3 orj2,
i 3, thenJ-G
is even.Proof.
From (A2.15), (A2.8),
(A2.9)
and(A2.10),
(A2.35)
gn(1
q-
el)
q-
e
q-
e3,(A2.36)
g=(1
q-
e)r
cos/ en,sinfl -k
car/asin./,(A2.37)
gla(1 q-
e)r/a
cosq-
e=r/asinfl
-t- ea.w
sinfl,(A2.38)
g:r/
q-
r/
sin’,
g. vr/cosfl,
ga r/aq-
r/
sin:ft.
Let
P,
Q,
R, S
denote nth degree homogeneous polynomials in r/,., r/a with coefficients that are regular functions of r/ r/.From (A2.14),
(A2.4)
and(3.11),
(A2.39)
J
sin/
(
q-
r/a)
1P(r/,.,
r/)}
hence
(A2.40)
From
(A2.11)-(A2.13)
we obtain bystraightforward calculation(A2.41)
e2r/3+
ear/2(r/22
.gf_732)Ql(r/2,
73),
(A2.42)
e:r/a-
ear/.(r/2
--
r/a)Rl(r/2,
806 D. S. AHLUWALIA, R. M. LEWIS AND ft. BOERSMA
Now
from(A2.35)-(A2.38)
wecompute G
=
cofactor (g) using2\--III
(A2.41)-(A2.43).
We
find, e.g., that(2
-l-
is even; hence from(A2.40),
J-1Gll
is even, etc.LEMMA
3.If
a isoddand biseven,then Va. Vbisodd.Proof.
From (19),
(A2.44)
aVa.Vb(J-1)(J-1G)
Oa ObBy
using(A2.21)
andLemma
2wefind" in case(i),
Oa Ob iseven, hence in case(ii),
Oa Ob isodd, hence in case(iii),
(aj_I)(j-1G)
Oa Ob(aj_)(j_lG)
Oa Ob is odd; isodd;
(A2.48)
y.(-Av,,_l
A-
q,)
a-l/a[/,
(2a)/sin (’/2),
(2r)
1/2COS
(’/2)1
has an expansion in nonnegativeintegral powers of
,
i.e.,
is regularinThus
Since
z
isregularinaneighborhoodoftheedge,z
iseven.From Lemma
1,is
odd,
and fromLemma
3,zV0.Vz
isodd. Furthermore,fromLemma
4,zzA0
isodd;
hencefrom(5.13)
we seethatzqisodd.We
shallprove byinduction that for each m the integral
(5.17)
exists andv
is odd. Theassertion is clearly true form -1because v_ v_2 q_ 0. If V m-is
odd,
itfollowsfromLemma
4that zAv_isodd;
hence from(A2.46)
wesee that
a,(,
.,
a)
isodd, where(A2.47)
a
o)/2y
(--
Ym-1
"2t"q.,,)
o-12y[
trAYm-1
"3t- aq,].LEMMA
4.If
a isodd, thenrAa
isodd.Proof.
Let
J-1Gi Oa/O
h.
Thenit is easilyseenfromLemma
2thath is odd andh and h are even.
It
follows thatOh/O
is odd and from(A2.18)
that zAa(rJ-)Oh/O
isodd.Proof of
Theorem 1.From (A2.39)
and(A2.4)
weseethat(A2.45)
(z-lJ)/
isevenand(-lJ)-l/
iseven.Since y sin
j
j/2
(j/2),
itfollows that --12(A2.46)
y is evenand2y-
iseven.Oa
DIFFRACTION BY A l?LANE SCREEN 807
1 1/2_.-1
where at 0. Thus the integral in
(5.17)
exists, andvm
y o/,,(A2.49)
We
seethato/misodd becauseam
isodd.It
followsfrom(A2.46)
thatvm
isodd. This completes the induction argument. Since
Vm(W_,
VR, q3) isodd,
we seefrom
(A2.3)
that(A2.50)
lim v,v(vl,
0,0)
0.Acknowledgment. The authors wish to acknowledge their indebtedness
to
J.
B.
KellerandD.
Ludwigforseveral helpful suggestionsmade duringthecourseofthiswork.
REFERENCES
[1] J. BOERSMAANDP. H. M. KERSTEN,Uniformasymptotic theoryofelectromagnetic
diffraction by a planescreen, to appear.
[2] R. N. BucHhLAND J. B.KELLER, Boundary layer problemsindiffraction theory,
Comm. PureAppl.Math., 13 (1960), pp. 85-114.
[3] J. B. KE, Thegeometrictheoryofdiffraction,J. Opt. Soc. Amer., 52 (1962), pp. 116-130.
[4] R. M. LEwIsANDJ. BOERSMA, Uniform asymplolic theory of edge diffraction, J.
MathematicalPhys.,to appear.
[5] R. M. LEwIsANDJ.B.KELLER,Asymptotic methodsfor parlialdiffercnlial equa-lions: thereducedwaveequalion and Maxwell’sequalions,Res. Rep.EM-194,
NewYork University,Ne York, 1964.
[6] A. SOMMERFELD, Optics, Lectures o- Theoretical Physics, vol. IV, Academic
Press, New York, 1964.
[7] J. A. STaATTON, Electromagnetic "l’heory,McGraw-ttill, NewYork,1964. [8] P. WOL,Diffractionofascalarwavebyaplane screen, thisJournal,14(1966),pp.
577-599.