Материал: Волновое уравнение (96

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

∆G (M , M

 

 

) = ∆

 

 

 

1

 

 

+ v (M )

 

 

 

 

 

 

 

0

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M0M

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

+ ∆v (M )

= −4πδ(M , M

 

 

).

 

 

= ∆

 

 

 

 

 

0

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M0M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∆G (M , M

 

0

) = −4πδ

(M , M

0

),

M Ω,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.15)

 

 

 

G

 

 

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∑

 

 

 

 

 

 

 

M Ω,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∆v = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

(2.16)

 

v

 

 

 

= −

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∑

 

 

 

 

 

 

 

 

rM0P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При таком выборе функции Грина, когда G (P, M0 ) = G

 

∑ = 0,

 

 

применение формулы (2.14) даст решение

 

 

 

 

 

 

 

 

 

u (M0 )

= −

 

1

w∫∫

u (P)

∂G (P, M0 )

dσP .

(2.17)

4π

 

 

∂n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∑

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив в (2.17) заданное значение искомого решения u (P)

на поверхности ∑,

получим решение задачи Дирихле:

 

 

u (M0 ) = −

 

1

 

w∫∫

F (P)

∂G (P, M0 )

dσP .

(2.18)

4π

 

 

∂n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∑

Здесь функция Грина определена как решение задачи (2.15).

2.8.Функция Грина в задаче Неймана

Вэтом случае получить функцию Грина так, чтобы выполня-

лось условие

∂G

 

= 0, не представляется возможным, так как для

 

 

∂n

 

∑

 

 

 

 

21

гармонической функции v (M ) получим

∂G

 

 

 

∂

 

1

 

 

 

 

 

 

 

 

= −

 

 

 

и

∂n

 

 

 

 

 

 

 

 

 

 

∂n r

 

 

 

 

 

 

 

∑

 

 

M0M

 

 

 

 

 

 

 

w∫∫

∂v dσP ≠ 0, что противоречит свойству гармонических функ-

∂n

 

 

 

 

 

 

 

 

 

 

∑

 

 

 

 

 

 

 

 

 

 

 

ций из разд. 2.4. Поэтому функцию Грина для задачи Неймана следует искать из решения уравнения

∆G (M , M0 ) = −4πδ(M , M0 ), M Ω,

 

∂G

 

 

 

4π

 

 

 

 

= −

,

 

 

 

 

 

∂n

 

 

S

 

 

∑

 

 

 

 

 

где S – площадь поверхности ∑.

Решение задачи Неймана

∆u (M ) = 0,

M (x, y, z) Ω,

 

∂u

 

 

 

 

 

 

 

= g (P),

P ∑,

 

∂n

 

 

 

 

∑

 

 

 

 

 

запишется в виде u (M0 ) = 41π w∫∫∑ g (P)G (P, M0 )dσP .

2.9. Физический смысл функции Грина

Из определения функции Грина следует, что в каждой точке M0 Ω она удовлетворяет в обобщенном смысле уравнению Пуассона

∆G (M , M0 ) = −4πδ(M , M0 ), M Ω,

и обращается в нуль на границе области ∑ в случае задачи Дирихле. Поэтому функцию Грина можно интерпретировать как кулонов потенциал, порождаемый внутри заземленной поверхности ∑ единичным зарядом, находящимся в точке M0 Ω.

22

Приведем примеры решения краевой задачи Дирихле с помощью функции Грина для полупространства и шара. Для других областей функция Грина может быть получена методом отраже-

ний [1, 5].

2.10. Примеры построения функции Грина для различных областей

1. Задача Дирихле для полупространства. Постановка крае-

вой задачи имеет вид

∆u ( x, y, z) = 0,

 

−∞ < x, y < +∞,

 

u ( x, y, 0) = F

(x, y),

 

 

z > 0.

 

 

 

 

 

 

 

 

В этой задаче поверхность ∑ представляет собой плоскость z = 0, которую можно считать замкнутой в бесконечности. Для

нахождения функции Грина введем точку M0* (x0 , y0 , − z0 ), со-

пряженную точке M0 ( x0 , y0 , z0 ), т. е. точки M0* и M0 будут симметричны относительно плоскости z = 0.

В качестве гармонической функции в пространстве z > 0 возь-

мем функцию v (x,

 

 

y, z) = −

 

1

 

 

, а функцию Грина возьмем в

r

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

виде

 

 

 

 

 

 

 

 

 

 

 

M0M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∆G (M , M0 ) =

 

 

1

 

−

1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rM0M

 

rM0*M

где r

=

 

 

 

( x − x

)2 +( y − y

)2

+( z − z

0

)2 ;

 

 

M0M

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

r *

=

 

 

 

( x − x0 )2 +( y − y0 )2

+( z + z0 )2 .

M0M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При этом G

 

∑ = G

 

z=0 = 0 . Найдем

∂G

на поверхности ∑:

 

 

 

 

∂n

 

 

 

∂G

 

∂G

 

 

 

 

 

 

 

 

 

 

 

 

2z0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂n

 

∑ = ∂z

 

z=0 = −

 

.

 

 

 

((x − x0 )2 +( y − y0 )2 + z02 )3 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23

Тогда решение задачи Дирихле для полупространства

u (x , y , z

0

) =

z0

+∞+∞

F (x, y)dx dy

.

 

 

0

0

 

2π−∞∫ −∞∫ ((x − x0 )2 +( y − y0 )2 + z02 )3 2

 

 

 

 

 

 

Интеграл в правой части этой формулы носит название инте-

грала Пуассона для полупространства.

2. Задача Дирихле для шара. Постановка краевой задачи имеет вид

 

 

 

∆u

(M ) = 0,

 

M ( x, y, z) Ω ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

u

 

∑ = F (P),

P( x, y, z) ∑R ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где ΩR – шар радиуса R с центром в начале координат;

∑R – по-

верхность этого шара.

 

 

 

 

 

 

 

 

 

 

Введем точку M0* вне шара,

сопряженную точке M0 внутри

шара,

такую,

что

 

r r*

= R2

(рис. 3),

где r =

x2

+ y2 + z2

;

 

 

 

 

 

0 0

 

 

 

0

0

0

 

0

 

r* = x*2

+ y*2

+ z*2

;

 

 

 

r

= ( x − x

)2 +( y − y

)2 +

( z − z

0

)2

;

0

0

0

0

 

 

 

 

 

M0M

 

0

0

 

 

 

 

r *

= (x − x0* )2 +

(y − y0* )2 +

(z − z0* )2 .

 

 

 

 

 

 

M0M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3

24

Совместим точку M ( x, y, z) с точкой P(x, y, z) ∑R (рис. 4).

Рис. 4

Из условия

r r* = R2

и подобия

треугольников ∆OPM

0

и

 

 

 

 

 

0 0

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∆OPM0*

следует r

*

=

rM0P . В качестве функции Грина рас-

 

смотрим функцию

M

0M

 

r0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∆G (M , M0 ) =

1

 

 

−

R

 

1

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rM0M

 

r0 rM0*M

 

 

 

 

где G

 

∑

= G (P, M

0

) = 0,

а v (M ) = −

R

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r0

rM0*M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем ∂G

=

∂G . Для этого введем некоторые обозначения и

 

 

 

∂n

 

∂r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rM0M

 

воспользуемся

теоремой

косинусов

для

нахождения

 

и

rM0*M :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

OM = r = x2 + y2 + z2 , OM

0

= r = x2 + y2

+ z2 , M

0

OM = ϕ,

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25

Источник: https://studfile.net/preview/16724467/