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

Учебное пособие 1920

.pdf
Скачиваний:
14
Добавлен:
30.04.2022
Размер:
2.9 Mб
Скачать

Лемма. Если функция (x) нечетная, то решение

u(t, x) удовлетворяет условию u(t,x)

x 0

0, если же функция

(x) четная, то u(t, x)

удовлетворяет условию

u(t.x)

 

x 0

0.

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В самом деле, если функция (x) нечетная, то

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

u(t,x)

 

 

 

 

 

(z)e 2a

t dz 0

 

 

 

 

 

 

 

 

 

 

 

x 0

 

2a

t

 

 

 

 

 

 

 

 

в силу нечетности подынтегральной функции. Если функция (x) четная, то

u(t,x)

x

u(t,x)

x

так как функция

 

1

 

 

 

1

 

 

 

 

 

 

(x z)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x z) (z)e 2a

2t dz,

 

 

 

 

2

 

 

2a t

 

 

 

a t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 0 2a31t t z (z)e

 

 

dz 0,

2a2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x (x) четная. Чтобы получить решения задач

теплопроводности полубесконечного стержня, построим продолжение функции (x) на всю числовую ось Ох нечетным или четным образом, то есть построим функции:

 

(x),

если x 0;

 

 

 

 

 

1x)

 

 

если x 0.

 

( x),

 

 

 

 

 

 

(x),

если x 0;

2

 

 

 

 

(x)

 

 

если x 0.

 

( x),

 

 

 

 

 

1( х) 1(х), 2( х) 2(х).

Решение задачи теплопроводности для бесконечного

стержня с начальным условием u(t,x)

t 0

1(x)

 

 

(х z)2

 

u2 (t,0) 2 (z)e

4a2t

dz

 

60

удовлетворяет при х 0 основному дифференциальному уравнению, начальному условию u(t,x) t 0 (x) и гранично-

му условию u(t,x) х 0 0, то есть функция u1(t,x), рассматри-

ваемая только при

х 0,

 

является решением. Точно так же

функция u

 

 

1

 

 

 

 

 

(х z)2

dz, рассматриваемая

 

(t,x)

 

 

 

 

(z)e 4a2t

 

 

 

 

 

 

2

 

2a t

 

 

2

 

 

 

только при х 0, является решением основного дифференциального уравнения и удовлетворяет начальному

u(t,x)

 

t 0 (x)

и граничному условию

u(t.x)

 

x 0

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

Решения u2 (t,x) можно преобразовать:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

(x z)2

 

 

 

 

 

 

 

 

 

 

 

(x z)2

 

 

 

 

 

 

 

 

u1(t,x)

 

 

 

 

 

 

 

 

1(z)e 4a2t

 

dz 1(z)e

4a2t dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

(x z1 )2

 

 

 

 

 

 

 

 

 

 

 

(x z)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( z )e

 

4a2t dz (z)e 4a2t

 

dz ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a t

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

0

 

 

 

 

 

 

(x z1 )2

 

 

 

 

 

 

 

 

 

 

(x z)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4a2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4a2t

 

 

 

 

 

 

u

2

(t,x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( z

 

)e

 

 

 

 

dz (z)e

 

 

 

 

 

 

 

 

dz .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a t

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если в первых интегралах сделать замену переменной

z1 z, то получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x z)2

 

 

 

 

 

 

 

 

(x z)2

 

 

 

 

u1(t,x)

 

 

 

 

 

 

(z)e 4a2t

 

 

dz (z)e 4a2t

dz ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x z)2

 

 

 

 

 

 

 

 

 

 

 

(x z)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4a2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4a

2t

 

 

 

 

 

 

 

u2(t,x)

 

 

 

 

 

 

 

 

 

 

 

 

 

(z)e

 

 

 

 

 

 

 

 

dz

(z)e

 

 

 

 

 

 

 

dz

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

(x z)2

 

 

 

 

 

(x z)2

 

 

 

 

 

 

 

 

u1(t,x)

 

 

 

 

