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

книги из ГПНТБ / Арсенин, В. Я. Методы математической физики и специальные функции учеб. пособие

.pdf
Скачиваний:
145
Добавлен:
22.10.2023
Размер:
13.99 Mб
Скачать

§ 3. Решение задачи о распространении тепла на бесконечной прямой (задачи Коши)

ина полупрямой

1. Теперь мы можем построить решение задачи Коши. Рассмотрим сначала однородное уравнение

а2ихх = щ,

(36)

и (х, 0) = ср(х),

(37)

где ср (х) — непрерывная функция.

(п. 4) и

Согласно рассуждениям § 1 настоящей главы

формуле (23) решением будет функция

 

СО

 

и(х, t)= ^ G (х — I, t) ф (|) of|.

(38)

— со

 

Надо лишь показать законность вычисления производных ихх и ut путем дифференцирования по этим переменным под знаком интеграла. Мы это сделаем позже.

Формулу (38) можно получить и из наглядных сооб­ ражений. Для этого воспользуемся температурной интер­ претацией задачи. На прямой t 0 возьмем отрезок длины d\, содержащий точку х = |. Количество тепла, выделив­ шегося в момент ^ = 0 на этом отрезке, равно срф (£)

Это количество тепла можно отнести к точке |. Таким

образом, мы будем

иметь точечный источник, в котором

мгновенно в момент

времени t = 0 в точке х = 1 выдели­

лось

количество тепла

dQ= сру (%) d%.

Температура на

бесконечной прямой

для ^ > 0 , обусловленная действием

этого

источника, равна

 

 

 

