Файл: Козырев, А. П. Теория тепловых и гидродинамических процессов в атомных энергетических установках учеб. пособие.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