Файл: Харитонов, В. В. Некоторые нелинейные задачи теплопроводности.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 30.10.2024

Просмотров: 57

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

ГЛАВА II

ТЕОРИЯ ОБОБЩЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ БЕССЕЛЯ НУЛЕВОГО ПОРЯДКА

Введение

Было показано, что различного рода задачи теплопровод­ ности с граничными условиями третьего рода при наличии теплообмена излучением сводятся к нелинейным дифферента альным уравнениям вида

dr2

4. _ L . cK i rl +

А [Г4 (/■) — Г4] = О

(2- 1)

г .

dr

 

 

 

 

где искомое решение Т(г)

и Т0 больше нуля и А>0,-

 

Поскольку получение точного

решения уравнения (2-1) в

аналитическом виде представляет

существенную

сложность

[23], в практических расчетах

прибегают к

всевозможным

приближениям,

в частности к уже

описанной

замене закона

Стефана—Больцмана законом Ньютона. В результате такой

замены нелинейное

дифференциальное уравнение (2-1)

сво­

дится к линейному уравнению

 

агт{г)

-------- ^

- ± А 11 Т ( г ) - Т 0] = 0 ,

(2- 2)

dr2

г

dr

 

которое решается с помощью известных функций Бесселя нуле­ вого порядка [24, 33, 34].

Полученная таким образом оценка, естественно, дает не 'достаточно точное представление о реальных процессах тепло­ проводности и оказывается тем грубее, чем больше разность Т — Т0. Поэтому как использование указанного приема полу­ чения оценки решения (2-1), так и использование других при­ ближенных решений не обеспечивает во многих случаях тре­ буемой точности. Возникает задача — найти пути к точному

решению уравнения (2-1).

_

 

Заменой переменных х = г У А (2-1) приводится к виду

.

X

+ Г4 (х) = ± Т\,

(2-3)

dx2

dx

 

29


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

и" (х) — — и' (х) ± и"1(х) = f (х),

(2-4)

X

 

к которым, в частности (при т= 1), относится и

уравнение

Бесселя нулевого порядка.

 

Попытаемся отыскать в виде бесконечных рядов и обосно­ вать довольно широкий класс решений уравнения (2-4), кото­ рое в дальнейшем именуется обобщенным дифференциальным уравнением Бесселя нулевого порядка [25].

Поскольку доказательство сходимости полученных рядов и нахождение радиуса их сходимости в общем виде — вопрос далеко не тривиальный, авторам удалось решить только часть этой задачи. Однако для многих приложений, с которыми при­ ходится сталкиваться в теплофизике, этого оказывается впол­ не достаточно. В равной степени это относится и к доказатель­ ству единственности рассматриваемого класса решений. Сле­ дует заметить, что уравнение (2-4), вообще говоря, относится к более широкому классу уравнений, а именно — к уравнени­ ям Эмдена—Фаулера .[23, 35—37]

— (х? — \ + х° ит= 0.

(2-5)

dx \ dx ) ~

 

Однако результаты, приведенные далее, в очень малой степе­ ни перекликаются с теми, которые следуют из теории уравне­ ний Эмдена—Фаулера, и претендуют, таким образом, на но­ визну.

§1. Интегрирование обобщенного дифференциального

уравнения Бесселя без возмущающего члена [25]

Предположим, что в дифференциальном уравнении вида

и" -!-иг ± ит ■= f (х)

(2-6)

т — целое число, большее нуля, и будем искать положитель­ ное решение уравнения без возмущающего члена

и” + — и’ + ит = 0,

(2-7)

X

причем такое, которое удовлетворяет условию и(0) =и0>0. Поскольку в точке х= 0 уравнение (2-7) имеет особую точку, то, как будет видно из дальнейшего, для решения уравнения

30


достаточно только одного условия. Вторым условием является ограниченность решения в этой точке.

Заметим, что уравнение

