Файл: Ufimtsev P. Fundamentals of the physical theory of diffraction (Wiley 2007)(348s) PEo .pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 28.06.2024

Просмотров: 941

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

or

us =

uh =

 

 

 

2.2 Transition to the Plane Wave Excitation 39

π

π

 

I0H0(1)(kr0)

 

ei 2 vl Jvl (kr)[cos vl ϕ0) − cos vl + ϕ0)], (2.32)

2

l=0

 

 

 

π

π

 

I0H0(1)(kr0)

 

εl ei 2 vl Jvl (kr)[cos vl ϕ0) + cos vl + ϕ0)].

2

l=0

 

 

 

(2.33)

To clarify these expressions, we consider the solution of the wave equation (2.1) in the free infinite homogeneous medium, without the scattering wedge. For this problem, it is convenient to use the new polar coordinates (ρ, φ) with the origin at the radiating source (r0, ϕ0). It is clear that due to the azimuthal symmetry of the problem, its solution is a function of only one variable, u = u(ρ). In addition, as the wave equation is the differential equation of second order, its solution, in general, is the sum of two fundamental solutions of the related homogeneous equation:

u(ρ) = c1H0(1)(kρ) + c2H0(2)(kρ),

(2.34)

with constants c1 and c2. We retain here only the first term,

 

u(ρ) = c1H0(1)(kρ),

(2.35)

because the second term, with H0(2)(kρ), does not satisfy the Sommerfeld radiation condition (2.12) and represents a nonphysical wave incoming from infinity. The constant c1 is found again with the Green theorem (2.15) applied to the circular region S of a small radius ε and with the center at ρ = 0 (Fig. 2.3).

For small values 1, the function u(kρ) and its normal derivative du/dn = du/dρ (at the boundary of the region S) are described by the asymptotic approximations

i2

du(ρ)

 

i2

 

 

u(ρ) c1

 

ln(kρ),

 

c1

 

,

with 1.

(2.36)

π

dρ

πρ

By substitution of these quantities into the Green theorem and taking the limit with ε → 0, we find c1 = I0/i4 and

1

 

u(ρ) = i4 I0H0(1)(kρ).

(2.37)

Figure 2.3 A circular region (0 ≤ ρ ε, 0 ≤ φ ≤ 2π ) with the radiating source at the center.

TEAM LinG


40 Chapter 2 Wedge Diffraction: Exact Solution and Asymptotics

This solution explains the physical meaning of the factor in front of the series in Equations (2.32) and (2.33). It actually represents the field of the incident wave on the edge of the wedge

u0 =

1

 

4i I0H0(1)(kr0).

(2.38)

When r0 → ∞ and I0 → ∞, this field can be interpreted as the plane wave traveling to the wedge from the direction ϕ = ϕ0:

uinc = u0eikr cosϕ0).

(2.39)

As a result, one can rewrite Equations (2.32) and (2.33) in the classical Sommerfeld form (Sommerfeld, 1935):

us = u0 · [u(kr, ϕ ϕ0) u(kr, ϕ + ϕ0)]

for uinc = u0eikr cosϕ0)

(2.40)

 

 

 

 

 

 

 

 

 

 

uh = u0 · [u(kr, ϕ ϕ0) + u(kr, ϕ + ϕ0)]

 

(2.41)

where

 

 

 

2π

 

 

 

 

 

 

 

u(kr, ψ ) =

α

εl ei(π/2)vl Jvl (kr) cos vl ψ .

(2.42)

 

 

l=0

 

 

Equations (2.40) to (2.42) finally determine the total field generated by the incident plane wave in the presence of the perfectly reflecting wedge.

2.3 CONVERSION OF THE SERIES SOLUTION TO THE SOMMERFELD INTEGRALS

In his work Sommerfeld (1935) presented the solution of the wedge diffraction problem in integral form and then transformed it into infinite series. Here, we will perform a reverse procedure, and convert the infinite series (2.40) and (2.41) into integrals. To do this, we use the following expression of the Bessel function (Sommerfeld, 1935):

 

 

