книги из ГПНТБ / Палагин, Э. Г. Основы гидромеханики учебное пособие для метеорологов
.pdfхарактер турбулентности определяется градиентом скорости и кри визной ее профиля, от которых зависит возможность п степень впхреобразованнй в потоке. Анализ размерностей приводит нас к простои формуле
dvx
I= X dz d1 v x '
dz-
где у. — эмпирическая константа.
Ее можно переписать несколько иначе:
d I d v . y '
dz Vdz
Если скорость потока не параллельна осп .v, то последнее вы ражение легко обобщается и на этот случаи, приобретая вид
IS.5)
d Г/rf-z',.
IF [ Г Д т
Формулы (18.4) или (18.5) в принципе позволяют замкнуть за дачу, ибо, с учетом (18.1), сможем для коэффициента турбулент ности записать выражение:
,. |
1т' |
da |
(18.6) |
|
dz |
~dz |
|
|
|
Здесь принято /„ = /.
Однако обе рассмотренные схемы являются с физической точки зрения весьма грубыми и, кроме того, справедливы только для слу чая нейтральной стратификации. Расширить сферу их применения и уточнить можно с помощью уравнения баланса турбулентной энергии, которое в данном кратком очерке мы рассматривать не будем.
§ 2. ТУРБУЛЕНТНОЕ ТЕЧЕНИЕ ПО ПОВЕРХНОСТИ ГЛАДКОЙ ПЛАСТИНЫ
Для иллюстрации применения идей теории пути смещения п выявления отличий между ламинарными и турбулентными тече ниями, обратимся к простейшему, но принципиально важному примеру. А именно, рассмотрим движение несжимаемой жидкости
22U
I
вблизи поверхности плоской гладкой стенки бесконечной протя
женности. |
будем предполагать |
отсутствие внешних сил |
При этом |
||
п считаем, что |
движение происходит |
без ускорения. Последнее |
означает, что если ось .v направить |
вдоль осредпенпоп скорости i'.v, |
|
го v,, — 0 п |
все производные по .v от осредненных величин равны |
|
пулю ^1<)»х |
<)р = 0\| . Поскольку |
речь оудет идти о движении |
в области, находящейся вблизи стенки, то в задаче справедливы приближения пограничного слоя. Это означает, что изменение на пряжений по // значительно превосходит изменение их по .v. Поэто му логично оставить только производные от напряжений по // п не учитывать производных по .v. В такой постановке уравнения лами нарного движения вырождаются в лишь из первого следует ра венство
Интегрируя это соотношение, получим o.v = C//+Cy
Таким образом, имеем линейный профиль скорости. Постоян ные интегрирования легко определить, зная, что vx i v о = 0 , н по
лагая известной величину напряжений на стенке тц = р сь!
dy |
у |
о |
Используя эти условия, будем иметь |
|
|
( 1 |
S |
. / ) |
В случае турбулентного движения вся система уравнений в на шем случае сводится к равенству
г) v'y г[,
Интегрирование этого равенства даст
--vx'v,( = С2.
При j/->0 —>0, а следовательно, v / v , / —.>0. Поэтому С2 должно быть равно то, т. е. в самом тонком, непосредственно прилегающем к стенке слое, движение должно быть ламинарным (ламинарный подслой). При удалении от нее турбулентная вязкость начинает доминировать п влиянием молекулярного трения можно прене бречь. Поэтому
—[I ул/у,/ = т0= const.
221
С учетом (18.6) последнее равенство может быть переписано в виде
или, с использованием гипотезы Прандтля относительно пути сме шения 1 — у.ц, приходим к уравнению
du |
_ J _ |
, f |
J _ |
dy |
■/. |
| |
у |
После интегрирования получим
« = - ^ - | / -о°-In г/Ч-Сз. |
(18.8) |
Сравнивая (18.7) и (18.8), мы убеждаемся в том, что в турбу лентном потоке имеет место совершенно другое распределение ско ростей, чем в ламинарном *. При этом, как и следовало ожидать, профиль скорости при наличии турбулентности внизу изменяется значительно быстрее, что связано с большей интенсивностью обме на количеством движения.
В заключение заметим, что постоянную С% следует определять из условия сопряжения турбулентного слоя с ламинарным под слоем. ибо логарифмически!) профиль не может быть проэкстраполирован вплоть до самой стенки.
* Разумеется, этот вывод справедлив не только по отношению к дампов за даче, во мосмт общий характер.
222
П Р И Л О Ж Е Н И Е A
ЭЛЕМЕНТЫ ТЕНЗОРНОГО ИСЧИСЛЕНИЯ
Объекты различной физической природы требуют при матема тическом описании различного числа компонент. Так скалярные величины типа температуры, плотности, давления могут быть пол ностью охарактеризованы только своими численными значениями. Другие являются векторами (перемещение, скорость, ускорение, сила, момент силы и т. д.). Для их определения надо указать не только численное значение величины, но и ее направление в про странстве. Однозначно заданы они могут быть совокупностью трех величин, например проекциями на оси какой-либо координатной системы. Примерами более сложных объектов являются тензор напряжений и тензор скоростей деформаций в жидкости, требую щие для своего описания девяти компонент. Свойства анизотроп ных тел определяются совокупностью 81 величины и т. д.
Удобно с целью унификации назвать скалярные величины тен зором нулевого ранга (3° — одна компонента), векторные — тензо ром первого ранга (З1— три компоненты), тензор второго ранга требует З2 величин, четвертого — З4 и т. д. Таким образом, прихо дим к понятию тензора п-го ранга, имеющего Зп компонент, и будем с этих позиций рассматривать операции над конкретными физическим и объектами.
Примером возникновения понятия тензоров может служить на пряженное состояние в жидкости. Напряжение есть сила внутрен него взаимодействия частиц жидкости, отнесенная к единице пло щади. Ее векторное описание в принципе невозможно, ибо помимо величины и направления в пространстве должна быть также известна ориентация площадки, к которой она приложена. Чтобы избавиться от последнего ограничения, напряжение в точке сле дует выразить через три величины, приложенные к площадкам строго фиксированным в пространстве. В качестве последних удобно выбрать координатные поверхности пли плоскости в случае декартовых координат. Тогда сила, приложенная к каждой из этих поверхностей, является уже вектором и может быть задана тремя компонентами, а всего их, естественно, будет девять (см. гл. VII).
Напряженно, |
как |
мы |
видели, |
записывается |
в виде матрп- |
||||||
иы (P/j): |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ри |
“С to |
Pvt |
|
(Л.1) |
||
|
|
|
|
|
Р-2, |
Р22 Р23 |
|
||||
|
|
|
|
|
Р-м Р*2 Рзз |
|
|
||||
Отметим, что для тензоров вообще удобна матричная запись. |
|||||||||||
При этом |
скалярная |
величина |
|
запишется |
просто (а). |
||||||
Ьслп «-=(«!, |
оо, аз) - |
вектор, |
то |
будем |
иметь |
.матрицу-столбец |
|||||
|
|
|
|
|
<h |
|
|
|
|
|
|
|
|
|
|
|
<К |
|
|
|
|
|
|
Для тензора |
второго ранга |
(<?;/ |
) |
(/, |
/ = 1 , 2 , |
3) |
матрица имеет |
||||
вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Щ| |
и12 |
«и |
I |
|
|
|
|
|
|
|
I |
Оч| |
a-in |
023 |
| |
• |
|
(Л.2 ) |
|
|
|
|
1 |
о31 |
a:V2 |
O33 |
J |
|
|
|
Cy.MMoii (разностью) тензоров является тензор, компоненты ко торого представляют собой сумму (разность) компонент слагае мых. Для тензора второго ранга, например, имеем
Ясно, что складывать (вычитать) можно лишь тензора одного ранга и в результате получаем тензор того же ранга. 1 1о сути дела речь идет о сложении или вычитании двух аналогичного вида матриц.
Умножение тензора на скалярную величину а сводится к умно жению на псе всех его компонент, т. е. а (а • J - - (а <7у ).
Внешним произведением тензоров называется новый тензор, ранг которого равен сумме рангов сомножителей. Его компоненты представляют собой всевозможные комбинации произведений ком понент сомножителей. Пример произведения тензоров первого
г- г — >
ранга дает диада a v (v -■ скорость)
Н |
] |
| 0 , 0 , |
V ]V '2 |
0 1c':i |
(А 3) |
|
|
1>2~0 1 |
|
V 2 V 3 |
|
г . . |
1 |
1 |
с'зс'з |
У3^3 |
|
В итоге получаем тензор второго ранга. Аналогично при умно жении тензоров второго ранга и,- , ft/,,,, получаем тензор четвертого
ранга аикт — аи ft,,,,,.
224
Существует тензор, называемый единичным, при умножении на который каждый тензор 2-го ранга переходит сам в себя. Он обозначается U, а его компоненты б,у , причем 6,у = 1 , если /—/
и бг/- = 0, если i ф /. Таким образом,
|
|
' |
( 1 |
0 |
0 ] |
|
|
|
(А.4) |
|
|
и — |
0 |
1 |
О |
I |
|
|
|
Легко удостовериться, что |
1.0 |
0 |
1 |
|
|
|
|||
|
|
|
|
|
|
|
|||
1 |
О О |
а \ 1 |
0 12 |
° 1 3 |
] |
Я | | |
Я 12 |
« 1 3 |
|
aU->au 8у = О |
1 |
О |
Й 21 |
а 22 |
^ 2 3 |
| = |
а 2 1 |
^ 2 2 |
« 2 3 |
О 0 |
1 |
Л31 |
а 32 |
О з з |
1 |
« 3 1 |
# 3 2 |
« 3 3 |
Тензор может быть симметричным или антисимметричным по паре индексов, если при их взаимной перестановке его компоненты или не меняются (симметричный) или изменяют свой знак на про тивоположный (антисимметричный).
Поэтому, если 5,- симметричный тензор, то
<?..— S-.
Легко убедиться, что это эквивалентно равенству компонент, расположенных симметрично относительно главной диагонали, так что справедливо следующее:
|
512 |
Si3 |
i |
f |
^11 |
S]2 |
^13 |
5 2] |
5 22 |
^23 |
>— |
j |
^i2 |
^22 |
•S23 |
5 3, |
>^32 |
5 33 |
|
^ |
5)3 |
5гз |
^зз |
Следовательно, у симметричного тензора имеется лишь шесть не зависимых компонент.
Для антисимметричного тензора Л/у должно выполняться усло вие A u = — Ajt .
Ясно, что все диагональные элементы антисимметричного тен зора равны нулю, при i—j имеем А н = — A it, т. е. величина равна
себе самой с обратным знаком, что означает Л,•,= (). Матрица для
антисимметричного тензора имеет вид |
|
|
|
О |
А 12 |
Л13 |
|
Д12 |
О А 2з |
(А.6) |
|
A i з |
А 2з |
О |
|
В этом случае имеется три независимых компоненты.
15 Зак. 112 |
225 |
Любой тензор |
второго |
ранга а ;- |
может |
быть |
разложен |
на |
|
сумму симметричного и антисимметричного тензоров, т. е. |
|
||||||
О / |
' су |
(Li; |
с ji) i о |
0 '^' |
(А ./) |
||
Обозначая Sn ~ |
-75- |
(си 4- сп) и |
(cti |
cjl) |
. видим, |
что |
|
Sj j S -; И A Ij |
Aj j . |
|
|
|
|
|
Все рассмотренные операции над тензорами отражают конкрет пые свойства величин, с которыми приходится сталкиваться в раз личных областях физики и, в частности, гидромеханике.
226
ПРИЛОЖЕНИЕ Б
ЗАПИСЬ УРАВНЕНИЙ ГИДРОМЕХАНИКИ В КРИВОЛИНЕЙНЫХ КООРДИНАТАХ
§ 1. ОРТОГОНАЛЬНЫЕ КООРДИНАТЫ
При аналитическом описании физических явлений, как пра вило, используется координатный метод. Выбор системы коорди нат носит, в известной мере, субъективный характер. Поэтому есте ственно поставить вопрос, во-первых о математической формули ровке какой-либо физической задачи, выраженной в самой общей инвариантной форме, независящей от системы координат; и, во-вторых, о переходе от обобщенной записи к конкретной с использованием той координатной сетки, которая наиболее целе сообразна при решении поставленного вопроса.
Ниже мы приведем основные сведения, касающиеся методов записи уравнений в инвариантной форме и способов их преобразо вания при переходе от одной координатной системы к другой. При этом основное внимание будет уделяться наиболее употребитель ным в практике ортогональным системам координат.
Наиболее часто используются декартовы косоугольные или прямоугольные координаты (х, у, z) = х ь х2, х3). Однако, наряду с ними, широко применяются криволинейные координаты qu q2, q3 *. Они однозначно связаны с декартовыми, т. е.
Q\ = q\ (A'i, х2, х3) ; q2= q 2 (*ь х2,'Х3) ; q3 = q3 (*,, *2, *з)
и обратно |
|
|
xi=xi (?,, |
q2, <73) ; х2= х 2 (qu q2, q3) ; х3= х 3 (qu q2, q3) . |
|
Естественно, |
это предполагает: что якобиан det j! |
„ отличен |
|
■дхк |
от нуля и от бесконечности.
* Переход к криволинейным координатам производится с целью упроще ния рассматриваемой задачи, ибо за счет их удачного выбора можно, например, упростить уравнения или уменьшить число аргументов. Так при наличии осевой симметрии целесообразно ввести цилиндрические координаты. Направляя ось z по оси симметрии, получаем двумерную задачу, так как рассматриваемое явле ние не зависит от угла поворота.
15* |
227 |
У сл о в и е |
c/i = |
c o n s t |
о п р е д е л я е т к о о р д и н а т н у ю |
|
п о в ер х н о ст ь . |
||||||||||||||||||
Я сн о , |
что |
к о о р д и н а т н ы е |
п о в ер х н о ст и , |
с о о т в е т с т в у ю щ и е |
о д н о й |
и той |
|||||||||||||||||
ж е к о о р д и н а т е , н е п е р е с е к а ю т с я м е ж д у с о б о й . |
Н а о б о р о т п о в е р х |
||||||||||||||||||||||
ности , о т в е ч а ю щ и е |
= |
c o n s t |
и |
|
f/. = c o n s l , |
п е р е с е к а я с ь |
о б р а з у ю т |
||||||||||||||||
к о о р д и н а т н у ю |
л и н и ю |
<7/, = |
c o n s l. |
|
К а ж д а я |
|
точка |
п р о с т р а н с т в а |
ф и к |
||||||||||||||
си р у е т с я |
как |
р е з у л ь т а т |
|
п ер ес еч е н и я |
грех |
к о о р д и н а т н ы х |
п о в е р х |
||||||||||||||||
н ост ей |
или |
д в у х |
к о о р д и н а т н ы х |
|
л и н и й . |
В к а ч еств е |
п р и м е р а |
у к а |
|||||||||||||||
ж е м , что |
в ц и л и н д р и ч еск о й |
с и с т е м е |
к о о р д и н а т |
(/?, |
cp, |
z) |
к о о р д и |
||||||||||||||||
н атны м и п о в е р х н о с т я м и |
я вл яю тся : |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
/ ? = c o n s t — к р у г о в о й ц и л и н д р , |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Ф — c o n st — п о л у п л о ск о с т ь , |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
2 = |
co n st — п л о с к о ст ь |
|
п е р п е н д и к у л я р н а я |
оси |
z, а |
к о о р д и н а т н ы е |
|||||||||||||||||
линии : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R — c o n st, |
ф — c o n s t — |
п р я м а я , |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
Ф = |
соп $ Т |
2 = |
c o n st — |
п р я м а я , |
|
|
|
|
|
|
|
|
|
||||||||
|
|
R = c o n st, |
2 — c o n st |
- - |
|
о к р у ж н о с т ь . |
|
|
|
|
|
|
|
||||||||||
Р а д и у с - в е к т о р |
л ю б о й |
|
точки |
п р о с т р а н с т в а |
в |
д е к а р т о в о й |
п р я м о |
||||||||||||||||
у го л ь н о й |
с и с т е м е |
к о о р д и н а т |
м о ж е т |
бы ть |
|
п р е д с т а в л е н |
как |
|
|
||||||||||||||
|
|
|
|
|
|
• ► |
|
-> |
|
|
-> |
-> |
|
-> |
|
|
|
|
|
(Б.1) |
|||
|
|
|
|
|
|
г = |
.\:1г | -)-л'2г2 -f- л'3г'з = |
|
2,-г\, |
|
|
|
|
|
|||||||||
г д е i — орты . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Р а с с т о я н и е |
м е ж д у д в у м я |
|
б е с к о н е ч н о |
|
б л и зк и м и |
точк а м и |
с о о т |
||||||||||||||||
в ет ст в ен н о за п и ш е т с я |
в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
dr=dxi -ij . |
|
|
|
|
|
|
|
|
( Б . 2) |
||||||
П е р е й д я |
к |
о р то г о н а л ь н о й |
|
к р и в о л и н ей н о й |
с и с т е м е |
к о о р д и н а т q\, |
|||||||||||||||||
q2, q3, |
|
|
|
|
|
|
|
|
|
|
|
- ► |
|
|
|
|
|
|
|
|
|
|
|
мы |
м о ж е м |
р а с с м а т р и в а т ь |
dr |
к ак |
д и а г о н а л ь |
эл е м е н т а р н о г о |
|||||||||||||||||
к р и в о л и н ей н о го п а р а л л е л е п и п е д а , |
о б р а з о в а н н о г о к о о р д и н а т н ы м и |
||||||||||||||||||||||
п о в ер х н о ст я м и . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Е с л и ds} |
( / = |
1, 2, |
3) |
есть |
д л и н ы |
р е б е р |
(рис. |
Б . 1 ), |
то |
|
|
||||||||||||
|
|
|
|
|
^ |
|
—у |
|
—> |
|
|
—^ |
|
|
—у |
|
|
|
|
|
|
||
|
|
|
|
|
d r===cls\6\-\'dsoC2~\~ds^6^ ==dSj*Cj9 |
|
|
|
|
(Б . 3) |
|||||||||||||
-г |
— орты |
р а с с м а т р и в а е м о й |
си стем ы |
к о о р д и н а т . |
|
|
|
|
|
||||||||||||||
г д е е j |
|
|
|
|
|
||||||||||||||||||
В ел и ч и н ы |
|
dsj |
м о ж н о |
за п и с а т ь |
ч ер ез |
к о о р д и н а т ы |
dqу , |
в в ед я |
|||||||||||||||
к оэф ф и ц и ен ты |
п р о п о р ц и о н а л ь н о с т и |
# у , |
н а зы в а е м ы е |
п а р а м е т р а м и |
|||||||||||||||||||
Л а м е . |
Т о г д а |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dsj = |
H j ■dqj. |
|
|
|
|
|
|
|
|
(Б .4 ) |
|||||
Т е п ер ь |
в м ес т о |
(Б . |
3) |
б у д е м |
им еть |
в ы р а ж е н и е |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
dr — H j dq} ■е} . |
|
|
|
|
|
|
|
|
, (Б .5) |
228
Д и ф ф е р е н ц и р у я (Б . 5) по к а ж д о й из к о о р д и н а т , п о л у ч и м
дг |
и |
-> |
|
|
~а— — Н г ег |
|
|
||
d(]j |
1 |
1 |
|
|
В о з в е д е м в к в а д р а т о б е части |
п о с л е д н е г о р а в е н с т в а . |
Э т о р а в н о |
||
си л ь н о с к а л я р н о м у у м н о ж е н и ю |
л е в о й |
и п р а в о й частей |
у р а в н е н и я |
|
н а с е б я , т. е. |
|
|
|
|
|
= |
Н) |
-(ej-ey). |
|
И л и , п о ск о л ь к у еу— орты , то {eyeу) = | ву| 2 = I .
П о э т о м у
В св о ю о ч е р е д ь , п р о д и ф ф е р е н ц и р о в а в (Б . 1), п р и х о д и м к с о о т н о ш е н и ю
дг __ дх{ г
229