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

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

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

И ^ Xty,^ dSl

1

JJLx,yJdS.

(5.20)

4xa2t

$ ( 0, 0,0)

 

s lx*y>z )

 

 

1

 

at

 

 

Из (5.19) и (5.20), учитывая,

что,

согласно

(5.16),

А^. л,

получим vtt = a 2Av. Так как функция

ф,па Si°,0,0) непре­

рывна вместе с производными первого порядка, то существует М >0 такое, что| Ф|$1(одо)<./И’и |grad^lSi(o,o,o) < А<Г.Тогда из (5.17)

и (5.18)

следует, что

*и|/—о = 0

и Vf= ф(£.’1.065^ ) + ' ^ JJ (gradф, h )d s — ty(х, у, z)

 

$(0,0,0)

при at-*-0. Тем самым доказано, что решением первой из задач (5.14) является функция v из (5.15), а решением задачи (5.13)— функция

“<*■»• *>~7г(3^7 Я*«• *

Я ♦«'”•;)rfs'

S(x>!/>*)

s<x>y>z)

at

at

 

(5.21)

В случае неоднородного уравнения с правой частью f(x, уу z, /), применяя метод Дюамеля, получим, что к выражению в правой части равенства (5.21) следует прибавить слагаемое

0

S(*,y,2)

 

a{t—т)

Делая в этом интеграле замену a(t—т )= г , преобразуем его к

виду

 

y = J ï4ла*гS 7

Я / (5, * • с> l~~г ) ^ =

о

s(x>y>z)

 

Г

" ï Sу(х,у,г)r f f î T

^ h ' ^ - T

at

 

Полученное интегральное представление решения задачи Коши

, х à

и(х, у , 2 , /)= • dt f-d* ЯтЙ1'• °*)+

V

1

+ * = * f i

* » * + d * Î J J

$(*,{/.*)

Vr(r*^>z)

at

at

(5.23)

называется формулой Кирхгофа.

Интегральное представление решения задачи Коши для вол­ нового уравнения относительно

 

 

и{х,

у, t)

 

 

 

 

 

 

 

till—a2(uxx-\-UyU) в области

 

 

 

{U , y)c=R, if> 0 } ,

 

 

 

a|*-o=<pU, У), й/|/-о= Ф(а:,

*/)

 

 

 

 

V U , y ) ^ R 2,

 

 

 

 

а6 получим из формулы Кирхгофа ме­

 

 

тодом спуска.

 

переменное

Рис. 5.2

 

Введем

фиктивное

 

z и положим

 

 

 

и(х, у , z) = a U , у , z,

О

V U ,

у , Z )œ R3 и V / > 0 ,

(5.24)

<pU, у) = у{х, y t z \

ф и , #)=<|>и, у,

z)

 

V ( x ,y ,z ) ΠR 3.

 

По формуле Кирхгофа

(5.21),

учитывая,

что в силу равенств

(5.24) интегралы по сфере S¥t’v,z)с центром в точке

(*, у, z)

рав­

ны интегралам по сфере S%t,y,z) с центром в точке (я, у , 0), имеем

a U ,

у , /):

д

/

1

j'J* <р(£> “Л, C)ûf5^-f-

 

(9/

 

4яд2/

 

 

 

 

 

$(*,0,0)

 

 

 

 

 

 

 

 

at

 

 

 

 

 

1

 

J j Ф«. ч, Orfs.

 

 

 

 

+

 

 

 

 

 

 

5(-Г,У,0)

 

 

 

 

 

 

 

 

at

 

 

 

 

Сведем поверхностные интегралы к двойным по кругу

в

плоскости Оху с центром

в точке

(х, у)

радиуса at (рис. 5.2).

Учитывая, что

 

da

 

у — угол единичного вектора внеш-

d s = ------, где

 

cos Y

 

 

 

 

 

 

ней нормали к сфере S /t,y,0)(

в точке Л4(£, TJ, £)

с осью z VL da

 

 

Ы Х , у ) -

cos Y =

V (Л^>2— (ç— JC)2— (-n— w)2

элемент площади круга K at

»что

— - — -— ---

-

a t

и что в силу равенств (5.24) подынтегральные

функ­

ции четны относительно £ и, следовательно, интегралы

по пол-

ной сфере Sat’y't) равны удвоенным интегралам по верхней полу­

сфере, получим

- ^

V.JJi.a.W - i t - x P - h - y ) 2

 

