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

книги из ГПНТБ / Смольский Б.М. Нестационарный теплообмен

.pdf
Скачиваний:
7
Добавлен:
23.10.2023
Размер:
4.63 Mб
Скачать

Приведенная зависимость может быть обобщена иа случай переменных теплофизических свойств [24].

Метод средней температуры [24]. В этом методе ис­ пользуется особенность нагрева (охлаждения) тела, со­ стоящая в том, что всегда имеется некоторая изотерми­ ческая поверхность в теле, температура на которой рав­ на в каждый момент времени средней температуре тела. При условии Fo^0,5 координата этой плоскости для пластины при граничных условиях второго рода опреде­ ляется из условия

(б2

- 34)/G ö2 =

0,

(2.25)

а решение уравнения

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

случае

сводится к выражению

 

 

Ч ь , т ) - / „ = - f -

О

(2.26)

 

А

 

Из сравнения условий (2.25) и (2.26) следует, что тем­ пература в выражении (2.26) является средней темпера­ турой тела.

В работе [24] показано, что формула (2.25) спра­ ведлива и при других граничных условиях. Тепловой по­ ток в этом методе определяется с помощью выражения (2.26) по измеренной температуре в точке я*.

Метод определения теплового потока по изменению температуры поверхности. Если стенку можно рас­

сматривать

как

полубесконечное

тело,

достаточно

для

определения

q

измерять температуру

на поверхности

[104]. Решение

уравнения теплопроводности в

 

этом

случае имеет вид

[104]

 

 

 

 

 

t(X,Т) =

 

 

Т

 

 

 

 

 

 

у —

Ф (Т*) ехр

4а (т — т*)

X

 

2 у па J .

 

 

 

 

 

 

о

di*

 

 

 

 

 

 

 

X

 

 

 

(2.27)

 

 

 

 

 

 

 

( т - т *)3/2

гдеср(т)—-найденная из опыта зависимость температуры поверхности от времени. Если найти производную темпера­ туры по координате' на поверхности и подставить ее в за­ кон теплопроводности, получим выражение для теплового потока на поверхности [104]

4. Зак. 1284

49

 

т

 

dt*

 

 

 

dcp (т*)

 

 

?(т,= w

lо

dt*

(т— т*)1/2

Определение коэффициента теплообмена при движе­

нии газа в трубе [13,

18, 93].

Коэффициент

теплообме­

на в нестационарных

условиях

определяется

так же,

как и в стационарных условиях,

следующим

выраже­

нием:

 

 

 

 

 

а,

 

 

 

(2.28)

В нестационарных случаях необходимо измерять все

или некоторые из величин в

правой части

выражения

(2.27) в зависимости от времени. При исследовании ста­ ционарного теплообмена достаточно определить эти ве­

личины для одного

какого-нибудь момента времени.

При рассмотрении

предыдущих методов предполага­

лось, что температура жидкости не зависит от времени и координаты вдоль направления по потоку, а q и ta не зависят от координаты. В общем случае все три пара­ метра могут зависеть от времени и координат.

Для определения температуры жидкости

в дюбом

сечении для любого момента времени

(х, т)

з работе

[18] решено одномерное уравнение энергии для средне­ калориметрической температуры потока при известных начальной температуре, расходе н зависимости теплово­ го потока на внутренней поверхности трубы q от х и т. Найденная численным методом решения среднекалори­ метрическая температура принималась равной искомой температуре жидкости tiK.

Для нахождения теплового потока и температуры на внутренней поверхности трубы находилось решение уравнения теплопроводности для стенки трубы. При этом предполагалось, что перетечек тепла по стенке вдоль трубы нет, источники тепла по толщине стенки распределены равномерно, тепловыделение определяется заданной величиной qv. Решение получено методом ко­ нечных разностей с помощью вычислительной машины. Используя найденные зависимости, можно было по из­ меренным значениям температуры наружной поверхно­ сти трубы и тепловыделению qv найти тепловой поток щ температуру на внутренней поверхности. В отличие от

50

первоначальной постановки задачи на этой стадии пред­ полагалось, что q и tn зависят от времени и не зависят от координаты. Таким относительно сложным расчетным

путем

находились величины, определяющие в соответ­

ствии с выражением (2.28) коэффициент

теплообмена.

Положительной

особенностью такой постановки за­

дачи

является

максимальное приближение к одно­

му из

реальных

случаев нестационарного

теплообмена.

Метод определения нестационарного теплового потока [105, 117]. В работаіх [105, 117] были проведены исследования нестационарного теплообмена цилиндри­ ческой вырезки из плоскопараллельной пластины и ша­ ра с потоком жидкости. Первичными эксперименталь­ ными результатами являлись зависимости температуры образца от времени.

Для нахождения тепловых потоков и коэффициентов

теплообмена' на поверхности образцов были

получены

решения уравнения

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

для

пластины

[105] и для шара [117].

 

неограниченной

Уравнение теплопроводности для

пластины (одномерная задача) имеет вид

 

 

 

dt

дЧ

 

(2.29)

 

 

дх

= а -----

 

 

 

дх2

 

 

В качестве граничных условий, исходя из опыта, бы­

ли использованы выражения

 

 

х = 0 ,

t(0,

т) =

А0 [1— ехр (— К0х)] ,

(2.30)

х = д ,

Ң8,

т) =

