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

книги из ГПНТБ / Козырев, А. П. Теория тепловых и гидродинамических процессов в атомных энергетических установках учеб. пособие

.pdf
Скачиваний:
13
Добавлен:
22.10.2023
Размер:
15.28 Mб
Скачать

ния. Для теплотехнической системы в целом получаемые при этом уравнения баланса, или интегральные уравне­ ния сохранения, не позволяют исследовать внутренние процессы. Детальное изучение внутренних процессов в системе становится возможным лишь при переходе к диф­ ференциальной форме уравнений сохранения, называемых в этом случае уравнениями переноса. Дифференциальные уравнения переноса, как правило, включают большее ко­ личество неизвестных, чем число уравнений. Это объяс­ няется тем, что используемые законы физики описывают макропроцессы и не учитывают дискретный характер стро­ ения среды и микроскопические особенности процессов. Для замыкания системы дифференциальных уравнений и последующего их интегрирования на помощь приходят за­ коны о дополнительных связях между переменными про­ цесса. Такие законы позволяют считать среду, где про­ текает процесс, сплошной и непрерывной и дают возмож­ ность применять статистические свойства среды (давле­ ние, температуру и т .д .) к ее дифференциальным эле­ ментам. Полученные из наблюдений за явлениями приро­ ды указанные законы иногда в литературе[55, 841 назы­

вают феноменологическими (от слова "феномен" - явление), феноменологическим называют и сам метод исследования процесса

Коэффициенты пропорциональности в дополнительных связях между искомыми величинами (феноменологически­ ми законами) находятся опытным путем и называются физическими параметрами вещества. Так, при рассмотре­ нии молекулярного (диффузионного) переноса энергии, количества движения и массы применяются классические законы: теплопроводности - Фурье, внутреннего трения в вязкой жидкости - Ньютона и диффузии массы - Фика.

1 0

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

В отличие от феноменологического метода статистиче­

ский метод, исходя

из

определенной дискретной структу­

ры сред,

применяет

аппарат математической статистики

и теорию

вероятностей

к законам движения и распреде­

ления энергии молекул

рассматриваемой системы. Такой

подход позволяет вскрыть сущность термодинамических явлений и в принципе получить все результаты феномено­ логической термодинамики, но гораздо более сложным путем.

До недавнего времени в науке о теплообмене господ­ ствовал эмпирический подход. За последние два деся­ тилетия сделаны большие успехи в развитии аналитичес­ ких методов исследования теплообмена, особенно в те­ ории конвекции однофазной жидкости и теплопроводно­ сти, где на долю эксперимента все чаще отводится про­ верка математической модели процесса. Однако в целом ряде других разделов учения о теплообмене решающая роль все еще принадлежит эксперименту.

§ 2. Уравнение переноса^теплоты в вещественней среде

Рассмотрим неравномерно нагретую вещественную среду (твердое тело или поток жидкости). В этой среде выделим произвольный объем у с помощью воображаемой не-

I I

подвижной и проницаемой контрольной поверхности F • В выделенном объеме действуют внутренние источники тепла. Интенсивность внутренних источников тепла обоз­

начим через

(J,v bt/ m3(количество

тепла, выделяемое

в единицу времени в единице объема).

Тогда во всем

объеме будет

выделяться количество

тепла,

равное

 

/ q y civ-

 

 

 

(1Л)

Если

союзником минус,

то

в

среде

действуют

стоки тепла. Часть тепла будет вытекать через поверх­

ность

F

 

путем теплопроводности. Поток тепла через

поверхность

F

будет равен

 

 

 

 

J " ( 4 - < U d F ’

( '- 2 )

 

q

 

F

 

 

где

-

вектор плотности теплового потока, направ­

 

 

 

ленный в сторону,

противоположную гради-

 

_

 

енту

температур?

 

 

8 ^ -

единичный вектор по внешней нормали п к

 

 

 

верхности.

 

Полное

изменение

количества тепла в выделенном объеме

в единицу

времени

равно

 

 

 

 

 

 

(1 .3)

v

где Q y - количество тепла, получаемое единичным объемом в единицу времени, вт/м3 .

На основе закона сохранения тепловой энергии будем иметь интегральное уравнение баланса

j Q v d V =

-

d F + J f y v d V . (1 .4)

V

F

V

Изменение теплосодержания объема зависит от интенсив­ ности процесса отвода тепла теплопроводностью и подвода тепла за счет внутренних источников.

Т2

Используя теорему Остроградского-Гаусса о связи между потоком вектора через замкнутую поверхность и дивиргенцией вектора, интеграл по поверхности F (1 .2) можно преобразовать в объемный

 

 

d F = fdlirfydV.

(1 .5)

F

'

v

можно переписать в ви-

С учетом (1 .5) уравнение (1 .4)

j a v d v + § d i v q , d v = f a Y dv.

(1вб)

V

 

V

V

 

Из равенства

(1 .6)

получаем дифференциальное

уравне­

ние переноса

тепла

 

 

 

 

й у

+ d ('if =

Оуу

( l.7 )

В соответствии с феноменологическим методом исполь­ зуем закон теплопроводности Фурье о дополнительной связи между переменными процесса: тепловым потоком

и температурой Ь • Согласно закону плотность теплового потока пропорциональна по величине и проти­ воположна по направлению градиенту температуры:

 

^ = - 1 у г а 4 Ъ ,

(1 .8 )

где Л

- коэффициент теплопроводности среды,

 

вт/м .град.

 

Закон

Фурье представляет собой простейшую форму об­

щего закона переноса потока энергии и является стро­

гим тогда, когда система

однородна во всех отношениях,

исключая наличие градиента температуры.

 

 

С учетом закона Фурье

уравнение (1 .7)

принимает

вид

a v =

d e lr ( j l c p a d t ) - » - •

(1 .9 )

 

 

 

 

 

13

Прменяя оператор Гамильтона, последнее выражение пе­ репишем в виде

 

 

 

a v = V

(ЛVt)+

tyv *

(I.IO )

 

В соответствии

с

правилами векторного анализа

 

 

v ( J l v t ) = ? A ‘ V t + A

 

 

Из (1.10) получаем

 

 

 

 

 

 

 

 

 

Qv = v ] [ v t + X v z t +

 

( i . i l )

В

декартовых

координатах

 

 

 

 

 

 

уД *

A =

ЭЛ г

ЭЛ -

+

ЭЛ

 

 

-г—

L + —

j

Эх

 

 

 

 

 

 

д х

д у °

 

 

 

v t

= агос/Ь —

дЬ

т

д Ь

-

 

at

 

 

 

 

 

 

Эх

0 + Эу ^

+ Э*

 

 

 

v H

lit

. э ч

 

a*t

 

 

 

 

3 acl

 

Э и z

*

д Х г

 

 

 

 

 

 

 

 

Тогда выражение ( I . I I )

в декартовых координатах име­

ет

вид

 

 

 

 

 

 

 

 

 

 

а

_ М _

at

ЗА

Э £

 

<LiL.

 

 

 

+ ЭЧ + 3 4 )

a v -acc

дх

 

 

 

Эх

Эх

 

Эзс? Эиг

9z V + 9v-

 

 

 

 

 

 

 

 

 

 

*

(I.I2)

 

Для определения

величины

Gy

используем

первое

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

сов закон

сохранения и превращения

энергии.

В соот­

ветствии с

э т м за малый промежуток

времени

dJv

приращение полной энергии единичного объема движу­ щейся среды равно подведенному теплу GL и внешней работе, совершаемой над средой:

14

p ( d v + A d - ^ 2) = a v ^ + A L v c/<T>

( 1 . 13 )

где Ly

- работа внешних сил над единицей объема сре­

U

 

ды га единицу времени#

 

-

внутренняя энергия единицы массы среды#

J\

-

тепловой эквивалент механической работы#

jD

-

плотность

среды#

 

ИХ

-

скорость

двниения среды.

 

Полная энергия потока складывается из внутренней

энергии

JDU и кинетической энергии

.

Внутренняя энергия среды связана с энтальпией из­

вестным термодинамическим соотношением

 

 

 

I = I/ + A pif

 

или,

в дифференциальной форме,

 

 

 

d V = d i - A d ( р у ) ,

(I.M )

где

р

- давление#

 

 

1Г

- удельный объем среды.

 

Для

идеального газа

 

 

 

di = с р

(1 .15)

где ср - изобарная удельная теплоемкость среды, ко­ торую в дальнейшем изложении будем прини­ мать постоянной величиной. Тогда (I .I4 )

можно записать в виде

d U = сР d t - A p d i f - A v d p .

(1.16)

Подставляя ( I . 16) в ( I . 13), после деления на

d

получим

 

15

Используя полученное выражение, уравнение

( I . I I )

за­

пишем в виде

 

 

 

 

 

г

vA ■v t ♦Av>t * nv*ft (I*

 

 

 

 

 

 

 

,

.

 

 

(I.I8)

Проанализируем уравнение ( I . 18). В это уравнение

входят

четыре неизвестные величины

Ь

,

р

,

иГ

и J3

. Следовательно, для общего

решения

задачи

о

теплообмене (наховдение поля температур

и тепловых по­

токов) в движущейся среде необходимо решить гидроди­

намическую задачу о распределении

 

скоростей

и давле­

ний, т . е . уравнение (I .I8 ) необходимо замкнуть диффе­ ренциальным уравнением сохранения количества движения, уравнением сохранения массы и уравнением состояния.

Последнее уравнение необходимо для связывания термо­

динамических

параметров

среды.

Уравнение (I.т я )

содер­

жит полные,

или

субстанциальные,, производные

от ве­

личин р ,

, t ,

и /л .

Субстанциальная

произ­

водная характерна для движущейся среды и обусловлена зависимостью рассматриваемой величины от координат и времени.

Полное изменение величины -f равно

dt

- #

rf,r + М - °,х * Щ 'с,г

и

*

или

 

°

 

 

ell

дф

д£ olx М du ,

dA

(т.тя)

Первый член правой части характеризует скорость изме­ нения рассматриваемой величины ^ во времени в той точке пространства, в которой элементарный объем на-

ходится

в данный момент времени, и поэтому называется

локальной производной.

Для движущейся среды производ­

ные

 

 

 

,

dj&r

являются

компонентами

вектора

скорости

 

иУ

вдоль

декартовых

координат,

т .о .

 

 

 

 

 

 

 

 

 

 

 

d$

ЪЪ

 

yj-

+ u r M +>,rM -

 

(1. 20)

 

d u

 

 

 

 

x doc

'by

*>

 

 

где

,

uf*

,

uf^

-

компоненты вектора

иУ .

Соотношение

(f.20) можно

записать в векторной форме

 

 

 

 

 

 

 

 

 

 

 

а -г1)

величина ufv<ft

характеризует изменение величины ^

из-за перемещения рассматриваемого элемента вещест­

венной

среды из

 

одной точки пространства в другую

и поэтому называется конвективной’составляющей пол­

ной производной.

 

 

 

 

 

 

 

 

 

Таким образом,

полное

изменение

субстанции

^

в единицу времени обусловлено локальным изменением

со

временем

 

 

и конвективным

переносом Ы v-f .

В

литературе

для

 

обозначения

субстанциальной

произ­

водной

применяют

 

специальный символ

 

 

 

С учетом- (I.2T ) уравнение переноса тепла в вещест­

венной среде (ТЛЯ) запишется

в векторной форме

di*(bpadt)+qY+A(Lv+

 

 

 

 

 

» ( Ф )

 

 

 

 

 

 

 

 

 

 

 

dec

или в проекциях

 

 

 

 

 

 

 

U . 22)

на прямоугольные оси координат:

 

 

 

W

 

<эл

Ъь

‘дк'Ъь

 

<fc-t

 

 

 

bij2-

Ъ%у + Ъх,

Ьх

ду

Ъ у +Ъ%

 

 

2,

зак.

 

 

Г "

 

 

 

 

17

 

 

 

 

 

 

 

 

 

 

Рассмотрим некоторые частные случаи уравнения пере­ носа тепла.В общем виде решить уравнение (1.22) сов­ местно с замыкающими дифференциальными уравнениями движения, сплошности, состояния при заданных началь­ ных и граничных условиях не представляется возможным. Однако ири определенных допущениях и анализе порядка величины отдельных членов можно значительно упростить уравнение (1.22) и привести его к виду, поддающемуся аналитическому решению.

потоках с большими скоростями необходимо учитывать изменение давления и кинетической энергии, вклад кото­ рых соизмерим с внутренней энергией. Однако при уме­ ренных скоростях потока среды можно пренебречь работой внешних сил и кинетической энергией, поскольку их вли­ яние на изменение теплосодержания (.энтальпии) потока мало. Работа в данном случае связана с изменением дав­ ления и внутренним трением двиз^тщейся среды. Пренебре­ жение работой внешних сил означает, что изменение дав­ ления и диссипация энергии трением дают относительно малый вклад в изменение внутренней энергии элементар­ ного объема текущей среды. В этом случае уравнение переноса тепла значительно упрощается и принимает вид

\

18

\

В целом ряде задач по теплообмену можно считать ко­ эффициент теплопроводности Л постоянной величиной, что не будет вносить большой погреиности в окончатель­ ный результат. Тогда уравнение (1 .24) принимает еще более простую форму:

А

- j> Ср

(1.25)

'р'сли среда неподвижная,

в частности твердое

тело,

то

перемещением частиц среды можно пренебречь

( и / -

о ),

и тогда конвективные члены в правой части уравнения

(1.25)

исчезают. Уравнение (1.25) в этом случае при­

нимает

вид

 

 

Л V z t + ^ v = p c p - f |r •

(1.26)

Уравнение (1,24) при отсутствии внутренних источников тепла, умеренных скоростях течения и постоянных физи­

ческих

свойствах

среды (

Л =

const

)

принимает вид

 

 

^

p v H

*

 

I

с1-27’

или

 

 

 

 

 

 

 

 

 

 

 

 

2>t

 

 

 

 

 

C tV zb

- c f r

 

(1.28)

где

а

-

- коэффициент температуропроводно­

 

 

Р с ,

сти.

 

 

 

 

 

 

 

линейным дифференииальниы

Уравнение

(1.28)

является

уравнением второго порядка в частнж првизводних па­ раболического типа и называется уравнением Фурье-Ос- троградского конвективного теплообмена.

Параболические уравнения теплопроводности еиясмваит

неооратимый процесс распрестраиения тенла. Эта следу­ ет из уравнения (1 .2 8 ), где замена времени иа-Т

1 9