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

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

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

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

Добавлен: 28.06.2024

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

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

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

94 Chapter 5 First-Order Diffraction at Strips and Polygonal Cylinders

under the conditions q(1 ± α) 1 and q(1 ± α0) 1. Here, symbol O(qm) is used to show the asymptotic behavior of quantities Qs,h, α0, q) when q → ∞. The

definition of this symbol has been given above in conjunction with Equation (2.104).

Also,

 

 

 

/

 

 

 

 

 

 

lim '1 − α2

 

 

 

 

Qs, α0) = 0,

with α → ±1,

(5.57)

∂α

and

 

 

 

 

 

 

Qh, ±1) = Qh(±1, α0) = 0.

(5.58)

Equations (5.57) and (5.58) are the consequences of the fact that the scattered field (5.46) and (5.47) and its normal derivatives at the plane x = 0 are continuous outside the strip. In other words, according to approximation (5.46) and (5.47), no scattering sources exist outside the strip surface.

We have calculated functions (5.46) and (5.47) with the well-established estimations (5.55) and (5.56) to demonstrate the accuracy of the PO and PTD approximations. The results are plotted in Figures 5.5 and 5.6 in the decibel scale as 10 logs,h/kl2). It is seen in these figures that the accuracy of the first-order PTD approximation is higher for the soft boundary condition. Only a small discrepancy with the exact TED curve is observed in the vicinity of the direction φ = 90, and the PTD and TED curves infact merge in the region 100< φ ≤ 270. The reason for the better approximation provided by the PTD for the soft boundary condition is the faster attenuation of the primary edge waves of the scattering sources js(1) compared to the similar waves of jh(1). In this case, the amplitude of these waves (at the location of the opposite edge of the strip) is of the order (2ka)−3/2 for the soft strip, and of the order (2ka)−1/2 for the hard strip.

It is also seen in Figures 5.5 and 5.6 that PO cannot provide a reasonable approximation for the scattered field in the vicinity of the minima of the directivity patterns,

Figure 5.5 Scattering at an acoustically soft strip (Ez -polarization of the electromagnetic field).

TEAM LinG


 

 

 

 

 

 

 

 

5.1 Diffraction at a Strip

95

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Figure 5.6 Scattering at an acoustically hard strip (Hz -polarization of the electromagnetic field).

where PO predicts wrong pure zeros. The incorrect values of the PTD approximation for the function σh in the vicinity of the directions φ = 90and φ = 270are caused by the fictitious scattering sources jh(1) distributed outside the strip surface, as has already been discussed at the end of Section 5.1.2. This shortcoming of the first-order PTD is eliminated in the next section.

5.1.4 First-Order PTD with Truncated Scattering Sources jh(1)

The relationship uh = Hz exists between the acoustic and electromagnetic fields for this

problem.

The geometry of the problem is shown in Figure 5.1 and the incident wave is given by Equation (5.1). To calculate the scattered field uh, we apply Equation (1.10) and utilize the following observations.

The symbol r was used in Equation (1.10) for the distance'between the integration and observation points. Now we replace it by ρ = x2 + (y η)2 + ζ 2 and retain the symbol r for the polar coordinate r = 'x2 + y2 of the observation point, assuming that x = r cos φ and y = r sin φ.

The surface of integration in Equation (1.10) covers both faces of the strip: S = S+ S+ (Fig. 5.7).

TEAM LinG


96 Chapter 5 First-Order Diffraction at Strips and Polygonal Cylinders

Figure 5.7

The integration surface S = S+ S+ in Equation (1.10) applied to the strip problem.

 

 

 

The

operator

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

in Equation (1.10)

acts

on

 

coordinates

of

the

 

 

 

 

 

 

 

 

∂n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

integration/source point. According to Equations (1.14) and (3.5), it can be

 

replaced by the operators acting on coordinates of

the

observation point.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Namely,

 

=

 

 

for the illuminated face Sof the strip and

 

 

= −

 

 

 

∂n

∂x

∂n+

∂x

 

for the shadowed face S+, where ∂/∂x is the differentiation with respect to the

 

observation point (outside the strip).

 

 

 

 

 

 

 

 

 

 

 

 

 

The integral in Equation (1.10) over the variable ζ applied to the strip problem

 

can be expressed through the Hankel function H(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

(k

'x

2

+ (y η)

2

) as shown

 

in Equation (3.7).

0

 

 

We calculate the scattered field in the far zone where r

ka2. For this zone

the Hankel function can be approximated by its asymptotic expression

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

H0(1) k(x2 + (y η)2

 

 

eik(rη sin φ)iπ/4.

(5.59)

π kr

We denote the scattering sources as jh(0) = u0Jh(0) and jh(1) = u0Jh(1).

Then the scattered field in the far zone can be represented in the following form:

uh(0)

uh(1)

 

eiπ/4

 

eikr

 

a

= u0

2

 

 

 

 

 

 

Jh(0), x = −0)eikη sin φ dη,

∂x

 

2π k

r

 

 

 

 

 

 

 

 

 

a

 

eiπ/4

 

eikr

 

a

= u0

2

 

 

 

 

 

a Jh(1), x = −0)eikη sin φ dη

∂x

 

2π k

 

r

 

a

 

 

 

 

 

a Jh(1), x = +0)eikη sin φ dη .

(5.60)

(5.61)

TEAM LinG


5.1 Diffraction at a Strip

97

Remember that the uniform scattering sources jh(0) face of the strip where x = −0.

