
книги из ГПНТБ / Картвелишвили Н.А. Потоки в недеформируемых руслах
.pdfсплошной среды не дают. Следовательно, чтобы замкнуть получен ную систему уравнений, т. е. привести число уравнений в соответ ствие с числом неизвестных, необходимо обратиться к физическим законам или гипотезам, не сводимым к основным законам меха ники. С одной стороны, это есть соотношения, определяющие из менчивость плотности р (простейшее из них p = const), индивиду альные для различных групп конкретных задач. С другой стороны, это связи между компонентами тензора Т и другими элементами потока. Чтобы установить такие связи, попробуем выяснить, чем определяется деформация элементарного жидкого объема при дви жении жидкости: очевидно, что напряжения должны быть связаны именно с деформациями, а не с движением жидкого объема как целого (например, с движением сосуда, в котором заключена жид кость, покоящаяся относительно стенок сосуда).
Рассмотрим тензор производной вектора скорости v по радиусувектору г. По формуле (3.18) имеем
|
|
|
|
( |
да |
ди |
даrt" ^ |
|
|
|
|
|
|
|
|
|
|
дх |
ду |
dz |
|
|
|
|
|
|
|
|
|
|
dv |
dv |
dv |
|
|
|
(6.1) |
|
|
|
|
сіт |
|
дх |
~w |
dz |
|
|
|
||
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
dw |
dw |
dw |
|
|
|
|
|
|
|
|
|
~дх~ ~dy~ |
~~dT |
|
|
|
|
|||
Представим этот тензор |
в |
виде |
суммы симметричного |
тензора |
||||||||
5 и антисимметричного тензора А: |
|
|
|
|
|
|
||||||
|
|
|
да |
|
|
ди . |
dv |
du |
|
dw \ |
1 |
|
|
|
~~дГ |
|
|
ду |
1 |
дх |
|
|
dx |
|
|
5 |
= |
dv |
ди |
|
|
dv |
' |
dv |
, |
dw |
(6.2) |
|
dx |
dy |
|
|
ду |
* |
|
' |
dy |
||||
|
|
|
|
д г |
|
|||||||
|
|
dw |
да |
|
|
dw |
, |
dv |
dw |
|
|
|
|
|
дх |
dz |
|
|
dy |
1 |
dz |
|
dz |
|
|
|
|
о |
|
|
da |
|
dv |
|
du |
|
dw |
|
|
|
|
|
dy |
|
dx |
|
dz |
|
dx |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
du |
dv |
|
|
|
|
|
dw |
|
dv |
(6.3) |
|
|
dy |
dx |
|
|
|
|
|
dy |
|
dz |
|
|
|
du |
dw |
|
dw |
|
dv |
|
0 |
|
|
|
|
|
~dz~ |
dx |
|
dy |
|
~~d~F |
|
|
) |
||
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Из |
выражения |
(6.3) и |
формул |
для скалярного |
произведения |
тензора на вектор и для векторного произведения векторов выте кает, что
(A, rfr)=K dv], (6.4)
где
1 |
(6.5) |
|
<!)=—rot v |
||
|
Как известно из кинематики твердого тела, вектор [и, йт] есть та часть перемещения точек тела, которая вызвана его вращением вокруг мгновенной оси с мгновенной угловой скоростью (о. Следо вательно, антисимметричная часть тензора dv/dr определяет вра
щение |
жидкой частицы как |
твердого тела, с угловой |
скоростью |
- ^ - rotv |
и никакого влияния |
на тензор напряжений Т |
оказывать |
не может. Поступательная составляющая перемещения определя ется скоростью v. Симметричная же часть 5 тензора dv/dr связана с деформацией частицы и называется тензором скоростей деформа ции.
Уточняя сделанное предположение, можно утверждать, что
тензор Т |
напряжений, вызываемых отклонением свойств жидкой |
среды от |
свойств идеальной жидкости, определяется тензором 5 |
скоростей |
деформации. Простейшая связь такого рода состоит |
в том, что компоненты тензора Т есть линейные однородные функ ции компонент тензора 5 с коэффициентами, не зависящими от вы бора осей координат. Жидкости, удовлетворяющие этому условию,
называются ньютоновскими жидкостями. |
|
Пусть xi, уи zi — главные оси, а ех, еу, ег—собственные |
значе |
ния тензора скоростей деформации. Тогда, согласно определению ньютоновской жидкости,
ІХІХ^О-ХЄХЛ-ЬФУЛ- <hez, (6.6)
где аи аг, аз — постоянные. Переименуем |
ось Х\ в ось уІ, ОСЬ у І — |
||
в ось ZI и ось zi — в ось ХА. Тогда получим из предыдущего |
выра |
||
жения |
|
|
|
и аналогично |
|
|
|
^ l Z , = a i e z - f a 2 ^ + a 3 |
e r |
(6.8) |
|
Произведем теперь преобразование |
координат |
|
|
•**=*i. У * = 2 ь |
г* = |
—Уі- |
|
Новые оси А'*, ул., z.# опять суть главные оси тензора скоростей деформации. Если Ui, vi, Доі — компоненты скорости в прежней си стеме, то компоненты в новой системе будут
|
«* = |
«ь |
v*=Wu |
™ * = |
—Щ. |
Далее |
|
|
|
|
|
* |
ди* |
dti\ |
* |
dv* |
dw\ |
Є х ~ |
дх* ~~ |
дхг |
~ Є х ' Є>~~ду~—~дїІ |
е » |
|
* |
dw |
di/j |
|
|
|
Так как коэффициенты а.\, аг, аз не зависят от выбора системы координат, то
* . * . *
и тогда, выражая е*, е*, е* через ех, еу, ez, получаем
Сравнивая последнее выражение с формулой |
(6.6), находим, |
||
что аг = азПолагая. |
|
|
|
fl[=X+2(A, |
а2 =а3 =--Х, |
|
|
перепишем (6.6) — (6.8) так: |
|
|
|
т*,*,=^ ( е , + е у + е г ) + |
2|хвЛ, |
|
|
*у.у, = Че* + |
Єу + Є,) + |
2цЄу, |
|
**л=Ме,+<?у+е*) + 2|«?г . |
(6.9) |
Рассмотрим теперь какой-либо недиагональный компонент тен зора напряжений, например xXlVl. Согласно принятому предполо жению должно быть
*хіУі=а4ех+а-0еу-\-абег> |
(6.10) |
где 04, й5, ав — постоянные. Произведем преобразование координат
* * = * і . У* = - У ь 2» = —г,.
Новые оси опять будут главными осями тензора скоростей де формации, при этом очевидно, что тх, yt = —T .t,y i > так как направ ление оси х'і при преобразовании сохраняется, а направление оси t/i меняется на противоположное. Далее
и, = |
и„ г>* = |
- ї > „ ^ |
= |
* |
ди* |
да, |
- ® „ ^ = |
- ^ = |
- ^ L - = e „ |
||||
* |
dv# |
dvi |
* |
dw+ |
dw\ |
|
Но в силу независимости 04, as, ав от выбора системы координат должно быть
Яг.у, =а4е*х+аьеу+а6е*.
На основании предыдущих равенств получается
— У , = а 4 е у + а 6 е 2 .
Сопоставляя этот результат с (6.10), приходим к выводу, что
т=0. А так как главные оси тензора скоростей деформации S
названы осями Xi, yi, zi произвольно, то из этого следует, что все недиагональные компоненты тензора напряжений Т в системе ко ординат OJKir/iZi равны нулю, т. е. что главные оси тензора Т со впадают с главными осями тензора 5. Так как давление р опреде ляется по формуле (5.13), первый инвариант тензора напряжений равен нулю
ххх + lyy + ^zz — хх>х, + 1 у,у, + 1 г Л = 0. |
(6.П) |
Замечая далее, что
ex + ey+ez = divv,
находим из (6.10) и (6.11)
(3X+2H.)divv=0, откуда К=-—2(.і/З. Теперь можно написать
|
1 |
0 |
0 |
|
Т= |
=-p.divv О |
1 |
0 -2(х5 |
(6.12) |
О 0 1
и, следовательно, в любой системе координат будет:
-=.г.с= |
1 - [лdiw + 2[х |
|
, |
|
|
|_p.divv + |
dv |
||||
дх |
УУ |
|
2|x-g^- , |
||||||||
|
|
|
|
' |
|
3 |
|
|
|
||
|
|
|
|
2 |
|
,. |
, 0 |
dw |
|
|
(6.13) |
|
|
* « = - - 3 - I i |
d |
l v v + 2 |
| , |
, " 5 r ' |
|
||||
|
|
|
|
||||||||
|
|
/ да |
, dv \ |
|
|
|
|
I |
dv , |
dw \ |
|
- . v y — v — ^\~dT+~dT')' |
|
|
|
""^—"^у—^\дТ^"1у~) |
|||||||
|
|
|
|
|
|
du |
, |
dw |
|
|
|
|
|
-zx— -xz-* |
V- \-JT |
|
+~dl |
|
|
|
|||
Постоянная |
и. называется коэффициентом |
вязкости. |
|||||||||
Если теперь |
написать |
уравнение |
(5.14) |
в |
проекциях на коор |
динатные оси и использовать полученные выражения для компо
нентов тензора |
Т, то мы получим следующие три уравнения, назы |
|||||||||
ваемые уравнениями |
Навье—Стокса: |
|
|
|
|
|
||||
du |
du |
du |
-ISO |
du |
|
dp |
+'(4- |
дйіч |
v |
|
dt |
дх |
dy |
~d~T |
p |
dx |
dx |
|
|
||
dv |
dv |
dv |
•w |
dv |
|
dp |
|
ddiv |
V j |
•Дг»), |
dt |
dx -v |
dy |
dz |
|
dy |
|
dy |
1 |
||
dw |
dw |
dw |
•w |
dw |
p |
dp |
|
d dlv v |
-f-Awj . |
|
dt |
dx |
dy |
~dz~ |
dz |
|
dz |
|
(6.14) |
||
|
|
|
|
|
|
|
|
|
|
|
|
Здесь X , Y, Z — проекции |
на оси координат вектора |
F ускоре |
|||||||
ния |
объемных |
сил; Д — оператор Лапласа; v = |Vp — кинематичес |
кий коэффициент вязкости, обязанный своим названием тому, что его размерность [Ь2Т_1] не содержит силы.
Для жидкости постоянной плотности три уравнения Навье— Стокса (6.14) и уравнение неразрывности (5.18) образуют замкну тую систему.
Следует заметить, что предположения Навье—Стокса, приво дящие к зависимости (6.12), могут быть уточнены с термодинами ческих позиций, и тогда мы получаем более точные выражения, учитывающие так называемую вторую вязкость (бароклинной) жидкости. Однако для задач, рассматриваемых в данной книге, этот круг вопросов не играет роли, и мы не будем его затрагивать.
|
Большинство |
жидкостей |
(например, |
вода, |
воздух, |
технические |
|||||||||||||||
масла) |
с очень высокой точностью |
могут считаться |
ньютоновскими |
||||||||||||||||||
жидкостями. Это значит, что коэффициент вязкости этих |
|
жидкостей |
|||||||||||||||||||
2Л |
|
|
|
|
|
|
|
|
|
|
практически не зависит от |
||||||||||
|
|
|
|
|
|
|
|
|
/ |
характеристик |
|
|
течения |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
/ |
(например, |
для |
воды |
он |
|||||||
|
|
|
|
|
|
|
|
|
/ |
изменяется только на 1 % |
|||||||||||
2,2 |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
при |
изменении |
давления |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
от 1 до 100 атм). Но он |
||||||||||
2,0 |
|
|
|
70sji . |
|
|
|
|
может |
довольно ощутимо |
|||||||||||
|
|
|
|
|
|
«г |
|
|
|
зависеть от |
других |
пара |
|||||||||
|
|
|
|
/ |
|
|
|
|
метров, например, от тем |
||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
- |
|
|
|
|
|
|
пературы |
жидкости, |
|
как |
|||||||||
1,8 |
1 |
|
|
/ |
|
|
|
|
|
|
видно |
на |
рис. 7, где |
даны |
|||||||
ft—^ |
• |
У * |
|
|
|
|
|
|
температурные |
|
кривые р, |
||||||||||
|
|
|
• |
105-о |
|
|
|
|
|
|
1-і, v |
для |
воды |
|
и |
воздуха |
|||||
1,6 |
|
/ |
|
|
|
|
|
|
|
при |
атмосферном |
давле |
|||||||||
\ |
|
|
|
|
|
|
|
|
|||||||||||||
|
/ т |
|
|
|
|
|
|
|
|
нии. |
|
|
|
|
|
|
|
|
|
||
|
\ |
/ |
|
|
|
|
|
|
|
|
Поскольку |
нас |
интере |
||||||||
1,4 |
/Л |
|
|
|
|
|
|
|
|
|
суют |
почти |
исключитель |
||||||||
|
|
|
|
|
|
|
|
|
но |
водные |
потоки, |
то |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
рассматривать |
|
неньюто |
||||||||
1,2 |
|
|
|
|
|
|
|
|
|
|
новские жидкости |
(назы |
|||||||||
|
|
|
Юр |
|
|
|
|
|
|
ваемые аномальными) |
мы |
||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
Ю'2р |
|
|
|
|
|
|
не будем. |
|
|
|
|
|
|
|
||||
1,0 |
|
|
|
|
|
|
|
|
|
Вязкость |
|
|
идеальных |
||||||||
|
|
|
|
|
|
|
|
|
|
жидкостей, |
не |
|
способных |
||||||||
|
|
|
|
|
|
|
|
|
|
|
воспринимать |
|
касатель |
||||||||
|
|
|
|
|
|
|
|
|
|
|
ные |
напряжения, |
равна |
||||||||
0,8 |
|
|
|
|
|
|
|
|
|
|
нулю |
(v = 0), и в этом слу |
|||||||||
|
|
|
|
106\1 |
и |
05ч |
|
|
|
чае |
уравнения |
|
Навье— |
||||||||
|
|
|
|
|
|
|
|
|
|
|
Стокса |
обращаются в ура |
|||||||||
0,6 |
|
|
|
|
|
|
|
|
|
|
внения |
Эйлера. |
Однако |
||||||||
|
|
|
|
|
|
|
|
|
|
|
было |
|
бы |
ошибкой |
|
рас |
|||||
|
|
|
|
|
|
|
|
|
|
|
сматривать уравнения Эй |
||||||||||
OA |
|
|
|
|
|
|
|
|
|
|
лера |
просто |
как |
частный |
|||||||
|
|
1 |
|
|
|
|
|
|
|
|
случай уравнений Навье— |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
Стокса, отвечающий |
част |
|||||||||
0,2 |
|
|
|
|
|
|
|
|
|
700°с |
ному |
значению коэффици- |
|||||||||
|
|
20 |
40 |
60 |
|
80 |
|
ента |
кинематической |
|
вяз |
||||||||||
Рис. |
7. Зависимость |
плотности |
р кг • с/м4 'и ко |
кости |
v = 0, |
ибо |
между |
||||||||||||||
теми |
и другими |
уравне |
|||||||||||||||||||
эффициентов динамической |
Л| |
кг - с/м 2 |
и |
кине |
|||||||||||||||||
матической v |
м2 /с |
вязкости |
воды |
(/) |
и |
воз |
ниями |
есть |
существенная |
||||||||||||
|
|
духа |
(2) |
от температуры. |
|
|
принципиальная |
разница. |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Вязкость, |
как |
|
бы |
мала |
||||||
она |
ни |
была, автоматически |
приводит |
к тому, |
что |
в |
потоке |
од |
нородной жидкости поле скоростей непрерывно: разрывы в этом поле имели бы следствием, согласно формулам (6.13), физически не возможные бесконечно большие напряжения.
Таким образом, в уравнениях Навье—Стокса уже имеется пред посылка непрерывности их решений по координатам. В уравнениях же Эйлера такой предпосылки не содержится, и эти уравнения, вообще говоря, допускают разрывные решения. Естественное тре бование непрерывности решений этих уравнений привносится извне и при этом оно далеко не всегда может быть выполнено. Например, условие ограниченности касательных напряжений требует, чтобы на неподвижных твердых стенках скорость жидкости равнялась нулю (прилипаймо потока), что вполне согласуется с уравнениями Навье—Стокса. Для идеальной же жидкости такое требование обычно приводит к совершенно абсурдным выводам, вроде того, что при его выполнении жидкость может только покоиться. В ре зультате сплошь и рядом оказывается, что при v->-0 решения уравнений Навье—Стокса не переходят в соответствующие решения уравнений Эйлера.
Для несжимаемой |
жидкости |
уравнения |
Навье—-Стокса в не |
|||||||||||
сколько более подробном виде записываются так: |
|
|
||||||||||||
ди |
, |
ди |
, |
|
|
ди |
, |
|
ди |
v |
|
1 |
dp |
, |
dt |
' |
- й —з |
r - v - з |
ду |
\-HD-A—=Х |
|
|
р |
dx |
— |
||||
дх |
' |
|
|
' |
|
dz |
|
|
|
|||||
|
|
|
. |
/ |
д*и |
|
d*u . |
d*u |
\ |
|
|
|
||
|
|
+ |
V |
\ |
dx* |
+ |
ду* |
" I " |
dz* |
j • |
|
|
|
|
dv . |
dv |
, |
|
|
dv |
, |
dv |
|
|
1 |
dp |
|
||
dt |
|
д х |
-г» |
|
д у |
^ |
dz |
|
|
P |
dy |
|
||
|
|
+ 4 |
|
d*v |
, |
d*v |
d*v |
|
|
|
|
|||
|
|
|
dx* |
~ |
dy* |
' |
dz* |
|
|
|
|
|||
dw |
|
. dw . |
|
<9те/ , |
|
дїе; |
^. |
1 |
dp |
, |
||||
й т - а , |
|
' - ^ ^ |
+ і ™" - ^^ -= "z |
— |
: p |
dz |
+ |
|||||||
|
|
+ 4 ^ + ^ + ^ ) ' |
|
|
( 6 Л 5 ) |
|||||||||
а уравнение неразрывности имеет вид |
|
|
|
|
|
|||||||||
|
|
|
|
du |
|
, |
dv |
. |
dw |
|
|
|
|
(6.16) |
|
|
|
|
дх |
|
' |
dy |
' |
dz |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Принципиальное различие между идеальной и ньютоновской жидкостями стирается только тогда, когда жидкость находится в равновесии. Для покоящейся жидкости v = 0 и уравнение (5.12) приводит к таким же уравнениям равновесия—уравнениям гид ростатики, что и уравнения (6.14)
= ^ |
< 6 Л 7 > |
Уравнение же неразрывности для покоящейся жидкости удов летворяется тождественно.
7. Термодинамическое уравнение
Изменения плотности однородной жидкости происходят от изме нений давления и температуры. Изменения этих величии связы ваются уравнением, получаемым из первого начала термодинамики и называемым термодинамическим уравнением. Первое начало мо жет быть записано так:
|
|
|
|
|
'2 , |
|
| J ( 5 + t/)rfl/ J |
- |
j j |
{3 + U)dV^ |
= -j§p(v-n)dSdt |
+ |
|
+ 1 |
і |
P( F |
• v)dVdt |
+ |
§ j Q*dVdt. |
(7.1) |
/, |
v |
|
|
|
t, v |
|
Здесь U и |
|
|
|
|
|
|
|
|
Э = - £ - ( я 2 |
+ г; 2 +да 2 ) |
(7.2) |
соответственно внутренняя (т. е. обусловленная молекулярным дви жением и деформацией сжатия или расширения) и кинетическая энергия единицы объема жидкости; Q* — полный приток тепла к единице объема жидкости в единицу времени, остальные обозна
чения прежние. Левая часть |
уравнения (7.1) есть изменение |
пол |
||||||||
ной энергии объема |
V за время от момента |
ti до момента t%, первые |
||||||||
два члена правой части — работа сил давления |
(знак минус |
поя |
||||||||
вился потому, что нормаль |
к поверхности |
5 направлена |
наружу |
|||||||
объема) и объемных |
сил за то же время, |
последний член |
правой |
|||||||
части — полный приток тепла к объему. |
|
|
|
|
||||||
Преобразуя |
в (7.1) поверхностный |
интеграл в объемный, |
пола |
|||||||
гая U — t, t2 |
= t+dt, |
группируя члены |
и сокращая |
на dt, получаем |
||||||
I |
{ 4г |
<Э + |
^ ) + d i v р у - р (F • v ) - Q*\ d V=0 |
|
|
|||||
или в силу произвольности объема V |
|
|
|
|
|
|||||
|
|
йЭ |
, |
dU |
divpv= P ( F • v) + Q*. |
|
(7.3) |
|||
|
|
dt |
1 |
dt |
|
|||||
|
|
|
|
|
|
|
|
Внутренняя энергия U и давление р зависят от абсолютной температуры жидкости в и ее удельного объема (объема единицы массы р _ 1 = 1/р)
U=U(Q, р-1 ), p=p{Q, р-1 ). (7.4)
Дифференцируя калорическое |
уравнение |
состояния [первое |
|||
уравнение (7.4)], получаем |
|
|
|
|
|
dU^-^dQ |
+-%L- |
|
dp-*=cvde-E-f-, |
(7.5) |
|
где |
|
|
_, |
|
|
dU |
г? |
|
dU |
dU |
|
дв |
' |
v |
|
д?-1 ~~ v |
д? |
46
Эти величины |
имеют четкий |
физический |
смысл: Су — теплоем |
||
кость жидкости |
при постоянном |
объеме, а |
Е—модуль |
упругости |
|
при постоянной температуре |
(отношение — dpfp есть объемное сжа |
||||
тие жидкости). В обычном |
для |
естественных потоков |
и потоков |
в гидротехнических сооружениях диапазоне изменения р, р, Т теп лоемкость су и модуль упругости Е можно считать постоянными.
Выражения (7.2) и (7.5) и термическое уравнение состояния [второе уравнение (7.4)] связывают все величины, кроме Q*, вхо дящие в термодинамическое уравнение (7.3) с механическими ве личинами, входящими в уравнения Навье—Стокса и уравнение не разрывности. Дальнейшие преобразования дифференциала внут ренней энергии dU строятся на этих зависимостях. Но поскольку эти преобразования различны для различных конкретных задач, мы их сейчас не будем рассматривать.
Приток тепла можно представить |
в виде Q* = Q + qn, где qn — |
приток тепла извне, Q — приток тепла |
за счет внутреннего трения |
вжидкости. Внутреннее трение входит через тензор напряжений Т
вуравнение движения (5.14). Но получить из этого выражения ве личину диссипирующейся и переходящей в тепло энергии нельзя, ибо роль напряжений в этом уравнении сводится к чисто силовым воздействиям. Умножив (5.14) скалярно на pv, получим работу напряжений в виде v • div Т, но это будет только та часть их работы, которая связана с передачей движения от слоя к слою. Полная же работа напряжений по поверхности 5, ограничивающей некоторый объем V, в единицу времени равна
0(рп-рп) |
• vdS. |
(7.6) |
Под интегралом стоит разность полного напряжения р п по пло щадке dS и вектора рп, ибо давление р необратимой работы заве домо не производит. Но
(Рп-/>п) • v=(n, W), |
||
где W — вектор с компонентами |
(хх, v), (xv, v) , (xz, v) , причем |
|
т ^ = |
Іх .м + |
Кі-у + kz*z> |
т у = |
І ^ух ~Ь Іхуу "Т~ к^уг> |
|
|
|
т * = К.с + К у + к ^ |
|
|
|
|||
и, следовательно, по формуле |
(2.27) |
|
|
|
|
||||
|
§(?а-рп) |
- vdS=j>{n, |
W)dS=\u\vVidV. |
|
|||||
|
S |
|
|
|
S |
V |
|
|
|
А так как |
|
|
|
|
|
|
|
|
|
d l v W |
v) , a ( V v > , d ( * z |
, v ) ( < h x |
. by |
• dxz |
\ , |
||||
U 1 V V V |
dx |
^ |
dy |
T |
d z |
—[ д х |
-Г д у |
-Г d z |
. vy-t- |
|
+ T - |
* -dT+b |
• -dT+x> |
' - ^ - = v |
. d i v r + Q , |
|
где
„ |
ди . |
ди , |
dv . |
I ду , ди \ , |
то |
|
|
|
|
|
|
|
|
|
<j> (Pn-/>n) • ^ ^ 5 = j ( v |
• div |
T+Q)dV. |
|
||||
Величина Q есть |
разность |
между |
полной |
работой напряжений |
||||
на единицу объема в единицу времени |
и работой |
напряжений, свя |
||||||
занной |
с передачей |
движения |
от слоя к слою. Следовательно, |
|||||
именно |
Q характеризует |
диссипацию |
энергии в капельной |
жидко |
||||
сти, т. |
е. дает то тепло, |
которое сообщается |
жидкости |
за счет |
внутреннего трения. Если воспользоваться формулами предыду щего параграфа, получим
dy |
' |
dz |
* |
~[ |
N |
S |
+ |
( |
7 |
. |
8 |
) |
' |
) |
dz |
~ |
дх |
|
|
|
|
|
Приток тепла извне (если пренебречь очень малой теплопровод ностью жидкости) состоит из тепла q, приносимого движущейся жидкостью, и притока тепла qe за счет поглощения излучения. Очевидно, что
<7=-div(<:0v) |
(7.9) |
где с — теплоемкость жидкости. Следовательно, |
|
? „ = - d i v ( C e v ) + ?t . |
(7.10) |
Таким образом, все члены термодинамического уравнения (7.3), кроме qi, расшифрованы. Если поглощение или изучение лучистой теплоты отсутствует (qi = 0), то уравнения Навье—Стокса, уравне ние неразрывности и термодинамическое уравнение образуют замк нутую систему. Баланс лучистой теплоты мы не будем рассмат ривать.
8.Простейшие закономерности для потоков идеальной жидкости
Напишем уравнение (5.12) в виде
4 r = F - T g r a d / 7 |
( 8 Л ) |
и умножим его скалярно на элементарное перемещение жидкой час тицы у^.вдоль траектории вектора скорости v
. vdt=F • vdt—- grad p • v dt
или
(8.2)
Это уравнение интегрируется, если выполнены три следующих условия:
1) объемные силы имеют потенциал, т. е. существует такая функция V, что
х _ |
дУ |
у——-^- |
Z — |
d V |
|
дх |
' |
ду ' |
dz ' |
2)движение установившееся, т. е. dv/dt = dp/dt = 0 и траектории вектора скорости совпадают с линиями тока;
3)жидкость баротропна.
При выполнении первых двух условий |
v |
Xdx + Ydy+Zdz=-dV, |
^Ldx+-^dy+-&.dz=dp |
суть полные дифференциалы, и уравнение (8.2) может быть пере писано так:
а при выполнении третьего условия
|
где Г — величина, постоянная для данной линии |
тока, но меняю |
||||
|
щаяся, вообще говоря, при переходе от одной линии тока к дру |
|||||
|
гой. Интеграл |
(8.3) уравнения |
(8.1) называется |
интегралом Бер- |
||
|
нулли. Если жидкость однородна и несжимаема, а объемные силы |
|||||
|
сводятся к силе тяжести, т. е. V=gz, то уравнение |
(8.3) приводится |
||||
|
к виду |
|
|
|
|
|
|
|
* + - f - + |
- £ = const, |
|
(8.4) |
|
|
где y = pg— объемный вес жидкости. |
|
|
|||
' |
Это уравнение достаточно хорошо известно, и мы привели его |
|||||
! |
здесь только |
для сравнения. Из |
уравнения |
(8.4) выводится так |
||
: |
называемое уравнение Бернулли |
для целого |
потока, являющееся |
|||
': |
основой всей |
гидравлики установившихся течений. В § 14 мы да |
||||
|
дим вывод уравнения, обобщающего уравнение Бернулли для це |
|||||
|
лого потока на неустановившиеся |
течения, с совершенно иных ис- |
||||
|
ходных позиций. Это дает возможность оценить сильные и слабые |
|||||
|
стороны традиционных концепций гидравлики. |
|
|
Перейдем к безвихревому движению. Так называется движение, характеризующееся тем, что в каждой точке поля скоростей v вы полняется условие rotv = 0. Это значит, что v есть потенциальный
4 Заказ № 428 |
49 |