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

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

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

где arxi , - проекции на ось сс пульсационных добавок скоростей, измеренных одновременно в точ

ках I и 2, расположенных на

расстоянии ц,

друг

от

дру

га.

Зависимость

коэффициента

корреляции^от

расстояния

а-

показана

на рис. 1.3.

Масштаб турбулентности

оп­

ределяется интегрированием

Я (^ ) от О до

<>=

 

©о

Рис. 1 .3 . Зависимость коэффициента корреляции от расстояния между точками

Отдельный вихрь в параллельном потоке условно зани­ мает область диаметром 2 L .

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

40

Уравнения движения осредненного турбулентной!

потока

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

слоями макроскопическими объемами

жидкости.

оси ос

Рассмотрим уравнение движения в

направлении

(1*49). Подставим в это уравнение

истинные значения

скоростей и давлений, выраженные в

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

с(1 .б7)

через их средние и пульсационные значения, и произведем

осреднение уравнения по времени.

Тогда получим уравнение

осредненного турбулентного потока по оси

х

 

 

 

d/Ь

ЭР +

v z-ur -

Р дх

 

 

)

 

р дх

 

 

 

 

 

 

± д _

 

- I J - l o w

х. W

г

 

 

Р э*

f

 

 

р d z ly

:

 

 

 

 

 

 

 

(1.69)

Аналогичным образом получаются уравнения осредненного

турбулентного потока в направлении осей

у,

и

%

Эти уравнения были получены Рейнольдсом.

 

 

 

Как видно,

получено

3 x 3 =

9 дополнительных

членов

вида

р-и/г-ик . Эти члены имеют размерность напряжений,

Зак.

t

 

 

 

 

 

4 1

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

Действительно, рассмотрим

простейший случай пульса-

ционно установившегося

потока между двумя

плоскими бесконечно длинными пластинами, осредненное

движение

которого параллельно оси

сс

,

а

скорость

является

функцией только

координаты

%

,

т .е .

и ^ = 0 ;

< ^ = 0 ;

Эагх

 

 

 

 

 

 

= 0 :

 

 

 

=0.

Тогда получим

Эх

 

 

 

 

 

 

 

 

 

 

 

 

 

0= X

j_

Эр

с Э У Я

1

Э_

 

 

 

Р

Эос

Эх1

Р

Эзь pur^ k i . 7 0 )

 

 

Осредненное уравнение в напряжениях в этом случае имеет вид

 

 

- _ 1 .Э р

_

 

 

 

 

Х ~ Р

Эх

f Э*

( I *7I)

где

-

осредненное

касательное напряжение. Сопо­

ставляя уравнения (1.70)

и (I.7 I), ползгчим

 

 

 

 

X

 

 

(1.72)

 

 

oL%

 

 

 

 

 

 

 

Как

видно из

полученного

выражения,

полное

касательное

напряжение в осредненном турбулентном потоке можно пред­ ставить как сумму вязких напряжений, определяемых по закону Ньютона, и дополнительных турбулентных напря­

жений!

V * Ъ + ТГ

(1.73)

 

 

42

Аналогично закону Ньютона о внутреннем трении мож­ но положить, что турбулентное напряжение пропорцио­ нально производной от осредненной скорости по нормали к направлению потока, т .е .

 

 

V J 4- d x

 

(1.74)

 

 

 

 

Тогда

 

ч d-uJi

 

, dur.

-

,

,

 

 

 

 

(1.75)

где

-

величина, называемая

коэффициентом

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

так как зависит

от степени

перемешивания (турбулент­

ности) потока и

координат.

Величина j\(%)

может

быть определена экспериментальным путем по найденному

распределению

-от.

от

%

. На разных расстоя­

ниях от стенки значение слагаемых в (1.73) различно,

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

рости небольшой,

вязкие напряжения малы по сравнении

с турбулентными,

j^ T» J4

.

Вблизи стенки

унт « ум ,

так как поперечные составляющие скорости пульсаций

имеют ничтожно малую величину.

Поскольку чг

не яв­

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

43

Осредненн.ые уравнения переноса тепла

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

=

t= t + 0; р = р+-р'; р -=р + р'.

(1.76)

Для несжимаемой жидкости при умеренных скоростях и

постоянных физических свойствах среды (

c o a s t )

уравнение переноса тепла фурье-Остроградского имеет

вид ( о

полагаем равным нулю)

 

 

 

 

 

(1.77)

Подставляя в это уравнение вместо актуальных значений скоростей и температур их значения в соответствии с (1.76) и производя осреднение, получаем

где

d-biT \>у '0 -

Уравнение (1.78) отличается от обычного уравнения Фурье-Остроградского наличием дополнительного члена

ср р ou,ir-и/В , имеющего размерность удельного теплового потока [ккал/м3* ч] и характеризующего перенос тепла

врезультате действия турбулентных пульсаций скорости

