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

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

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

 

 

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

т\

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.

в

м‘-°С

 

 

 

 

 

Воздух

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

+ и

- ^

 

(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

( -

X - j - W

St

X

- j - St 1 /б f

2Bi

ехр ÖM X — —

X

1 — 2Ві

и

erfc I b і А - ~

(1.29)

Эти решения можно обобщить на случай нелинейного условия сопряжения, используя приближенный метод решения нелинейной задачи теплообмена при ламинар­ ном движении жидкости в канале [124]. Особенностью

18

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