Физика вод суши by Винников С.Д., Викторова Н.В. (z-lib.org)
.pdf-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 T— dx dz dx, |
|
|
|
|
|
|||
|
|
dy |
|
Г |
5/ |
, Л |
|
|
|
|
|
я |
(6.5) |
||||
|
( |
m |
t н---- dy |
|
|
|||
e ;= |
Л |
. |
Qy |
. |
|
dxdzdx, |
||
|
I |
dy |
|
|
dy |
|
|
|
|
|
dt |
|
|
|
|
|
|
Q'5 = -X T— dxdydx |
|
|
|
|
|
|||
|
|
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
