Документ взят из кэша поисковой машины. Адрес оригинального документа : http://acoustics.phys.msu.ru/teachers/shanin_files/~shanin/papers/sepvar.pdf Дата изменения: Mon Mar 14 17:05:48 2011 Дата индексирования: Mon Oct 1 20:32:44 2012 Кодировка:

A generalization of the separation of variables method for some 2D diffraction problems
A.V. Shanin
Department of Physics, Acoustics division, Moscow State University, Leninskie Gory, Moscow 119992, Russia

Abstract A new method is proposed for solving diffraction problems having piecewise linear ideal boundaries. According to this method, the initial Helmholtz equation is replaced by a system of matrix equations of order 1. The coefficients of these equations are rational matrix functions of the coordinates. The properties of the coordinate matrix system are close to that of the ordinary differential equation, therefore the new method can be treated as a generalization of the separation of variables.

1

Intro duction

Separation of variables (SV) is the most powerful technique for finding exact solutions of the problems related to partial differential equations. However, this method has some weak points. Two of them are the following. First, this method is applicable to a very restricted set of geometries. All these geometries, say for the Helmholtz equation, are classified by group theory analysis; the detailed discussion of this matter can be found, for example, in [1]. Second, even a simple diffraction problem leads to a complicated Fourier series in special functions. As a rule, the structure of the field cannot be easily found from the series. A special technique, such as Watson's transform, is required to study the exact solution obtained using SV. Sadly, such a transform is not known for many interesting cases. In the present paper we propose a technique that can help to overcome the disadvantages mentioned above, i.e. it is applicable to a wider class of problems and the solution is not represented as a Fourier series. The main idea of the method is to reduce the diffraction problem, which involves the Helmholtz equation, boundary, edge and radiation conditions, to an ordinary differential equation (ODE) of a 1


rather simple form. The coefficients of the equation are rational functions of the coordinates. Here we describe the main idea of the SV generalization. Consider the simplest case of SV: the wave field in a rectangular area. Consider the Helmholtz equation
2 u + k0 u = 0

(1.1)

in a rectangle 0 < x < L, 0 < y < M with Dirichlet conditions at the boundaries, u = 0. The solution is the product m n 2 u = sin k x sin y , k= , = k0 - k 2 = , m, n Z . (1.2) L M Construct the vector U = (u1 , u2 , u3 , u4 )T of dimension 4 with the following components: u u 2u u1 = u, u2 = , u3 = , u4 = . x y x y Note that all components of the vector are solutions of the Helmholtz equation in the rectangle, but the boundary conditions for these solutions should be chosen differently. One can check, that the vector U obeys the differential equations U = XU, x with the matrix coefficients 01 -k 2 0 X= 0 0 00 0 0 0 -k 0 0 , 1 0 U = YU, y 0 0 Y= - 0 0 0 0 - 1 0 0 0 0 1 . 0 0 (1.3)

2

(1.4)

2

2