1

III

 

J(kr)

=

ei[kr cos β+vl π/2)] dβ,

(2.43)

 

 

 

 

vl

2π

I

 

where the integration contour is shown in Figure 2.4. Then the function u(kr, ψ ) can be represented as

 

1

III

 

 

 

 

 

 

 

u(kr, ψ ) =

2α

I

eikr cos β 1 +

eivl π +ψ ) + eivl π ψ )

dβ, (2.44)

 

 

 

l=1

l=1

 

TEAM LinG


2.3 Conversion of the Series Solution to the Sommerfeld Integrals 41

with vl = lπ/α. Here, the series are geometrical leads to

 

1

III

eikr cos β

1

 

u(kr, ψ ) =

 

 

 

 

2α

I

1 − ei πα π +ψ )

progressions. Their summation

1

dβ. (2.45)

1 − ei πα π ψ )

With a new variable β

= β π , this becomes

 

 

 

u(kr, ψ ) =

1

III

e

ikr cos β

 

1

1

dβ . (2.46)

 

2α

I

 

1 − ei πα +ψ )

1 − ei πα ψ )

Here, the integration contour is shifted by −π compared to the contour shown in Figure 2.4. According to the difference inside the brackets, the function u(kr, ψ ) can be represented as the sum of two integrals. In the integral related to the first term inside the brackets, we replace β by β. In the integral related to the second term inside the brackets, we change β by −β. As a result, we arrive at the Sommerfeld integral (Sommerfeld, 1935)

u(kr, ψ ) =

1

 

 

eikr cos β

 

dβ.

(2.47)

2α

C 1

− e

i πα

+

ψ )

 

 

 

 

 

 

The integration contour C, consisting of two branches, is shown in Figure 2.5. The integrand of u(kr, ψ ) possesses the first-order poles

βm = 2αm ψ , with m = 0, ±1, ±2, ±3, . . . .

(2.48)

Application of the Cauchy theorem to the integral over the closed contour CD (Fig. 2.5) results in

u(kr, ψ ) = v(kr, ψ ) + eikr cos ψ ,

with −π < ψ < π ,

 

u(kr, ψ ) = v(kr, ψ ), with π < ψ < 2α π ,

(2.49)

u(kr, ψ ) = v(kr, ψ ) + eikr cos(2αψ ),

with 2α π < ψ < 2α + π ,

 

Figure 2.4 Integration contour in Equation (2.43).

TEAM LinG



u(kr, ϕ ϕ0) = v(kr, ϕ ϕ0) + eikr cosϕ0) u(kr, ϕ + ϕ0) = v(kr, ϕ + ϕ0) + eikr cos+ϕ0)

42 Chapter 2 Wedge Diffraction: Exact Solution and Asymptotics

Figure 2.5 Integration contours in Equations (2.47) and (2.50).

where

v(kr, ψ ) =

1

 

 

eikr cos β

 

dβ.

(2.50)

2α

D 1

e

i πα

+

ψ )

 

 

 

 

 

 

The integration contour D consists of two branches (Fig. 2.5). In the integral over the left branch, we replace the variable β by ζ π , and in the integral over the right branch we put β = ζ + π . Then the function v(kr, ψ ) transforms into the integral over the contour D0 (Fig. 2.6)

π

sin

v(kr, ψ ) = i n

2π n

 

 

 

eikr cos ζ

 

 

 

 

 

 

 

 

dζ

(2.51)

D0 cos

π

cos

ζ + ψ

 

n

n

 

 

 

 

 

where n = α/π .

The physical interpretation of Equation (2.49) is the following. The function v(kr, ψ ) describes the diffracted part of the field, and the residues relate to the geometrical optics. This interpretation becomes clear if we consider functions u(kr, ϕ ϕ0) and u(kr, ϕ + ϕ0).

In the case 0 < ϕ0 < α π , when only one face (ϕ = 0) of the wedge is illuminated (Fig. 2.7), these functions are determined by

with 0 < ϕ < π ϕ0

(2.52)

TEAM LinG