Л6 [1— ехр(— /С6т)]

 

и начальное условие t (х, 0). = 0.

Применяя преобразования Лапласа, уравнение теплопро­ водности и граничные условия можно записать в виде

4=о = А I ехр (— S t)

[1— ехр (— Д » ] dx = А0 ^

!

\

А0К 0

S + К0 )

S (S + К 0) ’

т

_

Л Д б

 

S(S + K6) ’

4*

51

d-t ( X , S)

5

t{x, S ) = 0.

 

dx2

 

 

 

 

 

Решение дифференциального

уравнения

для изобра­

жения

 

 

 

AüK6sh

 

 

 

t =

 

 

 

S (S + K6) Sh (

j / 1 . 6 J

h

A0K0sh I

Ä

(&— x)

 

+ '

S (S -f K0) sh

( 1/

Переходя к оригиналу, получим следующее выражение для распределения температуры по толщине пластины:

б — X

t (.V, т) = Аа

sin

sin l/ 4

Ах —

ехр (— а д —

•К ѵ

 

exp (— К6т) + 2А0 ^ ( —1)" X

sin

i

f

 

 

п—1

 

V

 

 

 

 

X

sin

И-п

б — X

ехр

аК

 

■2А6 V

(-1 )«

Рп

sin

 

ехр

Л—1

 

 

 

б2

 

 

 

 

 

б — X

— ! ! ] (

- ' ) ' ’

4 , 1

- V

sln

X

 

л=1

 

 

 

 

 

52

х е хр | - ^ т )

- 2 Л

6 V

Wn

H l ) - aiil_ Kstf X

 

'

n=I

 

xsin

X

exp

(2.31)

1-1,1

 

 

 

62

Для полых шаров, которые так же, как и вырезки из пластины, нагревались в жидкости постоянной темпера­ туры, уравнение теплопроводности имеет вид

dt

дЧ

, 2

dt

(2.32)

---- — d

------дг 2

-|- —

------дг

дх

г

 

В качестве краевых условий в этом случае были исполь­ зованы выражения:

при г = гх

 

 

 

t (rv

х) =

А1 [1— ехр ‘(— а д ] ,

 

при г = г2

 

 

 

t (га>

т) =

А 2 [1— ехр (— К2Х)]

(2.33)

и начальное условие: при X — О

t(r, 0) = 0.

Используя подстановку Т = tr, задачу для шара можно свести к задаче о нахождении температурного поля пла­ стины:

дТ (г, т) _

д2Т (г; т)

дх

дг2

Начальная температура Т (0, г) = 0. Граничными условиями являются выражения

т К т) = ГЛ [Г— ех Р (— Я » ] ,

 

Т (га, т) = г2А 2 [1— ехр (— /Сат)] .

(2.34)

Такая задача рассмотрена выше. Ее решением явля­ ется выражение (2.31). Чтобы получить решение для полого шара, необходимо использовать указанную вы­ ше подстановку T = tr, а также следующие соотношения:x

x = r — rv 8*=rt — rv

53

Решение для температуры полого шара в любой момент времени и в любом сечении имеет вид

t (г, т) = А, Гі (г2г)

Н- А„ гЛГ~ гі)

X

'' (X гі)

г(Го — гд

 

in[ l

 

 

 

 

exp (— Kjx) A2

 

X

X

 

 

 

 

 

 

sin

[ l

/ 4

(»•— rl)

 

 

 

 

 

X

 

exp (— K2t) + 2 A1

X

sin

l /

:

 

('2 -

rl)

 

 

 

 

 

 

x

j

(

-

1)"

sm

пл

 

 

X

(2.35)

 

 

n=l

 

 

ПЛ

rn— r1

J

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

an-n

 

 

 

 

 

 

 

Xexp

 

 

“"T T “*

 

+ 2j43

У

д - ѵ

- і -

X

 

( r , - r d 2

 

 

 

n=l

 

П Л

 

 

 

 

 

 

 

 

 

 

 

 

xsm

пл

r — r1 '

exp

(r2 — ri)2

2AX

X

 

 

r2— r!

 

 

 

 

 

X

 

 

 

 

апл

sm

П Л

 

 

X

 

 

 

 

 

 

 

 

n= I

 

 

ап2л2— /Ci (r2— г1)“

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 „ 2

 

00

 

 

 

xexp

 

 

an‘n

- 2 A2

V

(—l)n X

 

(Г.-Гх)8

 

 

 

 

 

л=і

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

апл

sm mt

r — r,

X

 

X an2n2— К.г (r2— r,)

 

r 2 —

Г1

J

 

 

 

 

 

Xexp

( 'W

, ) 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

54

т , сск

5

10

15

5

10

15

5

10

15

Т а б л и ц а 2

 

-

 

 

 

 

сумм

 

 

со

\

 

 

 

8H I

всех

-0{

 

=п

 

 

 

 

 

2

 

 

 

Вклад

рас■

 

 

 

а

а

 

а .

UP

 

 

 

 

 

" 1-

 

 

в4-

 

 

 

/1 =

1

 

 

 

 

0,8100

1,7800

+ 0,8500

+ 1,8500 0,0900

11,5

0,9550

0,0550

0,1180

+ 0,0600

+ 0,1340 0,0210

22,5

0,0935

0,0033

0,0071

+ 0,0036

+ 0,0080 0,0012

31,2

0,0038

 

 

 

/і =

2

 

 

 

 

+ 0,0470

0,0960

0,0460

+ 0,0980

0,0030

11,5

0,0260

+ 0,0003

0,0007

0,0003

+ 0,0007

0,0015

22,5

0,0068

0,000

0,000

0,000

0,000

0,000

31,2

0,000

 

 

 

/і =

3

 

 

 

 

0,0024

0,0037

+ 0,0024

+ 0,0037

0,000

11,5

0,000

 

0,000

0,000

0,000

0,000

0,000

22,5

0,000

 

0,000

0,000

0,000

0,000

0,000

31,2

0,000

В уравнениях (2.31). и (2.35) были сделаны оценки вклада отдельных членов в значение температуры при различных значениях х для пластины и шара.

В табл. 2 приведены в качестве примера результаты расчета вклада сумм в значение температуры для по­ лого шара из латуни с толщиной стенки 25 мм.

При обработке экспериментальных результатов вклад сумм в уравнениях (2.31) и (2.35) оценивался для всех образцов. Было установлено, что он не превышал 1 % от значения температуры.

После дифференцирования уравнений (2.31) и (2.35) и подстановки их в закон Фурье были получены уравне­ ния для теплового потока:

на поверхности цилиндрической вырезки из плоско­

параллельной пластины

 

q — X

Ар____ ^6

exp (—/С0т)+

б

б

 

sin

 

 

 

55

+ л»V

 

 

c,g

 

К,

 

ехр (— К 6т) -Ь

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

ап2я 2

 

2АЬ X

 

 

 

 

" exp

X —

 

 

л==1

 

 

б2

 

б

 

 

 

 

 

 

 

 

 

 

 

 

exp

an2n2

 

 

 

 

 

 

s

б2

-

- у

1

У (—1)"Х

(2.36)

 

 

 

 

л=1

2я 2

 

а/і“Я

 

Л = І

 

 

 

 

 

 

2ААа

 

X

 

 

 

exp

^

2 тг2

х

+

 

X

 

 

 

б2

 

 

 

!- / С 0б2

 

 

 

 

 

 

 

 

 

2—2“ТТ “•

 

 

 

 

 

 

X У 1

 

п-л

 

exp

I —

 

 

шг2я 2 — /Слб2

 

 

 

п = \

 

 

 

б2

 

 

 

 

 

 

 

 

 

 

 

 

н а вн еш н ей п о вер х н о сти п олого ш а р а

 

 

 

q ~ — k

Ѵ і

 

Ѵ і

 

 

 

 

 

 

ГЛГ2 ~ Гі)

Г2ІГ2 ~ Гі)

 

 

 

 

 

л

*1

 

 

ехр (— а д

+

 

 

г2 Sin

(^ — гі)

 

 

 

 

 

J ^ j/'

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

 

 

ехр (— /Сат) X

 

 

 

 

 

 

 

 

 

 

X cos

 

 

 

ft

 

 

тс

 

 

 

 

 

 

 

 

 

 

 

sin

 

'

к ,

Оя — ri)

-

 

У ( - i ) nx

 

 

 

 

 

 

 

 

 

 

rz(r-2— ri)

 

 

 

 

 

 

 

2„2

 

2Aa

 

 

 

 

X exp

 

aii-n

 

 

X

(2.37)

 

 

 

 

 

 

 

56

X ^

an2я 2

2Л1л1ая2

 

exp

T

 

 

n=l

(r2- r t f

r2 (r2~ ф)

 

X

 

exp

(O — ri f

+

оп2я2K x (r2— fj)2

 

2Л2ая2

’v S

пг

 

 

(О — ^i)

й/г2я 2

/(,, (r2— г,)2

 

 

 

ЯЯ2Я2

 

 

X exp

l ) 2 T.

 

 

 

( 'W

 

Оценка вклада

отдельных

членов, аналогичная

оценке в случае формул (2.31) и (2.35), была сделана и для формул теплового потока (2.36) и (2.37). В резуль­ тате выполненных оценок (табл. 3) было установлено, что добавка в тепловой поток за счет слагаемых форму­ лы (2.37), в которые входят бесконечные суммы, в са­ мом худшем случае не превышала 4% и в дальнейших расчетах не учитывалась. Это позволило упростить вы­ ражения для теплового потока на поверхности теплооб­

мена.

Окончательные

расчетные

формулы имеют вид:

для пластины

 

 

 

 

 

 

 

 

 

ехр (— Л » +

+

Л ‘ Ѵ

4 - Ч - / 4 - * exp (— /С6т)Ь (2.38)

для полого шара

 

 

 

 

 

Яя = — х

 

 

^2ГХ

 

 

Г2(О — О)

Г2(Г2 — Гі)

 

 

 

 

 

 

 

Г ~ к

1__

 

'

.__ Г-

 

А ггг

V

 

ехр (— K jt) 4- А .2 л /

X

(Г2 ~ >\) '

57 ,

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

3

 

 

 

 

 

 

оо

 

1

00

 

 

 

 

1 , сек

 

2А,

1

У ...

'2А,

V ...

2А„

 

 

 

 

 

1

V ...

^ д о п '10 3'

V 10 °

g,

%

 

«=1

"

rs—r L

п=[

 

гг (г,—Г,) _i_J

 

Г2—Г|

/I

 

 

 

 

«=1

 

 

am/ м 2

ет/м-

 

 

3

5

10

15

3

5

10

15

 

 

11=

+ 0 , 4310- ІО3

0, 9420- Ю3

0,4800- ІО3

+ 0,1520 ІО3

- 0, 3250- ІО3

0, 1640-103

+ 0 , 0110- ІО3

0,0230- Ю3

0, 0220- ІО3

+ 0,0007 ІО3

- 0,0015- ІО3

0, 0007- LQ3

 

 

п = 2

+ 0, 0038- 10s

0, 0082- Ю 3

- 0, 0039- ІО3

+ 0,0001- ІО3

- 0, 0002- ІО3

0, 0001- 103

0,0000

0,0000

0,0000

0,0000

0,0000

0,0000

+ 0, 0700- Ю 3

7,05

125,0

 

+ 0, 3700- ІО3

3,64

113,0

 

+ 0, 0270- 103

0,15

86,5

+ 0, 0007. Ю3

0,02

71,0

 

+ 0, 0085- ІО3

0,021

125

0,0167

4- 0, 0002- Ю3

0,001

113

0,0005

0,0000

0,000

86

0,0000

0,0000

0,000

71

0,0000

Соседние файлы в папке книги из ГПНТБ