The situation is similar for each classical SV. Let and be the coordinates, in which the variables are separated. Then an eigenfunction can be written as a product u = T ( )S ( ), and the functions T and S obey differential equations of order 2: T ( ) - f1 ( )T ( ) - g1 ( )T ( ) = 0, S ( ) - f2 ( )S ( ) - g2 ( )S ( ) = 0, One can see that the indices here): u u, u, u, following equations are valid (we denote the derivatives by the 0 1 g1 ( ) f1 ( ) = 0 0 0 0 g 2 0 0 0 ( ) f 0 0 1 1 ( u u u u )
, , ,

,

(1.5)



1


u u u u




, , ,

0 0 1 0 0 0 = g2 ( ) 0 f2 ( ) 0 g2 ( ) 0 f



0 1 0 2 (

u u u ) u




, , ,

.

(1.6)

The equations (1.3) will now be called the coordinate equations. These equations formally are in partial derivatives, but their properties are similar to that of an ordinary differential equation. Namely, the solution of a partial differential equation can be found when its values are known on a contour in the (x, y )-plane. Conversely, the solution of the system (1.3) is specified by the value of the vector U at just one arbitrary point P . Any other point P can be connected with P by a line, and an ordinary differential equation describing the behaviour of U on this line can be constructed as a restriction of (1.3) on this line. The value U(P ) should be used as the initial conditions for this equation. One can see that the coefficients of the equations (1.5) and (1.6) have a very strict structure. Some coefficients should be equal to zero, and the other ones should depend on a single variable. The idea of generalization these equations is by allowing the unknown vector U and the coefficient matrices to have a more flexible structure. We demand only that the first component of the vector is the solution of the given diffraction problem, all other components can be chosen arbitrarily. The vector can have an arbitrary dimension, and the elements of the coefficient matrices can be arbitrary functions of the coordinates. Obviously, the structure of the coefficents should be much simpler than the structure of the solutions, otherwise the equations are senseless. Surprisingly, equations of the form (1.3) with relatively simple coefficients can be written for a wide class of diffraction problems. Of course, some price is paid for this generalization of the SV method. While the order of the ODEs in the classical SV applied to the Helmholtz equation is normally equal to 2, the order of the ODEs emerging in our method is usually higher than 2. The more complicated scatterer leads to equation of higher order. The result is a situation, in which the analytical solution of the ODE is not known. However, the equations obtained here belong to a well-studied class: they are confluent Fuchsian. A powerful theory has been developed for these equations, so the situation is not completely hopeless. Our paper has a strong relation to the works [2], [3], [4] devoted to the strip problem. The ODEs obtained in these works for the field on the scatterer or for the far-field diagram can be found as particular cases of the coordinate equations obtained here. The procedure proposed below is close to that of [3]. We should mention here also the work [5], where the structure of the eigenfunction series for the several strips problem was constructed. However, the relation of these series to 3


our method still remains questionable. The paper is organized as follows. In Section 2 the coordinate equations are derived for the problem of diffraction of a plane wave by two Dirichlet or Neumann strips. The coefficients and the initial conditions are known up to several numerical constants, which should be found either by solving the eigenvalue problem or using some other technique. In Section 3 some special issues related to the coordinate equations are discussed. The solvability of the equations and obeying the Helmholtz equation and the boundary conditions by the solutions are these issues. The equations for the directivities are also derived in this section. In Section 4 we represent the analytical solution of the coordinate equations for the case of diffraction by a half-plane and the numerical solution of the coordinate equation for the case of two strips.

2
2.1

Derivation of the co ordinate equations for the problem of diffraction by two strips
Problem formulation

Let the Helmholtz equation (1.1) be fulfilled in the (x, y )-plane. The time dependence is assumed to have the form e-it and is omitted everywhere. The incident plane wave 2 2 -ik x-i k0 -k y , k = k0 cos uin = e ( is the angle of incidence) illuminates the scatterer which is a set of two strips (see Fig. 1) with Dirichlet boundary conditions ud = 0. The problem is effectively two-dimensional. The (x, y )-cross-section of the strips is the set of the segments (a1 , a2 ), (a3 , a4 ) of the x-axis. All the techniques presented below can be obviously generalized to the case of an arbitrary number of the strips located in one plane and bearing the Dirichlet boundary conditions, The case of 2 strips is chosen just because it is the simplest case, for which the classical SV is not known. Note that our method gives nontrivial results even for the cases of one strip and a half-plane, for which the classical SV is applicable. The total field ud can be split into the incident and the scattered field: ud = uin + ud . sc The scattered field is symmetrical. The boundary conditions for ud are as follows. sc ud (x, ±0) = -e- sc
ik x

for 4

x (a1 , a2 ) (a3 , a4 ),

(2.1)


y ud (x, 0) = 0 sc

for

x (-, a1 ) (a2 , a3 ) (a4 , ).

(2.2)

We assume that the radiation condition is fulfilled for the scattered field. Moreover, Meixner's conditions are fulfilled in the vicinity of the edges of the scatterer, i.e. near the points (am , 0). In our case Meixner's conditions conditions demand that the field has the asymptotics u A + B (r
1/2

) + O(r)

near the edges. We formulate the problem of the diffraction by a set of Neumann strips as well. Let the geometry of the strips be the same as for the Dirichlet problem. Let the Neumann boundary conditions be fulfilled by the field un on the strips: y un = 0 for y = ±0, x (a1 , a2 ) (a3 , a4 ).

The total field, again, can be represented as the sum un = uin + un , sc un is the scattered field, which is antisymmetric and satisfies the following boundary sc conditions: y un (x, ±0) = i sc un (x, 0) = 0 sc
2 2 k0 - k e- ik x

for

x (a1 , a2 ) (a3 , a4 ),

(2.3) (2.4)

for

x (-, a1 ) (a2 , a3 ) (a4 , ).

The radiation and the edge conditions are fulfilled by un . sc

2.2

Auxiliary solutions for the strip and the slit problem

The derivation of the coordinate equations for ud,n is performed in two steps. Firstly, sc the auxiliary (oversingular) functions are introduced and the coordinate equations for them are derived. Secondly, the coordinate equations for the ud,n are derived sc using the equations for the auxiliary functions. Introduce the local cylindrical coordinates (m , m ) near the edges of the scatterer (see Fig. 2). Consider the diffraction problem with the same Dirichlet strips as described above, but with the excitation by a line source rather than by a plane wave. Let the source be located near one of the edges (am , 0). Namely, let v m , m = 1 . . . 4 be the limit of the solution of the inhomogeneous Helmholtz equation
2 v m + k0 v m = 1/2 -3/2

(m - ) (m - ),

(2.5)

5


Here is the Dirac's delta-function; 0. The strength of the source is equal to -1/2 Note that the dependence of strength of the source on is chosen such that there exists the non-zero limit of the solution. The boundary condition for v m have the form v m (x, ±0) = 0 for x (a1 , a2 ) (a3 , a4 ), (2.6) Obviously, the function v is valid:
m

is symmetrical, hence the following boundary condition for x (-, a1 ) (a2 , a3 ) (a4 , ). (2.7)

y v m (x, 0) = 0

Also, the radiation condition is fulfilled by v m . Meixner's conditions are assumed to be satisfied for each non-zero , but the behaviour of the limiting function v m near the edge (am , 0) appears to violate Meixner's condition just because the presence of the source near this edge. This behaviour can be studied by the standard technique. Taking the limit, we obtain that the function -1/2 behaves like v m m . Therefore, instead of taking the limit of the line-source problem, one can find v m as the solution having the oversingular behaviour near the edge am . In more details, the asymptotics of the solutions v m near the edges (an , 0) have the form m m,n - n 2Cn 1/2 n v m (n , n ) = - n 1/2 sin + n sin + O(3/2 ), (2.8) n 2 2
m where Cn are some unknown coefficients; m,n is the Kronecker's delta. We also introduce the auxiliary oversingular functions wm , m = 1, . . . 4 for the Newmann strips problem. The boundary conditions have the form

y wm (x, ±0) = 0 wm (x, 0) = 0 for

for

x (a1 , a2 ) (a3 , a4 ),

(2.9) (2.10)

x (-, a1 ) (a2 , a3 ) (a4 , ).

The radiation condition is fulfilled. The Meixner's conditions are satisfied for each non-zero . The functions wm are characterised by the oversingular terms in their asymptotics: m m,n n 2En 1/2 n + n cos + O(3/2 ), (2.11) wm (n , n ) = - -1/2 cos n n 2 2 are some constants. m m The constants Cn and En play an important role in the subsequent derivation of the coordinate equations. Namely, they are the numerical parameters of the coefficients X and Y. At this stage these parameters are unknown, but we assume that they can be calculated by some technique. E 6
m n


2.3

An imp ortant note

We suppose that k0 belongs neither to the spectrum of the Dirichlet nor of the Neumann strips problem. It means that the functions ud,n are defined uniqiely and sc the following statement is true. If some function u(x, y ) obeys the homogeneous Helmholtz equation, radiation condition at infinity, homogeneous boundary conditions (2.4), (2.3) or (2.9), (2.10) on the x-axis and Meixner's conditions at the edges, then u 0. This feature is very important for derivation of the coordinate equations. If k0 belongs to the spectrum, the method should be modified. Note that the auxiliary solutions v m and wm are allowed to be non-zero only because they violate Meixner's condition at one of the edges.

2.4

Derivation of the co ordinate equations for the auxiliary functions
U = (v 1 , v 2 , v 3 , v 4 , w 1 , w 2 , w 3 , w 4 )T . (2.12)

Construct the vector of the unknowns

Our purpose here is to construct the equations of the form (1.3) for this vector. We find it convenient also to introduce the vectors of dimension 4: V = (v 1 , v 2 , v 3 , v 4 )T , Construct the derivative v v
m m

W = (w1 , w2 , w3 , w4 )T .

= (x - am )

-y vm, y x

where m is treated as a global cylindrical coordinate. Obviously, v obeys the homogeneous Helmholtz equation (1.1) and the radiation condition. Moreover, detailed study shows that v satisfies the boundary conditions (2.9) and (2.10) (note that v m 2 2 2 is antisymmetrical and smooth enough on the strips, and that y v m = -(x + k0 )v m ). Meixner's conditions are violated by v . Using the asymptotics (2.8) one can obtain asymptotic estimations for v near the edge an : vm = (-1)m m,n + (an - am )C 2
m n

1 - n

1/2

cos

n + O( 2

1/2 n

).

(2.13)

Using these asymptotics, one can construct the combination v


= v +

(-1)m m w+ 2 7

4 m (an - am )Cn wn , n=1


which satisfies the Meixners's conditions. Besides, v satisfies the Helmholtz equation, radiation condition and the boundary conditions (2.9), (2.10). Due to Subsection 2.3, v 0. Therefore, (x - am ) -y v m = Rm W y x where Rm is the string with the elements
m Rn = 4 m Rn w n , n=1

m = 1 . . . 4,

(2.14)

(-1)m-1 2

m,n

m - (an - am )Cn ,

n = 1 . . . 4.

(2.15)

Thus, one equation expressing the spatial derivatives of v m through the components of the vector U is constructed. We need more such equations to obtain the separate expressions for the derivatives. For this we should find several combinations containing the derivatives of v m and wm with respect to x and y and having -1/2 the singularities at (an , 0) not stronger than n . One can check directly that the combinations (x - am ) -y wm , y x wm vm + (-1)m , y x vm wm - (-1)m y x

satisfy this condition. Using the technique described above, i.e. studying the asymptotics of these combinations, subtracting the appropriate sums of v n or wn and taking into account the uniqueness of the solution, we obtain the following representations: (x - am ) -y w m = Qm V y x
4

Qm v n , n
n=1 4 m Sn w n , n=1 4 m Tn v n , n=1

m = 1 . . . 4,

(2.16)

vm wm + (-1)m = Sm W y x wm vm - (-1)m = Tm V y x
m m m

m = 1 . . . 4,

(2.17)

m = 1 . . . 4,

(2.18)

where Q , S and T are the raws of dimension 4: Q
m n

m Sn m Tn

(-1)m m,n m - (an - am )En , 2 m m = -Cn - (-1)m+n En , m m = -En - (-1)m+n Cn = 8

(2.19) (2.20) (2.21)


for n = 1 . . . 4. The equations (2.14), (2.16), (2.17), (2.18) can be solved with repect to the xand y -derivatives of v m and wm for each m. The result can be written in the form of the coordinate equations U = XU, x where X and Y are the block matrices X= X1 X1 1 2 X2 X2 1 2 , Y=
1 1 Y1 Y2 2 2 Y1 Y2

U = YU, y

(2.22)

,

(2.23)

i The elements of the matrices Xi and Yj are given by the following formulae: j 1m 1 )n m m m (-1)m-1 (x - am )(Qm - (x - am )Tn ) y (Rn - (x - am )Sn ) n 1m , (X2 )n = - , (x - am )2 + y 2 (x - am )2 + y 2 m m m y (Qm - (x - am )Tn ) (-1)m (x - am )(Rn - (x - am )Sn ) n 2m - , (X2 )n = , (x - am )2 + y 2 (x - am )2 + y 2 m m m (x - am )Rn + y 2 Sn (-1)m-1 y (Qm - (x - am )Tn ) n , (Y21 )m = , n (x - am )2 + y 2 (x - am )2 + y 2 m m m (x - am )Qm + y 2 Tn (-1)m y (Rn - (x - am )Sn ) n , (Y22 )m = . n (x - am )2 + y 2 (x - am )2 + y 2

