Lumf_p1&2(2010)
.pdfУпражнение. Проверить, что (193) удовлетворяет уравнению (180) и соответствующему начальному условию.
Далее необходимо вернуться к исходной переменной t: = a2t и, подставляя в (193), получим
|
1 |
|
1 |
|
|
|
|
|
(x )2 |
|
||||||
|
|
|
|
|
|
Z |
f( )e |
|
|
|
|
|
||||
|
u(x; t) = |
2ap |
|
4a2t d |
(194) |
|||||||||||
|
t |
|||||||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||
Можно проверить, что функция |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
1 |
|
|
|
|
(x )2 |
|
|
|
||
|
|
|
|
|
|
|
|
|
4a2t |
|
(195) |
|||||
|
|
|
|
|
2ap |
|
||||||||||
|
|
G(x; ; t) = |
t |
e |
|
|
|
|
также является решением исходного уравнения и ее называют
фундаментальным решением уравнения теплопроводности.
Физическим тепловым импульсом называется начальное распре-
101
деление температуры |
(0;0 |
|
xj x0 |
|
0j> ": |
(196) |
|
f"(x) = |
; |
|
|||||
|
u |
x |
x |
< "; |
|
||
|
|
j |
|
|
j |
|
В этом случае решение задачи будет иметь вид
|
u |
x0+" |
(x )2 |
|
|||
|
2ap0 |
|
Z |
e |
|
|
|
u(x; t) = |
|
4a2t d |
(197) |
||||
t |
x0 "
и по теореме о среднем оно может быть записано следующим образом
|
|
|
|
|
~ 2 |
|
|
2"u0 |
|
|
(x ) |
|
|
|
|
4a2t |
(198) |
|||
|
2ap |
|
||||
u(x; t) = |
t |
e |
|
Точечный тепловой импульс соответствует " ! 0. Количество теплоты, переданное стержню, пропорционально произведению
2"u0
102
и при " ! 0 должно оставаться конечным. Полагая
2"u0 = 1
получаем, что u0 ! 1 при " ! 0. Т.о., точечный тепловой импульс может быть записан в виде -функции Дирака:
f(x) = (x x0):
Подставляя записанное в таком виде начальное условие в (194), получаем решение
|
1 |
|
|
(x x0)2 |
|
||
|
|
2 |
|
|
(199) |
||
u(x; t) = |
2ap |
|
e |
4a |
t ; |
||
t |
которое есть фундаментальное решение G(x; ; t) при = x0. Т.о., мы можем утверждать, что функция
|
|
1 |
|
|
|
(x )2 |
|
|
|
|||
|
|
|
|
2 |
(t t0); |
|
|
|||||
G(x; ; t |
|
t0) = |
|
|
|
|
e |
4a |
(200) |
|||
|
|
|||||||||||
|
|
|
|
|
||||||||
|
2ap (t t0) |
|
|
|
|
|
|
103
дает температуру в точке x в момент времени t, если в начальный момент времени t = t0 в точке возникает точечный тепловой импульс. Функция G(x; ; t t0) носит название функции влияния точечного источника для неограниченной области или функции Грина, с ее помощью решение задачи записывается в виде:
|
1 |
|
|
u(x; t) = |
Z |
f( )G(x; ; t) d |
(201) |
|
1 |
|
|
104
4.3Решение задачи о теплопроводности для конечного отрезка
Рассмотрим задачу о теплопроводности на отрезке: |
|
|
ut = a2uxx + g(x; t); |
(0 < x < l; t > 0) |
(202) |
Начальное условие |
|
(203) |
u(x; 0) = f(x) |
||
и однородные граничные условия |
|
|
u(0; t) = 0; |
u(l; t) = 0: |
(204) |
4.3.1 Однородная задача |
|
Рассмотрим сначала однородную задачу |
|
ut = a2uxx: |
(205) |
Будем искать решение в виде |
|
u(x; t) = X(x)T (t)
105
Подставляя в уравнение, получаем
1 T |
0 |
|
X00 |
||||
|
|
|
|
= |
|
|
= |
a2 T |
|
X |
В результате получаем два обыкновенных ДУ:
X00 + X = 0; T 0 + a2 T = 0:
(206)
(207)
(208)
Из граничных условий для u получаем граничные условия для
X:
X(0) = 0; X(l) = 0:
В результате для функции X(x) мы получили задачу о собственных значениях (задачу Штурма-Лиувилля):
X00 + X = 0; X(0) = 0; X(l) = 0: |
(209) |
Ранее было показано, что собственные значения этой задачи
n = n 2 (210) l
106
соответствующие собственным функциям |
|
|
|
|||
X (x) = sin |
n |
x = sin |
n |
x |
(211) |
|
|
||||||
n |
|
|
l |
|
||
Далее находим функцию T (t): |
|
|
|
|
||
|
|
|
|
|
|
|
Tn(t) = Cne a2 nt |
|
|
(212) |
Таким образом, мы нашли частные решения однородной задачи:
un(x; t) = Cne a2 nt sin |
n |
x |
(213) |
|
|||
|
l |
|
Общее решение нашей задачи запишем как суперпозицию част-
ных |
1 Cne |
n |
2a2t sin lnx |
|
|||
|
|
||||||
u(x; t) = |
l |
(214) |
|||||
|
X |
|
|
|
|
|
|
|
n=1 |
|
|
|
|
|
|
Из начального условия получаем |
|
||||||
f(x) = |
1 C sin nx |
(215) |
|||||
|
|
X |
|
n |
|
|
l |
||
n=1 |
||
|
107
Последнее выражение есть разложение функции f(x) в ряд Фурье по синусам на интервале (0; l). Для нахождения Cn домножим уравнение (215) на sin lmx и проинтегрируем:
l |
|
1 |
|
l |
|||||
|
|
|
|
|
|
|
|
||
Z |
f(x) sin |
m |
x dx = n=1 Cn Z |
sin |
n |
x sin |
m |
x dx (216) |
|
l |
l |
l |
|||||||
0 |
|
|
X |
0 |
|
|
|
|
|
С учетом формулы
sin sin = 12(cos( ) cos( + ))
получим для интеграла в правой части
Z |
l |
lnx sin l x dx = |
2 nml: |
|||||
sin |
||||||||
|
|
|
|
m |
1 |
|||
0 |
|
|
|
|
|
|
|
|
108
В результате для коэффициента Cn имеем
|
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
Cn = l Z |
f( ) sin ln d : |
|
|
|
|||||||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Подставим в решение найденное значение Cn: |
|
|
|
|||||||||||||||||||
|
1 |
2 |
|
2 |
|
l f( ) sin |
n |
d 3e |
n |
2a2t sin |
n |
x |
||||||||||
u(x; t) = |
|
l |
||||||||||||||||||||
|
|
|
|
|||||||||||||||||||
|
n=1 |
|
|
l |
Z |
l |
|
|
|
|
|
|
|
|
|
l |
||||||
|
X |
|
|
|
|
0 |
|
|
|
7 |
|
|
|
|
|
|
|
|
|
|||
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
4 |
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|||
Поменяем порядок суммирования и интегрирования |
|
|
||||||||||||||||||||
|
l |
2 |
2 |
|
|
|
1 e |
n |
2a2t sin |
n |
sin |
n |
x3f( ) d |
|||||||||
u(x; t) = |
Z |
|
|
|
l |
|||||||||||||||||
|
|
|
|
|
||||||||||||||||||
|
|
l |
n=1 |
|
|
|
l |
|
l |
|
|
|
||||||||||
|
0 |
6 |
|
|
X |
|
|
|
|
|
|
|
|
|
7 |
|
|
|||||
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
(217)
(218)
(219)
109
Введем функцию |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
G(x; ; t) = 2l |
1 e |
n |
2a2t sin ln |
sin lnx |
|
||||||
|
l |
(220) |
||||||||||
|
|
|
|
X |
|
|
|
|
|
|
||
|
|
|
|
n=1 |
|
|
|
|
функцию мгновенного точечного источника или функцию температурного влияния мгновенного точечного источника тепла. С ее использованием решение нашей задачи будет иметь вид
u(x; t) = Z |
l |
|
G(x; ; t)f( ) d |
(221) |
0
Покажем, что функция G(x; ; t) представляет собой распределение температуры в стержне в момент времени t, если в начальный момент температура равна нулю и в этот момент в точке x = мгновенно выделяется некоторое количество тепла, при том что на краях стержня поддерживается нулевая температура.
110