Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Методы математической физики и задачи гидроаэродинамики

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
14.42 Mб
Скачать

 

Д о к а з а т е л ь с т в о . Окружим лежащую в D точку х0 сфе­

рой Se с центром лго такого радиуса е,

что ограниченный

 

ею

шар Qe целиком принадлежит области D. Удалим

 

из области

D шар Qe (рис. 6.3). Для функций и{х)

и G(x, х0)

в получен­

ной области D \ Q e с границей ruS8 согласно

формуле Грина

(6.6)

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

[u(x)AxG(xt x0)—G(x, x0)ku(x)]dx—

 

 

 

 

 

 

D \Q ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d G ( x ,

Xp)

G(x, xQ) du (x)

 

 

(

6

. )

 

 

 

 

 

 

 

 

dnx

 

 

dax

 

 

 

 

22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где A X G (X , X 0) — оператор

Лап­

 

 

 

 

 

 

 

ласа

от

G(x,

хо)

по

переменной

 

 

 

 

 

 

 

х,

пх— внешняя

нормаль

к гра­

 

 

 

 

 

 

 

нице

области

D \ Q B,

которая

в

 

 

 

 

 

 

 

точках

х е Г

остается

внешней

 

 

 

 

 

 

 

нормалью к Г, а в точках S8

 

 

 

 

 

 

 

является

внутренней

нормалью

 

 

 

 

 

 

 

к

сфере Se

(рис.

6.3).

Левая

 

 

 

 

 

 

 

часть равенства (6.22) равна ну­

 

 

 

 

 

 

 

лю,

 

так

как

в

области

D \ Q

 

 

 

 

 

 

 

 

AxG(x, Х о ) = 0

и Аи(х)=0. Пре-«

 

Рис.

6.3

 

 

 

образуя

правую часть равенства,

 

 

 

 

получим

 

 

 

 

 

 

 

ds—j G(x, xQ) да (JC)

 

 

 

 

 

 

 

 

 

 

д О (х,

х 0)

 

 

 

 

 

 

 

 

 

 

 

дпх

 

 

 

г

дпх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dG(x, дг0)

 

 

 

 

 

(6.23)

 

 

 

 

 

 

 

дпх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как G(x,

x 0) | x t r = 0 ,

то

второй интеграл

в

правой

части

равенства (6.23)

равен

нулю.

Согласно

свойству

Г функции

Грина в силу непрерывности

 

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дпх

 

 

 

 

 

 

 

 

 

 

 

 

Иш f

G(x, x 0) - ^ d s —О,

 

 

 

 

 

 

 

 

 

 

 

t-*o

J

 

 

 

дпх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

st

 

 

 

 

 

 

 

 

 

 

 

a на основании свойства 2° в силу непрерывности

и(х)

имеем

 

 

 

 

Ига f

и(х)

 

 

ds=u(x„y

 

 

 

 

 

 

 

 

 

 

•н-о J

 

 

 

 

дпх

u

 

 

 

 

 

 

Переходя в равенстве (6.23) к пределу при в-»-0, учитывая изло­

женное выше, получаем

о = Г « ( * ) да(х' Хо)- d s + u ( x Q),

 

 

 

 

 

.)

дпх

 

 

 

 

 

 

 

откуда следует доказываемое

интегральное

представление

(6.21) гармонической функции и(х).

полупространства,

круга,

Построение функции

Грина для

шара. Пусть D — полупространство в # 3, плоскость 5 — его гра­

ница, Хоe D — произвольная точка.

Функцией

Грина

 

первого

 

рода для полупространства D бу­

 

дем называть функцию

 

 

 

 

 

G(x, х0)=£(Ху *0)+«(■*» л:0),

 

 

где

& (я, х0) — главное

фундамен­

 

тальное решение оператора

Лапла­

 

са, v { x y

Хо) — гармоническая в об­

 

ласти D и стремящаяся к нулю на

 

бесконечности

функция

такая,

что

 

G ( x , Хо)

на границе 5 тождествен­

тельно, для построения

но обращается

в

нуль.

 

Следова­

функции

Грина

G { x t хо)

надо

по­

строить функцию v ( x , хо)

так, чтобы

 

 

 

 

 

 

 

•о (ХуxQ)Ues = — <$(Ху х0)\xes=

— г——-------1

 

.

(6.24)

 

 

 

4л (ДГ—х0) Ues

 

 

 

Если обозначить через х0 точку, симметричную х0 относительно

плоскости Sy T O V

XS имеем

\xiXo\ = \xi—xQ\ (рис. 6.4).

Поэтому, полагая

v ( x ) = -

1

__— , получим, гармоническую

 

4л 1х

х 0 1

в D функцию, удовлетворяющую условию (6.24), а тогда будет справедливо равенство

^U,JC0)U6S= 0 .