§ G ( x - l ,

0 = <p(|)G (*-g,

t)dl.

И так для каждого отрезка длины d% прямой ^ = 0. Имея в виду замечание на стр. 138, естественно предположить, что температура, обусловленная действием всех таких от­ резков, т. е. обусловленная заданием начальной темпера­ туры и(х, 0) = ф(х), будет равна

СО

 

и{х, t)= $ ф (l)G( x - l, t)d\.

(39)

— СО

140

Если это верно, то функция (39) и будет решением задачи Коши (36)—(37). Чтобы убедиться в справедливости по­ следнего, достаточно доказать, что функция (39) удовлет­

воряет уравнению (36) для всех

о о < л ; < о о и t > О,

а также начальному условию (37).

 

Проверим сначала условие (37). Согласно формуле (39)

и учитывая

также, что G(x £,

0) = 8 (х — £), имеем

 

00

 

СО

 

и(х, 0 )=

$

<p(£)G(x —£, 0)dl=

\

cp(£)8 (x -£ )d £ = cp(x).

СО

 

— со

 

Последний интеграл равен ср (х) согласно основному свой­ ству 8-функции. Таким образом, функция (39) действи­ тельно удовлетворяет условию (37).

Чтобы установить, что функция (39) является реше­ нием уравнения (36), достаточно доказать, что эту функ­ цию молено дифференцировать по л: (дважды) и по t под знаком интеграла. Действительно, если

 

СО

ихх =

$ 4>{%)Gxx( x - l , t)cll

СО

И

со

И/= $ 4>(l)Gt ( x - t , t)dl,

СО

Т О

оо

агихх — щ = \ ср (£) {a2Gxx - G,} dl = О,

ОО

так

как

функция

G(x —£,

t) является

решением уравне­

ния

(36).

 

 

 

 

 

 

Очевидно, достаточно показать, что интеграл (39) схо­

дится,

а

интегралы

 

 

 

 

 

СО

 

 

СО

 

 

 

со

 

^

Ф (s) Gt dc,

$

ф (ё) G, dl

и

$ ф (ё) Gxx dl (**)

 

— СО

 

 

СО

 

 

— СО

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

в

области

ВЕ= {— о о -< х < о о ;

t ^ e } с произвольным е >

0 .

 

 

Для упрощения выкладок будем предполагать при этом,

что ф (х)

ограничена,

т. е.

[ ф (х) | «с; М .

В интеграле (39)

141

произведем замену

переменной

интегрирования:

а =

— (£ — х)!У~\аЧ. Тогда

(см. формулу (31), стр. 137)

 

 

 

 

СО

 

 

 

и (х,

^) = — .

{

ф (х-\-2ааУ t) е~а2da,

 

 

у

л

J

 

 

 

 

со

 

— со

 

 

 

 

 

 

 

 

 

| и(х, t)\s^ ~=

( 'го (х +

2аа ]/7 ); е~а‘da ^

 

| ' я

J

 

 

 

 

 

 

— т СО

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

у п

{ e - ^ d a = м.

(40)

 

 

 

 

J

 

— СО

Таким образом, интеграл (39) сходится, притом равно­ мерно, в области В1 и \ и \ ^ М . Если предположить до­ полнительно, что ф (х) непрерывна всюду, то из этого

следует также непрерывность функции

(39) в замкнутой

области В1*).

следует, что если

З а м е ч а н и е . Из неравенства (40)

начальные значения фх (.х) и ср2 (х) отличаются меньше чем на 8, т. е. | 9 i (х) — Ф-2(х) | <С 8 для всех х, то соответст­ вующие им решения задачи Коши ut (х, t) и и2 (х, t) также отличаются друг от друга меньше чем на е, т. е.

[ (х, t) « 2 (х, t) | < e.

Таким образом, решение задачи Коши непрерывно зави­ сит от начальных значений.

Рассмотрим теперь первый из интегралов (**):

СО

0 0

$ Ф(1 ) М | = -

 

. 5

^ G ( x - l , t)dl +

 

 

— со

 

— оо

 

 

 

 

 

 

 

со

 

 

 

 

 

+

I

t)dl-

(***)

 

 

 

 

— со

 

 

Первый

интеграл

заменой переменной а =

(£ — x)j\!4аЧ

сводится

к интегралу

 

 

 

СО

 

 

 

 

 

 

\

—г— ф(х + 2 a a V t)e~ a*da.

 

 

J

2

у n t

 

 

*) При дополнительном предположении об ограниченности функ­ ции ф (х). Если ф(х) кусочно-непрерывна, то функция (39) непре­ рывна всюду в В ъ кроме точек прямой / = 0, в которых ф(х) раз­ рывна.

142

Этот интеграл равномерно сходится в области Вг

с произвольным е >

0 , поскольку подынтегральная функ­

ция мажорируется в этой области

функцией

е~а\

интеграл от

которой

 

 

 

2 у ль

 

сходится. Второй интеграл (***) той

 

 

 

 

 

со

а2

 

же заменой

переменной

сводится к

интегралу

f

х

\

 

 

 

 

 

J

У л t

 

 

 

 

 

 

00

 

 

Х Ф ( х - \ - 2 а а

Y t ) e ~ al

d a .

Этот интеграл равномерно схо­

дится в области Вг с произвольным е >. О, поскольку под­ интегральная функция мажорируется в этой области функ­

цией Му:-~ а 2е ~ а \ интеграл от которой сходится. Анало­ ев я

гично поступаем с третьим интегралом (**). Таким образом, мы доказали, что формула (39) действительно дает реше­ ние задачи Коши (36)—(37). Этот результат верен и для функций ф ( х ) , неограниченно возрастающих при х-»-оо, например для таких, для которых существуют постоян­

ные b и N такие,

что | ф (*) |

N e bx.

 

 

2 .

Аналогично строятся решения задач:

и (0,

/) = 0,

a)

a2uxx = Ui,

и(х,

0) = ф(х)

(0 ^ * < о о ),

 

и(х,

 

00

Ф (£)<?*(*> 1; t ) d l \

 

 

 

/) =

$

 

(41)

 

 

 

О

 

 

 

 

 

б) а2ихх = щ, и(ху 0) =

ср(х) (0 ^ * < о о ) , их (0 , /) = 0 ,

 

 

 

00

 

 

 

 

 

 