(X

=

2 (X1 )m = n

(Y11 )m = n (Y12 )m = n

Thus, the coordinate equations for the auxiliary solutions are derived. Note the following features of the these equations. - The coefficients of the matrices are rational functions of x and y . The denominators have the form (x - am )2 + y 2 . m m - The coefficients contain several unknown constants Cn , En , not depending on x and y . These constants should be found using a special eigenvalue problem or some other technique. - The auxiliary functions can be found from the system (2.22) without any Fourier series. - The classical separation of variables cannot be performed for two strips.

9


2.5

Derivation of the co ordinate equations for the plane wave incidence
2Cn n = An + 1/2 sin + O(n ), 2 n 2En = Bn + 1/2 cos + O(n ), 2

Let the functions ud and un have the following asymptotics near the edges an : u u
d

(2.24) (2.25)

n

where An , Bn , Cn , En are some constants. Consider the following combination: (ud ) = sc ud sc + i k u x
d sc

This combination obeys the Helmholtz equation, and it satisfies homogeneous boundary conditions (2.6), (2.7). The function (ud ) satisfies the radiation condition, but sc it does not satisfy the Meixner's edge conditions. However, this can be fixed by calculating the asymptotics of (ud ) at the edges and subtracting an appropriate sc combination of the functions v n . After performing this procedure and applying the statement of Subsection 2.3, one can prove that ud sc = -ik ud - sc x
4

(-1)n Cn v n .
n=1

(2.26)

Analogously, one can obtain three other equations: un sc = -ik un + sc x
4

(-1)n En wn ,
n=1 4

(2.27)

ud sc 2 2 = -i k0 - k un - sc y un sc 2 2 = - i k 0 - k ud - sc y

Cn wn ,
n=1 4

(2.28)

En v n .
n=1

(2.29)

These equations together with (2.14), (2.16), (2.14), (2.16) can be rewritten in the matrix form. Namely, introduce the vector ~ U = (ud , un , v 1 , v 2 , v 3 , v 4 , w1 , w2 , w3 , w4 )T . sc sc 10


The coordinate equations for this vector have the form ~ U ~~ = XU, x ~ ~ where X and Y are the block matrices: ~ X= H1 = -ik 0 G1 =


~ U ~~ = YU, y

(2.30)

H1 G1 0X ,

, H2 =

~ X=

H2 G2 0Y
2

,
2

(2.31) ,

0 -ik

0 2 -i k0 - k

2 -i k0 - k 0

C1 -C2 C3 -C4 0 0 0 0 0 0 0 0 -E1 E2 -E3 E4 0 0 0 0 C1 C2 C3 C4 E1 E2 E3 E4 0 0 0 0 .

,

G2 = -

The matrices X and Y are given by (2.23). Thus, the coordinate equation involving the scattered field usc is constructed. This equation has dimension 10. Note that the scattered field is the solution of a single coordinate equation, and it is not represented by a Fourier series. ~ Note also that the vector U can be redefined as follows: ~ U = (ud , un , v 1 , v 2 , v 3 , v 4 , w1 , w2 , w3 , w4 )T , i.e. the total fields can be considered instead of the scattered ones. The new vector obeys the equations (2.30) with the same coefficients.

3
3.1

Sp ecial issues
Identities for the unknown constants

The coordinate equations derived in the previous section ivolve unknown constants, m m namely Cn , En , Cn and En . Here we show that these constants are connected by several algebraic identities. Obviously, due to the reciprocity theorem, the following identities are valid
n m Cn = Cm ,

E

m n

n = Em

for m = n.

(3.1)

Consider the next terms in the asymptotics (2.8), (2.11), namely the terms having 3/2 Їm Їm the order n . To describe these terms, a set of new constants, like Cn and En , can 11


m m be introduced. These constants can be expressed through the constants Cn and En using the coordinate equations. This can be done by substituting the asymptotics into the coordinate equations and comparing the terms containing different powers Їm Їm of n . Note that there are two different ways to express the values Cn and En . The first way is to use the coordinate equation for U/ x, and the second one is to use the equation for U/ y . Comparing the results for similar values, we obtain the following identities: 4 l [(-1)m Elm + (-1)l Clm ]Cn + l=1 4 l [(-1)m Clm + (-1)l Elm ]En + l=1

al - am l l [(-1)n Clm En - (-1)m Elm Cn ] an - am al - am l l [(-1)n Elm Cn - (-1)m Clm En ] an - am

= 0, (3.2) = 0, (3.3)

both for m, n = 1 . . . 4, m = n, and
4 2 k0

+
l=1

(-1)

n+l

l ((Cln )2 + (Eln )2 ) + 2Cln En = 0

(3.4)

for n = 1 . . . 4. Analogously, by calculating the next terms of the expansions (2.24) and (2.25) one can obtain the following relations for Cn and En :
4

i

2 k0

- k En + i(-1) k Cn +
m=1 4

2

n

m m Cm [En + (-1)m+n Cn ] = 0,

(3.5)

i

2 2 k0 - k Cn - i(-1)n k En + m=1

m m Em [Cn + (-1)m+n En ] = 0,

(3.6)

both for n = 1 . . . 4.

3.2

Compatibility, Helmholtz equation and b oundary conditions

Let us formulate the conditions sufficient for the compatibility of the coordinate equations. It is well-known, that compatibility follows from the identity ~ ~ 2U 2U = , x y y x 12


where the corresponding derivatives are calulated by applying the equations for the x- and y - derivatives in different orders. Using the equations (2.30), we obtain the equation ~ ~ X ~ ~ ~ Y ~~ ~ + YX U = + XY U x y Obviously, the sufficient condition for the compatibility is as follows: ~ ~ ~~ ~~ Y + YX = X + XY. x y (3.7)

~ All components of the vector U are solutions of the Helmholtz equation (1.1) by construction, therefore ~ ~ X ~ 2 Y ~ +X + + Y2 + k x x Thus, the condition
2 0

~ U=0

~ ~ 2 ~ ~ X + X2 + Y + Y2 + k0 = 0 (3.8) x x quarantees that all the components of the solution of the coordinate equations satisfy the Helmholtz equation. One can check that the conditions (3.7) and (3.8) are fulfilled provided the conditions (3.1)­(3.6) are obeyed by the coefficients. ~ Each component of the vector U is the solution of a diffraction problem, i.e. the Helmholtz equation and the boundary conditions on the line y = 0. The xaxis is split into 5 parts, namely 1 . . . 5 by the points a1 . . . a4 . Note that the following important property of the coordinate equations is valid. If a solution of these equations satisfy the boundary conditions at any point of a segment i , then this solution satisfies the boundary conditions at all other points of i automatically. This can be proved by restricting the coordinate equations onto the x-axis. This m m property is valid for an arbitrary set of the values Cn , En , Cn , En . Studying the asymptotics of the solutions of the coordinate equations by the standard methods, we conclude that among their solutions their exist some solutions obeying the radiation and the appropriate edge conditions. However, there is no guarantee that a single solution of the coordinate equations satisfies the boundary, edge and radiation conditions simultaneously. A proper choice of the constants C and E is needed to provide this. Therefore, we have a specific eigenvalue problem for defining the constants. This problem is very sophisticated, so below we use another technique (namely, an approximate one) for the practical calculation of the constants. 13


3.3

The equations for the directivities
x = r cos , y = r sin . (3.9)

Introduce the polar coordinates (r, ) by the formulae:

Consider the coordinate equations (2.22). Each of the components of the vector U can be represented using its directivity: v w
n n

^ = V n ()eik0 r r-1/2 + O(eik0 r r-3/2 ), ^ = W n ()eik0 r r-1/2 + O(eik0 r r-3/2 ).

(3.10) (3.11)

as r . Substitute these representations into (2.22). Consider the r-derivatives of the components, i.e the combination U = r cos + sin x y U.

On the one hand, this combination should have the asymptotics of the form vn ^ = ik0 V n eik0 r r-1/2 + O(eik0 r r-3/2 ), r wn ^ = ik0 W n eik0 r r-1/2 + O(eik0 r r-3/2 ). r On the other hand, the same asymptotics can be calculated using the coordinate equations: vm = r wm = r
4 4

(-1) cos
n=1

m

^ T V n + sin
m n n=1 4

m^ Sn W 4

n

ei

k0 r -1/2

r

+ O(e
-1/2

ik0 r -3/2

r

),
-3/2

(-1)

m-1

cos
n=1

^ S W n + sin
m n n=1

m^ Tn V n () eik0 r r

+ O(eik0 r r

).

Comparing these asymptotics, we obtain the following formulae linking the values ^ ^ V and W : 4 1 cos ^ m ^ ^ W m = (-1)m-1 V+ T mV n, (3.12) sin ik0 sin n=1 n or ^ V
m

cos ^ W = (-1) sin
m

m

1 + ik0 sin

4 m^ Sn W n . n=1

(3.13)

14


Construct the -derivative of the vector U, i.e. the combination U U U =x -y . y x The asymptotics of this derivative can again be found using two different methods: directly from (3.10), (3.11) or by applying the coordinate equations. Comparing the results, we obtain the following equations: ^ dV m ^ = ik0 am sin V d
4 m

+
n=1 4

m^ Rn W n ,

(3.14) (3.15)

^ dW m ^ = ik0 am sin W d

m

+
n=1

^ Qm V n , n

Taking into account the relations (3.12) and (3.13), we obtain the closed systems of order 4: 4 4 ^ ^ dV m dW m m^n Їm ^ (3.16) = Kn V , = Kn W n , d d n=1 n=1 where K
m n

= ik0 am

m,n

sin + (-1)

n- 1

cos m 1 Rn + sin ik0 sin

4 ml Rl Tn , l=1 4 l Qm S n . l

Ї K

m n

= ik0 am m,n

cos m 1 sin + (-1) Qn + sin ik0 sin
n

l=1

The equations (3.16), sub ject to the change of the variables, coincide with the equations obtained in [3] (for a single strip) and in [6]. Consider the coordinate equations (2.30). Introduce the directivities for the scattered field: u u
d sc n sc

^ = U d ()eik0 r r ^ = U n ()eik0 r r

-1/2 -1/2

+ O(eik0 r r + O(eik0 r r

-3/2 -3/2

), ).

These representations yield ud sc ^ = ik0 cos U d ()e x un sc ^ = ik0 cos U n ()e x
ik0 r -1/2

r r

+ O(eik0 r r + O(eik0 r r

-3/2

), ).

(3.17) (3.18)

ik0 r -1/2

-3/2

15


The same derivatives can be found using the coordinate equations. Taking the main terms and comparing the results with (3.17), (3.18), we obtain that ^ U ^ U
d

=- =

1 i(k0 cos + k )

4

^ (-1)n Cn V n ,
n=1 4

(3.19) (3.20)

n

1 i(k0 cos + k )

^ (-1)n En W n .
n=1

These equations correspond to the embedding formula introduced in [3] and developed in [7] and [6]. A detailed discussion of the embedding formulae in various diffraction problems in preparation [8].

4
4.1

Examples
Diffraction by a half-plane

Diffraction by an ideal half-plane is a particular case of diffraction by several strips considered above. The theory developed above is applicable to this case. Since the simple analytical solutions are known for the half-plane, it is possible to check the validity of the coordinte equations. Let the Dirichlet or Neumann screen occupy the half-line y = 0, x > 0, i.e. a1 = 0. Introduce the cylindrical coordinates r, by the formulae (3.9). The Sommerfeld solutions for the half-plane problems look like t1 1 - i -ik0 r cos(-) 1 + i 2 ud,n = e + ei /2 d 2 2
0

1-i - e 2 where is the angle of incidence; t1 = 2 - k0 r cos , 2 t2 = 2
ik0 r cos(+ )

1 + i + 2

t

2

ei
2 /2

d ,

(4.1)

0

+ k0 r cos . 2

The auxiliary solutions have the form eik0 r v 1 = - sin , 2 r 16 eik0 r w1 = - cos . 2 r (4.2)


One can check directly that the vector ~ U = (ud , un , v 1 , w1 )T obeys the coordinate equations (2.30), where H1 = -ik 0


0 -ik

,

H2 = ,

0 2 -i k0 - k G2 = -

2

2 -i k0 - k 0

2

,

G1 = k = k0 cos ,

C1 0 0 -E

1

0 C1 E1 0

,

, E1 = (i - 1) k0 cos . 2 2 m m The matrices X and Y can be defined by using the identities for Cn and En . Choosing the constants, corresponding to the outgoing waves, we obtain that C1 = (i - 1) k0 sin
1 R1 = -Q1 = 1/2, 1 1 1 S1 = T1 = ik0 ,

and X= Y=

1 2+y x 1 2+y x

2

-x/2 - ik0 x2 -y /2 + ik0 xy y /2 + ik0 xy -x/2 + ik0 x2 -y /2 - ik0 xy x/2 + ik0 y 2 -x/2 + ik0 y 2 -y /2 + ik0 xy

, .

2

4.2

Numerical calculations for two strips

For demonstration purposes we represent here the computations for diffraction by m two strips. The most difficult point is the calculation of the constants Cn , En , Cn , m En and the initial conditions for the coordinate equations. The proper way for this is to solve the eigenvalue problem, but here for the simplicity of the calculation we use the iteration procedure close to the one described in [9], [4], [10]. According to this procedure the diffraction process is considered as the series of elementary diffraction acts occuring on the edges of the scatterer. Each elementary act can be described by the Wiener-Hopf method [11], so the scattered field is represented in the form of diffraction series. We take the truncations of these series for the approximations of the fields. The dimensionless parameters were chosen as follows: k0 = 1 + 0.12i, a1 = -20, a2 = -8, 17 a 3 = 8, a4 = 20.


The reference point P (in which the initial conditions were taken) was chosen having the coordinates (0, 1). m m ~ The constants Cn , En , Cn , En and the initial conditions U(P ) were calculated ~ ~ approximately using the diffraction series. The coefficients X and Y were calculated according to the formulae (2.31), (2.23). The coordinate equations have been solved according to the following scheme. ~ First, the equation for U/ y was solved along the y -axis on the segment (0.5, 25). ~ Second, the equation for U/ x was solved on the lines y = n/2, n = 1 . . . 50, ~ parallel to the x-axis. The values of U at the points (0, n/2) found in the first step have been used as the initial values. The 3D plots of the functions - Re[ud ] and Im[un ] are shown in Fig. 3 and sc sc Fig. 4, respectively.

4.3

Concluding remarks

A generalization of the separation of variables has been constructed above. Let us discuss the advantages and disadvantages of the proposed method. The main advantages are the following: · The class of problem which can be treated by new method is wider than the class of problem admitting the classical SV. · Unlike in the classical SV, the solution obtained using our method does not contain a series. The price paid for these benefits is the following: · Generally, the equations obtained using our method have order higher than 2 and, therefore, cannot be solved in elementary functions. The more complicated geometry of the boundary lead to the system of higher order. · The coefficients of the equations emerging in our method contain the complicated combinations of the coordinates, which generally cannot be split into a pair containing a single coordinate each. That is why analytically solving the coordinate system on the whole plane seems questionable. · The coefficients of the coordinate equations involve numerical parameters, which should be calculated somehow. A complicated procedure is needed to calculate these parameters with a good accuracy. It is not easy to predict whether the benefits of the new method will overweight its disadvantages. However, we suspect that the system of the form (1.3) possesses 18


some strong properties that simplify the solution significantly. Besides, the system (1.3) can be treated as a new exact property of the diffracted field, and some specific tasks can be solved using this knowledge. In our examples the elements of the matrices X and Y are rational functions of the coordinates. Therefore, the corresponding scattered wave fields belong to a very narrow class: they are the solutions of 2D confluent Fuchsian systems. This makes these fields a sub ject of a powerful mathematical theory. As an example, the analytic continuation of the wave field into the area of complex x and y can be constructed, and the monodromy group of the field can be calculated. We hope that this new link between applied and fundamental mathematics will be a good donation to the theory of diffraction.

Acknowledgements
The current work was supported by the RFBR grant NN 00-15-96530 and by the program "Universities of Russia". Author is grateful to Prof. V.M. Babich and Dr. V.P. Smyshlyaev for the valuable discussions. Author is also very grateful to Dr. R.V.Craster for help in preparation of this paper.

19


References
[1] Miller, W. Symmetry and separation of variables. Addison-Wesley, 1977. [2] Latta, G.E., "The solution of a class of integral equations", J. Rat. Mech., V.5, N.5, pp. 821­834 (1956). [3] Williams, M.H., "Diffraction by a finite strip", Quart. Journ. of Mech. and Appl. Math., V.35, N.1, pp.103­124 (1982). [4] Shanin A.V. "Three theorems concerning diffraction by a strip or a slit" Quart. Journ. of Mech. and Appl. Math., V.54, N.1, pp.107­137 (2001). [5] Shinbrot M., "The solution of some integral equations of Wiener-Hopf type", Quart. Appl. Math. V.28, N.1, pp. 15­36 (1970). [6] Shanin A.V. "Diffraction of a plane wave by two ideal strips" Submitted to QJMAM. [7] Biggs N. R. T., Porter, D., Stirling, D.S.G. "Wave diffraction through a perforated breakwater", Quart. Journ. of Mech. and Appl. Math., V.53, N.3, pp.375­ 391 (2000). [8] Craster R. V., Shanin A. V., Doubravsky E. M. "Embedding formulae in diffraction theory", in preparation. [9] Jones, D.S., Acoustic and electromagnetic waves, Oxford: Clarendon press, 1986. [10] Shanin A.V. "To the problem of diffraction by a slit. Some properties of the Schwarzchild's series", Zapiski seminarov POMI V.275 (2001) (in Russian). [11] Noble, B., Methods based on the Wiener-Hopf technique, Pergamon Press, 1958.

20


Figure 1: Geometry of the problem

21


Figure 2: Local cylindrical coordinates

22


Figure 3: 3D plot for - Re[ud ] sc

23


Figure 4: 3D plot for Im[un ] sc

24


Fig.1: Fig.2: Fig.3: Fig.4:

Captions for the figures Geometry of the problem Local cylindrical coordinates 3D plot for - Re[ud ] sc 3D plot for Im[un ] sc

25