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

Физика вод суши by Винников С.Д., Викторова Н.В. (z-lib.org)

.pdf
Скачиваний:
23
Добавлен:
19.02.2024
Размер:
9.8 Mб
Скачать

-X,

dt

(5.84)

тdz

 

Пользуясь здесь тем же приемом, что и при выводе критерия Фурье, введем следующие соотношения:

Х„ „ :

^Т,М

^Х,Ат,Н> К

mfcH,

 

 

 

(5.85)

zM= miz

L

KP, u =mL L4»»’ Pm=wpPh; *м = щ (н-

Относя уравнение (5.84) к натуре, а затем к модели, мо­

жем написать:

W р L —— = X —

'тнг'нА'кр,н ,HQz

 

-0

^ МРм4~ » —

dL

= Х„

^ ' м Эх.

dz.,

т , н

_^Н

(5.86)

& и

 

 

'Sz„

(5.87)

 

 

Заменив в последнем уравнении все величины через их вы­ ражения по (5.85), найдем:

mwm mL m,

P 4

щтпх

qz

Qt

 

 

н ^

= mx Xwh^ l

- ^ A , h

. (5.88)

H “ H K p ,H о

Л „ M «H л

 

dxH

dzH-o

dztI +o

 

При условии, что mx - m x^- mx, уравнение (5.88) упроща­

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

.2

 

= . а ,

 

 

mWmpmL<™l

^нРнАф.н

-X ,

(5.89)

mtmxmx

м,н dz„

 

 

Решая совместно уравнения (5.89) и (5.86), получаем:

 

mwmPmr Щ / Щ mxmx =1 ■

(5.90)

Отсюда, согласно соотношениям (5.85),

 

161

- 'мтА >

1»ЬХ» = idem .

(5.91)

^ mPmW m ^ .Р А ^ А

 

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

ла, что и натура, и при одинаковой влажности, т. е. при

mw = 1,

mL = 1, тр = 1, тх = 1., уравнение (5.90) значительно упрощается

и приобретает вид

 

rnf/mt mx = 1.

(5.92)

Кроме этого требования на модели должно быть выполнено требование критерия Фурье, согласно которому при принятых

выше условиях

 

m jm 2 = 1.

(5.93)

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

Забегая несколько вперед, отметим, что если выполняется те­ пловое моделирование подвижной среды - потока в реке, водохра- нилище-охладителе, конвективных потоков в водоеме и воздухе, окружающем исследуемый объект, то для вывода соответствующе­ го критерия подобия должны воспользоваться уравнением энергии (см. главу 6, формулы (6.10) или (6.13)). Тогда, по аналогии с крите­ рием Fo (5.75), получим критерий Пекле:

 

=

Ре,

(5.94)

° н

« М

а

 

где vM, /м, ам и vH, /н, ап - соответственно скорость, характер­

ный размер потока и коэффициент физической (турбулентной) температуропроводности на модели и в натуре.

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

Re = — .

(5.95)

v

 

162

Если теперь воспользоваться полученными критериями по­ добия Bi, Ре и Re, то можно составить комплексные критерии, ко­ торые хотя нового ничего не создают, но иногда более удобны в исследованиях, например,

Ре

у

(5.96)

— = — = Рг - критерий Прандтля

Re

аТ

 

и

 

 

 

Bi = У1 (Re, Рг) = ,

(5.97)

откуда

 

 

 

a = y / 1(Re,Pr).

(5.98)

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

критерий Грасгофа:

Щ ^ = Ог,

(5.99)

v

где At - характерный температурный напор (разность температу­ ры поверхности воды и воздуха на удалении), (3, - коэффициент объемного расширения жидкости (воздуха). Если же свободная конвекция вызвана разностью плотности в двух точках жидкости на одной вертикали, обусловленной как различием температуры, так и другой какой-либо причиной, то получим критерий Архиме­ да, используемый в качестве критерия устойчивости частиц жид­ кости (воздуха) в неоднородной среде,

g£z^ 4 = Ar>

(5л°°)

Ро

v

 

Р-Ро

~

-

где g - —— - ускорение подъемной силы, действующей на тело

Ро

 

 

(частицы жидкости) с плотностью р; р0 -

плотность окружающей

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

163

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

а /

Bi = — = / 2(Gr,Pr) = N u. (5.101)

А

Произведение критерия Gr на критерий Рг называют крите­ рием Рэлея, используемый в качестве критерия термической неус­ тойчивости в среде:

Ra = Gr-Pr = ^ ^ //- ,

(5.102)

 

va

 

 

где At - разность температур на границах горизонтального слоя

 

воды толщиной /; обычно принимают / = Н, т. е. равной глубине

 

водоема.

 

 

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

 

нике, помимо рассмотренных применяют также другие критерии

I

моделирования, соответствующие определенным задачам. С ними

|

можно познакомиться в специальной литературе [16, 22, 35, 59].

 

Рассмотрим пример использования критериальных зависи-

!

мостей для решения практической задачи. Пусть нам требуется,

 

например, рассчитать теплопотери с водоема в условиях свобод­

 

ной конвекции в атмосфере над ним (при штиле и ветре до 2 м/с).

 

Для этого воспользуемся известным в теории теплопередачи соот­

 

ношением между критериями Нуссельта и Рэлея

 

 

Nu = ,4Ra1/3,

(5.103)

 

гдеЛ = 0,14.

 

 

После подстановки соответствующих выражений для Nu и

 

Ra и некоторых (опущенных здесь) выкладок сотрудниками Физи­

 

ки атмосферы и океана Г.С. Голицыным и А.А. Грачёвым было найдено, что в режиме свободной конвекции в приводном слое воздуха:

1) явный поток теплоты (п. 3.4)

2) скрытый поток теплоты (п. 3.6)

 

