
книги из ГПНТБ / Картвелишвили Н.А. Потоки в недеформируемых руслах
.pdfНетрудно проверить, что
_ g r a d v 2 = ( v , V ) v + [ v , r o t v ] .
Используя эту формулу и полагая, что внешние силы имеют потенциал, т. е. что F = —grad У, а жидкость имеет постоянную плотность, приведем уравнение (11.1) к виду
- | f + g r a d r = [ v , r o t v ] . |
(11.2) |
Здесь
Для установившегося движения уравнение (11.2) имеет вид
g r a d r = j v , rotv] . |
(11.3) |
Так как векторное произведение двух векторов есть вектор, нор
мальный к обоим векторам, то из |
(11.3) следует, |
что вектор |
grad Г перпендикулярен к плоскости, |
в которой лежат |
векторы v |
и rotv. Это значит, что линия тока вектора скорости, совпадающая при установившемся движении с его траекторией, и вихревая ли
ния, т. е. траектория rotv, лежат |
на поверхности |
Г = const, т. е. Г== |
||||||
= const вдоль линии тока; этот |
вывод |
был |
получен |
в § |
8 иным |
|||
путем. Но с переходом от одной |
траектории |
к |
другой |
Г, |
конечно, |
|||
в общем случае изменяется. |
Чтобы |
Г |
было |
постоянно |
во всей |
|||
массе жидкости, как видно из |
|
(11.3), |
должно |
быть |
[v, rotv] = 0, |
|||
т. е. rotv = 0, что соответствует |
безвихревому |
движению, либо век |
торы v и rot v должны иметь одинаковое направление. Этот послед ний случай — единственное вихревое движение с одинаковым зна чением Г для всей массы жидкости. Такое движение называется однородно-винтовым. Это движение в прямоугольном канале впер вые было исследовано И. С. Громеко [20], значительное внимание уделено ему А. Я. Миловичем [41]. М. В. Потапов [47] обратил внимание на то, что однородно-винтовое движение характеризуется резко неравномерным распределением скоростей по сечению по тока, очень далеким от реально наблюдаемых распределений, ко торые в призматических руслах более или менее близки к равно мерным (поступательная скорость одинакова во всех точках попе речного сечения). Отступление от равномерного распределения есть следствие главным образом гидравлических сопротивлений, кото рые приводят к тому, что максимальная в сечении скорость превы шает среднюю только на 20—30%, а не в 2,5 раза, как в одно родно-винтовом движении. В связи с этим М. В. Потапов рассмат ривает иное, так называемое продольно-винтовое движение, с вих ревыми линиями, параллельными дну и стенкам канала, которое значительно ближе к движениям реальных жидкостей, чем однород но-винтовое.
Уравнение (11.3) в проекциях на координатные оси записыва ется так:
|
— 5Р=2( - іто - С<а), |
(11.4) |
|
|
(11.5) |
|
|
(11.6) |
Здесь s, |
т], £— проекции на оси х, у, z компонентов |
угловой |
скорости — |
rot v: |
|
Присоединим |
к написанным |
уравнениям |
р и с |
^ |
Поперечная |
|||
также уравнение неразрывности |
|
циркуляция |
в |
прямо- |
||||
ди |
, до |
, dw |
п |
/ 1 1 0 Ч |
Уольном канале |
(один |
||
дх |
ду |
dz |
=0. |
(11.8) |
|
винт, |
т=п=1). |
Рассмотрим прямоугольный канал с горизонтальным дном; на
правим оси координат как показано на |
рис. 12. Так как течение |
||||
не изменяется вдоль оси х, |
то |
|
|
|
|
дТ |
ди |
ди |
dw |
^ |
п л п . |
дх |
дх |
дх |
дх |
=0, |
(,11. У) |
|
|||||
а так как вихревые линии |
параллельны |
оси х, то т] = £=0. |
Урав |
нение (П.4) при этом удовлетворяется тождественно, уравнения (11.5) и (11.6) приобретают вид:
дТ |
дТ |
(11.10) |
|
dv |
dz |
||
|
|||
а из уравнений (11.7) и (11.9) следует |
|
||
— = 0 |
да |
=0. |
|
|
|||
dz |
и' |
|
|
Таким образом, и = const, т. е. продольная составляющая ско |
рости |
одинакова во всех точках потока. Уравнение же неразрыв |
|||
ности |
(11.9) приобретает вид |
|
dz -0. |
|
|
ду |
' |
(11.11) |
|
|
d v |
|
d w |
|
Продифференцируем |
первое уравнение (11.10) по z, |
а второе |
по у и вычтем второе |
уравнение из первого. Учитывая |
при этом |
уравнение |
(11.11), получаем |
|
|
|
|
или |
|
|
|
|
|
где |
|
(W, gradS)=0, |
|
|
(11.13) |
|
W^Vrf+w1 |
|
|
(11.14) |
|
|
|
|
|
||
— скорость поперечной |
циркуляции. |
|
|
|
|
Первое |
уравнение |
(11.7) и уравнения (11.11) |
и (11.12) |
обра |
|
зуют систему трех уравнений с тремя неизвестными v, |
w, |
|. Эта |
|||
система нелинейна, так как нелинейно уравнение |
(11.12) |
и в квад |
ратурах она в общем виде неразрешима. По М. В. Потапову, мы
ограничимся |
частным случаем, в |
котором |
квадратура |
возможна. |
|||||||||
К этому случаю мы придем, |
написав уравнение (11.12) |
в виде |
|||||||||||
|
|
|
V |
|
|
|
w |
|
|
|
|
|
|
|
|
|
ді |
|
~~ |
~~дТ~ • |
|
|
|
|
|
||
|
|
|
дг |
|
|
|
ду |
|
|
|
|
|
|
Обозначая каждое из этих отношений через К, получаем |
|||||||||||||
|
|
|
|
|
w = |
- X - g - . |
|
|
(11.15) |
||||
Частный |
случай, |
который |
мы |
имеем |
в |
виду, |
есть |
случай ?» = |
|||||
= const>0. |
Полагая |
А, = 2/А2, |
уравнения |
(11.15) |
запишем так: |
||||||||
|
|
дї |
|
k2 |
|
di |
|
k2 |
|
|
|
/111 r\ |
|
|
~W=~ |
— W> |
-dr=Tv- |
|
|
|
|
( l u 6 > |
|||||
Дифференцируя первое уравнение |
(11.16) по у, |
а второе по z |
|||||||||||
и складывая, получаем, |
учитывая |
первое |
уравнение |
(11.7), |
|||||||||
|
_дЧ_ ,_дЧ_ |
|
k2 (I dv |
|
dw \ |
, 9 t |
|
|
|||||
|
ду2 |
• |
дг2 |
|
|
|
|
|
|
|
|
|
|
или |
|
|
ду* ~' |
|
' - ^ = 0 . |
|
|
|
|
(11.17) |
|||
|
|
|
дг2 |
|
|
|
|
||||||
|
|
|
д Ч |
|
д Ч |
|
|
|
|
|
|
|
|
Таким образом, задача свелась к системе трех линейных урав |
|||||||||||||
нений (11.16) и (11.17). |
|
|
|
|
|
|
|
|
|
|
|
||
Положим |
% = YZ, |
где Y зависит |
только |
от у, |
a |
Z |
только от z. |
||||||
Подставляя это значение | в уравнение |
(11.17), получаем |
||||||||||||
|
|
|
- V |
+ |
4 - + # = 0 . |
|
|
|
(11.18) |
||||
Отсюда следует, что |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
^у-^-о2, |
|
|
^ |
= |
- р |
2 |
, |
|
|
(11.19) |
где а |
н (З — вещественные |
постоянные, |
удовлетворяющие |
условию |
||||||||
|
|
|
|
а2 + |
р2=АЗ_ |
|
|
|
|
(11.20) |
||
Из уравнений (11.19) |
|
|
|
|
|
|
|
|
|
|
||
|
|
Y = |
С[ cos ay -(- С2 |
sin ay, |
|
|
||||||
|
|
Z = D , |
cos p z + D 2 |
sin pz, |
|
|
||||||
где Сі, Co, Di, A>-—произвольные |
постоянные. Найденное |
решение |
||||||||||
должно |
удовлетворять очевидным |
краевым |
условиям: |
и = 0 при |
||||||||
у = 0 |
и |
у = Ь; до = 0 |
при |
z = 0, |
z = /i. |
Подставляя |
\=YZ |
|||||
в (11.16), находим из этих условий |
|
|
|
|
|
|
||||||
|
|
С , = Д = 0 , |
sina6=0, |
|
sinp/z=0, |
|
|
|||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
* = • — |
, |
Р = |
— . |
|
|
|
||||
где /п и п — произвольные |
|
целые |
числа, |
положительные |
или отри |
|||||||
цательные. Следовательно, |
|
|
|
|
|
|
|
|
|
|||
|
|
i j = = C s i n ^ ^ s i n - ^ , |
|
C = O D 2 , |
|
(11.21) |
||||||
|
|
2Сля . |
|
/пяу |
|
ляг |
|
,., w n m |
||||
|
|
v = = ~ M ~ s m |
|
Ъ C 0 S _ A — ' |
|
( 1 L 2 2 > |
||||||
|
|
™ = - П Р Ї - С 0 8 |
- Т - 8 Ш |
Т - |
|
О 1 - 2 3 ) |
||||||
|
|
£2 |
|
- 2 ( 4 + 4 ) - |
|
|
(П.24) |
|||||
Если вместо выражений |
(П.19) |
взять |
|
|
|
|
||||||
|
|
Г II |
|
|
"711 |
|
|
|
|
|
||
|
|
Г |
|
' |
|
Z |
- = р |
или |
|
|
||
|
|
Y" |
|
., |
|
Z" |
|
|
-Р2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(условие (11.20) в обоих случаях может быть удовлетворено), или если рассматривать случай К= const < 0, т. е. k2 < 0, то функции Y и Z (или по крайней мере одна из них) будут выражаться не три гонометрическими, а гиперболическими функциями, и тогда крае вые условия будут удовлетворяться только при С = 0, т. е. при отсутствии вихрей. Поэтому указанные комбинации интереса не представляют.
Кинетическая энергия поперечной циркуляции на единицу длины потока составляет
А ь |
.,2 |
Р ^ rfy dz =—cj— bit,
О О
где v0 — среднеквадратическая |
поперечная |
скорость. Подставляя |
|
в последнее выражение формулы |
(11.22) и (11.23) и выполняя ин |
||
тегрирование, находим |
|
C2bh |
|
а _ |
|
|
|
— |
Р |
2/г2 |
|
и, следовательно, |
|
|
|
где е = Э/рЫг — удельная энергия |
поперечной |
циркуляции, т. е. ки |
нетическая энергия поперечной циркуляции, отнесенная к единице массы жидкости.
Уравнения линий тока |
(2.2) |
в данном случае будут |
|
|
|
|
||||||||
dx |
|
dy |
|
|
|
|
|
|
dz |
|
|
|
- J J |
|
и |
ІСш |
тку |
nnz |
|
|
2Ст% |
/тису |
, |
пкг |
• |
\ |
• ' |
||
|
г ^ т — s i n |
т — C O S — г — |
|
|
— Г о ? — c |
o s — |
|
г — |
|
|
|
|||
|
k2h |
b |
h |
|
|
k2b |
|
b |
|
h |
|
|
|
|
|
Последнее равенство можно записать в виде |
|
|
|
|
|
||||||||
|
т |
тку . |
гшг |
, |
, |
п |
|
тъу |
nnz |
, |
„ |
|
|
|
или |
так: |
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, уравнение проекций линий тока |
на |
плоскость |
||||||||||||
yOz есть |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s |
i n ^ s i n ^ |
£ . |
= A . |
|
|
|
|
(П.27) |
||||
|
Здесь А — постоянная, |
которая |
может |
иметь |
значения |
от |
—1 |
|||||||
до 1, при этом, если m и п положительны, то |
|
|
|
|
|
|||||||||
|
А = 1 п р и у = ( 2 / + 4 ) А |
|
г = ( 2 ^ + ^ ) 4 , |
|
|
|
||||||||
|
Л = - 1 п р и у = ( 2 / + 4 ) ^ , |
г = ( 2 А + 4 - ) 4 , |
|
|
||||||||||
|
Л = - 1 П р и у - ( 2 / + 4 - ) А , |
z = ( 2 A + 4 - ) 4 , |
|
|
||||||||||
и при любых целых тип |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
А—0 |
при у—i |
— , |
z=k |
— . |
|
|
|
|
|
В этих формулах i n k — произвольные положительные целые числа. В потоке имеется \тп\ винтов; на рис. 12—14 это показано соответственно для т=1, п=\\ т=2, я = 1; т = 3, п = 2. При изме нении знака т или п (но не т и п одновременно) направления
вращения этих винтов и направления вихревых линий, показанных прямыми стрелками, изменяются на противоположные.
Чтобы получить уравнение проекций линий тока на плоскость
хОу, |
обратимся |
к первому |
равенству |
(11.26), которое представим |
||||||||||||
в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ішу |
|
|
nizz |
|
|
dy |
g |
a№h |
|
|
|
|
|
|
|
Sin—JT-COS |
Л |
= |
|
в.dx |
' |
2ъпс ' |
|
|
|
(11.28) |
||||
Возводя |
в квадрат |
(11.27) |
и |
(11.28) и |
складывая, |
получаем |
||||||||||
уравнение, не содержащее z, |
|
|
|
|
|
|
|
|
|
|
||||||
|
sin-- тку |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставляя |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
cos- |
тку |
] Л |
—A* sin с? |
|
|
|
|
|
|
|
|
|
|||
в это уравнение, |
получаем |
|
|
|
|
Рис. |
13. |
Поперечная |
цирку |
|||||||
|
dx2-- |
|
£>2 rf?2 |
|
|
|
|
ляция в |
прямоугольном |
ка |
||||||
|
1 — |
(1 — |
.42) |
Sin2 =р |
|
|
|
нале |
(два винта, |
т—2, |
п= |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
= |
!)• |
|
|
|
|
|
D |
иЫгЧ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2т&\тп\С |
- |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Отсюда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-x0=D |
|
Йф2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
VI—К2 |
Sin2 <}/ |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
где |
/ С 2 = 1 — Л 2 ; |
F — символ |
эллипти |
|
|
|
|
|
|
|
||||||
ческого интеграла первого рода. Урав |
|
|
|
|
|
|
|
|||||||||
нение (1.28) равносильно |
следующему |
Рис. |
14. |
Поперечная |
циркуля |
|||||||||||
|
телу -Ksn(K, |
|
|
|
|
|
|
ция |
в |
прямоугольном |
канале |
|||||
cos |
X ~ D |
X a |
) , |
|
(П.29) |
(шесть |
винтов, т=3, |
|
п=2). |
|||||||
где |
sn — символ |
так |
называемого |
эллиптического |
синуса |
|
Якоби. |
|||||||||
Это |
уравнение определяет |
сложную периодическую |
по |
х |
кривую. |
В заключение необходимо подчеркнуть, что мы рассмотрели только возможные течения с поперечной циркуляцией в прямо угольном канале. На вопрос о том, какие значения т и п будут иметь место в действительности, будет ли действительное циркуля ционное течение представляться формулами (11.22) и (11.23) или двойными рядами из членов вида (11.22) и (11.23), в которых сум мирование будет производиться по т и п и т. д., изложенная тео рия ответа не дает. Выяснение всех этих обстоятельств требует анализа конкретных причин возникновения циркуляции.
Глава IV
Г И Д Р А В Л И Ч Е С К А Я И Д Е А Л И З А Ц И Я Т У Р Б У Л Е Н Т Н О Г О Т Е Ч Е Н И Я
12. Уравнения Рейнольдса
Математическая модель идеальной жидкости, рассматривав шаяся в предыдущей главе, позволяет с удовлетворительной точно стью решить довольно широкий круг задач о движении жидкости, охватывающий те случаи, в которых влияние сил сопротивлений несущественно для решения поставленного вопроса. Однако су
ществуют |
и имеют первостепенное |
значение для |
ряда |
областей |
|
науки, в первую очередь |
для гидрологии и гидротехники, |
задачи, |
|||
в которых |
пренебрегать |
эффектом |
сопротивлений |
принципиально |
невозможно, т. е. нужно рассматривать не идеальную, а реальную жидкость, описываемую уравнениями Навье—Стокса.
Экспериментально установлено существование ламинарной и турбулентной форм движения реальной жидкости. При детермини рованных краевых и начальных условиях, скорости и давления в ла минарном потоке есть вполне детерминированные функции коорди нат и времени, теоретически определяемые как решения уравнений Навье—Стокса. Ламинарное течение возможно при сравнительно
малых |
скоростях потока, |
точнее |
при |
малых |
числах |
Рейнольдса |
||
Re = VL/v, где |
V-—характерная |
(например, |
максимальная в сече |
|||||
нии) |
скорость |
потока, |
L — характерный |
линейный |
поперечный |
|||
размер потока |
(например, диаметр |
трубы), |
v — кинематический |
коэффициент вязкости, фигурирующий в уравнениях Навье— Стокса.
В гидрологической и гидротехнической практике приходится сталкиваться почти исключительно с турбулентным движением жидкости. Турбулентное движение характеризуется тем, что при детерминированных краевых и начальных условиях скорости и дав ления в потоке не детерминированные, а случайные функции коор динат и времени, детерминированными же функциями оказываются только математические ожидания, т. е. средние значения скоростей и давлений. Случайные функции, выражающие мгновенные скоро сти и давления (точнее, распределения вероятностей этих функ ций), так же как и их математические ожидания (называемые обычно локальными усредненными скоростями и давлениями), не могут быть найдены из уравнений Навье—Стокса, по крайней мере при современных знаниях о свойствах решений этих уравнений. К этому очень тонкому вопросу мы вернемся ниже, после того как установим уравнения турбулентного потока.
Используя модель идеальной жидкости, мы попросту пренебре гаем флуктуациямп скоростей и давлений вокруг их математичес ких ожиданий, равно как и вязкими напряжениями в жидкости, а вместе с тем пренебрегаем и теми дополнительными членами, ко торые появляются в результате флуктуации и влияют на сопротив ления.
Чтобы получить уравнения турбулентного потока, рассмотрим уравнение (5.14)
- I f - + ( v • grad) v = F - - f grad p + -j - div T. |
(12.1) |
Изучая несжимаемую жидкость, характеризующуюся тем, что div v = 0, можно написать
|
|
(v |
• grad)v=( v • grad) v-f-v div v. |
(12.2) |
|||||||||
Но, с другой стороны, написав проекции |
(12.2) |
на координатные |
|||||||||||
оси Ox, Оу, Oz, можно убедиться, что |
|
|
|
|
|
||||||||
|
|
(v |
• grad) v-f-v div v = d i v ( v v ) , |
|
(12.3) |
||||||||
где тензор |
(vv) |
есть диадное произведение |
вектора |
v на самого |
|||||||||
себя. Теперь уравнение |
|
(12.1) |
можно записать так: |
|
|||||||||
|
d v |
' - d i v ( v v ) = F — - g r a d / 7 + |
— div Т. |
(12.4) |
|||||||||
Рассматривая турбулентное течение, положим |
|
|
|||||||||||
|
|
|
|
|
Р=Р+Р', |
v = v |
+ v \ |
|
|
(12.5) |
|||
где р и v — математические |
ожидания |
давления |
и скорости; р' и |
||||||||||
v' — флуктуации |
этих |
функций. Математическое |
ожидание можно |
||||||||||
мыслить как результат усреднения |
р и v по реализациям процесса. |
||||||||||||
Пусть |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рі=Рі(х, |
У> z . |
0. |
|
|
|
|
|
|
|
|
|
|
|
Vi = V ; ( A ' , у , Z, |
t) |
|
|
|
|
||
— давление |
и скорость |
|
в момент t |
в точке |
(х, у, |
z) при і-том вос |
|||||||
произведении исследуемого |
движения |
(при |
і-той |
реализации). |
|||||||||
В силу случайности |
процесса |
ріфр,, Vi=^Vj, |
если |
£ = / . По опреде |
|||||||||
лению |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
п |
|
|
|
п |
|
|
|
|
Р= |
l |
i m |
|
1Г 2 |
А> |
v , = |
lim 4" |
2 |
v,., |
(12.6) |
||
|
|
|
л->• оо |
1 — \ |
|
я - » - с с |
; = = 1 |
|
причем предел понимается, конечно, не в обычном, а в
вероятностном смысле. Очевидно, что
p'=v' |
= {vv') |
= |
0, |
( v ' v ' ) ^ O , |
|||
dv |
dv |
|
• |
- j |
— |
, |
— |
~~дТ~~д~Г |
• |
grad /7 = grad/>, |
|||||
div /==div T, |
|
divv = |
divv, |
||||
dv |
dv |
dv |
|
dv |
|
dv |
dv |
dx |
dx ' |
dy |
|
dy |
' |
dz |
dz |
Подставив p и v яз |
(12.5) в (12.4) и взяв |
операцию |
матема |
|
тического ожидания от обеих частей уравнения, получим |
|
|||
d t |
, div(vv) = |
F-grad/? + ^ d i v | 7 , - P |
( v ' v ' ) ) . |
(12.7) |
Строго говоря, намследовало бы также учесть флуктуацию |
||||
плотности, |
положив р = р + р', чего мы, однако, |
делать не |
будем, |
так как в дальнейшем будет рассматриваться почти исключительно
однородная жидкость, для которой p = const. Учитывая, |
что divv = |
||
= 0, |
и используя (12.3) и (12.2), можно записать (12.7) |
в |
виде |
- |
| f + (v • g r a d ) v = F - ^ g r a d ^ + - ^ d i v { r - P ( v 7 v T ) } . |
(12.8) |
Уравнение (12.7) носит название уравнения усредненного дви жения. Компоненты тензора вязких напряжений Т определяются формулами (6.13). Вследствие линейности этих формул компоненты
усредненного тензора Т выражаются через компоненты и, v, w усредненной скорости так же, как компоненты актуального тен зора Т через компоненты и, v, w актуальной скорости.
Актуальное турбулентное движение не может быть установив шимся. Когда говорят об установившемся турбулентном движе нии, имеют в виду движение усредненное, а именно движение, ха рактеризующееся тем, что р и v не зависят от t. В этом случае флуктуации давления р' и скорости v ' есть стационарные случай ные процессы, в силу эргодичности которых определения р И V по формулам (12.6) совпадают со следующими определениями:
|
т |
|
|
т |
|
р= 11т |
[pdt, |
v = lim - |
L |
f v dt. |
(12.9) |
T - V CO |
J |
T CD |
~ I |
J |
|
B этих формулах интеграл берется по одной реализации. За метим, между прочим, что большинство гидрометрических прибо ров, в частности все приборы, построенные на принципе трубки Пито, а также вертушки, пружинные и жидкостные манометры измеряют усредненные скорости и давления именно по формулам (12.9).
Диада (v'v') в уравнении (12.8), являющаяся симметричным тензором
(v'v')= |
(12.10) |
называется тензором пульсации. Тензор —p(v'v') часто называют тензором турбулентных напряжений. Этот тензор добавляется к тензору вязких напряжений и по отношению к усредненному дви жению действительно играет роль тензора дополнительных на пряжений. Обозначим через Т% тензор
|
|
7-, = 7 - - p ( v V ) |
|
(12.11) |
||
и запишем уравнение |
(12.7) в виде |
|
|
|
|
|
|
^ ' . dlv(vv) = F — 1р |
- grad/>+' - £1 - dlv |
7V |
(12.12) |
||
|
dt |
|
|
|
|
|
Уравнение |
(12.12) |
формально совпадает с (12.4) и может |
трак |
|||
товаться не |
только |
как уравнение |
усредненного |
движения, |
если |
под р и v понимать усредненные давления р и скорость v, а под тензором Т* — тензор (12.11), но и как уравнение актуального дви
жения, если под р, v, |
понимать |
актуальные |
давления, скорость |
и тензор вязких напряжений. |
|
|
|
Векторное уравнение |
(12.12) |
эквивалентно |
трем скалярным. |
В ортогональной прямолинейной системе координат эти три урав
нения и уравнение неразрывности |
divv записываются так: |
||||||
да |
да* |
|
dav |
|
daw |
-X-± |
dp |
dt |
dx |
|
dy |
|
dz |
dx |
|
|
|
p |
|||||
|
|
|
dz- |
|
dzr |
, d-r |
(12.13) |
|
|
|
dx |
|
dy |
dz |
|
|
|
|
|
|
|||
dv |
dva |
|
dv* |
|
dvw |
|
|
dt |
dx |
|
dy |
|
dz |
|
dy |
|
|
|
dz yx |
|
d' yy |
dz J |
(12.14) |
|
|
|
dx |
|
dy |
|
|
dw |
dwu |
, |
dwu |
|
dw* |
=z |
dp |
dt |
dx |
|
dy |
|
д. |
dz |
|
|
|
|
|||||
|
J_ |
|
|
dzzy |
dz, |
(12.15) |
|
|
+ P |
I |
dx |
|
dy |
dz |
|
|
|
da |
. |
dv |
dw |
0, |
(12.16) |
|
|
dx |
1 |
dy |
1 dz |
где Xxx, txy, и т. д. — компоненты тензора Т%