«(*. ».« - 4

 

1

ГГ

at

 

______ Ф(£, ri)dldi\_______

(5.25)

2яа

JJ

V(aOa- ( 6 - ^ ) 2- ( 4 - ÿ ) 2

 

 

К(*>У)

 

 

В случае неоднородного уравнения с правой частью f(x, у, t) согласно методу Дюамеля к выражению (5.25) следует приба­ вить слагаемое

J =

1

W f

_________ /(€.ч,0<№1

(5.26)

2па

V û '2 ( t — T ) 2 — ( £ — J f ) 2 — ( î | —

t j ) 2

X(ÏJy)>

Полученное согласно (5.25), (5.26) интегральное представление решения задачи Коши для волнового уравнения относительно

и(х, у, t) называется формулой Пуассона.

§ 5.2. ЭНЕРГЕТИЧЕСКОЕ НЕРАВЕНСТВО. ЕДИНСТВЕННОСТЬ РЕШЕНИЯ ЗАДАЧИ КОШИ. УСТОЙЧИВОСТЬ РЕШЕНИЯ ЗАДАЧИ КОШИ

Рассмотрим заДачу Коши иц=й2Аи в полупространстве

Я++1= { (* , t)\x<=Rn, / > 0},

 

a|* -o= ?U ), и*|/-о=?(*) Y x ^ R n.

(5.27)

Введем энергетическое неравенство, рассматривая для нагляд­ ности свободные колебания мембраны. Пусть в положении равно­ весия мембрана находится в плоскости Оху. Выделим в мембра­ не круг Q радиуса R с центром в,точке (х, у). Пусть в момент времени / = 0 вне этого круга и на его границе возникли возму­ щения. Пусть скорость распространения волн от возмущений рав­ на а и пусть и(х, у, /) — отклонения точек (х, у) в момент t от

положения равновесия. В момент времени t—x, 0 < т < — , волна

CL

от возмущений вызовет колебания точек круга Q, отстоящих от границы круга на расстояние, не превосходящее ах, а невозму­

щенными останутся

точки

концентрического круга QT радиуса

R—ax. Изобразим круг Qx в плоскости /= т

(рис. 5.3).

Полная энергия

Е(х)

невозмущенной

части

QT равна (см.

§3.4)

 

 

 

 

Е ( t ) = -J- JJ

[ря?+1Л gradjw] dxdy.

(5.28)

 

а

 

 

 

Из физических соображений можно заключить (строгое матема­

тическое доказательство мы опускаем), что для V t,

0 < т <R}a,

£ ( т ) < £ ( 0 ) .

(5.29)

Это неравенство называется энергетическим неравенством. Оно имеет место и для решений задачи Коши (5.27) в случае пф2.

Т е о р е м а е д и н с т в е н н о с т и .

Два решения задачи

Ко­

 

 

ши (5.27), удовлетворяющие од­

 

 

ним и тем же начальным услови­

 

 

ям, совпадают для\i

t > 0.

 

 

 

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

 

 

и2— два

решения

задачи

Ко­

 

 

ши (5.27). Тогда и = щ и2— ре­

 

 

шение задачи Коши

 

 

 

 

ult= a 2ÙM>

и|/-0'=О, Я/|/2о=0.

 

 

Фиксируем

произвольное

т>0.

 

 

В плоскости t = х возьмем произ­

 

 

вольный круг QT и построим ко­

рИС4 5 з

 

нус с осью, параллельной оси /,

 

образующими, составляющими с

 

 

осью

t

угол

a = a r c t g а, и

осно­

и £=0. соответственно (рис.

ваниями Qx и Q в плоскости t=%

5.3). Согласно энергетическому не­

равенству (5.29) и равенству

(5.28)

имеем: Уt, 0< t< R /a ,

 

0 < £ ( T) < £ ( 0 ) = Y

jj*

(ри?+ |Г|gradîtt)|/=0d x d y —0

 

 

Q

 

 

 

 

 

 

ИЛИ

 

 

 

 

 

 

 

J j (рИ/+

|Г| gradua) dxdy.=0>

 

 

откуда в силу непрерывности подынтегральной функции следует, чтой*=0, grad* и = 0. Поэтому ws=const внутри построенного ко­

нуса. На основании конуса Q по условию йн=0. Следовательно, и==0 внутри конуса. В силу произвольности т и Qt получим, что

йен0 и, следовательно, и \= и 2всюду в полупространстве /?++1.

Из теоремы единственности следует, что всякое решение за­ дачи Коши в случае я = 1 , 2, 3 может быть записано соответст­ венно по формулам Даламбера, Пуассона, Кирхгофа.

Т е о р е м а ( у с т о й ч и в о с т и ) !

Пусть ui(x,t ) и u2(xf t) —

решения задачи Коши ,с начальными

данными соответственно