и(х,

/) =

$ G** (х,

1\ Оф ‘X)dl-

 

(42)

 

 

 

0

 

 

 

 

 

 

Функция Грина задачи

 

 

 

 

 

a2uxx-{-f(x,

t) = u(,

(43)

 

 

и (х,

0) = ф (х),

 

(44)

 

 

 

и (0 ,

0 = Ф(0 ,

 

(45)

 

 

 

 

х, / 3 =0,

 

 

определяется следующим образом.

Функцией Грина G* (х, t) задачи (43)—(45) называется решение задачи

a2G%x= Gf,

G* (х, I; 0) = 6 (* -£ ),

G*(0, |; 0 = 0, I, х, t > О,

143

непрерывное всюду в замкнутой области

=^ЗгО},

кроме точки (|, 0).

Для нахождения

этой функции

Грина воспользуемся

следующим свойством решения задачи Коши.

Для решения задачи Коши

 

 

 

 

 

а2ихх = щ,

 

(46)

 

 

 

и(х, 0) = ф(х)

(47)

с нечетной функцией ф (х) выполняется тождество

 

 

и (0 , 0 = 0 ,

f> 0 .

Действительно,

решение задачи

(46)—(47) имеет вид

 

 

 

СО

 

 

и ( х ,

0 =

\ G ( x — l ,

0

ф (Ю dl,

где G (х — |, 0

________р - (х -| ) У (4 аЧ)

Полагая здесь л" = 0,

VЫаЧ

 

 

 

 

получим

00

,) = 7 Ш S- 5,л,я,ф © ^ = 0 ,

так как подынтегральное выражение есть нечетная функция. Решение задачи Коши

a?G$x= Gf,

G*(x, g; 0) = 6 ( х - Ъ ) - 6 ( х + 1), l > 0 ,

непрерывное всюду в области Ви кроме точек (— |, 0) и (I, 0), согласно упомянутому свойству и будет искомой функцией Грина. Очевидно,

G*

(x,l\ t) = :

1

(

(■* —1

— ехр

(* + £ )2 ]

 

У Ы а Ч

Г ХР

АаЧ

 

4 аЧ

4.

Функция

Грина задачи

 

 

 

 

a2uxx + f(x, t) = uf,

(48)

 

 

 

и {х,

0) = ф (х),

 

(49)

 

 

 

их (0,

0 = Ф(0

 

(50)

определяется следующим образом, 144

Функцией Грина G** (х, g;

t) задачи (48)—(50) назы­

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

 

a2G*x = Gf,

G**(x, g; 0) =

6 ( x - g ) ,

GT (0 , g; 0 = 0 ,

x, g, t > 0 ,

непрерывное всюду в замкнутой

области Вг, кроме точки

&0 ).

Для нахождения ее воспользуемся следующим свойст­ вом решения задачи Коши (46)—(47).

Для

решения

задачи

Коши

(46)—(47) с четной

функцией

ср (х) выполняется тождество

 

их (0 , 0 = 0 ,

* > 0 .

Действительно,

дифференцируя функцию

 

 

СО

 

 

 

и (х,

0 = 5

G(x —g,

0 ф (g) dl

СО

по переменной х и полагая затем х = 0 , получим

со

“'< 0'

так как подинтегральное выражение есть нечетная функ­ ция.

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

 

G**{x, g,

aaG £=*G r,

 

 

 

 

0) = 6 (x -g ) + 6 (x + g),

g > 0 ,

 

непрерывное

всюду в области

Blt

кроме точек (— g, 0) и

(g, 0), согласно упомянутому

свойству и будет искомой

функцией Грина. Очевидно,

 

 

 

 

G** (х,

^

 

ехр

Ш

+ ехр

О Д Н О

 

 

У аЧ

 

ш

J

5.

Доказательство того,

что функции (41) и (42) яв­

ляются решениями задач а) и б), проводится совершенно

аналогично предыдущему.

 

 

 

 

З а м е ч а н и е

1, Из формулы (39) следует, что тепло

распространяется вдоль стержня мгновенно. Действительно,

пусть начальная температура ф (х) положительна на конеч­

ном

отрезке

(хх, х2) бесконечного стержня

и равна

нулю

вне

(хг,

х2).

Тогда температура

произвольной точки х

145

стержня равна

и(х, t)= \ (p(l)G(x-l, t) d t

Xi

Очевидно, при сколь угодно малых t > О эта функция положительна для любого х. К такому выводу мы пришли вследствие неточности физических предпосылок, которыми мы пользовались при постановке задачи Коши (например, при написании уравнения (36)).

З а м е ч а н и е 2. Полученную формулу (39) можно рассматривать как свертку (см. Дополнение) фундаменталь­ ного решения

Glx, Л = _ Д ==е-^/(4аЧ)

К4яа2^

сначальной функцией ср(х), т. е.

и(х, t) = G (х, t) * ср (х).

Если в этой формуле в качестве начальной функции ф (х)

брать произвольную финитную обобщенную функцию, то

и (х, t) будет также решением. задачи Коши. 6 . Рассмотрим несколько примеров.

П р и м е р 1. Решить задачу Коши

 

 

2

U(,

/ m

0) =

/ ч ( “1. * < 0,

 

a 2u xx —

и(х,

ер (х ) =

<

х :> 0.

По формуле (38)

 

 

 

 

(. и2,

 

 

 

 

 

 

и (х,

t) =

 

 

 

 

 

 

 

-

J( p©G( JC- E,

t ) d l = Ul I

G ( x - l ,

+

G ( x - 1 , t)dl.

Производя замену переменной ос = (х —£)/|/4а2^, получим

 

 

x)YАа.4

 

 

— со

 

 

и (х,

t) = — —

f

е~“3d a ------

 

,1

 

 

 

V л

J

 

V Л

 

 

 

 

СО

 

 

x/YAaH

 

 

 

 

 

/ 0

xfV^a2i\

f

2

 

 

 

« i f f ,

f

 

 

 

 

 

«3 /

f

 

 

 

V n [ ) ^

}

) ' |Лт

\

 

 

 

 

 

 

 

'x/Y4.

 

 

 

 

 

 

 

«l +

«2

Ф

 

 

 

 

 

 

%

'

 

 

 

 

 

 

V ЬаЧ

m

П р и м е р 2. Решить задачу Коши

“\ x x = ut, и(х,0) = Ае~х\

По формуле (38) имеем

u( x, t ) = A ] e - ^ G (x —l, t)dl-

00

Произведя замену переменной а = (| —х )//4 а + получим

и (х, t) = - 4 -

X

Уп

3

е- ( х + 2а а Г ( У е- а * d a _

Ле

- A x a a V t

— (4а П +

1 ) а 2 j a = = :

 

Vп

 

 

 

 

— со

 

 

 

 

 

 

 

 

 

lx*a*t

со

f

2axY t ,

r^-r~.—

,

 

 

 

Лрх г --------------

(•

■■■■--------- НУ 1+4 0 ^0

) d a _

 

 

 

А

, е 1 + 4 « ^

^ е \ K l + 4 аЧ

 

 

 

 

)

Я

 

 

 

 

 

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

фор-

муле

 

Qcix~\ft

______

получим

 

 

 

у--------—+

У~1 + 4 аЧ « = р,

 

 

 

 

У

1 + 4а2/

 

 

 

 

 

 

 

и (х,

t) -

j L«-*V(i+ 4o*/)___У

 

 

Р2£ф =

 

 

 

 

У я

 

V 1 + 4 а-И

 

 

 

 

 

 

 

 

 

 

 

А

-*»/(!+4а»0

 

 

 

 

 

 

 

/1 + 4 а2*

 

 

7. Обратимся к неоднородному уравнению. Решение задачи Коши

a2uxx + f(x,t) = uh

(51)

и (х, 0) = ф (*)

(52)

будем искать в виде суммы двух функций и (х, t) = v (х, t) + -tw(x, t), являющихся решениями следующих задач:

v. a2vxx = vt,

v(x, 0 ) =

ф (х );

w: a2wxx + f(x,

t) = wt,

(51д)

w (x,

0) = 0 .

(53)

Функция v(x, i) дается формулой (39). Решение задачи (51j), (53) можно получить методом, описанным в гл. III,

147

§ 6 , п.1. Это решение равно интегралу t

w {х, t) = ^Ul (х, t, г) d%,

о

где Щ (х, t, х) есть решение однородной задачи Коши

аЧЦхх = Щь Щ\,_х= Пх, т).

Согласно § 1

настоящей

главы оно

равно

 

Щ (X, t, т) =

СО

G(_x— l , t

— T)f (|, т) d\.

 

5

 

Следовательно,

 

■—СО

 

 

СО

 

 

 

 

t

 

 

 

w{x,

t) — ^

^

G (х — £, t — т )/(|, т) d%dr.

(54)

0 — со

Эту функцию можно также получить, если воспользоваться температурной интерпретацией функции Грина и урав­ нения (5 К).

В этом уравнении cpf(x, t) есть плотность тепловых

источников

в

единицу времени *). Следовательно,

на

отрезке

длины

d\,

содержащем точку £, за промежуток

времени

(т,

т +

с/т)

выделится

количество тепла dQ—

— cp/(£,

x)d%dx. Если это количество тепла считать выде­

лившимся в

точке

| мгновенно

в момент времени т,

то

температура, обусловленная действием этого источника, будет равна

=x ) G ( x - l , t — x)dldx.

Имея в виду замечание на стр. 138, естественно предполо­ жить, что температура, обусловленная действием всех таких источников (распределенных по всей прямой) в тече­ ние промежутка времени от 0 до /, будет равна

t

СО

 

w(x, 0 = $

5 f{l,x)G(x — \ , t — x)d%dx.

О — оо

 

Если, в частности, тепловой источник действует лишь

в точке | 0, но изменяется со временем,

то функция / (х, t)

в уравнении (51) будет иметь вид

 

/(X, 0 = /( О б ( * - ы .

 

*) См. вывод уравнения теплопроводности

(стр. 30).

148

Тогда решение задачи (51х), (53) с такой неоднородностью будет иметь вид

t

СО

 

w( x, *) = $

S бЙ - 1 о )/(т )0( дг -|,

/-T)d§rfT =

О — со

t

 

 

 

=

\ f (т) G (x —g0, t - x )c h .

 

 

0

Мы здесь воспользовались свойством б-функции.

З а м е ч а н и е 3.

Решение задачи

Коши

для неодно­

родного уравнения

с нулевыми

начальными

значениями

а2ихх + /(*,/) = «,,

и (х,

0) = 0

 

также можно записать в виде свертки (по двум перемен­ ным!) фундаментального решения G(x, t) с функцией f(x, t):

t

00

и(х, t) = G(x, t)*f.(x, 0 = $

$ G(x — l, t — x)f (g, x)d\dx.

0 — со

З а м е ч а н и е 4. Решение задачи (5К), (53) на полу­ прямой с краевым условием w (0 , t) = 0 (или wx (0 , i) = 0) строится аналогично и дается формулой

 

t СО

w{x, t) = \ \ f (1, т) G* (х, l , t —т) di dx

 

о о

t

о о

или оу= §

$ fG* * d\ dx

оо

§4. Решение задачи о распространении тепла в трехмерном (двумерном) пространстве

Теперь обратимся к рассмотрению задачи Коши для уравнения теплопроводности в двумерном и трехмерном пространствах.

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

a2Au = ut.

(55)

О п р е д е л е н и е . Функцией Грина G(M,M0; i) задачи Коши для уравнения (55) называется такое его решение, которое:

а) удовлетворяет начальному условию

и(М, 0) = 8 (Л4, М0);

(56)

149

Соседние файлы в папке книги из ГПНТБ