Р*Г*

1/3

4 / 3

(5.105)

а = л А ,р вд?

 

 

 

W2

где ср - теплоемкость воздуха при постоянном давлении; рв -

плотность воздуха; а - коэффициент теплового расширения воз­

духа; kt и кЕ - коэффициенты теплопроводности воздуха и диф­

фузии водяного пара; v - кинематический коэффициент вязкости воздуха; т « 0,075; Во - число Боуэна (9.7); Ьи - удельная тепло­

та испарения; At я Aq - разность значений температур и удель­

ной влажности у поверхности воды и на удалении, т. е. на верхней границе пограничного слоя атмосферы. Метеорологи и гидрологи за верхнюю границу приводного пограничного слоя атмосферы принимают высоту, равную 2 м, океанологи - 1 0 м; это высоты, на которых ведутся стандартные метеонаблюдения.

165

ГИ Д Р О Т Е Р М И Ч Е С К И Й Р А С Ч Е Т В О Д О Е М О В

ИВ О Д О Т О К О В

До сих пор мы рассматривали задачи, связанные в основном с изучением распределения теплоты в неподвижных средах.

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

6.1. Д ифференциальное уравнение теплопроводности турбулентного потока

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

ражающего частный случай за­

 

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

 

энергии. Оно может быть запи­

 

сано для произвольного объема

 

воды изучаемого водотока.

 

Итак, выделим в пределах

 

водотока в системе декартовых Q,+Q

 

координат х, у, z

элементарный

*■ х

параллелепипед

с

гранями dx,

 

dy, dz (рис. 6.1). Рассмотрим его

 

тепловой баланс. Через грани

 

параллелепипеда

теплота будет

 

распространяться двумя путями:

1)вместе с водными масРис­ . 6.1. Схема к выводу диффе­

сами, пронизывающими грани параллелепипеда со скоростями vx, v , vz - молярный перенос;

ренциального уравнения тепло­ проводности потока жидкости.

166

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

Уравнение теплового баланса для выделенного элементарно­ го объема жидкости в этом случае будет иметь следующий вид:

й + б г + бз + 6 4 + 0 5 + 6 б + 6 1 + 6 2 + 6 3 + 6 4 + 6 5 + 6 б = 6 7 > (6 -1 )

где 6 1 > 6 2 j 6 з и Т-Д- - количество теплоты, обусловленное скоро­ стью потока жидкости через соответствующие грани в направле­

нии осей х, у, z за время dx, a <2i, 6 2

> 6 3 и Т-Д- ~ количество теп­

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

теплопроводностью потока

через эти же грани и за то же время dx.

В том случае когда потоки теплоты, проходящие через грани параллелепипеда, взаимно не компенсируются, т. е. в него входит теплоты больше, чем выходит, или, наоборот, будет наблюдаться изменение энтальпии рассматриваемого объема dxdydz, которое

в уравнении (6 .1 ) обозначено через Qj. Определим составляющие уравнения (6.1).

Количество теплоты, поступившее в параллелепипед через

грань dy dz молярным путем за время dx, оценим по формуле

 

Qx= ср vxt dydzdx,

(6.2)

где с и р - удельная теплоемкость и плотность жидкости; vx- про­ екция скорости на ось х; р vxdy dz - расход жидкости через грань

параллелепипеда dy dz; t - температура жидкости, проходящей че­ рез грань dy dz.

Количество же теплоты, выходящее из элементарного па­ раллелепипеда через противоположную грань, отстоящую от пер­ вой на расстоянии dx,

 

dv.

dt

(6.3)

Q2 -~ cp vx н---- -d x

t-\---- dx dydzdx,

 

дх

dx

 

где dvx/dx и dt/dx -

изменение скорости и температуры жидкости

внутри выделенного

