![](/user_photo/_userpic.png)
книги из ГПНТБ / Палагин, Э. Г. Основы гидромеханики учебное пособие для метеорологов
.pdfКроме коэффициента вязкости р, большое применение в теории вязкой жидкости находит так называемый кинематический коэф фициент вязкости v. Под этим понимают отношение коэффициента вязкости жидкости к ее плотности:
Как известно, ц характеризует собой величину касательного напряжения вязкости, создающегося при данном поле скорости. Пели поставить вопрос о том, какое ускорение может вызвать это напряжение, будучи приложенным к одной из гранен единичного кубика жидкости, то для получения ответа следует р > разделить
на плотность жидкости. В результате ускорение оказывается рав-
у |
dv |
dv |
„ |
образом, кинематический коэффициент |
ным - |
-г" —- ■' -г- . |
1акпм |
||
о |
on |
|
|
|
вязкости характеризует собой ускорение частиц жидкости, вызы ваемое напряжением вязкости.
V см 2/сек
700
720
740
760
780
|
|
[|J-] |
СМ1 |
|
|
|
|
|
|р| |
сек |
|
|
|
Приведем значения v для воздуха, |
взятого |
при |
давлении |
|||
760 мл и при различных температурах: |
|
20 |
30 |
|||
Т [°С1 . . . |
20 |
— 10 |
о |
10 |
||
•у [ем*1сек] . . . |
н,11Л |
0,121 |
0,130 |
0,139 |
0,149 |
0,170 |
100
Как видно нз этих данных, с ростом температуры v возрастаем
причем эта закономерность имеет |
место у всех газов. Напротив, |
|
с ростом давления v убывает Ма |
рис. 8. 2 показана зависимость |
|
кинематического |
коэффициента вязкости воздуха от температуры |
|
п давления, а на |
рис. 8.3— зависимость v от высоты над уровнем |
|
моря в условиях (Международной стандартной атмосферы. |
||
Из формулы |
(8.4) имеем следующее выражение для кинемати |
|
ческого коэффициента вязкости: |
|
Отсюда и следует, что с возрастанием температуры, вызыва ющем увеличение с, v возрастает, а с возрастанием давления, вы зывающем уменьшение /, v убывает.
Рис. 8.3
Приведем значения v для воды, взятой при различных темпе ратурах:
7 ГС] . . . |
О |
К) |
Д) |
40 |
•' \см'Цсек\ . . . |
0,0181) |
0,0130 |
0,0101 |
0,0060 |
Сравнение с данными для воздуха обнаруживает, что у воды v значительно меньше, чем у воздуха, хотя р у воды и больше. Это связано, очевидно, с тем, что плотность воздуха в несколько сот
101
раз меньше плотности воды, и кинематический эффект от прило жения напряжения вязкости в случае воздуха оказывается боль шим, хотя сами напряжения здесь являются меньшими.
Убывание v с возрастанием температуры объясняется темп же. причинами, что и соответствующее изменение р, поскольку влияние температуры па плотность жидкости, очевидно, является весьма малым в сравнении с ее влиянием на коэффициент вязкости.
§ 2. УРАВНЕНИЕ НАВЬЕ — СТОКСА
Ввиду малой изменчивости р в интересующем нас диапазоне температур (для большинства метеорологических задач -50°С-ь -j-50°C), в дальнейшем будем полагать p = const, что по зволяет выносить ее из-под знака дифференцирования.
1. |
Запись |
составляющих |
силы ---div А |
Подставляя |
вначале |
||||||||||||
формулы (8.3) в выражение для |
|
Р |
|
|
|
|
|
|
|||||||||
проекции силы трения на ось х. |
|||||||||||||||||
получим ее запись |
через компоненты |
вектора скорости. |
Так, |
||||||||||||||
|
|
|
|
|
|
|
д Vх |
|
о |
in (1 iv v j = |
|
|
|||||
|
дх |
|
дх |
— р -г '2ц |
дх |
|
"3 |
|
|
||||||||
|
|
|
|
|
|
: |
> |
|
|
|
|||||||
|
|
|
|
|
|
d2v x |
|
|
|
|
О |
|
|
|
|||
|
|
|
|
|
+ 2^ |
|
- |
|
3 |
|
‘ |
|
div -v. |
|
|
|
|
|
|
|
|
|
|
|
|
д х |
|
|
|
|
|||||
Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L |
|
дР*х= |
1 |
дР |
I |
о.. д2А |
2 |
|
д div у |
|||||||
|
о |
|
дх |
|
|
р |
дх |
|
|
|
дх2 |
3 |
|
дх |
|
||
Далее: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
д_Рух_ |
|
|
i d v |
|
|
(ПщЛ _ ^ |
|
|
д-vy \ |
||||||
|
~р |
|
ду |
~~ * ду |
I |
дх |
|
' |
ду |
/ |
' |
I ду2 |
|
' дх ду) |
|||
|
1 |
д 1 \ х |
|
д |
( d v z |
i |
д v x\ _ |
( д2 vx |
^ |
d2vz \ |
|||||||
|
~j |
|
dz~~~^dz |
\ 'дх |
|
dz |
l ~ |
' \ |
dz2 |
|
dx dz) |
||||||
Объединяя результаты получим |
|
|
|
|
|
|
|
||||||||||
1дРхх |
|
д Руv и дР “ ) - |
|
1 |
op |
|
; d- vx |
i |
d2vx |
dz2 / |
|||||||
|
|
p |
|
дх |
J \ дх2 ~ |
dy2 |
|||||||||||
\ |
дх ' |
|
ду |
|
■ dz |
1 |
|
|
|||||||||
|
о |
(дгу |
, |
|
|
|
|
|
2 |
|
д |
|
|
1 |
др_ |
||
+ |
v дх |
\ |
дх |
1 |
ду |
|
az ) |
— -- v |
-— |
div v — ----- |
дх |
||||||
|
|
3 |
|
dx |
|
|
p |
||||||||||
|
|
|
|
|
+ |
'4 t ’>+ |
Т |
|
Ж |
|
d iv ,;’ |
|
|
(8.5) |
|||
|
|
|
|
|
|
|
|
|
|
где Avx — оператор Лапласа.
102
Проведя аналогичные выкладки, получи,vi
1 |
, а р ху , 0 Руу |
j- |
o P ,v\ |
|||
и \ |
дх |
dy |
dz |
|||
|
||||||
1 |
/ dPxl |
dPvz |
|
dP„' |
||
о |
\ |
dx |
dy |
|
dz |
1 |
dp |
vA v v |
\- |
л. |
Tx; |
flivu; |
(S.5'J |
|
о |
dv |
|||||||
|
|
|
|
|
|
|||
|
dp |
|
|
v |
C) |
cliv v. |
(8.5") |
|
|
— Н\'ЛУгН—-,y |
—y |
||||||
|
dz |
|
|
3 |
dz |
|
|
2. Различные формы представления уравнения Навье — Стокса. Подставляя (8.5), (8.5'), (8.5") в (7.21), получим уравнения дина мики вязкой жидкости в проекции на осп координат. Именно:
dvх |
dvx |
dvr |
|
dvx |
|
1 dp |
|
|
|
|
||
W + v |
Th |
+ V'J~di +Vz~dF==~ T ^ x |
+ x + vV v+ |
|
||||||||
|
|
|
|
a |
cliv t1, |
|
|
|
|
|
|
|
|
|
3 |
|
dx |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||
dvu |
dvi |
dv |
dVy |
1 |
dp |
+ ) |
+ v.Ac1,,-b |
|
|
|||
Tt |
dx |
v |
dz |
о |
dy |
|
( 8. 6) |
|||||
|
|
|
d |
cliv v, |
|
|
|
|
|
|
||
|
|
- y |
- y |
|
|
|
|
|
|
|
||
|
|
3 |
dy |
|
|
|
|
|
|
|
|
|
dv. |
dv, |
dv, |
|
dot |
|
1 |
dp |
* P'/. -!-vAtij-f |
|
|||
y ,----l-C’.v ~5-----1" ; 'i/ ~z— + |
c'j:" 3 — |
|
n |
dz |
|
|||||||
di |
dx |
■ dy |
|
dz |
|
|
|
|
|
|
||
|
|
v |
|
d |
|
|
|
|
|
|
|
|
|
|
+ T |
|
' aT |
‘Hv ”■ |
|
|
|
|
|
|
|
В тензорных обозначениях: |
|
|
|
|
|
|
|
|
||||
dvj |
|
dv, |
|
dp |
+ F, + V |
d-v, |
|
d-v, |
(8.7) |
|||
~dl |
Vj |
~ d x ,~ ~ |
|
dXi |
dxjdxj |
3 |
dxjdx; |
|||||
|
|
|
|
Умножая уравнения (8.6) соответственно на i, j и k и склады вая результаты, получим уравнение движения в векторной форме
dv
(v ■v)t' = grad p + F-\- yAl4- — orpaddivL'. (8.8)
dt
В итоге имеем запись второго закона Ньютона для жидких сред. Сила инерции dv равна сумме внешних массовых сил гра-
диента давления — grad/; и вязкости v | Дц-f-— graddivf J .
Напомним, что все они отнесены к единице массы.
103
—>
Для случая несжимаемой жидкости (div t>= 0) они упрощают ся и вместо (8.8) будем иметь
+ {v • v) v = ---- — grad p — F -f Ar. |
(8.9) |
|
( 1 ! |
'j |
|
Аналогичным образом |
перерабатываются (8.6) и |
(8.7). |
В заключение заметим, что наличие вязких членов в уравне ниях движения обуславливает необратимость процессов, проте кающих в жидкости при ее движении. Действительно, если мы по пытаемся вернуть систему в исходное положение, подставив вместо
/, .V, у, z, vx, vu, v: значения —t, —х, — у, —z, —vx, —ty —vz, то получающиеся при этом уравнения не будут тождественны систе ме (8.6), ибо перед вязкими членами появится знак минус.
ГЛАВА IX
УРАВНЕНИЕ ПРИТОКА ТЕПЛА
§ 1. ВНУТРЕННЯЯ ЭНЕРГИЯ СОВЕРШЕННОГО ГАЗА
Внутреняя энергия есть функция состояния, которая для слу чая однокомпонентной среды однозначно определяется двумя тер
модинамическими параметрами, например и— и(а,Т), где а = ------
Р
удельный объем.
Рассмотрим вид этой функции для совершенного газа. Как известно, эта модель предполагает, что молекулы совершенного газа имеют сферическую форму и взаимодействие между ними сво дится к столкновениям, подчиняющимся закону абсолютно упру гого удара. Определяющим фактором в этом процессе является кинетическая энергия молекул и, следовательно, температура. Поскольку других энергетических переходов между молекулами нет, расстояние между ними а, следовательно, и удельный объем на величине внутренней энергии совокупности молекул не отражается.
Ввиду того, что кинетическая, а, следовательно, для рассматри ваемого случая и внутренняя энергия, зависит только от темпера
ди
дТ объемной теплоемкости или, иначе, изохорной теплоемкостью.
Кинетическую энергию можно представить себе состоящей из энер гии поступательного, вращательного и колебательного движений молекул, в соответствии с чем и величина сг разлагается на три составляющих: cv= c uU-j- cvB-|- cvK. Как показывают эксперименты,
для газов, близких по своим свойствам к совершенным, из трех слагаемых лишь величина, связанная с колебаниями (с„к), зависит
104
от температуры, да и то, если 7>2000°К. К тому же она мала в сравнении с первыми двумя слагаемыми, так что ею, а с ней и
зависимостью |
cv от |
температуры, можно пренебречь, т. е. |
cv ~ с,,и + cva, |
причем |
с-о =const. Поэтому du = cvdT. |
Возьмем интеграл от обеих частей этого уравнения в пределах
ох 7 = 0 до 7 = 7, чему соответствует и— 0, и= и. Тогда |
получим |
u — cvT. |
(9.1) |
Теперь ясно, что изменение внутренней энергии эквивалентно изменению количества тепла в единице массы и связано с измене ниями ее температуры.
Атмосфера по своим свойствам, п это подтверждается много численными экспериментами, является совершенной смесыо совер шенных газов. В силу этого выведенное соотношение (9.1) приме нимо к ней в полной мере. В дальнейшем мы будем использовать именно это выражение для внутренней энергии тем более, что оно применимо и ко многим другим средам.
Молекулярный поток внутренней энергии для случая u — cvT представляет собой поток тепла, который будет определяться эмпирическим уравнением Фурье (о нем говорилось в вводной части):
Приток тепла:
div ■ |
,. |
(■ |
дГ |
(9.2) |
div |
лт |
-л- |
||
|
|
|
ап |
|
§ 2. ДИССИПАТИВНАЯ ФУНКЦИЯ. УРАВНЕНИЕ ТЕПЛОПРОВОДНОСТИ ДЛЯ ВЯЗКОЙ ж и д к о с т и
1. Диссипация в жидкости. Прежде всего, опираясь на вывед ные в предыдущей главе формулы (8.3), получим явный вид функ ции Е. Как установлено (см. 7.24), она представляет собой сумму скалярных произведений вида
|
—> . |
Е = Рх |
д о ) |
dz /’ |
Раскрывая его, получим:
/Г: |
dvx |
Р |
dv„ |
|
dvz |
|
dvu |
dvx |
• дх |
—- -U Р |
------ р Р |
Ж + |
~djf) + |
||||
|
vv |
ду |
' « |
п |
ху |
|||
|
+ Рх, |
dvz |
|
|
|
/ dvz |
dvy |
|
|
д Г |
^ |
~dz I ~ |
Fy* \. 'ду |
dz |
|
105
Ё промсжуточнь1х выкладках здесь использовалось свойство взаимности касательных напряжений, выражаемое полуденными
выше равенствами: Р.\-ц = Рцх. Рхг— Ргх, Рцг=Р |
|
|
|
|||||
Подставляя в последнее |
выражение формулы (8.3), после про |
|||||||
стейших преобразований получим |
|
|
|
|
|
|
||
И |
о |
|
dvx v |
, |
1dv,A- |
j |
di': |
|
it (div v) 2+ - Р |
дх ] |
|
\ ду / |
1‘ \ |
дг |
|
||
|
|
|
|
|
||||
|
d o . v \ - |
idi': |
ор.-у |
. |
1 Ch’z |
i да,, |
г |
(9.d) |
|
7й/ I |
’ ' l v |
с)г ) |
|
1 ду ! д.г ) . |
|
||
Из (9.3) |
мы убеждаемся |
в том, |
что механическая энергия |
пере |
ходит в тепло за счет следующих факторов. Во-первых, иод дей
ствием сил давления (член ■fullv v) происходит изменение объе
ма, что ведет при его уменьшении (divc<;0) к возрастанию тем
пературы. а при возрастании (divy>0) к ее понижению. Во-вто рых, из-за наличия сил вязкости, действие которых проявляется
двояким образом, при изменении объема при растяжении |
(сжатии) |
жидких линий внутренняя энергия уменьшается |
(член — |
о |
|
— Li(divu)2); при изменении длины жидких линий п перекосе силы
О
грения совершают работу, которая полностью переходит в тепло (два последних члена).
Величина
носит название диссипативной функции. Еще раз подчеркнем, что она представляет собой отнесенную к единице объема и единице времени часть кинетической энергии потока, затраченную на пре одоление сил вязкости и переходящую в тепло *.
Очевидно, что в случае несжимаемой жидкости выражения (9.3)
и (9.4) упрощаются за счет того, |
что уходят члены с div и. Если |
жидкость идеальна (р = 0), то |
D — 0, а Е — —р div и. Наконец, |
для идеальной несжимаемой жидкости справедливо D= E —0.
* Здесь вновь имеет место эффект необратимости процессов за счет действия сил вязкости, ибо кинетическая энергия переходит п тепло, а обратное не имеет места.
106
2. Запись уравнений теплопроводности применительно к ра личиым моделям жидких сред. Из всего сказанного выше выте кает, что уравнение притока тепла для вязкой сжимаемой жидкости должно иметь вид (см. (7.29), (9.1), (9.2), (9.3), (9.4)
dT |
■ |
, |
ОТ \ |
|
?CvW |
+р <Jiv Il=,liv |
( /"г |
~дп ) ^ Г) ■ |
19-5> |
Для несжимаемой жидкости
где D\ = 2 р (ei2—)—г--22—{- кз2)-|-(.1 (912—(—022—1~0в2) -
Мы видим, что в этом случае работа давления не оказывает влияния на тепловой баланс, ибо изменения объема в данном слу чае ие происходит.
Для идеальной жидкости будем иметь вместо (9.5) и (9.6) соот
ветственно: |
|
|
|
|
|
dT |
|
|
I дТ \ |
(9.7) |
|
['cv - J f |
+Р div l' =:div I M |
F |
|||
|
dT |
/ |
dT \ |
|
(9.8) |
[jC„ —77- — div |
( Хт |
~дп)ш |
|
||
‘ - |
dt |
|
|
||
неподвижной среды вместо уравнения |
(9.5) и (9.6) |
будем |
|||
иметь |
|
|
|
|
|
г‘С., |
ОТ |
/ |
дТ \ |
|
(9.9) |
- 57- — div |
|
||||
|
at. |
(/-т |
ЮГ Г |
|
|
Вэтом случае уравнение теплопроводности жидкости совпадает
суравнением теплопроводности твердого тела. Если же при этом
и= const, то получаем классическое уравнение , теплопровод
ности — уравнение Фурье
|
ОТ |
(9.10) |
|
— = аЧ Т , |
|
где а2= |
Лт |
|
-------коэффициент температуропроводности; |
|
§3. ПЕРВОЕ НАЧАЛО ТЕРМОДИНАМИКИ
I.Различные формы представления. По сути дела уравнение (9.5) представляет собой форму записи 1-го начала термодина мики. Обозначив
|
|
D |
|
|
можем записать |
|
|
|
|
Lt l -и JL |
dQ |
|
||
cliv v — dt |
19.11) |
|||
dt ' |
о |
|||
Поскольку из уравнения |
неразрывности следует, |
что |
1
P
dp dt ’
do |
d , 1 \ |
dT |
19.1*2) |
|
где а — удельный объем.
Подставив этот результат в последнее уравнение п помножив обе его части на dt, получим первое начало в классической форме записи
cvdT-\-p da=dQ. |
(9.13) |
Таким образом, изменение внутренней энергии слагается из ра боты внешних сил и притока извне.
Как известно, уравнение состояния совершенных газов имеет вид (формула Клапейрона)
p = p R T , |
(9.14) |
где R -—газовая постоянная. |
|
Уравнение состояния атмосферы имеет такой же вид. |
можно |
Воспользовавшись (9.14), уравнение первого начала |
модифицировать. Для этого сначала, в соответствии со вторым равенством соотношения (9.12), перепишем его в виде
р d[
cvd T --------------= d Q .
О о
Подставляя р из (9.14) иделя обе части получаемого выраже ния па Г, будем иметь
dT |
d[j |
dQ |
(9.15) |
Lv |
^ rj |
~ ~ • |
![](/html/65386/283/html_B1T5k7iRfa.sGto/htmlconvd-jpSMzr110x1.jpg)
Прологарифмировав и продифференцировав (9.1-1), можем по лучить
do |
dp |
|
dT |
|
|
~ Т ^ ~ Т |
|
~ Т " |
|
|
|
Используя последнее равенство, |
запишем |
(9.15) |
в форме |
||
d l |
_ П |
dp |
dQ |
|
(9.16) |
_ |
_С _ __ _ |
|
|||
Т |
|
р |
~ Т ■ |
|
|
Величина cl,= cr-\-R называется изобарной теплоемкостью, т. е. |
|||||
теплоемкостью при постоянном давлении. с |
> с г, |
ибо сс количе |
ственно равно затратам тепла на изменение температуры одного грамма массы жидкости на один градус при постоянном объеме, а с,, при постоянном давлении. Но в последнем случае для сохра нения давления внутри объема жидкость должна расширяться, совершая при этом работу на преодоление сил внешнего давления. Поэтому и необходимо для того, чтобы изменить ее температуру на один градус затрачивать соответственно большее количество энергии.
Подчеркнем еще раз, что все выведенные выше формулы суть различные случаи и формы первого начала термодинамики.
2. Политропические процессы. Адиабатический случай. Если предположить, что приток тепла к частице пропорционален темпе
ратуре, т. |
е. dQ = cc/T, |
то такие процессы |
называются |
политропи- |
|||||
ческими, |
а величина |
с — политропической теплоемкостью. |
В этом |
||||||
случае вместо |
(9.15) |
и |
(9.16) |
будем соответственно иметь: |
|
||||
|
|
|
|
|
d T |
rfp |
|
|
|
|
|
|
|
( с „ - Е ) — = R — • |
|
|
|||
|
|
|
|
(.*■/> О |
dT |
|
|
|
|
|
|
|
|
т* |
|
|
|
|
|
Интегрируя |
эти уравнения |
в пределах |
от То до Т, |
от |
ро до р |
||||
и от ро до р, получим: |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
(9.17) |
|
|
|
|
|
|
|
|
|
(9.18) |
Важным частным случаем политропического является адиаба тический процесс, протекающий в частице без теплообмена со
109