Итак,

G(XyX0) =

I

 

I x Xo I

*ol ]■

 

Пусть D — круг радиуса R с центром в точке О, Г — его граница. Пусть Хо— произвольная точка круга D. Для нахож­ дения функции Грина первого рода в круге D надо построить

гармоническую в D функцию о(х, х0)

такую, чтобы

функция

G(Ху Хо) =<§Г(Ху Xo)+v(x, Хо)

удовлетворяла

условию

О [xt Хо) 1*6$ = 0 , т. е. чтобы

 

 

 

*о(Ху л'0) |,es = 2 U , л:0) |,65= —

 

Ь ------!-----

 

*0| *е$ *

Обозначим

через

*о* точку,

инверсную точке Хо относительно

окружности

Г, так что Хо* расположена

на луче Ох0 и рр*=

= R 2, где

р, р* — соответственно

расстояния

точек хо и х0* от

точки О (рис. 6.5).

 

 

 

 

 

точка

окружности

Г.

Пусть

теперь

х\ — произвольная

В &OX\XQ и кОх\Хо*

стороны, образующие угол с вершиной

в

точке О,

пропорциональны: — = - ^ - .Следовательно,

эти тре-

угольники подобны и

 

R

 

р*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_£i_—£о_[__1L

 

 

 

 

(6.25)

 

 

 

 

 

 

 

xi ~ x Q

R

 

 

 

 

 

 

Тогда,

полагая

v {x ,x Q)—

In PU

R

 

 

получаем, что

 

 

 

функция и(х, Хо) — гармоническая по х в круге D и

 

 

О (х, л0) Uer=

——• In

 

I

 

In-

 

 

 

= 0.

 

I * — I

?\х —х0

 

Итак,

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

G(x, х0)

in—

!—

 

 

 

 

(6.26)

 

 

 

 

 

 

I

 

— *0 I

 

 

 

 

 

 

Пусть D — шар

радиуса R с центром

в точке

О, Г

его

граница,

хо — произвольная

 

точка

 

 

 

 

 

 

шара D. Для построения функции

 

 

 

 

 

 

Грина

первого

рода

в области

D

 

 

 

 

 

 

надо

найти

гармоническую

в

D

 

 

 

 

 

 

функцию v(x, х0) такую, чтобы вы­

 

 

 

 

 

 

полнялось условие

 

'

 

 

 

 

 

 

 

 

 

G(x,x0)jxer=

1#(х,х0)-1-ъ(х, х0)]хбг=0

 

 

 

 

 

 

т. е. чтобы

 

 

 

 

 

 

 

 

 

1

Рис. 6.5

 

*уи,л:0)и ег= —^(х,х0)\хет=

 

 

.сбг

 

 

 

 

 

 

 

 

через х0* точку

 

 

4я | .V — лг0 |

 

Обозначая

на луче Ох0 такую,

что рр*=/?2,

где р и р* — соответственно расстояния точек х0 и А'о* до точки

О, как и в предыдущем случае для любой точки .VIœT,

полу­

чаем равенство (6.25).

 

ft

 

 

Полагая v (л:, х0) ■

 

имеем v(x, Хо) — гармони-

4л | лг — JTSI р

 

 

 

ческая по х в шаре D функция,

9 I х Ь 5г]

 

 

