Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Палагин, Э. Г. Основы гидромеханики учебное пособие для метеорологов

.pdf
Скачиваний:
14
Добавлен:
20.10.2023
Размер:
8.24 Mб
Скачать

характер турбулентности определяется градиентом скорости и кри­ визной ее профиля, от которых зависит возможность п степень впхреобразованнй в потоке. Анализ размерностей приводит нас к простои формуле

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

Соседние файлы в папке книги из ГПНТБ