«PiM» Ф1 (х) и Ф2.М» фа(*). Пусть Т > 0 — произвольное число.

f i
$(0,0,0) 1

Тогда для V

е > О з б (е, Т) такое, что из неравенств

 

l? iC * )

<p2(-*0l<!8)

14*1 (-^) — Ф2 (-^)1

S, V X œ R 1 ,

(5 .3 0 )

следует неравенство

 

 

 

 

М * . *) — и2(х, 0 | < е ,

Y X œ R1 и

V/,

0 < * < 7 \

(5.31)

Д о к а з а т е л ь с т в о . Представляя

решения и{ и и2 с

помо­

щью формулы Даламбера, имеем

 

 

 

l«iU. О —й2и*, 01 <

—«О —?2U —^01+

 

«г+а/

+ ÿ l? i( - * + « 0 - ïï( .* + a O I + ; g -

j № i ( » - > b « ) | r f ? < Ÿ 8+

 

x —at

“b "7Г ®H— &*2л^^ S (1 -|-7’)<[ e

2

2a

 

П р и 8 < 8 ( e , T) = e / ( T - J - 1 ) .

 

2. Пусть щ (х, у, г, i) и и2 (х,

Т е о р е м а ( у с т о й ч и в о с т и )

уу z, 0 — решения задачи Коши с начальными условиями соот­

ветственно у\(х, у, z) o|)i

(х, у, z)

и ср2(х,

у, z),

ty2(x, у, z). Пусть

5Г>0 — произвольное число. ToedaV s > 0

3ô(e, f ) > 0 такое, что

из неравенств

 

 

 

 

 

|?1 ?2!<&> |grad(<PÏ —<р2)|<8,

|<!>1—ф2|< 8 ,

V(x> tjy Z)œ R*

следует неравенство

 

 

 

 

(5.32)

 

 

vO

0 < / < 7 \

(5.33)

|й1 —и2|< е V U , у, z)<=R* и

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

Положим и= щ —и2. Тогда функция и

будет решением задачи Коши с начальными, условиями ф1—<р2 и

—^2. Записывая решение и по формуле Кирхгофа (5.23), пре­

образуя интегралы по сфере 5 if,ÿ,z)B интегралы по сфере

*Si°,0,0>

и дифференцируя по / [см. равенства

 

(5.16) — (5.18)], получим

 

V

J

*rfs=

 

s(*M

 

= 1 7 И 'tdSl + ^

(grad?,

 

f f ♦ * ,.

(5.34)

$ ( 0,0,0)

1

где n — единичный вектор внешней нормали к учитывая условия (5.32), находим

$ ( 0,0,0)

Sa^x,v,z). Отсюда,

и(х,

я* 0 ^

5 4 я - | - 8 4

я = 8 [ 1 + ^ ( а Ч"^)]<Се

 

 

 

при 8 < 8 ( е , Т ) =

[ +Г(е- + |) . Теорема

доказана.

Так как формула Пуассона является следствием формулы Кирхгофа, то эта теорема будет справедлива и для решения за­ дачи Коши (5.27) при л = 2 .

З а м е ч а н и е . Следует обратить внимание на тот факт, что в условии (5.32) теоремы 2 в отличие от условия (5.30) теоремы 1 мы потребовали малость модуля |grad(<pi—<рг|. Это не связано с методом доказательства, а соответствует физической сущности волнового процесса. Действительно, из (5.34) при ф = 0 , | ф|< 0 ,

2й |grad<p|^ гй > 6 и i l > - — имеем

\и(х, у, г, t)\^atk—

Это означает, что расположенное на расстоянии г от точки М (я, у, г) возмущение ф, малое по модулю, но с большим модулем градиента, распространяющееся со скоростью а, через время f = = г /а вызывает в точке М «всплеск», величина которого превос­

ходит -£-r|grad<p|.

§ 5.3. ЗАДАЧИ КОШИ ДЛЯ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ. ФУНКЦИЯ ИСТОЧНИКА. ИНТЕГРАЛ ПУАССОНА

Рассмотрим задачу Коши для уравнения теплопроводности

и ,= а 2Ди в области

{(я, /)|д :е /? я, /> 0 } ,

(5.35)

tt|/-o=<p(jc)

V x ^ . R n%

 

ф(х)~— ограниченная функция.

Для простоты рассмотрим случай л = 1. Аналогичные рассуж­ дения справедливы и для произвольного п. Убедимся в том, что

1) функция

 

и(Х' / ) = 1 ^ Ь г

)

? а )

 

(5.Э6)

приУ R* n V t > 0

удовлетворяет уравнению

теплопроводно­

сти uta2Uxx=0;

 

 

 

 

 

имеет

место ра­

2)

во всякой точке XQ непрерывности ф(*)

венство

 

 

 

 

 

 

 

 

 

 

<p(*o)=lim

Ж*, О.

 

 

Положим

 

 

,/) - ( х 0, +0)

 

 

 

.

 

_

(j—б)‘

 

 

 

 

_

!

« • '

п р и

 

0(х, t, » =

 

_

е

 

2а У ль

 

 

 

(5.37)

 

 

 

О

 

 

при *< 0 .

 

Тогда равенство (5.36) можно записать в виде +**•

U(x, t ) = j (р&)G(X, t , \)d\.

(5.38)

Интеграл в выражении (5.36) называется интегралом Пуас­ сона, а его ядро G(x, t, £), определяемое равенством (5.37),—

функцией Грина или функцией источника.

Подставляя в уравнение теплопроводности производные

с , =

1

г

 

_L. Д2(Х — Ç)2

,

 

(.c-Ü *

 

J

р

e

4a1/

 

2 Ÿ n

L

2 (Д2/)3/2

i

4 ( a > /) 3/2

 

 

I

Г

1

| .

( x - W

.

P

 

< * -£ )*

о X X '

1

 

4a*/

2 К я

L

2

1

4{аЩза

I

e

 

 

 

 

убеждаемся, что Gta2Gxx= 0.

По условию функция ф(х) ограничена, т. е. 3Â f>0 такое, что

|<p(jt) I < M v хщЯК

Далее, V

х,

\х\<А, nVt,

 

 

 

имеем

\х—Ъ\2> \1—А \й приУ£>Л и

 

 

 

 

 

 

 

 

 

«

(*-£)»

 

.

(е- л

 

 

 

 

 

1 _

ç 4а1/ ^

 

1

е 4а5/,

 

 

 

 

2a

Y n t ‘

^

2а У nti

 

 

 

 

 

Тогда функции <p (|) G (х, /,

£), ф(îj) Gt (х, t, £) и

tfÇ)0«(*>

О

 

 

 

 

 

 

 

/1/1

— i

1

(е-л)*

мажорируются

функцией

F&y

м

0

2

4а»/, »

2а Уnti

 

 

 

 

 

 

 

 

 

 

 

 

 

для которой J F{\)d\ сходится,

и, следовательно,

 

интегралы

 

А

 

 

 

 

 

 

 

 

 

 

• о

в

^

 

 

 

 

 

 

 

 

 

j<р(?/Д,

J <pG^,

f yGxxd\

равномерно сходятся по (х, t),

\х\^.

л

а

л

 

 

 

 

 

 

 

 

 

г^А и

 

 

Аналогично, при тех же х, t устанавливается

 

 

 

 

 

 

 

—А

 

 

—А

 

равномерная сходимость

интегралов

J

<pGd£,

 

J

<?Gtdlt

- А

 

 

 

 

 

 

— м

 

 

— оо

 

 

 

 

 

 

 

 

 

 

 

 

^<t>Gxxdb. Отсюда следует, что интеграл Пуассона можно диф-

^ • в о

ференцировать по параметру под знаком интеграла и для него справедливо равенство

+•»

ut—a?uxx= J <р(0[<?/(*, A l)—a2Gxx{x, t, S)]<ft=0.

— оо

Тем самым доказано, что интеграл Пуассона является решени­ ем уравнения теплопроводности.

Докажем теперь, что если функция ф(*) непрерывна в точке *о, то '

?U*o)=Hm и(х, t).

(x,t)-+(xо, +0) 1

Фиксируем е> 0 и построим окрестность а точки

(х0 , 0) такую,

ЧТО | и (х, t) — ф (А'о) | < 8 V

{x,t)Œa.

 

Полагая

Р— X

= z y получим

 

*

_

 

+ -

 

 

~

+~

~

 

 

 

 

 

 

+“

(JT-D*

 

 

 

 

 

— Г e - * ' d z = s l

 

 

 

 

/ Я

J

и

 

 

i~ от

+ во

 

 

 

 

 

И * , t) ср(л'0)| =

\ ) d \ — <рС*о) J G ( x y /, £)^£

+ °*

 

+ °°

 

__

 

 

<

j

l ? t t ) - < p U o ) |û ( * ,< . О di = j

\f(x+2aVtz)

 

 

 

 

 

- I V

N

 

 

 

 

 

-

?(x0) | - J ^ e - ^ = j + J

+ J

(5.39)

 

 

 

—O,

—TV

N

 

 

Так как

J

erz'dz сходится и | ф (* ) | <M V R1, то

м о ж н о .

 

—••

 

 

 

—IV

"

выбрать столь большое N, чтобы каждый из интегралов

J

и J

 

 

 

 

 

 

—во

N

 

 

 

'

 

 

 

N

был меньше е / 3 .

При выбранном N рассмотрим интеграл

^ .

-IV

Так как функция ф(х) непрерывна в точке х0, то 3ô(e) > 0 такое, что |ф (а + 2 a V tz ) —ф(*о) | < е/3 при .|* + 2 Y^lz—лго| < ô . Мы удо­ влетворим последнему неравенству при всяком 2 е ( —Ny N), ес­

ли выберем такую окрестность а точки

(лг0; 0), чтобы V

(я, t) œ <J

выполнялись

неравенства

,\х—*о|< 0/2

и 2aY tN < àf2

(или t<.

 

 

N

N

 

 

 

^ (law )2)'

Т°гда будет

f <C ~

j* e~z*dz<^

Таким обра-

 

 

- N

— N

 

 

 

30M , из (5.39) получаем

И * . *)—*C*o)l<

Функция G{x, t, I) не случайно называется функцией источ­ ника. Она определяет в области |*|< оо, £>0, поле температур,

создаваемое единичным точечным источником, сосредоточенным

в точке 5.

Это можно объяснить таким образом. Пусть функция (ре (sJ удовлетворяет следующим условиям:

1)при |И1>е>0;

2) <р,Ш>0 V5, 151 <оо;

3)| <p.«)rfS=l.

Эту функцию можно рассматривать как плотность распределе­ ния источников тепла, выделяющих суммарно единицу количест­ ва тепла.

Пусть в монотонно стремится к нулю. Естественно при этом рассматривать lim<pe(£) как единичный точечный источник тепла, сосредоточенный k точке 5= 0 . При (р=<рв темпе­ ратурное поле определяется интегралом Пуассона

 

1

+-

(х-Е)»

 

 

и,{х,

] ?.(5)е 4а,/

ûf5 =

*) = 2a Y ni

 

Î

(х-Е)'

 

 

 

 

Л

 

 

 

Ш < «

 

 

 

 

Учитывая, что (|* | —е)2< \х| | 2^ (\х\ + 5)2

Прй | е | < е и ч т о

фе(5) удовлетворяет условиям 2 и 3, получаем

 

 

1 = г е

(ln+е)1

 

1= - е

(Ш-О*

4аи < « , ( * , * ) <

4а4/

2a Ynt

 

 

2a Ynt

 

 

откуда следует, что

 

 

 

 

 

1

—J£ll

 

t,

 

 

' ) = 1 7 ? 5 Г е 2“’'

 

0).

В заключение заметим, что решением задачи Коши

щ = а 2Ьл в /?++1= {л ;е/?я, *>()},

и |/ - о = < р ( л : )

является интеграл Пуассона

U(x, < > - J ,«)< ?(* . ,,

где

 

e

_!£=!!!

G(x, /, Ç)=

ЛаЧ при

(2а ■/яt)n

 

 

0

при /><0

и

i * - s i J= i2- 1

ГЛАВА 6

КРАЕВЫЕ ЗАДАЧИ ДЛЯ УРАВНЕНИИ ЛАПЛАСА И ПУАССОНА

§6.1. ГАРМОНИЧЕСКИЕ ФУНКЦИИ. ПРИНЦИП МАКСИМУМА. ФОРМУЛЫ ГРИНА. ЕДИНСТВЕННОСТЬ И УСТОЙЧИВОСТЬ

'РЕШЕНИИ КРАЕВЫХ ЗАДАЧ

Пусть DczR2— ограниченная область с границей Г.

Гармонической в области D функцией называется дважды не­ прерывно дифференцируемая функция, удовлетворяющая в этой области уравнению Лапласа Aw=0.

П

где г =

Например, и = const; и==2

л-1

 

— гармонические функции.

 

Как известно из теории функций комплексного переменного, гармоническими функциями двух переменных являются действи­ тельная и мнимая части аналитической функции. Так, взяв ана­ литические функции ez, In z, zn, где z —x+iy, получим соответст­ венно следующие гармонические функции:

и(х, y).=Reez= e xcosy; v(x,

ÿ ) = I m e z= e x sin у;

u(xt y ) = R e \n z = I n \z\

In(x2+ y2);

 

Z

v(x, v )= In iln 2 = a r g 2 = a rctg — + 0Jt» a==l ?’ x '> ®t

* l 1, * < 0 ;

u(x, y ) = R e z 2= x 2—y2; u(x, y ) = R e z n= r ncos n<?; v (x f y ) —\m zn—rns\nti<çt где r —\z\, <?=argz.