I Г

1

 

G (ХуЛд)—

I |

— -*o I

(6.27)

Интегральное представление решения задачи Дирихле в кру­ ге и шаре. Пусть D — круг радиуса R с центром в точке О, Г — его граница. Рассмотрим внутреннюю задачу Дирихле: найти функцию U œ C2(D)Ç\C(D) такую, что

Аи=0 в

D ,

 

 

 

 

(6.28)

 

 

 

 

 

 

 

Пусть существует решение

этой задачи

и пусть также K G

œ C2(D). Тогда

по теореме

об Инте­

лу, тральном

представлении

гармониче­

ской

функции

с

помощью

функции

Грина имеем: V XD

 

 

 

 

 

u(xQ) = — f f ix )

dG- (-x '*°-- dst

 

 

 

 

J

 

 

dnx

 

 

 

 

 

г

 

 

•*

 

 

 

 

 

 

 

 

 

(6.29)

где

fix— единичный

вектор

\

 

внешней

нормали

к Г в точке х (рис. 6.6),

а

G(x,

Хо)— определяется,

равенством

(6.26). Пусть точка Хо* инверсна точ­

ке хо относительно

окружности

Г.

Обозначим через р и р* соответственно расстояния точек х0 и хо*

до точки О. Положим г = \ х х0\,

г * = Хо*\.

Найдем выра-

жение для

dQOc.xo) Имеем

 

 

 

дпх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d G ( x ,x о)

 

д

М

Л

----- 1п-^---- In— ] ] =

дпА

 

дпх

1

[

 

 

1

Г 1

дг

1

дг* I

 

 

г

дпх

г*

 

 

 

 

 

1

-(gradr*,îïx)1 -

1

I

 

 

г*

 

 

 

 

\

 

 

где a = (g ra d r,

nx)t

a*=(gradr*, пх) (здесь

учтено, что

| grad г| = 1 ).

 

Л

 

 

 

 

У?2 — Г2 _

р2

 

 

 

 

 

 

Из АОххо и АОхо* имеем

c o s a = -------— и cosa* =

- f г *2 — р*2

Подставляя эти выражения в выражение для

2Rr*

 

 

 

 

 

 

 

 

 

дО_ и учитывая,

что

р р * = # 2

и что, согласно (6.25),

дпх

г * = — г, получаем

р

д а ( Х , Х о) _

1

П

Ф + Г 2 - р2

р

/ ? 2 ( / ? г / р ) 2 - ( / ? 2 / р ) 2

 

 

2 я

I г

2/?г

7?г

2# (/? г/р )

 

 

 

 

 

 

/?2-р2

 

 

 

 

 

 

(6.30)

 

 

 

 

 

 

2Л/?Г2

Из ДО**0

имеем

г2= /? 2+ р 2—2/?pcos6,

где 0 — угол, между

Охо и Ох (рис.

6.6). Вводя полярную систему координат с по­

люсом О,

обозначая

через р, <р и R, ф полярные координаты

точек х0 и х, получаем, что cos 0= cos (ср—ф), а дифференциал длины дуги окружности Г равен ds—Rd\\>. Тогда в соответст­

вии с равенствами (6.29) и (6.30) имеем

интегральное

пред­

ставление решения задачи Дирихле в виде

 

 

и( р,<р)=-

 

 

 

/?2 — р2

 

 

 

 

j / «

>*)

а д .

(6.31)

2nR

 

 

 

/?2

р'2 _ 2 /? COS ( f — 40

 

Входящий в равенство

(6.31) интеграл называется интегралом

Пуассона,, а его ядро

_________ /?2-р2_________

 

 

Р (X, •£())—

 

 

 

 

 