объема вдоль оси х. Знак минус в этом урав-

нении свидетельствует о том, что Q2 - уходящее из элементарного

параллелепипеда количество теплоты.

Для остальных граней параллелепипеда будем соответствен­ но иметь:

Q3 - ср v t dx dz dx,

dv„ . Y

Q a = ~CP vy + -^ ~ dy dy

Qs =cpv2t dx dy dx,

dv Qb = ~CP v. + —- dz

z dz

Qt

Л

 

t +~~dy

dxdzdx,

 

. дУ .

(6.4)

 

 

dt ,

dxdydx.

 

t н---- dz

 

dz

Другие шесть слагаемых уравнения (6.1) (Q[ , Q'2, Q3, Q'4,

Qs > Qc, )= обусловленные турбулентной теплопроводностью, опре­ делим по следующим формулам:

Q[ = -Х тdy dz dx, dx

 

 

 

d

 

dt

 

 

 

 

.

дК

t +— dx

 

 

02 =

л

 

dx

 

dy dz dx,

Хт+-—-dx

 

dx

 

 

 

dx

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

Q3 = -X Tdx dz dx,

 

 

 

 

 

 

 

dy

 

Г

5/

, Л

 

 

 

 

я

(6.5)

 

(

m

t н---- dy

 

 

e ;=

Л

.

Qy

.

 

dxdzdx,

 

I

dy

 

 

dy

 

 

 

 

 

dt

 

 

 

 

 

 

Q'5 = -X Tdxdydx

 

 

 

 

 

 

 

dz

 

 

dt

 

 

 

 

 

 

d

 

 

 

 

 

 

 

t н---- dz

 

 

 

i

8k

^

.

dz

J

dx dy dx,

 

А,т н------ dz

 

dz

 

 

 

dz

 

 

 

 

 

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

Изменение энтальпии рассматриваемого объема Q1 опреде­

лим по формуле

168

 

 

 

 

Q1 = cp— d xd yd zd x.

 

 

 

(6.6)

 

 

 

 

 

дт

 

 

 

 

 

 

 

 

 

Решая совместно уравнения (6.1) -

(6.6), получаем [28]:

 

 

/

dt

dt

dt

dt

 

 

f t

d2t

d2t

\

 

 

ср — + — + и — + u7 —

 

 

dx1

dry

dz

+

 

 

5т:

dx

dy

dz

 

 

 

(6.7)

 

 

 

 

r dkdt

dk

dt

dk

dt ^

 

 

 

 

 

+

 

 

 

 

 

 

---—---+ ----T---+ ----3L---

 

 

 

 

 

 

 

dx dx

dy

dy

dz

 

dz

 

 

 

 

При совместном решении уравнений (6.1) - (6.6) учтено ус­

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

 

 

 

 

 

 

 

dvx/dx +dvy/dy +d vjd z = 0

 

 

(6.8)

и

отброшены

слагаемые

 

 

 

dx d yd z,

dvy & j j

2 j

 

dx

—- — dxdy

dz ,

 

 

 

 

 

 

dx

 

 

dy

dy

 

dv, dt dxdydz

,

а также

dK

d2t dx2dy d z ,

dK

d2t dxdy2 d z ,

dz

dz

 

 

 

dx

dx

 

 

dy

dy2

 

dk^

d2t

,

, 2

 

 

 

 

 

 

 

 

у

 

— - — jd x dydz

из-за их малости по сравнению с другими. Урав-

dz

dz

 

 

 

 

 

 

 

 

 

 

 

 

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

При постоянном значении коэффициента турбулентной теп­ лопроводности А.т для всего потока уравнение (6.7) примет вид:

dt

dt

dt

dt

к т/ d2t

d2t

Э2Л

— + ux — + u v — + u

= —

dy2

(6.9)

дт

dx

dy

dz

cp dx2

dz2

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

169

Учитывая, что левая часть уравнения (6.9) - полная произ­ водная от температуры по времени, его можно представить в виде

dt

dx ~ а' дх2

ду2

dz2

(6.10)

 

или

 

 

 

~

= aTV

2t ,

(6.11)

dx

 

 

 

где а т=А,т/(ср) - коэффициент турбулентной (конвективной) тем­ пературопроводности.

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

dt

d2t

d2t

д21Л+ W/(cp),

(6.12)

dx

удх2

ду1

dz1 j

 

где W - интенсивность внутреннего источника (количество теплоты, которое выделяется или поглощается единицей объема жидкости).

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

— = aS/2t ,

(6.13)

dx

 

где а = V(cp) ~ коэффициент температуропроводности жидкости.

В случае установившегося температурного режима водного потока температура в каждой точке его остается неизменной во времени (dt/dт = 0) и меняется лишь по направлениям х, у, z,

а уравнение (6.9) принимает следующий вид:

170