![](/user_photo/_userpic.png)
книги из ГПНТБ / Смольский Б.М. Нестационарный теплообмен
.pdfПриведенная зависимость может быть обобщена иа случай переменных теплофизических свойств [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
![](/html/65386/283/html_pl2nArjsp7.6VUK/htmlconvd-BVHByI56x1.jpg)
Решение для температуры полого шара в любой момент времени и в любом сечении имеет вид
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я2— K 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 |