2я/? [/?2 +

р2 —2/?р cos (<р —4')]

 

— ядром Пуассона.

 

 

представление

реше­

Аналогично

получается интегральное

ния внутренней задачи Дирихле

(6.28) в шаре DczR3 с грани­

цей Г. Для решения этой задачи также имеет место равенство (6.29), где G(xf х0) находится по формуле (6.27). Выполняя в плоскости Оххо те же построения, что и выше, получаем

àG ( х , х 0)

____ д т J

/ J

____Д \ 1

дпх

дпх [ 4я

{ г

рr* J J

1

/ 1

R п

*\

 

Я2- Р 2

[ -----cos а ----------- cos а*

)

-----------— =

\ г2

рг*2

 

4я/?г2

_______ /?2- Р2________ = Р ( х ,х 0)

(6.32)

 

4JtR (/?2 +

р2 _ 2/?р cos 8)3' 2

 

.Тогда из (6.29) и (6.32)

следует

искомое интегральное

пред­

ставление решения, называемое интегралом Пуассона: YXQ^D

(

и(х )0= f / (х)Р(х, x0)ds.

(6.33)

Ядро интеграла Р(х, Хо)

называется ядром Пуассона.

так­

В случае шара DczRn решение задачи Дирихле (6.27)

же представимо

в виде

интеграла

Пуассона, ядро которого

P {x ,x Q)

_______ R2-

P2________

(6.32')

 

*

“яД (Я 2 -ЬР2 -2 Я р со 8 0)л/2

где (ùn— площадь единичной сферы в Rn.

З а м е ч а н и е . Мы получили представление решения задачи Дирихле в круге или шаре в виде интеграла Пуассона в пред­ положении, что решение существует. Можно доказать, что для

любой непрерывной функции / на границе Г круга

или ша­

ра интеграл Пуассона (6.32') является решением

задачи Ди­

рихле.

 

§ 6.3. СВОЙСТВА ГАРМОНИЧЕСКИХ ФУНКЦИЙ

 

В ОГРАНИЧЕННЫХ ОБЛАСТЯХ.

 

О РАЗРЕШИМОСТИ КРАЕВЫХ ЗАДАЧ

 

ДЛЯ УРАВНЕНИЙ ЛАПЛАСА И ПУАССОНА

 

1°. Интеграл от нормальной производной гармонической функции по любой замкнутой поверхности, принадлежащей вмес­ те с ограниченной ею областью области гармоничности, равен нулю.

Пусть и(х) — гармоническая в ограниченной области Q c ^

функция и Г — произвольная

гладкая замкнутая

поверхность,

принадлежащая области Q

вместе с ограниченной ею

об­

ластью D. По формуле Остроградского— Гаусса, получаем

 

J

= J budx= 0 .

(6.34)

ГD

2°. Гармоническая функция имеет производные любого по- рядка.

Рассмотрим

произвольную

точку х0еЙ . Построим шар D

радиуса R, целиком принадлежащий области Q и содержащий

точку х0. Пусть

Г — граница

шара D0. Согласно равенству

(6.33), гармоническая функция и(х) в точке х0 представима интегралом Пуассона. Этот интеграл — собственный, зависящий от координат точки х0 как от параметров. Ядро Пуассона име­ ет непрерывные производные по параметрам любого порядка. Поэтому интеграл Пуассона можно дифференцировать по пара­ метрам под знаком интеграла любое число раз, откуда следует существование в точке Хо производных и(х) любого по­ рядка.

(первая теорема о среднем). Пусть и(х) — гармоническая

в области QczRn функция и Q произвольная

точка. Пусть

D шар радиуса R с центром х и границей Г и /)с Й . Тогда

и (л:) = ~ ^ и (х ) ds>

(6.35)

г

 

где S площадь сферы Г.

Действительно, подставляя в выражение для ядра и инте­ грала Пуассона из равенств (6.32') и (6.33) х0—х и следова­ тельно, р = 0 , получаем

р ( х , х ) =

и “ W = T f “ W * -

 

 

Г

 

(вторая теорема о среднем). При том же условии, что и

в первой теореме о среднем, справедливо равенство

 

u ( x ) = ~ [ u ( x ) d x ,

(6.36)

 

О

 

где Vобъем шара D.

Пусть г, 0<r<R, произвольно. Пусть

Д о к а з а т е л ь с т в о .

Гг— сфера радиуса г с центром в точке х и Sr— ее площадь. По первой теореме о среднем, согласно (6.35), имеем

и (х)= -^~ j* u(x)dsn

где dsr— элемент площади сферы Гг. Умножая это равенство на Sr и интегрируя по г на отрезке [О, R], получаем

 

 

R

(6.37)

и (х ) j* Sfdr—J j* u(x)dsfdr.

0

 

0 rf

 

Учитывая, что dsTd r= dx — элемент объема шара D,

имеем

Г Srd r =

j* f

dsrd r = f d x = V ,

 

о

о vr

b

 

R

 

 

 

J J u(x)ds,dr— ^u(x)dx.

 

о rr

 

b

 

Подставив эти выражения в (6.37) и разделив полученное вы­ ражение на У, получим требуемое равенство (6.36).

(об оценке производных первого порядка гармонической функции). При том же условии, что и в теоремах о среднем V / = 1,...,п, справедливо неравенство

да

..'С — шах|й(л:)1.

(6.38)

дх I

х - х

R

г

 

Д о к а з а т е л ь с т в о .

Согласно

свойству 2°

гармоническая

функция бесконечно дифференцируема, поэтому

А-^ -(Л и )= 0 .

dxt

Следовательно, производные гармонической функции также яв­ ляются гармоническими функциями, и для них по второй тео реме о среднем справедливо равенство

да

да

~dx, I ■1| •••) я*

dxi

D дхI

 

Рассматривая ------

как div(0,..., О, и, 0,..., 0), где и i-я ком-

дхI

 

 

понента вектора, и применяя формулу Остроградского— Гаус­ са, получаем

да

и cos df ds,

дхi Х~Х ,

где т — угол внешней

нормали

к прверхности Г с осью Oxi.

Тогда

 

 

 

 

да

x-7

•< — max \ti\S.

дхi

V

г

1

Так как

j d x = ^ S , d r = ^ D О О

(здесь ©я— площадь единичной сферы в Rn), то

S_

i

да

 

шах \и\,

V

I

àxi х-х

R

г

чтр требовалось доказать.

(теорема об устранимой особенности). Пусть и(х) гармо­ ническая функция в проколотой окрестности точки Хо и и (х) =в = о (S’(х, хо)) при х-+-х0, где S’(х, х0) — фундаментальное реше­ ние оператора Лапласа. Тогда и(х) можно доопределить в точ­ ке хо так, что она будет гармонической в полной окрестности точки Хо.

Д о к а з а т е л ь с т в о . Пусть шар D радиуса R с центром Хо, имеющий границу Г, содержится в рассматриваемой окрест­

ности точки Хо. Пусть v(x) — решение

задачи Дирихле

 

 

Д-и = 0

в D,

 

 

 

 

 

тзlr=#Jr*

 

 

 

 

Рассмотрим

разность w —u—u. Имеем

ш |г = 0 . Так

как

V ограничена в D, а по условию

теоремы

и (х )= о (ёГ(х, хо))

при х-*-хо, то w = o (&(х, Хо))

при х-+#о.

Это

означает,

что

V > 0 з г ( е ) > 0

такое, что V г<г(в)

на сфере

Sr : \х—Хо\ = г

гЗ (AT^VQ) \Sf< w (x) |5f < гЗ (A * , * 0)|v

Фиксируем точку х* и выберем г<г'(е) так, чтобы х* находи­ лась в полом шаре Qe между сферами Г и S. Тогда для функ­ ции w, гармонической в Q8, согласно принципу максимума спра­ ведливо неравенство

Устремляя е к 0, получаем до(а'* )= 0, а в силу произвольности х* имеем до=0 в проколотом шаре D \ { x о}, а, значит, u = v в проколотом шаре D \ { x 0}.

Доопределяя функцию и(х) в точке х0 значением о(х0), по­ лучаем, что в полном шаре D функция и(х), как совпадающая

сv(x), является гармонической.

Оразрешимости внутренних краевых задач. Теорема об устранимой особой точке позволяет установить, что внутренняя задача Дирихле может не иметь решения, если область имеет

проколы. Рассмотрим

пример. Пусть

D={(x,

y)\Q<x2+ y 2<

<1} — единичный круг

на плоскости

с проколотым центром.

Тогда граница

Г состоит

из окружности у— {(х, у)\х2+ у 2= \

и центра 0(0,

0) круга. Рассмотрим задачу Дирихле: найти

функцию U œ C2 (D) Ç)C {D)

такую, что

 

 

 

 

 

àa—0 в D

 

 

 

 

1

в точке 0(0,0),

 

 

 

0 в точках окружности у.

Если бы решение и{х,

у) существовало, то как непрерывное

в ограниченной замкнутой области Д

оно было

бы ограничен­

ным. Но тогда

начало

координат О являлось

бы устранимой

особой точкой решения

и решение й(х, у) можно было бы до­

определить в начале коорЛинат так, чтобы оно в полном круге

ф=:{(л;, у) \х2+ у 2< \ ) являлось

гармонической

функцией. Со­

гласно принципу максимума для

и в Q в этом

случае имели

бы место равенства max w=max и—0, что противоречит тому,

1

что и(0, 0) = 1.

Существуют различные достаточные условия разрешимости внутренней задачи Дирихле для уравнения Лапласа. Так, для областей с гладкой или кусочно-гладкой границей без компо­ нент, размерность которых не превышает п—2, задача Дирихле разрешима при любой непрерывной функции f в краевом усло­ вии. Условия разрешимости для областей с негладкой грани­ цей формулируются более сложно и здесь не приводятся.

Внутренняя задача Неймана для уравнения Лапласа (Пуас­ сона) состоит в нахождении в ограниченной области D реше­ ния уравнения, непрерывно дифференцируемого в замкнутой области и удовлетворяющего краевому условию второго рода на границе Г области D:

да

(6.39) '

дп

Она разрешима не при любой функции р в краевом условии (6.39). Действительно, согласно свойству 1° гармонических

функций, предполагая,

что и е С 2(5 ), имеем

d s = 0, от-

куда для функции р из

(6.39) получаем

г

 

 

f prfs=0.

(6.40)

 

Г

 

Равенство (6.40) является, таким образом, необходимым условием разрешимости внутренней задачи Неймана для урав­ нения Лапласа Д «= 0 .

Вслучае задачи Неймана для уравнения Пуассона

Ди=/ в £),

да

интегрируя уравнение Пуассона по области D и применяя фор­ мулу Остроградского — Гаусса, получаем

f d x = J

badx —J ~ ^ d s — j* pds,

D

T

Г

Таким образом, равенство

 

 

£ f d x — J prfs

(6.41)

является необходимым условием разрешимости внутренней за­ дачи Неймана для уравнения Пуассона.

$ 6.4. ГАРМОНИЧЕСКИЕ ФУНКЦИИ В НЕОГРАНИЧЕННЫХ ОБЛАСТЯХ, РЕГУЛЯРНЫЕ НА БЕСКОНЕЧНОСТИ. ЕДИНСТВЕННОСТЬ ВНЕШНИХ КРАЕВЫХ ЗАДАЧ

ДЛЯ УРАВНЕНИЯ ЛАПЛАСА

Внешние краевые задачи ставятся для неограниченных об­ ластей D a R nt содержащих внешность шара достаточно боль­ шого радиуса р. Чтобы постановка краевых задач для таких