книги из ГПНТБ / Харитонов, В. В. Некоторые нелинейные задачи теплопроводности

Рис. 10. Зависимость J™ — а и

— б

от т при различных г:

/ —1,0, 2—0,9,

3—0,8,

4—0,7

чения функций возрастают, т.

е.

I0(z1) < I o2 (zj <

/о ( z . X l t (zj,

J0(Zl) < Jl (Zl) <

Jo (Z,) < Jo (Zt).

-1. Зак. 1208

49

Jo (г) = 1

V2

I

т

I z

\ 4

m(3m—2)

/ г

21

(2!)2У3

VІ Т2

І)

(ЗШ!)22

І Т

j_

m (18m2 —29m+ 12)

I _z_ \s_

(4!)2

I

2 j

m (l80ms — 487m2 + 452m—144) / z u °

(3-1)

(51)2

VT

/0' (z) = 1 +

m

m (3m—2)

/ г

2 j

'

(2!)2

\ 2

J '

(3!)2

+

+

in (18/n2 — 29m + 12)

I

z ' 8

(4!)2

\

2

m (180m3 —487m2 + 452m —144)

/+ _ '10

(3-2)

(5!)2

[ 2

Аналогичный вид принимает и решение уравнения (2-6), когда

в правой части ) (х) =

b0 = const [26]:

и± (х, т, Ь0) =

1 +

(1 -г Ь0) I — j +

т(1+Ь0)

+

(2!)2

m (1 + b0) [т + 2 (m— 1) (1 + £0)]

(З!)2

\ 2 )

т(1 + Ьй) [т 2+ 11т ( т —1) (1+ Ь0) + 6 (т —1) (т —-2)(1+Ь0)2] ѵ

^

(4!)2

Jо" (г)

при г < 1 и т < 4 можно приближенно вычислять, исполь­

зуя первые четыре члена соответствующего бесконечного сте­

пенного ряда (3-1).

Оценка погрешностей такой

аппроксимации

приведена в табл.

Л.

Т а б л и ц а

А

Z

in

1

0,8

0,6

0,4

4

<0,006

<0,001

<0,0001

<0,00301

3

<0,002

<0,0004

<0,00004

<0,00001

2

<0,0004

„0,00007

<0,00001

<0,00001

Далее с ошибкой не более 0,002

при

т О 4

и |й0| -<

-С 0,1 ряд (3-3) аппроксимируется выражениями

«+ (х,

т,

Ь0) л? Jo (х) — Ь0 [Jöl(х) — і] ,

(3-6)

ІХ,

т,

Ь0) « /от (*) -і- Ь0 [ /5* (X) - 1 ] .

(3-7)

Кроме того, при Xl— ö0) Äi0 решение

уравнения (3-4). можно

приближенно вычислять по формуле

и+(х, т,

Ь0) » 1 + ——

[•/„ [Vпг лг) —і] ,

(3-8)

т

а при (1 + Ьо) « 0

решение уравнения (3—5)—по формуле

и_ (х, т, Ьа) ж Н — *

■[/0 (V т х) — і] .

(3-9)

т

Кроме того, в работе [25]

показано, что при /?г^4 и x=sC0,5

значения обобщенных функций

можно,

с достаточной

для

инженерных приложений точностью, вычислить

по выраже­

ниям вида

_

J'o (х) яй J Ü T Z L

_j_

W

ni 4 s

(3.10)

w

m

m

7&"(s)^

/n -

.1 - +

