книги из ГПНТБ / Арсенин, В. Я. Методы математической физики и специальные функции учеб. пособие
.pdf§ 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** (х, |
^ |
|
ехр |
Ш |
+ ехр |
О Д Н О |
||
|
|
У 4лаЧ |
|
ш |
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
