u |
|
|
|
|
u |
|
||||
|
|
|
a |
|
2 |
|
|
; |
||
|
|
|
2 |
|
|
|
||||
|
|
|
|
|
|
|
|
|||
t |
|
|
|
x |
2 |
|
||||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
||
|
u |
|
|
|
f |
x , 0 x ; |
||||
|
t |
0 |
||||||||
|
|
1 |
|
|
|
|
||||
|
u |
|
|
|
|
|
|
|
||
|
|
|
0, 0 t . |
|||||||
x |
x 0 |
|||||||||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
158
Эта ситуация соответствует случаю, когда левый конец стержня теплоизолирован.
Доопределим функцию f1 x |
чѐтным образом относи- |
|||
тельно начала координат, введя новую функцию |
F2 x : |
|||
f |
x , если x 0, |
|
||
|
1 |
|
|
|
F2 x |
f |
x , если x |
0. |
|
|
|
|||
1 |
|
|
|
|
|
|
|
|
|
Записывая формулу Пуассона и совершая преобразования, аналогичные предыдущим, получим:
|
|
|
u x, t |
|
|
1 |
|
|
|
|
|
|
e |
x 2 |
|
||||
|
|
|
|
|
|
|
|
|
F |
4a2t |
d |
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
||
|
|
|
|
|
|
2a |
t 0 |
x |
|
|
|
|
|||||||
|
|
|
|
|
x |
2 |
|
|
2 |
|
|
||||||||
|
1 |
|
|
1 |
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
e |
|
4a |
t |
|
e |
4a |
t |
|
d . |
|
||||||
|
|
|
f |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
2a |
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
159
Проверим, что эта функция удовлетворяет граничному усло-
вию. Найдѐм производную |
u |
: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
x |
|
|
|||||
|
u |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
1 |
|
|
1 |
|
|
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|||
|
|
|
|
|
|
|
|
x e |
|
|
|
|
|
|
x e |
|
|
|
d . |
|||||||
|
|
|
|
|
|
f |
|
|
|
4a |
t |
4a |
t |
|
||||||||||||
|
x |
|
4a |
|
t |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
3 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Отсюда видно, что если |
x 0, |
то |
u |
|
x 0 |
0. |
Согласно форму- |
|||||||||||||||||||
|
||||||||||||||||||||||||||
x |
|
|||||||||||||||||||||||||
|
159 |
граничное условие u |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
ле |
x 0 |
0 будет выполняться, ес- |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|||
ли мгновенные точечные источники тепла одинаковых знаков
96
разместить симметрично относительно
Поэтому в силу чѐтности функции |
f1 x |
начала координат.
производная |
u |
в |
|
x |
|||
|
|
нуле всегда равна нулю.
Итак, решение начально-краевой задачи даѐтся форму-
лой 159 .
3.5. Уравнения эллиптического типа
Задачи, приводящие к эллиптическим уравнениям. Постановка краевых задач.
При изучении стационарных явлений разной физической направленности (колебания, теплопроводность, диффузия и др.), как правило, получают уравнения эллиптического класса. Самое известное уравнение данного класса – уравне-
ние Лапласа u 0.
Зависимость |
u |
именуется гармонической в области |
T |
при условии еѐ непрерывности в данной области одновременно со своими производными до второго порядка и выполнения для них уравнения Лапласа.
При исследовании свойств гармонических зависимостей получены разные математические способы, которые можно с успехом использовать и для уравнени гиперболического и параболического классов.
Рассмотрим стационарное тепловое поле. При изучении уравнений параболического типа получено, что для температуры нестационарного теплового поля справедливо дифференциальное уравнение теплопроводности
|
|
k |
|
ut a2 u a2 |
|
|
. |
|
|||
|
|
c |
|
При условии, что явление стационарно, устанавливается распределение температуры u x, y, z , постоянное во времени и, значит, для него выполняется уравнение Лапласа
97
u 0.
160
Если присутствуют источники тепла, то приходим к уравнению
|
|
u f , f |
F |
, |
|
|
|
k |
|
||
|
|
|
|
|
|
здесь |
F |
– плотность источников тепла, |
k |
||
теплопроводности. Неоднородное уравнение
161
– коэффициент
Лапласа |
161 |
обычно именуют уравнением Пуассона.
Выделим объѐм |
T , |
который ограничен некоторой по- |
верхностью . |
Задача стационарного распределения темпера- |
туры u x, y, z |
внутри данного объѐма T ставится так: опре- |
делить зависимость u x, y, z , для которой внутри T выпол- |
|
няется уравнение
|
|
|
|
|
u f x, y, z |
|
162 |
||
и граничное условие одного из трѐх типов: |
|
|
|||||||
|
а) |
u f1 |
на |
(первая краевая задача), |
|
|
|||
|
б) |
u |
|
f2 на |
(вторая краевая задача), |
|
|||
|
n |
|
|||||||
|
|
|
|
|
|
|
|||
|
в) |
u |
h |
u f3 0 на (третья краевая задача), |
|||||
|
n |
||||||||
|
|
|
|
|
|
|
|
||
здесь |
f1 |
, f2 , |
f3 |
, h – известные зависимости, |
u |
– производ- |
|||
n |
|||||||||
|
|
|
|
|
|
|
|
||
ная по внешней нормали к поверхности .
Замечание. Ясно, что стационарное распределение температуры устанавливается только тогда, когда суммарный поток тепла сквозь границу области нулевой. Таким образом,
для зависимости |
f2 |
необходимо выполнение ещѐ одного |
условия: f2d 0.
98
Первую краевую задачу для уравнения Лапласа обычно именуют задачей Дирихле, а вторую краевую задачу – задачей
Неймана. |
|
|
При условии, что находится решение в области |
T0 , |
|
внутренней (или внешней) относительно поверхности |
, |
со- |
ответственная задача называется внутренней (или внешней)
краевой задачей.
Задача Дирихле в круге.
Предположим, что имеется круг с радиусом R и центром, являющимся полюсом O полярной системы координат. Найдѐм зависимость u r, , которая является гармониче-
ской в круге и на его окружности для неѐ справедливо требо-
вание
u |
r R |
f , |
|
|
f
– известная зависимость, являю-
щаяся непрерывной на границе круга. Для неизвестной зависимости требуется выполнения в круге уравнения Лапласа
r |
2 2u |
r |
u |
|
2u |
0. |
163 |
r2 |
r |
2 |
|
||||
|
|
|
|
|
Для нахождения решения данной задачи воспользуемся методом разделения переменных. Предположим, что решение уравнения, удовлетворяющее заданному граничному усло-
вию, находится в форме u Q r T t . |
В этом случае имеем |
|||||||||||||||||||||||||||||||||
|
|
|
|
r |
|
Q |
|
r T |
rQ r |
T Q r T |
|
|
0. |
|||||||||||||||||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Производим разделение аргументов: |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
Q |
|
r |
|
T |
|
|
T |
|
θ |
|
|
|
r2Q |
|
r |
|
|
rQ |
|
r |
|
|
; |
|||||||
|
|
|
|
|
|
|
|
|
T |
|
|
|
r |
|
Q |
|
r rQ |
r |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
Q r |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Полагая обе части |
последнего соотношения равными кон- |
|||||||||||||||||||||||||||||||||
станте k |
2 |
, |
|
будем иметь пару обыкновенных дифференци- |
||||||||||||||||||||||||||||||
|
|
|||||||||||||||||||||||||||||||||
альных уравнений: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
T |
|
|
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|||||||
|
k T |
0, r Q |
|
r rQ |
|
r k Q r 0. |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
99 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Следовательно, для |
k 0 получим |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
T A B , |
|
|
|
|
|
|
164 |
|||||||
|
|
|
|
Q |
|
r |
|
C D ln r. |
|
|
|
|
|
165 |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Если же |
k 0, |
то |
|
|
|
|
|
|
|
|
|
|
166 |
||||||
|
|
|
|
|
Acos k B sin k , |
|
|
|
|
|
|||||||||
|
|
|
T |
|
|
|
|
|
|
|
|
|
|||||||
а решение |
второго |
уравнения будем определять в форме |
|||||||||||||||||
Q r r |
|
, |
отсюда |
имеем |
r |
m m 1 r |
m 2 |
rmr |
m 1 |
k |
|
r |
|
0 |
|||||
|
m |
|
|
|
|
|
|
|
|
2 |
|
|
|
2 |
|
m |
|
||
или |
r |
|
m |
|
k |
|
0, |
||
|
m |
2 |
2 |
||||||
|
|
|
|
|
|
|
|||
отсюда |
m k. |
|
|
|
|
|
|
|
|
Таким образом, |
|
|
|
|
|
|
|
|
|
|
Q r Cr |
k |
Dr |
k |
. |
|
|
||
|
|
|
|
|
|||||
167
Замечание. Зависимость u r, |
по аргументу явля- |
||
ется периодической с периодом |
2 , |
поскольку для одно- |
|
значной функции величины u r, и |
u r, 2 совпадают. |
||
Следовательно, из соотношения 163 |
вытекает, что B 0, а в |
||
соотношении
1, 2, ..., k 0 .
166
k
может принимать одно из значений
Далее, в соотношениях
165
и |
167 |
необходимо вы-
полнение условия D 0, иначе зависимость будет терпеть
разрыв в точке |
r |
Таким образом, |
мы |
дифференциального
0 и не являться гармонической в круге. имеем бесконечное семейство решений
уравнения 163 , удовлетворяющих за-
данному граничному условию, непрерывных в круге, которые (несколько изменив обозначения) можно представить в форме
100