(z) e 4a2t

e 4a2t

 

 

dz,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

 

 

 

t 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

61

 

 

1

 

 

 

 

 

(x z)2

 

 

(x z)2

 

u2(t,x)

 

 

 

 

(z) e

 

4a2t e

 

4a2t dz.

 

 

 

 

 

 

 

2a

 

t

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

1.13.Решение задачи теплопроводности

вконечном стержне

Пусть в стержне длины l задано начальное распределение температур. Требуется найти температуру точек стержня при t>0. Рассмотрим случай, когда на концах поддерживается нулевая температура, т. е. рассмотрим следующую задачу. Найти функцию u(t,x), удовлетворяющую при 0<x<l, t>0 дифференциальному уравнению

 

u(t,x)

2

2u(t,x)

 

 

 

 

 

a

 

 

,

 

 

 

t2

 

t

 

 

начальному условию u(t,x)

 

x 0

(x) и граничным условиям

 

u(t,0)=u(t,l)=0.

Найдем сначала ненулевые решения u(t,x) нашего уравнения, удовлетворяющие только граничному условию. Эти решения будем искать в виде u(t,x)=T(t)X(x), где X(0)=X(l)=0. Дифференцируя функцию по переменным t и x дважды, и подставляя результаты дифференцирования в основное уравнение, и поделив его на T(t)X(x), получим

T'(t) X''(x) . a2T(t) X(x)

Отсюда следует, что

T'(t) a2T(t) 0; X (x) X (x) 0; X(0) X(l) 0.

Задача

имеет ненулевые решения только

при

λ=λk=

k2

2

(k 1,2,...). Решения имеют вид Xk(x)=sin

k x

.

l

2

 

l

 

 

 

 

Тогда из уравнения для T(t) получим

62

 

 

(ak )2 t

 

 

 

T (t) A e l2 .

k

k

Решение поставленной задачи будем искать в виде ряда

 

(ak )2 t

u(t,x) Ake

 

 

sin

k x

,

 

l2

 

 

k 1

 

 

 

l

коэффициенты Ak которого подбираются таким образом, что-

бы выполнялось начальное условие. При t=0, получим

 

n x

 

 

2

l

 

n x

 

(x) An sin

l; An

 

0

(x)sin

dx.

 

l

 

n 1

l

 

 

l

Рассмотрим теперь случай, когда на концах отрезка [0,1] поддерживается постоянная температура и . В этом случае требуется найти решение основного уравнения, удовлетворяющее начальному условию и граничным условиям u(t,0) =u0, u(t,l)=u1.

Эта задача сводится к предыдущей с помощью замены неизвестной функции

 

v(t,x) u(t,x) [u

u1 u0

x].

 

 

 

 

 

 

 

 

 

0

 

l

 

 

 

 

 

2v

 

2u

 

 

Так как

v

 

u

 

,

то функция v(t,x) является

 

 

 

,

 

 

 

 

t

 

t

t2

 

x2

 

 

 

 

решением. Функция v(x) удовлетворяет условию

v(t,x)

 

t 0

(x) [ u ,

u1 u0

x] (x).

 

 

 

 

 

 

 

 

 

 

0

 

 

 

l

 

 

 

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k x

 

 

 

 

2

l

 

k x

 

v(0,x)

Ak

sin

.

Ak

 

(x)sin

dx.

 

l

 

k 1

 

 

l

 

 

 

0

 

l

Тогда

63

 

 

u1 u0

 

(ak )2 t

u(t,x) u0

 

x Ake

l2

l

 

 

k 1

 

 

sin k x, l

где Ak вычисляется по вышеприведенным формулам. Аналогично решаются задачи теплопроводности конечного стержня при граничных условиях других типов.

1.14.Решение задачи теплопроводности

воднородном шаре