In the far zone approximation one can set

 

eikr

ik cos φ

eikr

 

 

 

 

 

∂x

 

r

r

exist only on the illuminated

. (5.62)

According to the definition (1.31), one has Jh(0), x = −0) = 2eikη sin φ0 .

The nonuniform scattering source jh(1) is found by utilizing the exact solution for the half-plane diffraction problem presented in Section 2.5. According to this solution

The function jh(1) is anti-symmetrical,

jh(1), x = +0) = −jh(1), x = −0) and consequently

Jh(1), x = +0) = −Jh(1), x = −0);

– The function Jh(1), x = −0) is determined as

 

 

 

 

Jh(1), x = −0) = 2

0v )k(a η),

 

π

φ0

* eika sin φ0

 

 

 

 

 

 

 

2

 

 

+ v )k(a + η),

π

+ φ0

* eika sin φ0

1 .

(5.63)

 

 

2

 

In view of above observations, the functions (5.60) and (5.61) can be approximated by

uh(0)

uh(1)

=

u i2 cos φ

sin[ka(sin φ0 − sin φ)]

 

ei(kr+π/4)

,

 

 

 

 

 

 

 

 

 

 

 

0

 

 

sin φ0 − sin φ

2π kr

 

 

 

 

 

ei(kr+π/4) a

 

 

π

 

= u02ik cos φ

 

 

 

a 0v )k(a η),

 

φ0

*

2

2π kr

 

 

 

 

π

 

 

 

 

 

 

 

 

+ v)k(a +

η),

 

 

+ φ0* eika sin φ0

1 eikη sin φ dη.

 

 

2

(5.64)

eika sin φ0

(5.65)

The asymptotic expression (5.64) was found here by the direct integration of the uniform component of the scattering sources. It is totally identical to the expression (5.11) derived by the summation of two edge waves.

The field (5.65) is calculated by integration by parts.

We provide here the final results in terms of the directivity patterns:

ei(kr+π/4)

uh(0,1) = u0 (h0,1), φ0) √ (5.66)

2π kr

TEAM LinG


98

 

Chapter 5

 

First-Order Diffraction at Strips and Polygonal Cylinders

 

 

 

 

 

 

 

 

 

 

 

 

 

 

where h(0) is given by Equation (5.13) and

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

 

 

cos φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ei3π/4

 

 

 

 

 

 

 

 

 

 

2ka(1+sin φ0)

 

2

 

 

 

h

 

, φ0) =

 

 

 

 

 

 

 

 

i4 sin(kaξ ) + 2

 

 

 

eikaξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eit

 

 

dt

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt − 2 √π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− eikaξ

0

 

 

 

 

 

 

 

 

 

 

 

 

eit

 

 

 

1+ sin φ eikaξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ka(1

 

 

sin φ0)

 

2

 

 

 

 

 

 

 

 

 

ei3π/4

 

 

 

 

1

 

sin φ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×

0

 

 

 

 

 

 

 

 

 

 

 

 

eit

dt

 

1sin φ eikaξ

0

 

 

 

 

 

 

 

 

 

 

eit

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ka(1

 

 

 

sin φ)

 

 

2

 

 

 

 

 

 

 

1

 

 

sin φ0

 

 

 

 

 

 

2ka(1

 

 

 

sin φ)

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.67)

 

 

where ξ = sin φ0 − sin φ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

The directivity pattern of the total scattered field uhtr = uh(0) + uh(1)

 

(found by

truncation of the nonuniform sources) equals htr = h(0) + h(1):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ei3π/4

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tr

 

 

cos φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ikaξ

 

 

2ka(1+sin φ0)

 

 

 

 

it2

 

 

 

 

 

h

, φ0) =

 

 

 

 

 

 

i2 sin(kaξ ) + 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

dt

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt − 2 √π

 

 

 

 

 

 

 

 

 

 

− eikaξ

0

 

 

 

 

 

 

 

 

 

 

 

 

eit

 

1+ sin φ eikaξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ka(1

 

 

 

sin φ0)

2

 

 

 

 

 

 

 

 

 

ei3π/4

 

 

 

 

1

 

sin φ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×

0

 

 

 

 

 

 

 

 

 

 

 

 

 

eit

 

dt

 

 

1sin φ eikaξ

 

0

 

 

 

 

 

 

 

 

 

 

eit

 

 

 

 

 

 

 

 

 

+

sin φ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ka(1

 

 

 

 

 

2

 

 

 

 

 

 

1

 

 

sin φ0

 

 

 

 

 

 

2ka(1

 

 

sin φ)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.68)

In the ray region (where ka(1 ± sin φ0)

 

 

 

1 and ka(1 ± sin φ)

1), the asymptotic

expression of this function is given by

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

trh , φ0) = g(1)eika(sin φ0sin φ) + g(2)eika(sin φ0sin φ)

− cos φ

 

 

eika(sin φ0+sin φ)

+

eika(sin φ0+sin φ)

 

(1 + sin φ)

 

(1 − sin φ)

 

1 + sin φ0

1 − sin φ0

×

ei(2ka+π/4)

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

(5.69)

 

2π ka

 

 

 

 

where functions g(1) and g(2) are determined in Equations (5.31) and (5.33). These functions relate to the primary edge waves arising at the ends of the strip. The terms in the second line of Equation (5.69) represent the radiation of the primary edge waves when they reach the opposite end of the strip. These terms can be interpreted

TEAM LinG