и" -I-и' — ит = 0

(2-8)

z

 

получается из (2-7) заменой переменных z^ ix ,

поэтому перей­

ти от решения уравнения (2-7) к аналогичному решению урав­ нения (2-8) будет несложно.

Представим искомое решение дифференциального уравне­

ния (2-7)

в виде бесконечного ряда

 

со

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

 

 

 

И (х)

=

2

 

°кХк = ы° +

2 йкХ!і

 

 

(2-9)

 

 

 

 

Ä=0

 

 

 

k—\

 

 

 

 

и попытаемся определить неизвестные коэффициенты ак.

 

Для упрощения последующих выкладок будем

считать, что

и0 = 1 .' Это нисколько

не

уменьшает

общности

рассуждений,

так как уравнение

(2-7)

заменой и = и0Ѵ и z =

х ] / и'*~1

приво­

дится к уравнению

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V" -)— — V' -j- V"1= 0

 

 

 

(2-10)

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

с начальным условием V (0) =

Ѵ0=

1.

 

 

 

 

 

Представив,

таким образом, решение (2-7) в виде ряда

 

 

 

 

и (х) =

1 + 2

акхк,

 

 

 

(2-11)

будем иметь

 

 

 

 

 

 

к—1

 

 

 

 

 

оо

 

 

 

 

 

 

оа

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и' (х) — V

/га^Д“1;

и" (х) =

^

k (k — 1) akxk~2

(2-12)

 

 

ft=l

 

 

 

 

 

 

к=2

 

 

 

 

и, согласно формуле бинома Ньютона,

 

 

 

 

 

 

 

1+ Сш

00

 

 

п _і_ Г2

00

 

00

 

 

 

 

=

2

 

 

 

 

2

а^ + * "

 

 

1

a,iX

1 V-Wt 1

 

 

 

 

 

 

 

00

 

п—1

2=i

 

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

anXn 2

djX’ • • • 2

asxs =

2

« /

(2-13)

 

 

Л=1

 

 

/=1

 

S=1

 

ft=0

 

 

а 0 1> ак

 

 

 

к—1

 

 

k—2 /г—І—

 

 

-|- Ст2

апак-п

! ^"12

0,12

 

 

 

 

 

 

 

я=1

 

 

п=1

 

/=1

 

 

 

к—т+1 к—ш+2—п

 

 

к—1—n—j-----

 

 

 

 

• • • + О,!

2

ап

2

 

аі'

' '

2

 

asak—n—j------ s •

(2-14)

 

п=1

 

/=1

 

 

 

s=l

 

 

 

 

 

31


Действительно, из .бинома Ньютона

 

 

 

(1 + 2

а-х'1)"! = [ + с гп 2 а»хП+ с "‘

2

а* 2

+

П=1

П=1

П=1

/=1

 

со оо

СО СО

+ 2

^

2

öj ’ 2

ЯгЛ:1+ / + л + ■

•• +

Cm 2

йп2 ^ ’

 

 

П = 1

/ = 1

І = 1

 

 

 

/ ! = 1

/ = Г

 

 

 

 

...

2 « s2

a^

+s+’"+/+"

 

(а)

 

 

 

 

S = 1

< = 1

 

 

 

 

 

где Cm— число сочетаний из т по /.

 

 

сте­

Найдем выражения для коэффициента при А-ой (/г > 1)

пени X в (а).

Коэффициент

при

хк во

втором

слагаемом

(а),

очевидно,

равен Сmahl . Для того чтобы

найти аналогичный

ко­

эффициент из третьего слагаемого, нужно воспользоваться усло­

вием

(б\

j + n = k,

из.которого следует, что каждому определенному л соответст­ вует единственное j — k п. Кроме того, поскольку / ^ 1 , п может, следуя условию (б), принимать только значения от 1 до (k—1). Таким образом, коэффициент при хк в третьем слагаемом (а) будет равен

k — 1

Cm ^ anah-n•

п= 1

 

Из условия

 

І -|—j -}- Д — k

(в)

легко найти коэффициент при k-oii степени х в четвертом сла­ гаемом. Действительно, если мы зафиксируем п, то / и і могут

принимать только те значения,

которые

ведут

к

равенству

i+ j = k п, т. е. каждому

/

 

соответствует

единственное

i=.k п — /. Далее,

поскольку

 

1, то j

 

может

изменяться

только от 1 до

(k п — \), а п

(так как и

/ ^ 1

и

1) —

только от 1 до

\k — 2), откуда следует, что коэффициент при

 

 

 

 

 

 

