Файл: Харитонов, В. В. Некоторые нелинейные задачи теплопроводности.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 |