Пусть в шаре x2 +y2задано начальное распределение температур, и пусть в точках сферы x2 +y2 =R2 поддерживается нулевая температура. Мы рассмотрим частный случай этой задачи, когда начальная температура зависит лишь от расстояния точки от центра шара и не зависит от угловых координат этой точки. Поэтому и при t>0 температура в точках шара зависит лишь от расстояния этой точки от центра шара. По этой причине задачу целесообразно решать в сферической системе координат ρ, θ, φ. Как известно, оператор Лапласа в сферической системе координат имеет вид

u

2u

 

1 2u

 

2 u

ctg u

1

 

 

2u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

2

2

 

2

 

 

2

 

2 sin2

 

 

 

 

 

 

 

 

 

2.

Рассматриваемая задача заключается в следующем. Найти функцию u(t,ρ), удовлетворяющую при t>0

0 ≤ < уравнению

u(t, )

a2[

2u(t, )

 

2

 

u(t, )

]

t

2

 

 

 

 

 

 

начальному условию u(t, )t 0 g( ) и граничным условиям

u(t,0)≠ ∞ , ( , ) = 0. Первое условие говорит о том, что температура в центре шара не может быть бесконечной. Для решения этой задачи введем вспомогательную функцию

64

v(t,ρ)=ρu(t,ρ). Тогда u(t,ρ)=v(t,ρ)/ρ,

 

 

u

 

1

 

v

,

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

u

1 v 1

v,

 

2u

 

 

1 2v

 

 

 

 

2 v 2

v,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

2

 

 

 

 

2

3

 

 

 

 

 

 

 

 

 

2u

 

 

2 u

 

 

1 2v

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя результаты дифференцирования в основное уравнение, приходим после сокращения на 1/ρ к следующей задаче.

Найти функцию v(t,ρ), удовлетворяющую при t>0, 0<ρ<R дифференциальному уравнению

 

v(t, )

a2

2v(t, )

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

2

 

 

 

 

 

 

 

 

начальному условию v(0,ρ)=ρg(ρ)

и

 

 

граничному условию

v(t,0)=v(t,R)=0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение этой задачи имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ak )2 t

 

 

 

 

k

 

 

v(t, ) Ake

 

 

 

l2

 

sin

,

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

R

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 R

 

 

 

 

k

 

 

 

 

 

 

Ak

 

 

g( )sin

 

 

 

 

 

d .

R

 

R

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, решение задачи теплопроводности в од-

нородном шаре имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ak )2t

 

 

 

 

 

 

k

 

 

 

1

 

 

 

 

 

 

 

 

u(t, )

Ake l2

 

 

 

sin

,

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

R

где Ak вычисляется по формуле, приведенной выше.

65

1.15. Задачи, приводящие к уравнению Лапласа

Рассмотрим задачу об определении установившейся температуры в однородном твердом теле. Ранее было установлено, что распространение тепла в изотропном однородном теле (при отсутствии источников тепла) имеет вид

u

 

 

 

 

2

u

 

 

2

u

 

 

2

u

 

a

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

2

.

t

 

 

x

 

y

 

z

 

 

 

 

 

 

 

 

 

 

Пусть температура в каждой точке (х,у,z) тела установилась (она не меняется с течением времени). Тогда производная по времени равна нулю, и уравнение будет иметь следующий вид:

2u

 

2u

 

2u

0.

x2

y2

z2

 

 

 

Для определения u(х,y,z) не надо задавать начальное распределение температуры (начальное условие), а достаточно задать одно граничное условие, не зависящее от времени.

Задача определения решения уравнения Лапласа по его значениям на границе нашей области называется задачей Дирихле, а задача определения решения уравнения Лапласа,

удовлетворяющего граничному условию du S 0, называется dn

задачей Неймана.

Потенциальное движение несжимаемой жидкости

Рассмотрим установившееся движение несжимаемой жидкости и движение жидкости невихревое (потенциальное), т. е. вектор v=grad . Плотность ρ постоянна, и из уравнения неразрывности имеем div =0.

2 2 2

Так как div(grad ) 0 или x2 y2 z2 0,

т. е. потенциал скорости удовлетворяет уравнению Лапласа.

