Файл: Ufimtsev P. Fundamentals of the physical theory of diffraction (Wiley 2007)(348s) PEo .pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 28.06.2024
Просмотров: 945
Скачиваний: 0
48 Chapter 2 Wedge Diffraction: Exact Solution and Asymptotics
term by term, Pauli obtained the asymptotic expansion of function v(kr, ψ ) argument kr. The first term of this expansion is determined by
|
|
√ |
|
|
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin |
|
|
|
1 + cos ψ |
|
∞ |
e− |
krs2 |
|
|||||||
|
|
|
|
|
|
|
|||||||||||||
|
|
|
n |
|
|
|
|
||||||||||||
v(kr, ψ ) |
= |
2 |
|
|
|
|
ei(kr−π/4) |
|
ds. |
||||||||||
π |
n |
|
|
|
|
|
0 |
s2 − s02 |
|||||||||||
|
|
|
cos π |
− |
cos ψ |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
n |
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Here the integral can be represented as
|
|
|
|
∞ e−krs2 |
|
|
|
|
|
2 |
∞ |
∞ |
2 |
|
2 |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
ds = e−krs0 |
ds |
e−(s |
−s0 )t dt. |
|
|
|
|
|
|||||||||
|
|
|
|
0 s |
2 |
|
2 |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
− s0 |
|
|
|
|
|
0 |
kr |
|
|
|
|
|
|
|
|
|
||||||
By changing the order of integration we obtain |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
∞ e−krs2 |
|
|
|
|
|
∞ |
|
dt |
|
∞ |
|
|
√ |
|
|
∞ |
|
|
|
2t |
|||||||
|
|
|
|
2 |
2 |
2 |
|
π |
2 |
i |
s0 |
|
|||||||||||||||
|
|
|
|
ds = e−krs0 |
es0 t |
√ |
|
|
|
|
e−x |
dx = |
|
|
|
e−krs0 |
e |
| |
|
| |
|
||||||
0 s |
2 |
2 |
|
|
2 |
|
|
|
|||||||||||||||||||
t |
0 |
|
|
|
|||||||||||||||||||||||
|
− s0 |
|
√ |
|
|
kr |
|
|
|
|
|
|
|
|
kr |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
π |
e−ikr|s0| |
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
= |
|
|
|
|
|
√ |
|
|
|
eiq dq. |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
| |
0| |
|
|
|
|
|
kr |
|s0| |
|
|
|
|
|
|
|
|
|
|
|
for large
(2.71)
(2.72)
dt
√
t
(2.73)
As a result, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
π |
cos |
|
ψ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
2 |
|
|
sin |
|
|
|
|
|
|
|
|
|
|
|
|
|
ikr cos ψ e−i π4 |
∞ |
|
|
|
|
|
|
||||||||||||||||||
|
|
n |
2 |
|
|
|
|
|
|
iq2 |
|
||||||||||||||||||||||||||||||||
v(kr, ψ ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e |
dq |
(2.74) |
|
= n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ψ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
√π |
|
|
|
|
|
ψ |
|
|||||||||||||||||||
|
|
|
|
π |
− |
|
|
|
|
|
|
|
|
cos |
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
√2kr |
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
cos |
n |
cos |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
or |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
π |
|
ψ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
2 |
|
|
sin |
|
|
cos |
|
|
|
|
|
|
e−ikr cos ψ |
e−i π4 |
|
∞ cos ψ2 |
|
eiq2 dq. |
|
||||||||||||||||||||||
v(kr, ψ ) |
|
n |
2 |
|
|
|
|
(2.75) |
|||||||||||||||||||||||||||||||||||
= n cos |
π |
|
|
|
ψ |
|
√ |
|
· |
√ |
|
cos ψ2 |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
π |
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
− cos |
|
|
|
|
|
2kr |
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
n |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
This expression represents the slightly modified first term in the Pauli asymptotic expansion. The next term is of order (kr)−1/2 near the boundaries ϕ = π ± ϕ0 and it is of order (kr)−3/2 away from them.
The upper limit of the Fresnel integral in Equation (2.75) should be read as sgn(cos ψ/2)∞. It always equals infinity but changes its sign when the observation point intersects the geometrical optics boundaries (ψ = ϕ ± ϕ0 = π ). Here, the function v(kr, ψ ) undergoes the discontinuity and in this way it ensures the continuity of the function u(kr, ψ ) and therefore the continuity of the total field. Indeed, by using the formula
∞ |
|
√ |
|
i π |
|
iq2 |
π |
|
|||
e |
dq = |
|
|
e 4 |
(2.76) |
|
|
02
TEAM LinG
|
|
|
|
|
|
|
|
|
|
|
|
|
2.5 The Pauli Asymptotics 49 |
|||
one can show that |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
eikr , |
|
|
|
|
|
1 |
eikr |
|
|||
|
v(kr, π + 0) = |
|
|
v(kr, π − 0) = − |
|
(2.77) |
||||||||||
|
|
2 |
2 |
|||||||||||||
and |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
eikr. |
|
|
|
|
|||
|
|
|
u(kr, π ± 0) = |
|
|
|
|
|
(2.78) |
|||||||
|
|
2 |
|
|
|
|||||||||||
Also, with the help of the asymptotic approximations, |
|
|||||||||||||||
p |
iq2 |
eip2 |
−p |
iq2 |
|
|
eip2 |
|
|
|
|
|||||
e |
dq |
|
, |
|
|
e |
dq |
− |
|
, |
with p 1, |
(2.79) |
||||
2ip |
|
|
2ip |
|||||||||||||
∞ |
|
|
|
−∞ |
|
|
|
|
|
|
|
|
|
|
it is easy to verify that the Pauli expression√ (2.75) converts to the Sommerfeld asymptotics (2.60) under the condition kr|cos ψ/2| 1.
As shown in Section 5.5 of Ufimtsev (2003), the Pauli asymptotics (2.75) can be considered as a “stenographic form” of the more physically meaningful expression
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
π |
|
|
|
|
|
|
|
|
|||
|
= |
|
2 |
|
|
− |
|
+ |
|
|
|
|
|
|
ψ |
− ψ π |
√2π kr |
|||||||||||
|
|
|
|
|
|
π |
|
|
||||||||||||||||||||
v(kr, ψ ) |
|
V |
|
|
kr |
(ψ |
|
π ) |
eikr |
|
|
n sin |
n |
1 |
ei kr+ π4 |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
cos |
|
|
− cos |
|
|
|
− |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
n |
n |
|
|
(2.80) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
with
V (τ ) = e−iτ |
2 e−i π4 |
∞ sgn τ 2 |
|
|||
|
√ |
|
|
eiq dq, |
(2.81) |
|
|
|
π |
τ |
|
which follows from the solution of the parabolic equation. Here, the first term V [√kr/2(ψ − π )]eikr describes the transverse diffusion of the wave field in the
vicinity of the geometrical optics boundaries and does not depend on the reflective properties of the wedge faces. The second term in Equation (2.80) can be interpreted as the diffraction background.
It is of interest that in the particular case when the angle α = 2π and the wedge transforms into the half-plane, the Pauli asymptotics (2.75) transforms to the function
|
e−iπ/4 |
√ |
|
cos ψ2 |
|
|
|||
v(kr, ψ ) = e−ikr cos ψ |
2kr |
2 |
|
||||||
√ |
|
|
|
|
|
ψ |
eiq dq |
(2.82) |
|
π |
∞ cos |
||||||||
|
|
|
|
2 |
|
|
and provides the exact (!) solution to the half-plane diffraction problem. Indeed, in this case, n = 2 and Equation (2.51) becomes
v(kr, ψ ) = − |
i |
|
eikr cos ζ |
|
|||
|
D0 |
|
|
|
dζ , |
(2.83) |
|
4π |
cos |
ψ + ζ |
|
||||
|
|
|
|
|
|||
|
|
2 |
|
|
|
||
|
|
|
|
|
|
|
TEAM LinG |
50 Chapter 2 Wedge Diffraction: Exact Solution and Asymptotics
which can be converted to the Fresnel integral. To do this, let us separate the contour D0 (see Fig. 2.6) into two parts at the point ζ = 0. Summation of the integrals over these parts of the integration contour leads to the expression
v(kr, ψ ) |
|
i |
π2 −i∞ eikr cos ζ |
|
1 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
dζ |
|||||||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|||||||||||
|
= − |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ cos |
|
|
|
|||||||||||
|
4π |
0 |
|
|
|
|
|
|
|
ψ + |
ζ |
|
ψ − ζ |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
cos |
|
|
|
|
|
dζ . |
|
|||||||||||||
|
|
|
|
|
|
|
π2 −i∞ eikr cos ζ |
|
cos |
|
ζ |
|
|
|||||||||||||||||
|
|
i |
|
ψ |
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
cos |
|
2 |
|
|
(2.84) |
|||||||||||||||||||||||
|
= − π |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
0 |
|
|
|
|
cos ψ + cos ζ |
|
|||||||||||||||||||||||
We then introduce the integration variable s = √ |
2ei π4 |
sin ζ2 and apply the procedure |
||||||||||||||||||||||||||||
outlined in Equations (2.68) to (2.71). As a result we obtain |
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
√ |
|
|
|
|
|
∞ e−krs2 |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
2 |
|
ψ |
π |
|
|
|
|
|
||||||||||||||||
|
v(kr, ψ ) = − |
|
cos |
|
ei(kr− 4 ) |
|
|
|
|
|
|
|
|
ds. |
(2.85) |
|||||||||||||||
|
π |
2 |
|
|
|
2 |
|
|
|
2 |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
s |
|
− s0 |
|
|
|
|
|
With the help of Equation (2.73) this expression transforms to Equation (2.82). The latter, together with Equations (2.39) to (2.41) and (2.52) to (2.54) provides the exact solution to the half-plane diffraction problem.
Thus, the Pauli asymptotics (2.75) possesses valuable properties. It is simple. It provides the exact solution to the half-plane diffraction problem. It describes both the transverse diffusion of the wave field near the geometrical optics boundaries and the diffracted rays away from these boundaries. However, it is not free from certain drawbacks. These drawbacks are as follows:
•The total field us,h determined with the Pauli asymptotics (2.75) exactly satisfies the boundary conditions (2.2) and (2.3) on the face ϕ = 0. However,
on face√ ϕ = α, these boundary conditions are satisfied only asymptotically, |cos ψ/2| 1 and the Pauli asymptotics converts to the Sommerfeld
expression (2.60).
• The Pauli asymptotics (2.75) provides correct values for the wave field in the direction of the shadow boundary (ϕ = π + ϕ0) and in the direction ϕ = π − ϕ0 of the plane wave reflected from the face ϕ = 0. However, it fails at the direction ϕ = 2α − π − ϕ0 of the plane wave reflected from the face ϕ = α (Fig. 2.8). It predicts a wrong infinite value for the field in this direction.
One can suggest the following remedy to diminish these drawbacks, to some extent. The asymptotics (2.75) should be used only in the region 0 ≤ ϕ ≤ α/2. In order to calculate the field in the rest of the region α/2 ≤ ϕ ≤ α, it is necessary to introduce new polar coordinates with the angle ϕ measured from the face ϕ = α and then to apply the expression (2.75) in the region 0 < ϕ ≤ α/2. In this way, one can obtain correct values for the field at the boundary of the plane wave reflected from
TEAM LinG
2.6 Uniform Asymptotics: Extension of the Pauli Technique 51
the face ϕ = α and satisfy the boundary conditions on this face, but at the expense of the field discontinuity in the direction ϕ = α/2.
The mentioned discontinuity of the field at ϕ = α/2 is manifestation of the fact that asymptotics (2.75) does not satisfy the fundamental physical principle. It is not invariant with respect to choice of the coordinate system. Indeed, if we choose the polar coordinates ϕ and ϕ0 measured from the face ϕ = α, the Pauli asymptotics leads to the relationships
vPauli(kr, ϕ |
|
− |
ϕ |
|
|
) |
|
= |
vPauli |
kr, α |
− |
ϕ |
− |
(α |
− |
ϕ |
) |
] = |
vPauli(kr, ϕ |
− |
ϕ |
) |
(2.86) |
||||||||
|
|
|
0 |
|
|
|
[ |
|
|
|
|
|
|
|
0 |
|
|
|
0 |
|
|
||||||||||
and |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
vPauli(kr, ϕ |
+ |
ϕ |
|
|
) |
= |
vPauli(kr, 2α |
− |
ϕ |
− |
ϕ |
|
) |
|
|
vPauli(kr, ϕ |
+ |
ϕ |
). |
|
(2.87) |
||||||||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
0 |
|
= |
|
0 |
|
|
|
The last inequality indicates that, strictly speaking, the Pauli asymptotics does not satisfy the invariance principle. However, it satisfies this principle approximately, when it transforms into the Sommerfeld ray asymptotics.
In the next section, we derive new asymptotics applicable in all scattering directions (0 < ϕ < α).
2.6 UNIFORM ASYMPTOTICS: EXTENSION OF THE PAULI TECHNIQUE
Here we derive asymptotic expressions under the condition that the incident wave does not undergo double and higher-order multiple reflections at faces of the wedge. This condition is always realized for convex wedges (π < α ≤ 2π ) and also for the concave wedges/horns (π/2 < α < π ), but only for certain directions of the incident wave. However, the theory developed below can be easily extended for any narrow horns (0 < α < π/2) with multiple reflections.
Now we return to Equation (2.66) and we observe that only two poles, |
|
|||||||||||
|
|
|
ψ |
|
|
|
|
ψ |
|
|
||
s1 = √2eiπ/4 cos |
and |
s2 = −√2eiπ/4 cos α − |
, |
(2.88) |
||||||||
|
|
|
||||||||||
2 |
2 |
can approach the saddle point s = 0 when ψ = ϕ ± ϕ0 → π or ψ = ϕ + ϕ0 → 2α − π . The pole s1 approaches the saddle point when the direction of observation ϕ tends to the shadow boundary ϕ = π + ϕ0 or to the boundary ϕ = π − ϕ0 of the wave reflected from the face ϕ = 0 (Fig. 2.7). The pole s2 approaches the saddle point when the direction ϕ tends to the boundary ϕ = 2α − π − ϕ0 of the wave reflected from the face ϕ = α (Fig. 2.8). All other poles in (2.66) can be ignored as they are aside the integration contour and never reach the saddle point in the absence of multiple reflections.
Taking these observations into account we multiply and divide the integrand in
Equation (2.66) by the factor |
|
(cos ζ + cos ψ )[cos ζ + cos(2α − ψ )] = −(s2 − s12)(s2 − s22) |
(2.89) |
TEAM LinG