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