7о ( У ^ 4 .

•

(3-11)

m

m

4*

51

~ТѴ~ -г — ■--- ----А {ТА— То) =0,

dr-

г

dr

где А = г0еДгн, заменой переменных Ѳ =

Т/Т,,, ,ѵ = r/rH приво­

дится к безразмерному выражению

d-Q

j ___ de

—В (ѲА— Ѳо) = 0,

rf.v2

X

dx

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

2 =

х У В примет вид (2-4) (здесь

B = ArlT%

0" (2)

J_

0 ' (г) -

04 (2) =

_

04о .

(3.12)

2

При Ѳ0 = 0 (Т0 =0)

получаем уравнение

типа (2-8)

и решение

в виде ряда

(2-11), определенного как

обобщенная

функция

П (г).

_

вш/см-грш3\

е<С 1 н к порядка

Поскольку с0 =

5,67-Ю 12

1 вт!см-град,

то при га <

ІО3

см и Т„ Т

ІО3 град В= АгІТІ sC

-< 1, что при .V= /; /у, < 1 дает z — х У В

,< 1.

В общем случае для г ^ 1 необходимо

удовлетворить усло­

я

J0(x) =

1 ('

cos (.vsin Ѳ) de,

(3-13)

Л М =

JX „

exp (л:cos Ѳ) de,

(3-14)

■о

выражения (3-10) и (3-11) можно записать как

Jo (х) — ------------ 1---------(

cos (j/m xsinö))

de, (3-15)

m

mn

J

о

/ІЧ Ф

in—1

j exp (. in Xcos e) de, (3-16)

in

пт

/ 0(х):

exp X

1 +

1

1-9

1-9-25

У 2ях

1 !8 а'

2! (8 а-)2

3!

( 8 * ) 3

(3-19)

ях

U cos I X---- — 1 -f Fsin ( х —

Я

1

(3-20)

где

1-9

1-9-25-49

U = 1 —

2! (8 а ) 2

+

4! (8 а ) 4

V =

1

1-9-25

,

8 а

3! ( 8 а )3

Поэтому в области 0 < а < 0,5

и 0 < arg / т х <

я

будут спра­

ведливы формулы вида

Іо (X) :

III —1

exp {у m a )

1

in

1 +

8

Y m X

t

n 4

У

2 л а

1 !

+

1-9

+

•

П9-25

(3-21)

2! (8 1 m a ) 2

3! (8 ~\/~inx)3

1

(х)

т — 1

L

1

/

2

2

я

U cos У т X—

tI

т '

4

\

ях

1

т 4

4

Кsin

I' У т х ---- ,

(3-22)

где

1-9

,

1-9-25-49

и = 1

2! (8

I

m x f

'

4 !(8 ]/т л ')4

V =

1

1-9-25

m X

3! (8 Y m x f

8 I

J o’ (z) и /о (z) однородного обобщенного уравнения

Бесселя ну­

левого порядка

и" (z) Н— —и' (z) ± ит (z) = 0,

(3-23)

и" (х) 4 —

и' (л-) 4

и"' (,ѵ) = Ъ0

(3-24)

X

И

и” (х) 4 —

и' (х) — и"1(Д-) = Ь0

(3-24')

X

при граничных условиях и(0) = 1 и

для т = 1, 2,

3, 4

и при изменении аргумента от 0 до 1.

рассчитаны для

0,1

При этом таблицы решений уравнения (3-24)

а

таблицы

решений

уравнения

(3-24') — для

—0,8 < Ь 0< — 0,1

с шагом АЬа=0,1. Кроме того,

в таблицах

приведены значения первых

производных Jön (г)

и /о"! (z)

обоб­

щенных функций Бесселя.

ІО-5 во всем диапазоне

[0,

Для достижения точности порядка

1] изменения

аргумента при расчете были

использованы

11 членов соответствующих рядов.

Этого количества членов

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

10~6 при изме­

Бесселя (3-24) и (3-24'). Функции Jß' (z), J0m(г), Го (г)

и Іо” (г)

затабулированы с шагом Az =0,001,

а функции

и+(х,

т, Ь0) и

и_(х, т,

Ь0)—с шагом Az = 0,01.

задания аргумента

больше

В тех

случаях,

когда точность

шага таблиц, для

нахождения значений J™(z),

Jo”{z),

1о1(z) и

Напомним,

что суть

линейной интерполяции

состоит в

предположении

пропорциональности

приращения

функции

приращению аргумента.

Так, если заданное значение аргу­

мента z лежит между z0 и Z i = 2 0+ / i (

h — шаг таблиц), кото­

F(z) = F(za)

А-

п

Если две последующие разности

Ao = E(z0+ /i)—F(z0) и

Ai = F(zo-\-2h)—F(zc + h)

отличаются

в последнем знаке-не

более чем на 4.единицы,

ошибка линейной интерполяции ие

формулу

квадратичной

интерполяции

по Бесселю:

F (z) = F (z0) +

сД0 -

(Ах - А_,),

где с =

; Д_х =

F (z0) — F(z0 — h).

!.

Л ы к о в

А. В. Теория теплопроводности. «Высшая школа»,

М„ 1907.

2.

Л ы к о в

А. В. Тепломассообмен. «Энергия». А\., 1972.

М„ 1962.

3.

К у т а т е л а д з е

С. С. Основы теории теплообмена. Машгиз,

4.

М и х е е в

М. А.

Основы теплопередачи. Госэнергомздат, М.,

1956.

5.

Ш а ш к о в

А. Г.,

А б р а м е н к о

Т. II. Теплопроводность

газовых

смесей. «Энергия», М„ 1970.

Т. Н. Свойства переноса

газов п

6 . Ш а ш к о в

А. Г.,

А б р а м е н к о

жидкостей. «Наука и техника», Минск, 1973.

л о в

В. П. Методы

определения

теплопроводности и

температуропро­

водности. «Энергия». М., 1973.

8 . К а р е л оу

Г. С. Теория теплопроводности. ОГНЗ, М„ 1947.

9.

К а р е л оу

Г. С.,

Ег е р

Д. Теплопроводность твердых тел. «Наука»,

10.

М„ 1964.

Э. Р„

Д р е н

к

Р. М. Теория

теплообмена. Госэнергонздат,

Э к к е р т

_

М„ 1961.

Л. Л.,

Фр а й м а н

Ю. Е.

Теп.пофнзпческие

свойства

11.

В а с и л ь е в

12.

плохих проводников тепла. «Наука и техника», Минск, 1967.

по­

В а с и л ь е в

Л. Л„

Т а н а е в а

С. А. Теплофнзнческие свойства

13.

ристых материалов. «Наука и техника». Минск, 1971,

свойства

Н о в и ч е н о к

Л. Н.,

Шу л ь м а н

3.

П. Теплофизическпе

14.

полимеров. «Наука и техника», Минск. 1971.

П л а т у н о в

Е. С. Теплофнзнческие измерения в монотонном режиме.

15.

«Энергия», Л.. 1973.

В. С. Скоростной

метод

определения

теплофизн-

В о л ь к е н ш т е п н

16.

ческих характеристик материалов. «Энергия». Л.. 1971.

14,

№

5.

873,

М аи т у л е й ко

С. РІ.,

Х а р и т о н о в

В. В.

ИФЖ,

17.

1968.

А. Б. Автореф. каид. днсс. ІІТМО АН БССР, Минск,

В ер ж ни ск а я

1964.

Г. М. Автореф. канд. дисс. ІІТМО АН БССР, Минск. 1968.

18. В о л о х о в

19.

М ан т у л е и к о

С. И.,

Х а р и т о н о в

В. В.

Теплофизика

высоких

температур, 9, ДР 2, 373, 1971.

20.

Х а р и т о н о в В. В. ИФЖ,

14, № 2 , 347. 1968.

21.

X а р и т о н о в

В. В. Изв. АН

БССР,

серия

физико-энергетическая,

№ 3, 90, 1969.

22.

С м и р н о в

В. И. К\рс

высшей

математики,

т. 3. Гостехнздат,

М„

1951.

'

сл

Т а б л и ц а 1

ПРИЛОЖЕНИЕ