к—2

k—n—l

 

 

 

хк в четвертом слагаемом (а)

равен Ст^’

а п

2

а^к-п-і-

 

 

 

 

 

 

 

Я=1

 

/= 1

 

 

 

 

Рассуждая аналогичным образом и далее, нетрудно найти

коэффициенты при хк в каждом слагаемом (а),

а с ними и все

выражение для коэффициента при k-он

(k ^ m )

степени х:

 

 

 

к—1

 

 

ft—2

 

ft—л—1

 

 

 

 

& к — Ст ßfe + Cm

 

 

+

Cm

 

 

 

 

^ j ^ h - n - i

+

 

• • •

 

 

n = l

 

 

л=1

 

 

;=1

 

 

 

 

 

ft—m+1

ft—m + 2 — n

 

 

k-l—n—j-

 

 

 

 

 

 

 

a,

2

aj

■■

V

 

 

üflh-n-j-.

- t

(Г)

 

 

jieesi

 

 

 

Л = 1

/=i

 

 

 

/=1

 

 

 

 

 

 

 

32


Заметим, что если k < т, то

последние

т k -j-1

слагае­

мые в (г) равны нулю.

и (2-13) в (2-7), получаем

 

Подставляя далее (2-12)

 

2 k (k - 1 ) ahxk~'+

2

k a ^ + 2

« кл*+1 = 0 ,'

(2-15)

к=2

к—1

/;=0

 

 

откуда методом неопределенных коэффициентов нетрудно определить постоянные а/£:

°і == 0; а2 =

; ан= — ^ ~ .

(2-16)

Кроме того, поскольку в каждое слагаемое коэффициентаап 12-14) входит сомножителем или ап, или всевозможные двой­ ные ßißj, тройные aißjßi и так далее произведения коэффици­ ентов исходного ряда (2-11), у которых сумма индексов оста­ ется постоянной и равной п, приходим к следующему выводу. Если п — число нечетное, то в каждое слагаемое коэффициен­ та ап входит сомножителем по крайней мере один коэффици­ ент a k (k ^n ) с нечетным индексом k, а это в свою очередь с учетом (2-16) приводит к тому, что

£4

0; о5 =

3*

 

Стаз -\-2Стві&2 “I- Сто.

(2-17)

52

и так.далее, т. е. коэффициенты а* ряда (2-11) с нечетным ин­ дексом k равны нулю:

 

 

fl2„+1 =

0,

А=0,

1, 2,

3, ...

 

(2-18)

Таким образом, искомый ряд имеет вид

 

 

 

 

и (X) = 1 -

- 1 X2 +

У ]

a2hx^,

 

(2-19).

где

 

 

/с.—2

 

 

 

h—m

 

 

 

 

 

1

 

 

! C,na2h_2 — Ст 2

а2па2/і-2п-2

 

2 ai

агк

 

• ■♦ "1" ^/;i

1

/1=1

 

 

 

ti=1

 

 

 

 

 

 

к—т+1—п

Ä—2

 

2п—...—2^

1

j

X

2

 

 

 

S—1

 

(2k f

 

/=і

 

 

 

 

 

 

 

k = 2 , 3, 4,

. . . ,

 

 

 

3. Зэк: 1208

33