Файл: Ufimtsev P. Fundamentals of the physical theory of diffraction (Wiley 2007)(348s) PEo .pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 28.06.2024
Просмотров: 939
Скачиваний: 0
2.3 Conversion of the Series Solution to the Sommerfeld Integrals 43
Figure 2.6 Integration contour in Equation (2.51).
u(kr, ϕ |
ϕ0) |
v(kr, ϕ |
ϕ0) |
+ |
e |
− |
ikr cos(ϕ |
− |
ϕ0) |
with π − ϕ0 < ϕ < π + ϕ0 |
|||
u(kr, ϕ |
− ϕ0) |
= v(kr, ϕ |
− ϕ0) |
|
|
|
|
||||||
|
+ |
= |
|
+ |
|
|
|
|
|
|
|
|
(2.53) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u(kr, ϕ ϕ0) v(kr, ϕ |
ϕ0) |
|
with π + ϕ0 < ϕ < α. |
(2.54) |
|||||||||
u(kr, ϕ |
− ϕ0) |
= v(kr, ϕ |
− |
ϕ0) |
|
||||||||
|
|
+ |
= |
|
+ |
|
|
|
|
|
|
|
|
In Equations (2.52) and (2.53), |
the |
term |
e−ikr cos(ϕ−ϕ0) determines |
the incident |
|||||||||
plane wave, which exists only |
in |
the illuminated |
region, 0 < ϕ < π + ϕ0, and |
||||||||||
the term e−ikr cos(ϕ+ϕ0) |
relates to the reflected plane wave existing in the region |
0 < ϕ < π − ϕ0. In agreement with the geometrical optics, Equation (2.54) does not contain either the incident or reflected plane waves, because the region π + ϕ0 < ϕ < α is shadowed by the wedge. The boundaries of the incident and reflected plane waves are shown in Figure 2.7.
Figure 2.7 The incident plane wave propagates from the direction ϕ = ϕ0. The line ϕ = π − ϕ0 is the boundary of the reflected wave, and the line ϕ = π + ϕ0 is the shadow boundary.
TEAM LinG
44 Chapter 2 Wedge Diffraction: Exact Solution and Asymptotics
Figure 2.8 The incident plane wave illuminates both faces of the wedge. The line ϕ = π − ϕ0 is the boundary of the wave reflected from the face ϕ = 0, and the line ϕ = 2α − π − ϕ0 is the boundary of the wave reflected from the face ϕ = α.
Figure 2.8 illustrates the situation when both faces of the wedge are illuminated. In this case α − π < ϕ0 < π , and functions u(kr, ϕ − ϕ0), u(kr, ϕ + ϕ0) are determined by
u(kr, ϕ − ϕ0) = v(kr, ϕ − ϕ0) + e−ikr cos(ϕ−ϕ0) u(kr, ϕ + ϕ0) = v(kr, ϕ + ϕ0) + e−ikr cos(ϕ+ϕ0)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.55) |
u(kr, ϕ |
ϕ0) v(kr, ϕ |
|
ϕ0) |
|
e |
− |
ikr cos(ϕ |
− |
ϕ0) |
with π − ϕ0 < ϕ < 2α − π − ϕ0 |
||||
u(kr, ϕ |
− ϕ0)= v(kr, ϕ − |
ϕ0)+ |
|
|
|
|||||||||
|
|
+ |
= |
+ |
|
|
|
|
|
|
|
|
(2.56) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
and |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u(kr, ϕ |
− |
ϕ0) |
v(kr, ϕ |
− |
ϕ0) |
e−ikr cos(ϕ−ϕ0) |
|
with 2α − π − ϕ0 < ϕ < α |
||||||
u(kr, ϕ |
ϕ0) |
= v(kr, ϕ |
ϕ0) + e−ikr cos(2α−ϕ−ϕ0) |
|||||||||||
|
+ |
|
= |
+ |
|
+ |
|
|
|
|
|
|
|
|
(2.57) The term e−ikr cos(2α−ϕ−ϕ0) describes the plane wave reflected from the face ϕ = α.
2.4 THE SOMMERFELD RAY ASYMPTOTICS
The relationships us = Ez and uh = Hz exist between the acoustic and electromagnetic
edge-diffracted rays.
A simple asymptotic expression for the function v(kr, ψ ) with kr 1 can be found by the steepest descent method (Copson, 1965; Murray, 1984). With this purpose we replace the integration variable in Equation (2.51) by
√ |
|
|
iπ/4 |
|
ζ |
|
|
|
|
|
|
||||
s = |
2e |
|
sin |
|
. |
(2.58) |
|
|
2 |
TEAM LinG
2.4 The Sommerfeld Ray Asymptotics 45
Then s2 = i(1 − cos ζ ) and
|
|
sin |
π |
∞ |
|||||
|
|
|
|
|
|
||||
v(kr, ψ ) |
n |
ei(kr+π/4) |
|||||||
|
|
|
|||||||
= nπ |
√ |
|
|
||||||
|
2 |
|
−∞ |
||||||
|
|
|
|
|
|
|
|
|
|
n |
− |
e−krs2 ds |
|
2 |
|
(2.59) |
||
|
|
n |
|
|||||||
|
cos |
π |
|
cos |
ψ + ζ |
|
cos |
ζ |
|
|
|
|
|
|
|
|
where n = α/π .
Here, s = 0 is the saddle point. Indeed, when the point s moves from the saddle point along the imaginary axis, the function exp(−krs2) increases most rapidly. In contrast, this function decreases most rapidly when the point s moves away from the saddle point along the real axis. Because of that, the vicinity of the saddle point provides the major contribution to the integral when kr 1. According to the steepest descent method, the slowly varying factor of the integrand is expanded into the Taylor power series near the saddle point, and then it is integrated term by term. If the integrand expansion is convergent only in the vicinity of the saddle point, the series obtained after integration will be semiconvergent, that is, asymptotic. Retaining the first term in this series for the function v(kr, ψ ), we obtain
|
|
sin |
π |
|
|
ei(kr+π/4) |
|
|
∞ |
|
|
|
|
1 |
|
sin |
π |
|
|
ei(kr+π/4) |
|
||||||||
v(kr, ψ ) |
|
n |
|
|
|
e−krs2 ds |
|
|
n |
n |
|
. |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
nπ |
√ |
|
|
π |
|
ψ |
|
= |
|
|
π |
|
|
ψ |
√ |
|
|
||||||||||||
|
2 |
|
− cos |
−∞ |
|
|
|
− cos |
|
2π kr |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
cos |
|
|
|
|
|
cos |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
n |
n |
|
|
n |
n |
|
|
|
|
|
(2.60)
The next terms of the asymptotic series for function v(kr, ψ ) are small quantities of order (kr√)−3/2 and higher. The asymptotic expression (2.60) is valid under the condition kr|cos(ψ/2)| 1 and describes cylindrical waves diverging from the edge, that is, the edge waves.
According to Equations (2.40), (2.49), and (2.60), the wave diffracted at the edge of the acoustically soft wedge is determined as
d |
= u0[v(kr, ϕ − ϕ0) − v(kr, ϕ + ϕ0)] u0 f (ϕ, ϕ0, α) |
ei(kr+π/4) |
|
|
||
us |
√ |
|
|
, |
(2.61) |
|
2π kr |
||||||
|
|
|
|
|
|
|
where |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin |
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
1 |
|
|
|
|
|||
f (ϕ, ϕ0, α) |
|
n |
|
|
|
|
|
|
|
. (2.62) |
|||||||||||
|
|
|
|
|
n − |
|
n |
|
|
n |
− |
|
n |
|
|
||||||
|
= |
|
|
|
|
|
|
− cos |
|
|
|
|
|
||||||||
|
n |
|
|
|
|
π |
|
|
ϕ − ϕ0 |
|
π |
|
|
ϕ + |
|
|
|
|
|||
|
|
|
|
|
|
cos |
|
|
cos |
|
|
|
|
cos |
|
|
|
|
|
TEAM LinG
46 Chapter 2 Wedge Diffraction: Exact Solution and Asymptotics
Equations (2.41), (2.49), and (2.60) determine the wave arising at the edge of the acoustically hard wedge
d |
= u0[v(kr, ϕ − ϕ0) + v(kr, ϕ + ϕ0)] u0g(ϕ, ϕ0, α) |
ei(kr+π/4) |
|
|
||
uh |
√ |
|
|
, |
(2.63) |
|
2π kr |
||||||
|
|
|
|
|
|
|
where |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. (2.64) |
|
|
|
sin |
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
1 |
|
|
|||
g(ϕ, ϕ0 |
, α) |
|
n |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
ϕ − ϕ0 |
|
|
|
|
|
|
ϕ + ϕ0 |
|
||||
|
|
= |
n |
|
|
|
cos |
π |
− |
cos |
|
+ cos |
π |
− |
cos |
|
|
|||
|
|
|
|
|
|
|
n |
|
n |
|
|
n |
|
n |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Functions f and g describe the directivity patterns of the edge waves. The asymptotic expressions (2.60) to (2.64) were introduced by Sommerfeld (1935) and are well known. It is easy to verify that they satisfy the boundary conditions (2.2) and (2.3).
One should mention that the two-dimensional edge waves (2.61) and (2.63) can be interpreted as continuous sets of edge diffracted rays. They arise due to diffraction, but propagate from the edge in the first asymptotic approximation as ordinary rays in accordance with geometrical optics laws for 2-D fields. As shown in the work of Pelosi et al. (1998), the edge diffracted rays had already been visually observed already by Newton, although he did not use such a terminology. The term “diffracted ray” was introduced by Kalashnikov (1912), who was also the first to present an objective experimental proof of the existence of edge diffracted rays by recording them on a photographic plate. Theoretically, their existence was established first by Rubinowicz (1924) and later on by many other researchers. Keller (1962) formulated the concept of diffracted rays in a general form.
Here it is also pertinent to remind one about Sommerfeld’s warning against the too formal ray interpretation of diffraction phenomena. He wrote that shining diffraction points on edges do not exist in reality and they are just optical illusions: “Das ist naturlich eine optische Tauschung” (Sommerfeld, 1896, p. 369). He explained that such seemingly shining edge points are the result of our perception, or in Sommerfeld’s words, the result of “analytical continuation” of diffracted rays by our eyes.
Because of the ray structure of edge waves, the asymptotic expressions (2.61) to (2.64) derived above can be called the ray asymptotics, as emphasized in the title of the present section. These asymptotics have an essential drawback. They are not valid near the shadow boundary (ϕ ≈ π + ϕ0) and near the boundaries of reflected plane waves ( ϕ ≈ π − ϕ0, ϕ ≈ 2α − π − ϕ0) . The mathematical reason for this drawback
is given in the following. Two poles |
|
|
|
|
|
|
|
|
|
||||||||
|
√ |
|
|
π − ψ |
, |
|
√ |
|
|
|
|
|
π + ψ |
|
(2.65) |
||
s |
2eiπ/4 sin |
s |
2eiπ/4 sin |
α |
− |
||||||||||||
|
2 |
||||||||||||||||
|
1 = |
|
|
2 |
|
2 |
= |
|
|
|
|
||||||
of the |
integrand |
in Equation |
(2.59) |
approach the |
saddle point s = 0 |
when |
ψ = ϕ ± ϕ0 → π and ψ = ϕ + ϕ0 → 2α − π . In this case, the Taylor expansion
TEAM LinG
2.5 The Pauli Asymptotics 47
for the integrand becomes meaningless because its terms tend to infinity. There is a physical background behind this mathematics. The vicinity of boundaries of incident and reflected waves is the region of the effective transverse diffusion where the field cannot be described in terms of diffracted rays and has a more complex structure. This phenomenon is considered in detail in Section 5.5 of Ufimtsev (2003).
2.5 THE PAULI ASYMPTOTICS
In 1938, Pauli suggested the asymptotic expansion for the function v(kr, ψ ) that is valid at the geometrical optics boundaries ϕ = π ± ϕ0 and transforms to the Sommerfeld asymptotics away from these boundaries (Pauli, 1938). In this section we provide the derivation for the first term of the Pauli expansion. Usually, for engineering analysis, only the first terms in asymptotic expansions are of practical value. Higher-order terms commonly are not utilized, because they are smaller in magnitude, and are quite complicated to evaluate. Besides, the high-order terms can occur beyond the frames of validity of idealized mathematical models used for description of real physical phenomena. That is why we focus here on the first asymptotic term.
According to Equation (2.59),
|
|
sin |
π |
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ds |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
v(kr, ψ ) |
= |
|
|
|
ei(kr+π/4) |
|
|
|
|
|
|
e−krs |
|
|
|
|
|
|
|
|
|
|
|
, |
(2.66) |
||||||||||||||
|
|
√ |
|
|
|
|
|
|
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−∞ |
|
− |
|
ψ |
ζ |
|
|
|
|
|
|
ζ |
|
||||||||||||
|
|
nπ |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
n |
|
2 |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cos |
|
|
|
cos |
|
+ |
|
|
|
|
cos |
|
|
|
|
|||||||||
where |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s = |
√ |
|
|
iπ/4 |
|
|
ζ |
|
|
|
|
|
|
|
2 |
= i(1 − cos ζ ), |
|
|
n = |
α |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
2e |
|
|
|
|
sin |
2 |
|
, |
|
|
|
|
s |
|
|
|
π |
. |
(2.67) |
|||||||||||||||||||
Let us multiply and divide the integrand by |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
cos ψ + cos ζ = i(s |
2 |
|
2 |
|
|
|
2 |
|
|
|
|
|
2 ψ |
|
||||||||||||||||||||||||
|
|
− s0) |
|
with s0 |
= 2i cos |
|
|
|
. |
|
|
|
(2.68) |
||||||||||||||||||||||||||
|
|
|
|
2 |
|
|
|
||||||||||||||||||||||||||||||||
Then |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin |
π |
|
|
|
|
∞ |
|
|
|
e−krs |
2 |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
v(kr, ψ ) |
|
|
n |
ei(kr−π/4) |
f (s, ψ ) |
|
ds, |
(2.69) |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
= nπ √ |
|
|
|
s2 − s02 |
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
−∞ |
|
|
|
|
|
|
|
|
|
|
|
where the poles s = ±s0 = ±√ |
2ei π4 |
cos ψ2 |
are outside the integration contour and |
||||||||||||
approach it at the saddle point s = 0 when ψ → π . The function |
|||||||||||||||
f (s, ψ ) |
= |
|
|
|
|
cos ψ + cos ζ |
|
|
|
|
(2.70) |
||||
|
|
|
|
|
|
2 |
|
||||||||
|
|
n − |
|
ψ |
n |
ζ |
|
|
|||||||
|
|
|
cos |
|
|
|
cos |
|
+ |
|
|
cos |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
does not have a pole at the saddle point s = 0 (ζ = 0) when ψ = ϕ ± ϕ0 → π . Therefore it can be expanded into a regular Taylor series. By integrating this series
TEAM LinG