книги из ГПНТБ / Козырев, А. П. Теория тепловых и гидродинамических процессов в атомных энергетических установках учеб. пособие
.pdfния. Для теплотехнической системы в целом получаемые при этом уравнения баланса, или интегральные уравне ния сохранения, не позволяют исследовать внутренние процессы. Детальное изучение внутренних процессов в системе становится возможным лишь при переходе к диф ференциальной форме уравнений сохранения, называемых в этом случае уравнениями переноса. Дифференциальные уравнения переноса, как правило, включают большее ко личество неизвестных, чем число уравнений. Это объяс няется тем, что используемые законы физики описывают макропроцессы и не учитывают дискретный характер стро ения среды и микроскопические особенности процессов. Для замыкания системы дифференциальных уравнений и последующего их интегрирования на помощь приходят за коны о дополнительных связях между переменными про цесса. Такие законы позволяют считать среду, где про текает процесс, сплошной и непрерывной и дают возмож ность применять статистические свойства среды (давле ние, температуру и т .д .) к ее дифференциальным эле ментам. Полученные из наблюдений за явлениями приро ды указанные законы иногда в литературе[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, |
зак. |
7д |
|
|
Г " |
|
|
|
|
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