![](/user_photo/_userpic.png)
книги из ГПНТБ / Смольский Б.М. Нестационарный теплообмен
.pdf
|
|
1. ТЕПЛООБМЕН В КАНАЛЕ |
||||
на |
Первые решения задач нестационарного теплообме |
|||||
при |
течении |
в канале |
относятся |
к 1959—1960 гг. |
||
[1, |
4, 6]. |
Были рассмотрены задачи о нестационарном |
||||
теплообмене в трубе |
при |
стабилизированном течении |
||||
жидкости и на начальном |
термическом |
участке трубы |
||||
[4, 6]. В отличие от работы |
[2] рассматривалась не со |
|||||
пряженная задача, |
а |
задавались условия на стенке. |
Предполагалось, что в некоторый момент времени тем пература стенки или тепловой поток менялись до нового постоянного значения. Такая постановка задачи соот ветствует случаю тонкой стенки из температуропровод ного материала. Кроме того, предполагалось, что физи ческие свойства жидкости постоянны, диссипация энер гии в потоке отсутствует, теплопроводность вдоль оси трубы несущественна. В некоторых работах, в частности в работе [8], принималось, что скорость потока постоян на по сечению. Авторы пренебрегали также теплоем костью стенки. При решении задачи в некоторых случа ях использовалось уравнение энергии в интегральной форме.
Задача о теплообмене в круглой трубе при скачкооб разном изменении температуры стенки во времени в без размерной форме может быть сформулирована следую щим образом [1]:
где
В момент времени т —0 температура стенки изменяется ступенчато от ta до tn:
Т (х, 1, Fo) = 1..
9
От записанного в общем виде выражения для темпе ратуры требуется, чтобы оно удовлетворяло уравнению энергии (1.1) в интегральной форме.
Ненззеетная функция cpn(x, Fo) находится из полу ченного уравнения методом характеристик. Безразмер ная плотность теплового потока на стенке
Qifо |
|
|
e-VuFo; |
р 0 ^ |
|
- |
-2W |
_е2- |
_• (1.2) |
||
M t u - Q |
|||||
|
Л J |
е Еп' , |
F o > a„x |
||
|
|
п=іО |
|||
9гв |
|
|
|
|
|
•Ѵйі ^оУ |
|
|
|
|
Рис. 1. Плотность теплового потока на стенке |
круглой трубы |
при |
скачкообразном изменении температуры стенки: |
1 — х=0,03; |
2 — |
0,04; 3 — 0,07; 4 — 0,1; 5 — 0,15; 5 — 0,20; 7 — 0,30; 5 — 0,40 |
|
В случае скачка температуры стенки плотность теп лового потока на стенке зависит от координаты вдоль трубы xjd и критерия Fo (рис. 1). Величина qn с ростом Fo уменьшается, приближаясь к постоянному значению, зависящему от xfd. Расчет, приведенный в работе [1],
10
показывает, что при определенных значениях x/d и Ре новое стационарное состояние наступит для трансфор маторного масла через 18 сек, для воды через 8 сек и для воздуха через 0,05 сек.
Используя линейность уравнения энергии, с помощью принципа Дюамеля можно получить выражение для теплового потока на стенке при произвольном изменении
температуры |
стенки Ѳп |
(Fo) [1]: |
|
|
Чпго |
(Fo, |
х) = 2 |
2Вп {еХР(“ е^ (Ѳп)ро-апГ |
- |
I |
||||
|
Fo |
|
|
|
|
|
Ynexpl—T»(Fo—Fo*)]0„(Fo*)dFo*l— |
(1.3) |
Fo- 0j1a-
CO Fo
—2Vnß n j* exp [— yn (Fo — Fo*)] Ѳп (Fo*) dFo*.
n^N 0
По сравнению с рассмотренной выше задачей случай теплообмена при изменении на стенке теплового потока во времени ближе к реальным условиям теплообмена в трубе. Уравнение энергии в интегральной форме на на чальном термическом участке круглой трубы имеет вид
|
д |
|
|
д |
— |
I О'о — У) ( * — * о ) dlJ |
+ |
|
I ( Д — У ) (t — to) dlJ = |
. |
о |
|
|
о |
|
= |
го |
<7п |
(1.4) |
|
рСр |
|||
|
|
|
|
Принимается квадратная зависимость для распределе ния скорости и температуры в пограничном слое
-^ - = |
12 J L _ |
ф |
и t — tо = |
?ПА |
1 — £ |
(1.5) |
2 1 |
т\ |
2К |
||||
w |
'( г0 |
|
Д |
|
После подстановки (1.5) в (1.4) полученное уравне ние для безразмерной толщины теплового' пограничного
11
слоя А решается методом характеристик [1]. Изменение температуры стенки определяется выражением
|
= |
(1.6> |
а А находится для малых Fo из уравнения |
|
|
— А2-----— K3 = Fo. |
(1.7) |
|
6 |
24 |
|
Изменение температуры стенки при произвольном из менении во времени теплового потока находится с по мощью принципа Дюамеля для нестационарного ре жима теплообмена (F o ^ F o CT):
tn(x, Fo) — t0= |
A foct (<7niFo=0 |
-i-A F o ^ o .d F o , .. |
|||
|
|
AK |
|
öFo* |
|
|
|
|
|
|
( 1.8) |
В случае ступенчатого изменения теплового потока в |
|||||
2 раза по сравнению |
с первоначальным |
при реальных |
|||
значениях |
параметров |
потока жидкости |
в трубе время |
||
протекания |
процесса |
нестационарного |
теплообмена до- |
установления нового состояния составляет около 90 сек. При периодическом изменении во времени градиента давления (пульсирующий поток) также происходит про цесс нестационарного теплообмена, степень отличия ин тенсивное^ которого от интенсивности стационарного теплообмена зависит от частоты пульсаций [1], причем с уменьшением частоты процесс становится квазистационарным. Влияние числа Рг таково, что если при опре деленной частоте и других постоянныхпараметрах при значении Рг<СІ процесс теплообмена является квазистационарным, то при достаточно больших значениях числа Рг интенсивность нестационарного теплообмена существенно превосходит стационарный тепловой поток. Разница в величинах растет приблизительно пропор
ционально амплитуде колебаний.
Для случая изменения скорости теплового потока во времени результаты [8] отличаются на 25% от соот ветствующих величин для стационарного процесса. При постоянной скорости и изменении теплового потока ис
12
пользование постояных значений коэффициента тепло обмена давало хорошее совпадение с полученными ре зультатами.
Используя допущения, перечисленные выше, автор' работ [9, 10] свел задачу к решению уравнения тепло-' проводностп и получил зависимости для Nuu для ряда конкретных случаев.
Результаты анализа процесса нестационарного тепло обмена при линейном изменении во времени температу ры стенки [10] свидетельствуют о том, что значение нестационарного критерия NuH в несколько раз отлича ется в первый период времени от значения критерия Nuct для стационарных условий.
Впредыдущих исследованиях нестационарная тепло проводность в стенке заменялась заданием условий на. внутренней поверхности трубы. Это ограничивает об ласть применения полученных результатов случаем тон ких теплопроводных стенок.
Вработе [21] рассмотрено изменение температуры жидкости на входе в трубу с учетом термического со противления стенки. Труба предполагается теплоизоли рованной снаружи; расход, теплофизические константы жидкости и стенки постоянны. Не учитывается радиаль ная и осевая теплопроводность в жидкости, течение при нимается стержневым. Вводится средняя температура
стенки 7 в «параболическом приближении». Коэффи циент теплообмена, отнесенный к средней температуре стенки, принимается не зависящим от температуры.
Учитывая указанные допущения, получаем математи ческую формулировку задачи [21]
d t m |
ву |
|
,к |
— h ( t |
t m), |
( 1 . 9 ) |
+ |
|
|||||
ÖT |
|
|
o x |
|
|
|
|
|
|
|
|
|
. ( 1 . 1 0 ) ' |
і ж = |
f |
(j) |
п р и А- = |
0 , |
( 1 . 1 1 ) |
|
г д е h = a"F/f. KC m у ж , |
Я |
= |
а Т / / тС тѴт, |
«•■=,( |
|
|
|
|
|
|
|
, а + |
3 « ' ) ■ |
13
Введя переменные I = hx/w, т) = Я (т— x/w), получаем уравнения
діш = t t},., |
(1.12) |
|
dl |
|
|
Ö11 |
и - |
(1.13) |
|
|
Используя преобразования Лапласа, находим реше ние для температуры жидкости и стенки при произволь-- ном законе изменения температуры на входе жидкости в трубу.
Из общего выражения следуют решения для линей ного, экспоненщіалы-юго, ступенчатого и синусоидально го изменения температуры жидкости в сечении .ѵ=0. В частности, в случае скачка температуры на входе тем пература жидкости и стенки равна соответственно
іж= Аг'жехр (— £) [срх + |
ехр(— rj)/0 (2 і grj )], |
(1.14) |
|
Т= Д#жФхехр(—S), |
|
(1.15) |
|
|
00 |
|
|
Фх = exp (— ё) — exp(— 11) 2 |
К (2 / |
ІЛ ) |
|
|
п=0 |
|
|
при 11/1 > 1 , |
|
(1.16) |
|
оо |
|
|
|
Фх = ех р (— |
Іп[2 Ѵ%ц) |
при ii/g < l. |
(1.17) |
л = 1
При изменении условий на стенке канала по экспо ненциальному закону число Nun во времени неограни ченно возрастает.
В работе [22] рассмотрено течение несжимаемой жидкости в трубе при изменении во времени и вдоль оси ее температуры. Принято, что отношение разности ра диусов трубы к внутреннему радису мало, скорость жид кости постоянна во времени и по сечению, теплоемкостью единицы объема жидкости по сравнению с теплоемкостью единицы объема стенки можно пренебречь. Последнее допущение, как видно из табл. 1, справедливо, например, для сочетания воздуха с нержавеющей сталью. Оно про тивоположно допущению о пренебрежимо малой тепло емкости стенки, сделанному в рассмотренных выше раі ботах.
14
|
|
|
|
|
|
Т а б л и ц а 1 |
|
|
О |
О |
|
|
|
|
|
Вещество |
о |
|
•ч |
а* |
ср, |
срвещ |
с^вещ |
f\j |
|
||||||
|
1C |
5 |
|
X |
кдж |
срвозд |
Фвод |
|
- £ |
6. |
в |
м‘-°С |
|||
|
5С |
|
|
|
|
||
Воздух |
1,005 0,0245 |
1,29 |
0,068 |
1,255 |
1,00 |
0,00031 |
|
Вода |
4,186 0,6060 ■1,00 |
0,0005 |
4186 |
3230 |
1,00 |
||
Медь |
0,385 |
375 |
8950 |
0,40 |
3438 |
2650 |
0,82 |
Сталь |
0,502 |
16,1 |
7900 |
0,0145 |
3980 |
3060 |
0,95 |
1Х18Н9Т |
В отличие от предыдущей задачи учитывается изме нение температуры стенки . по радиусу tT{x, г, т). Так же, как в работе [21], коэффициент теплообмена а прини мается постоянным. Снаружи трубы происходит тепло обмен с известным коэффициентом теплообмена а\.
Уравнения энергии для потока жидкости и теплопро водности имеют соответственно вид
Gcp%*- |
рсг |
д(т |
* |
дх1 |
|
|
2 а |
|
||
р дх |
|
|
дх |
|
|
А |
|
|||
дН,т |
^ |
J_ |
д/т |
|
d2ty |
_ |
рср |
dt^ |
(1.19) |
|
|
Зга |
' |
г |
дг |
' |
дх2 |
~ |
Т " ’ |
дх |
|
|
|
|||||||||
С помощью новых переменных |
|
|
|
|||||||
ТJ m= |
|
|
|
|
— ta |
J |
Г |
2ax |
||
U |
|
|
|
|
RGcp |
|||||
^BX |
|
|
У |
— |
’ |
|
||||
где tB— температура среды |
вне |
трубы, |
/вх — температура) |
|||||||
жидкости на входе в трубу. Уравнения (1.18), (1.19) |
при |
|||||||||
водятся к безразмерному виду. Упростив уравнения, |
с уче |
|||||||||
том сделанных допущений |
получим |
|
|
|
||||||
|
|
а2т т |
, _і_ |
_ |
а г т |
_ |
дТт |
|
( 1. 20) |
|
|
|
аг]2 |
г] |
|
ÖTi |
|
öFo |
|
||
|
|
|
|
|
|
|||||
|
|
|
|
дТ.., 4 -Т |
— Т |
|
( 1-21) |
|||
|
|
|
|
д\ |
I |
1 Ж |
1 |
т* |
|
|
|
|
|
|
|
|
|
|
|
|
В начальный момент времени температура жидкости на входе в трубу скачком меняется от температуры сре ды tBдо значения tBX.
Граничные условия определяются выражениями
|
|
7\к(0, Fo) = l, |
|||
|
|
д Т \ |
— Тк T-г,,1=ь |
||
Ra |
V дг\ |
||||
Д1=І |
( 1. 22) |
||||
X |
I |
дТ, |
|
||
|
= 7 ,т,г1=1+о- |
||||
Ra, |
\ |
— 1 |
|||
<3і] |
11= Н-а |
|
|||
Используя преобразования |
Лапласа, находим вна |
чале решение уравнения (1.20). При нахождении обрат ных преобразований используется условие а<СІ.
Решения для внутренней поверхности трубы и по тока жидкости имеют соответственно вид
Т
Тт,ъ=1 = L'e-<A+ ^ - Lt / 0 (2 У Ц х ) + ß 'j e-w+i)T-t6 х
|
о |
|
X / 0 (2 |/T F ) drc, |
(1.23) |
|
Тж= в-И-И)х-і6/0 (2 У Щ ] |
+ в j |
х |
X / 0 (2 V 'IF ) |
о |
(1.24) |
rfT, |
||
где L, L', В, В', А — постоянные [22]. |
сводят |
|
В работе показано, что полученные решения |
ся к известным зависимостям, полученным в предполо жении бесконечно большой теплопроводности стенки.
Расчеты показали, что для нестационарных условий теплообмена применительно к гиперзвуковой аэродина мической трубе [22] в стенке трубы из нержавеющей стали возникает значительный перепад температуры. Скорости изменения температуры на наружной и внут ренней поверхностях существенно различны. Результаты работы подтверждают необходимость учета влияния на теплообмен процесса нестационарной теплопроводности в стенке, т. е. толщины (даже если, относительная тол щина невелика) и теплофизических характеристик ма териала стенки.
Вывод о необходимости учета конструктивных раз-
16
меров и теплофизических свойств стенки при расчете теп лообмена содержится также в работе [121]. В ней рас смотрена сопряженная задача стационарного теплооб мена твэлов реакторов. Решена система уравнений, включающая уравнение теплопроводности для стержня делящегося материала (с учетом тепловыделения), уравнение теплопроводности для оболочки и уравнение энергии для потока теплоносителя. В результате получе ны выражения для распределения температуры и теп ловых потоков, в которые входят как характеристики потока жидкости, так и характеристики материала деля щегося' стержня и оболочки. Для режимов расхолажи вания и выхода на номинальную мощность, о которых упоминается в работе, необходимо было учесть времен ные члены уравнений теплопроводности и энергии.
Вслучае существенных перепадов температуры по толщине стенки трубы необходимо при постановке зада чи нестационарного теплообмена учитывать уравнение теплопроводности для стенки [123].
Рассмотренная авторами работы [123] система урав нений описывает изменение скорости, давления и темпе ратуры газа вдоль канала и температуры стенки в ра диальном и, осевом направлении. На вХЪде в канал мо гут меняться скорость, давление и температура газа.
Вкачестве условия сопряжения между потоком и стенкой используется одна из критериальных зависи мостей для коэффициента теплообмена при течении в трубе. Коэффициент теплообмена, таким образом, явля
ется функцией параметров потока -газа и температурй стенки.
Для решения задачи авторы предлагают численный
метод, основанный на замене частных |
производных по |
|||||
координатам разностями величин, в результате чего за |
||||||
дача сводится к системе обыкновенных дифференциаль |
||||||
ных уравнений. |
постоянных |
газодинамических |
парамет |
|||
Для случая |
||||||
ров во времени |
и ■постоянного |
коэффициента |
теплооб |
|||
мена по длине трубы получено |
аналитическое решение |
|||||
задачи. |
|
задача |
формулируется сле |
|||
При таких допущениях |
||||||
дующим образом: |
|
|
|
_____________ |
||
|
= |
. _L |
|
Гис. |
||
дх |
дг |
учио- |
||||
\ д г |
^ |
г |
■}■. .-с-те .а С-с с |
|||
2. Зак.' 1284 |
|
|
|
|
|
"ішгр |
|
|
|
|
|
Ч: ІТА |
О . XT? B A J |
+ и |
- ^ |
|
4а |
(1.26) |
|
c ppd [ ' ж - У . |
|||
дх |
дх |
|
|
|
a(t„ |
« = - |
4 |
f |
.(1.27) |
|
|
где /т — температура стенки трубы; £ж — температура жид кости.
При т = 0 tn. = ta = t0. В начальный момент времени происходит ступенчатое изменение температуры жидкости на входе в трубу от t0 до /Жо. Используя безразмерные переменные
гр |
__ |
' Т |
' _____ |
, |
л |
м* --- |
) і |
IT! |
|
|
|
|
Т |
|
^Жо
уравнения (1.25) — (1.27) приводятся к безразмерному виду. Для решения задачи используется преобразование Лапласа. Обратные преобразования для начальной ста дии процесса дают [123]
Тт (т, X, |
г) |
= (Г~ Гп| аВі Г Тж(х* |
X ) X |
|
|
|
|
/•/•о |
J |
|
|
|
|
|
о |
|
|
X ехр [Ь2(т— т*)] erfc [ö ]/"т — т* J dx -f |
X |
(1.28) |
|||
|
Т |
|
dx* |
|
|
X I |
|
(т*. X) |
|
|
|
|
|
|
|
||
|
|
]/"п (т — X*) |
|
|
|
Г ж (г, X) = exp |
( - |
4х |
X - j - W |
St |
X |
- j - St 1 /б f |
2Bi |
ехр ÖM X — — |
X |
|
1 — 2Ві |
и |
erfc I b і А - ~
(1.29)
Эти решения можно обобщить на случай нелинейного условия сопряжения, используя приближенный метод решения нелинейной задачи теплообмена при ламинар ном движении жидкости в канале [124]. Особенностью
18