66

1.16. Общий вид уравнения эллиптического типа

Уравнение

n

2u

 

n

u

 

a

 

 

b

 

cu f ,

x x

 

x

, 1

 

1

 

 

где aij ,bi ,c , f - функции, заданные в области V, являются в этой области эллиптическим типом, если квадратичная форма

n

a q q, 1ˆ

сохраняет в области V знак и не равна нулю.

Число n будет называться числом измерений области V. Будем рассматривать трехмерные области (n=3), но ре-

зультаты применимы для плоских (n=2) и для многомерных (n>3) областей.

Предполагается, что функции aij ,bi ,c , f непрерывны и функции aij , а также функции

n ai

ei bi 1 x

будут иметь непрерывные первые производные. Тогда наше уравнение можно преобразовать к виду

n

 

 

u

n

u

 

 

a

e

cu f .

x

x

 

, 1

 

1

x

Выражение, стоящее в левой части вышеприведенного уравнения, обозначим через u.

Основное уравнение будет записано в виде u f .

Основные граничные задачи

1. Первая граничная задача (задача Дирихле): u f , ко-

гда x V FV ;

u ,

когда x FV .

67

2. Вторая граничная задача (задача Неймана): u

f ,

когда x V FV;

 

du

u ,

когда x FV,

где

 

 

 

dv

 

 

коэффициент не равен нулю на поверхности FV;

a d - это дифференцирование по направлению нор- dv

мали к FV.

3. Смешанная (третья) граничная задача: u f , когда

x V FV;

 

du

u ψ,

когда x FV ,

 

 

 

dv

 

где коэффициент , не обращается на всей поверхности FV в нуль тождественно, но равен нулю на части FV.

Эти граничные задачи будем называть внутренними или внешними, если они ставятся внутри или вне конечной замкнутой поверхности FV.

Сопряженные граничные задачи

Рассмотрим граничное условие задачи Неймана:

 

du

u ,

когда x FV

( 0).

 

 

dv

 

 

Его можно привести к виду

 

du

gu ,

когда x FV,

 

 

dv

 

где g и - известные функции.

Запишем его в виде

Pu , когда x FV.

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

Qu a

du

(g b)u ,

когда x FV,

 

 

dv

 

где дифференциальное выражение Q определено вторым из равенств, a - некоторая функция, которая определена на FV, назовем сопряженным граничному условию.

Назовем далее сопряженными граничные задачи:

68

u f , когда x V FV ;

Pu ,

когда

x FV;

 

N f ,

когда x V FV ;

~

когда

x FV,

 

Q ,

 

где uи N -

сопряженные дифференциальные выражения, и

Pu и Q - соответствующие

им выражения;

~

~

f , f и

, -

функции, определенные соответственно в изучаемой области V

и на ее границе FV.

 

 

 

 

Для граничной задачи Дирихле:

 

 

 

u f , когда x V FV ;

u ,

когда x FV .

 

Тогда сопряженной к ней назовем задачу:

 

 

N f ,

когда x V FV ;

~

когда

x FV .

 

,

 

Если функции f и тождественно равны нулю, то задачу будем называть однородной. Сопряженная задача называет-

ся однородной, если тождественно равны нулю

~

~

f

и .

1.17. Фундаментальные решения. Функция Грина

Фундаментальным решением уравнения

 

L 0 назы-

вают функцию Леви L( ,x) , которая при x

удовлетворяет

этому уравнению по координатам х или и зависит от коор-

динат другой точки, как от параметров. В выражениях L

или xL индексы ξ, x – это переменные, по которым произво-

дится дифференцирование. Выражения u и xu понимают

как xu(x) и u( ).

Рассмотрим задачу Дирихле: Найти решение уравнения

xu f , когда x V FV ; u(x) , когда x FV , где f и - непрерывные функции.

Пусть решение u(х) этой задачи и функция Леви L( ,x)

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

области V вместе со своими первыми производными. Применим к функции u(х) формулу Грина — Стокса и получим

69