итемпературы. Можно показать что плотность теплового потока [ккал/м3* ч] , обусловленного турбулентным

переносом тепла, равна

^ = - СрР и г '0 .

(1.80)

По аналогии с законом Фурье, вводя в рассмотрение ко­ эффициент турбулентной теплопроводности _ДТ , запииеы

 

 

 

о Л 1 '

 

 

( х . м )

Для плоского потока,

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

градиенты

скоростей и температур много больше продольных,

имеем

 

 

= -с оигге=-

 

д ±

 

(1.82)

 

 

■>т д%

 

 

%

 

 

 

Тогда

коэффициент турбулентной теплопроводности

опре­

деляется выражением

 

 

 

 

 

 

 

,

 

 

 

 

(1.83)

Соответственно

 

д1/дъ '

равная

коэффициен­

величина, численно

ту турбулентной теплопроводности

]\т в некоторой

точке,

деленному

на

объемную теплоемкость

жидкости

р с 0

в этой точке, называется

коэффициентом тур­

булентной температуропроводности.

45

 

 

(1.84)

Следует отметить, что уравнение

( I .8 I ) не является за­

коном, поэтому величины

, а т

не являются физичес­

кими параметрами среды, а зависят от тех же характери­ стик, что и J^T(0T) •

Уравнение, выражающее количество тепла, переносимое в жидкости за счет процессов теплопроводности и турбу­ лентного перемешивания, будет записано в виде

(1.85)

Из уравнений (1.74) и 1(.82) видна тесная связь между процессом турбулентного переноса тепла и количе­ ства движения.

§ б. Условия однозначности для процессов теплопереноса

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

Геометрические условия определяют, форму и размеры твердого тела, в котором протекает процесс, располояе-

46

нив поверхности нагрева в потоке

жидкости,

границы по­

тока жидкости.

 

 

 

Физические условия задают значения теплофизических

параметров

среды О

, ср , р

г м ),

в которой

протекает

процесс.

 

J

 

Совокупность временных и граничных условий называет­ ся краевыми условиями.

Временные краевые условия называются начальными ус­ ловиями и определяют условия протекания процесса во времени (стационарность или нестационарность, началь­ ное распределение температур и т .д .) . Во многих зада­ чах теплопроводности принимается равномерное распре­ деление температуры в начальный момент времени, т .е . при ъ = 0 -fc = c o rc s t. Начальное условие мо­ жет быть задано некоторой функцией, определяющей за­ кон распределения температуры в теле в начальный мо­ мент времени, т .е . при b~t0(х, f, х,0).

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

Для стационарного процесса задаются только гранич­ ными условиями, поскольку процесс от времени не за­ висит.

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

4 7

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

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

 

V V

•(1-ж>

На практике и в математической физике рассматриваются

три типа граничных условий.

 

I.

Граничное условие первого рода задает

распределе­

ние температуры по поверхности тела в любой момент вре­

мени,

т .е .

 

Требуется

найти функцию

i (

ос, и^Х/с) ,

удовлетворяю­

щую внутри области основному уравнению

и принимающую

на границе

заданное значение

f (т;,эс( l^

x) . В матема­

тической физике такая задача называется первой краевой задачей, а в случае стационарного процесса (уравнение

процесса - уравнение Лапласа

v H

= 0)

имеем зада­

чу

Дирихле. На практике такие задачи встречаются редко

и,

как

правило,

при интенсивной

теплоотдаче

на поверхно­

сти

( о б - о»)

, когда температура стенки

близка к тем­

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

поддержания постоянной температуры

стенки.

 

 

2.

Граничное условие второго

рода (втоцая краевая

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

ИВ

(1. 88)

Решением уравнения будет функция » Удов­ летворяющая внутри области основному уравнению, нормаль­ ная производная которой на границе удовлетворяет урав­ нению (1 .88).

В простейшем случае плотность теплового потока при­

нимается

постоянной (a

^ c o r ts t)

.

Такой случай

теплообмена встречается

на

практике

в

задаче нагрева

системы

внешним источником

при помощи излучения, когда

температура тела значительно меньше температуры излу­ чающей поверхности (топка котла, высокотемпературная печь, солнечная батарея и т .д . ) .

Для уравнения Лапласа вторая краевая задача называет­ ся задачей Неймана.

3. Граничное условие третьего рода характеризует закон конвективного теплообмена между поверхностью тела и потоком жидкости (третья краевая задача). В со­ ответствии с законом Ньютона количество тепла ^ , пе­ редаваемое в единицу времени через единицу площади по­

верхности

тела, прямо пропорционально разности темпера­

тур

поверхности тела и жидкости:

 

 

(1.89)

где

оС

- коэффициент пропорциональности, называемый

коэффициентом теплоотдачи и характеризующий интенсив­ ность переноса тепла от поверхности нагрева к потоку яидкости, вт/м2» град.

С учетом уравнения Фурье (1 .8 ) граничное условие (1.89) принимает вид Ч зак. 7д