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

книги из ГПНТБ / Картвелишвили Н.А. Потоки в недеформируемых руслах

.pdf
Скачиваний:
7
Добавлен:
25.10.2023
Размер:
12.62 Mб
Скачать

Далее

г„

- й / 2 '

г*

- й / 2

_ [ Л [ „ ( 4 . ) _ . ( _ 4 ) ] _ 0 .

так как вектор скорости на смоченной поверхности равен нулю. Здесь b — ширина русла на высоте z над осью х, а г^. — значение zo при у = 0. И, наконец,

 

В/2

 

Н

 

 

В/2

 

l J d ^ d F

= 1 а

У \

 

J b l - d

z =

1

wHdy.

 

-В/2

 

 

 

 

 

-В/2

 

Используя теперь

выражения

(13.3)

и

(13.17) и учитывая, что

dF = BdH, получаем

dUF

 

,

dF п

 

 

 

 

 

 

 

( ] 3 2 4 )

 

дх

 

 

dt

=0.

 

 

 

 

 

 

 

 

 

Учитывая это уравнение и принимая, как обычно в гидравлике, что полный корректив количества движения а = а' + а" постоянен, нетрудно найти

Г"

1

д

\udF-

дх

и2

dF\

= 1

dU

 

 

dt

F

 

 

 

£*

dt

 

 

 

 

 

 

 

 

 

 

 

 

U

 

dU

 

L

JL

dF

 

 

 

 

g*

 

дх

 

 

F

dt

 

 

Наконец,

в качестве

последнего

 

предположения

отбросим

в уравнении

(13.14)

два последних члена в фигурных

скобках, со­

держащих квадраты вертикальных составляющих скорости. Осно­

ванием для этого

упрощения

является

малость вертикальных со­

ставляющих,

которая

имеет

следствием

то, что отбрасываемые ве­

личины имеют второй порядок

малости. Теперь уравнение (13.14)

окончательно может

быть записано так:

 

 

 

 

 

dH

 

dL-i

..k

, 1

dU

4

 

dU

 

dx

 

dx

 

 

dt

1

g*

dx

 

 

U_

dF

1

 

d*

- F*

 

 

, dH

+

g*

F

dt

g*F

 

dtdx

В

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dH

 

 

dH

 

 

(13.25)

 

дх*

 

 

dx

 

 

dt

 

 

 

 

 

 

 

 

 

 

Уравнения

(13.24)

и (13.25)

и есть уравнения турбулентного по­

тока в одномерной

идеализации.

 

 

 

 

 

 

Вернемся к предположениям, которые оказались

необходимыми

для вывода уравнения (13.25):

 

 

 

 

 

 

 

1) русло имеет вертикальную плоскость

симметрии;

 

2) корректив количества движения а постоянен. Кроме того, хотя в процессе вывода уравнения (13.25) нам не пришлось делать никаких предположений о коэффициентах .6^, Вг, рз, 84, но для ре­ шения этого уравнения такие предположения необходимы. Прос­ тейшее из возможных и пока единственное предлагавшееся пред­ положение такого рода состоит в том, что

3) коэффициенты B'j, В2, Вз, В4 постоянны.

Второе и третье предположения часто встречаются в гидравли­ ческой литературе. Первое предположение в этой литературе не упоминается, но как видно из нашего вывода, оно довольно су­ щественно: во всяком случае, применение уравнения (13.25) к рус­ лам, не удовлетворяющим этому предположению, вызываемое не­ избежностью, есть определенная натяжка.

В гидравлической литературе большое значение придается двум гипотезам, которые мы не использовали:

а) давления распределяются по глубине потока по гидростати­

ческому закону, т. е. в правой части выражения

(13.5) исчезающе

малы все члены, кроме g$(H 2);

 

б) течение обладает свойством медленной

изменяемости как

в пространстве, так и во времени.

 

Эти гипотезы вызываются не существом дела, а следующими

обстоятельствами: 1)

принятым в гидравлике

способом

вывода

уравнений движения,

который сводится к тому,

что законы

меха­

ники применяются к потоку в целом или (что теперь уже устарело) к так называемой элементарной струйке (само это понятие не отве­ чает физической сущности турбулентного потока); 2) гипотеза мед­ ленной изменяемости течения порождается также и тем, что ее реа­ лизация есть необходимое условие допустимости определения гидравлического уклона по формуле (13.21), а входящей в нее про­ пускной способности русла К—-по обычным формулам гидравлики. Эти формулы, полученные из опытов над равномерным течением, нельзя, конечно, экстраполировать на течения, сильно отличаю­ щиеся от равномерного.

Заметим, что при обычном для гидравлики способе вывода урав­ нений движения не учитываются силы, связанные с искривлением линий тока локальных скоростей, отражаемые в (13.25) членом D и членами, содержащими коэффициенты В. И лишь после того, как уравнения движения таким путем получены, в них вводятся соот­ ветствующие поправки. Это, конечно, не строгий прием. Следует, однако, иметь в виду, что силы, связанные с искривлением линий тока, во многих задачах несущественны.

Оценку коэффициентов р можно получить, рассматривая поток в бесконечно широком прямоугольном русле. Если принять, что горизонтальные и вертикальные составляющие скорости распре­ делены по вертикали по степенному закону, т. е.

/ z \i/ft'

k + l r 7 f z у / А

/ z \>п

где Um — скорость на поверхности; U — средняя скорость на вер­ тикали; и т — постоянные, то

О

 

(от + 2) +

1 >

 

4

'

' '

 

 

1

1

 

 

 

 

 

 

1

 

 

 

 

 

 

от +

2

'

 

 

о

 

л

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

0

 

А

 

* (ш -г 2) + 2 '

б

й

*

 

+

2) +

1 '

 

 

/г (от + 2)

+ 1

 

 

(коэффициент (5j будет играть роль в главе V I ) . Для

распределения

скоростей и по глубине в гидравлике распространен так называе­

мый закон

одной седьмой

(к = 7). Однако в потоках с

большим

уклоном

и

быстрым течением, которые будут

рассматриваться

в главе V I , распределение

скоростей по вертикали близко

к равно­

мерному

(k — oo). Коэффициенты |3, значения

которых

даются

в табл. 1, отвечают так называемым ондуляциям, т. е. вторичным волнам, накладывающимся на основные волны. Анализ волновых движений в § 9 позволяет отметить, что эти вторичные возмущения не распространяются на большую глубину, т. е. что значения т должны быть достаточно большими, по-видимому, /?іЗ>5, т. е. ко­ эффициенты (3 должны быть достаточно малыми числами.

В заключение необходимо остановиться еще на следующем. Усредняя уравнения Рейнольдса по площади сечения потока с ве-

сом единица, мы

получили

уравнение

(13.25). Можно

выполнить

такое усреднение

и с другой

весовой функцией, например и. Тогда

нужно было

бы, подставив

р

и рн

в уравнение

(12.13)

из уравне­

ний (13.5) и (13.6), умножить получившееся уравнение

на и и про­

изводить усреднение по F после такого

умножения. В

результате

вместо уравнения

(13.8) мы получили бы выражение

 

 

 

 

 

 

 

(

"

 

 

 

 

 

F

 

 

 

 

да ,

ди*

, duv ,

д2

Н

...

Н

,„

д-Ф \

С

. д* Г

 

 

 

 

 

 

 

 

(13.26)

Производя с уравнением (13.26) те же операции, которые были сделаны с уравнением (13.8), придем к уравнению, которое будет отличаться от уравнения (13.25) только коэффициентами. В част-

Т а б л и ц а 1

 

Коэффициенты ондуляций

 

 

 

 

т.

 

 

>

10

20

 

 

 

 

р;

0,318

0,140

0,082

0,045

Pi

0,680

0,299

0,176

0,097

0,333

0,143

0,084

0,045

Р2

Рз

0,043

0,020

0,012

0,006

Р4

0,363

0,160

0,094

0,052

 

 

k = 10

 

 

р;

0,323

0,141

0,083

0,045

Pi

0,667

0,294

0,173

0,091

h

0,333

0,143

0,084

0,045

Рз

0,031

0,014

0,008

0,005

Р4

0,344

0,153

0,090

0,045

 

 

k — оо

 

 

РІ

0,333

0,143

0,084

0,045

Pi

0,666

0,286

0,167

0,091

Р2

0,333

0,143

0,084

0,045

Рз

0

0

0

0

Р4

0,333

0,143

0,084

0,045

мости, производная dU/dt войдет в это уравнение с множителем а',

а субстанциональная

производная UdU/dx будет множиться не на

ос, а на коэффициент

Кориолиса

 

а * = " ? 7 7 Г f u * d F >

 

F

дополненный поправкой, связанной с пульсацией скоростей. Все эти коэффициенты, вообще говоря, переменны, но в рамках гидравличе­ ской идеализации их принимают постоянными, как и коэффици­ енты уравнения (13.8). В этом, собственно, и заключается смысл гидравлической идеализации: учесть изменчивость этих коэффици­ ентов можно только с переходом к трехмерной модели потока, а если такой переход будет сделан, то вводить исследуемые коэф­ фициенты не надо.

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

Если учесть их изменчивость, то при усреднении с любыми весами мы будем получать всегда один и тот же результат: усреднение с произвольным весом есть чисто математическая операция, не свя­ занная с какими-либо упрощающими предположениями, а такая операция не может давать результаты, зависящие от произвольной весовой функции. Вследствие этого рассуждения о том, какой ве­ совой функции следует отдать предпочтение, могут иметь только одни аспект: при какой весовой функции предположение о посто­ янстве коэффициентов, получающихся за счет перехода к одномер­ ной (или к двумерной, как в § 15) идеализации, дает наименьшую погрешность. Очевидно, что ответ на этот вопрос не однозначен, а зависит от того, о какой конкретной задаче идет речь: в разных задачах оптимальные весовые функции будут разными. Поиски оп­ тимальных весовых функций для задач разных классов могут соста­ вить предмет специального исследования. Однако опыт показывает, что простейшая весовая функция — постоянная, равная единице — дает для всех задач, к которым применима одномерная идеализа­ ция, результаты, если не самые лучшие, то по крайней мере хоро­ шие. Поэтому здесь и в дальнейшем за весовую функцию прини­ мается единица.

14.Уравнение напорного турбулентного потока

водномерной идеализации

Рассмотрим течение в трубе переменного сечения, имеющей две плоскости симметрии, из которых одна вертикальна. Прямую, по которой пересекаются плоскости симметрии, примем за ось х и на­ правим ее по течению. Ось zi разместим в вертикальной плоскости симметрии, направив ее вниз. Тогда ось у будет горизонтальна, J = g-sintp, У = 0, Z = g-cos\p, где ip — угол между осью х и горизон­ тальной прямой в плоскости xOzi. Уравнения (12.13) и (12.15) на­ пишутся так:

ди

,

ди2

,

duv

,

duw

. ,

1

dp

Іл..

1 Ч

-дГ + ^Г+-ТГ+—=£sm6—г

 

 

 

(14.1)

dw

,

dwu

,

dwv

,

dw2

,

1

dp

,

n 4

-ЬТ+-1£Г+—+—=£COs6—r

 

 

 

(H.2)

Здесь формулы для сої и соз те же, что и в предыдущем пара­ графе. Уравнение (12.14) в данном случае не представляет инте­ реса по тем же причинам, что и в случае открытого потока.

Из уравнения (14.2) находим, полагая, что на оси симметрии р = ро, w = 0,

р, |

,

С dw ,h С dwu ,„

ґ

о

о

z,

z,

или, так как z\ не зависит от t, х,

у, то

 

 

 

z,

г,

рр Рор

1 • gzx cos ф

~д~Г~ J ^ ^—5т] t<y/J^"

 

•г,

 

 

 

• J" таг;

—- w2 j (о3 dfj.

(14-3)

оо

Подставляя это значение р в (14.1) и интегрируя результат по площади сечения F потока, получаем после деления на F

1

С

да

,,-, .

1

['

да2

,,-, .

1

f

daw ,

.

1

дра

,

F

J

dt

1

/•

.'

дх

'

F

.)

dz,

6

о

дх

1

 

 

 

 

J ^ W E T f * * +4- J d / 7

- ^ г 1

 

 

 

l ^ - s - K ^ - r l " ^ . -

<14-4>

Из условий симметрии, а также учитывая, что на смоченной по­ верхности вектор скорости равен нулю и что F не зависит от t, можно найти, проделав преобразования интегралов, аналогичные преобразованиям предыдущего параграфа, что все интегралы в пра­ вой части уравнения (14.4), кроме интегралов, содержащих сої и «з, и последний интеграл в левой части равны нулю. Так как

J ~dT d F ==~dTJ

u d F =

= — d T ~ = F ~ dt

F

F

 

то уравнение (14.4) можно записать так:

+ - И * / ? т И , » з * .

04.5)

Член

расшифровывается совершенно аналогично тому, как это было сде­ лано в предыдущем параграфе, в результате чего получим из (14.5)

 

 

• — S T -

I

с"

S

= Sintb

• — — { - ] * .

(14.6)

 

 

g

dt

1 gF

 

дх

-{

дх ' J

\

• t

Здесь

у — объемный вес жидкости. Но

UF = Q и в

напорных

потоках

не зависит от х. Полагая

а = const, получаем

 

 

 

 

 

1

_

д (aU*F)

=

gQ

dU

U

ди

 

 

 

 

 

gF

 

дх

 

gF

дх

g

дх

'

 

 

Учитывая,

что sini|) = —dzo/dx,

где го — отметка

оси

трубы,

и

отбрасывая индекс нуль

при z и

р, представим

уравнение (14.6)

в окончательном виде

 

 

 

 

 

 

 

 

 

 

 

 

Т ' Т + ^ + і + ^ Н * -

 

 

0 4 . 7 ,

Здесь /* — гидравлический уклон. В случае установившегося те­

чения

(dU/dt = 0)

при отсутствии гидравлических

сопротивлений

(/* = 0)

это уравнение

интегрируется по х,

в результате

получаем

 

 

 

 

 

z+ - 7 - + - TTr - =const .

 

 

 

(14.8)

Уравнение

(14.8)

внешне

совпадает с

уравнением

Бериулли

(8.4). Но между этими уравнениями есть глубокая разница: урав­ нение (8.4) относится к линии тока идеальной жидкости, уравне­ ние же (14.8) ко всему потоку реальной жидкости при условиях, когда сопротивлениями можно пренебречь.

Предположения о симметрии потока, сделанные при выводе уравнения (14.7), существенны: без этих предположений соответ­ ствующие интегралы в уравнении (14.4) не обратятся в нуль и уравнение (14.7) не будет иметь места. Поэтому применение этого уравнения к напорному потоку, не удовлетворяющему условиям симметрии, должно рассматриваться как приближение.

15.Уравнения открытого турбулентного потока

вдвумерной идеализации

Двумерная идеализация приемлема, когда глубины потока нич­ тожны по сравнению с главными радиусами кривизны свободной поверхности и поверхности дна, например для построения планов течений в искусственных и естественных водоемах, расчета водо­ сливных поверхностей двоякой кривизны, проектирования виражей на быстротоках и т. д.

Предполагая, что это требование выполнено, и учитывая урав­

нения (12.17)

и (12.18), примем за координатную поверхность

(xi,

Хг) некоторую

гладкую поверхность так, чтобы глубины потока

по

нормалям к этой поверхности мало отличались от глубин по нор­ малям к свободной поверхности и к поверхности дна. За координат-

иые линии л'з примем нормали к поверхности (Л'І, А'г), а за коорди­ наты А'з расстояния по этим нормалям. Тогда коэффициент Лямэ Ьз будет равен единице, а коэффициенты L , и L 2 можно считать не­ зависящими от А-3. Учитывая все это, уравнения (12.17) и (12.18) можно записать в развернутом виде:

dt

 

 

dx,

 

 

 

 

 

 

 

 

dx2

 

 

 

}+

 

 

дх3

 

 

 

P

іг\(^-—)

•41

--дх7+

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

чі

 

1

^ 2 1. . p

 

1

dx\

 

 

(15.1)

 

 

 

 

 

Li

 

 

 

 

 

 

 

 

 

 

dx,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dvr,

1

 

 

Li \v2vJ

4 \

 

 

 

 

LXL2

\ v\

41

 

dt

L \ L

i

dx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx3

 

 

 

T23

1

 

 

 

 

^22

1

dxo

 

 

 

 

 

 

 

^2

 

 

 

 

P

/-2

 

 

 

 

 

 

 

 

 

дЦ

 

•Fo-

 

1

 

_P

 

(15.2)

 

 

 

 

 

 

Z.1

dx.

 

 

І 2

dx2

I p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

Z,

0

u

 

x3l

 

 

 

 

L x [v3v2

 

 

dt

LXL2

\

dx

 

[V,V3 l

 

 

 

 

dx2

 

 

 

 

 

 

 

 

dx3

 

 

T 33

 

•-F,

 

 

d

P_

 

 

(15.3)

 

 

 

 

 

 

p

 

 

 

dxz

I p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dL2v,

 

dLiv2

 

^ L

 

^ - = 0 .

 

 

(15.4)

 

 

 

 

dx,

 

 

dx<>

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исключим из этих уравнений давление р. Интегрируя уравне­ ние (15.3) по А'з, получаем

_Р_

Р

 

 

 

- 2 - к

Зи1~

^31

 

 

 

 

 

 

 

 

VoV,

 

 

4-'°«.

•1»з-

^33

 

(15.5)

 

 

 

Здесь р*, А*, v#,

—значения

р, хз, из, т 3 3

на поверхности по­

тока, £ — переменная

интегрирования.

 

 

 

Используя формулу дифференцирования

под знаком интеграла

и учитывая очевидное соотношение

 

 

 

дх±

 

dx%

dx*

(15.6)

дх.

L i

дхо

dt

 

7 Заказ Hi 428

97

в котором a j и

— значения Vi и и2

на поверхности потока, вы­

ражение (15.5)

можно привести к виду

л"»

 

Р

 

 

 

 

'*

/ ? , ( Я - х 3 + г ) + 4 -

1

*s ^ + T 7 Z r X

р

Р

' я ^ ' - л я Т ^ П

37- I

ч-'я"?"?--; ,

 

X

i f

 

 

 

\

д

f

 

 

 

- u § + W , .

(15.7)

 

- \ L 2 \

vzvxdi

 

 

 

U . J - w f t

 

 

 

J

 

 

 

 

 

 

 

 

х3

 

 

 

\ х3

 

 

 

 

 

 

Здесь z— значение

Хз на дне потока,

Н=х^

— z — глубина по­

тока,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v 33

L,L2

 

дх,

 

хъ

 

дх<>

L, J т3 2 д'£

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т зі

 

дх*

х32

дхл

 

 

 

(15.8)

 

 

 

 

 

 

 

дх,

 

дх2

'

 

 

 

 

 

 

 

 

 

 

 

 

 

хзі

и т з 2

— значения тзі

и Тз2 на свободной

поверхности.

 

 

Подставим р из (15.7) в (15.1) и проинтегрируем результат по

глубине потока Я. Полагая,

как и в случае

одномерной идеализа­

ции, что на дне потока

вектор скорости равен нулю, а стало быть и

1>з = 0, получаем

— f

 

 

 

і

f

5 V 2 . aLr3

 

 

d

v i

dL

Lv\

 

 

 

 

1

 

Г*

 

 

 

 

 

 

 

 

 

^ 2

І

дх,

> Л з ~ г "727-J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1^2

 

 

 

 

 

 

 

 

 

 

 

 

а я

• ) - Т Г І ^ з - щ ^ - І

« 8 ^ -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CUT]

І і / . 2

 

 

 

 

 

1

2

1 L,]\v2d^

 

 

| +

- z r j

- ^ - ^ а — Е Г - Г - ^ - Л С . .

(15.9)

 

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Г*,.

rf„

1

f

dL,L2x„

,

/ - i l 2

cbr- j

.22ахъ

— pL

j

^

а х , -

 

 

 

 

1^2

г

 

 

 

\

 

 

 

 

 

^\^0 2

T+ — значение Тіз на дне потока.

Снова используя формулу дифференцирования под знаком ин­ теграла, а также формулу (15.6), приведем уравнение (15.9) к виду

- я [ Л + ґ , ( / , + ^ - ^ ) ] + « , + ^ Г ^ , +

^~dT

 

 

 

 

 

 

L.\

V{V2

dx.

 

 

dLx

^ - j

v\dx3

 

 

 

L4

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x*

i'*

 

 

 

 

 

1

 

 

 

 

 

 

 

 

52

 

 

 

 

 

 

 

j

v\dx3

 

 

 

 

j

j

vzd%dx3-\-

 

 

 

L\L2

 

dx\

 

 

 

dt

dxx

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

Л'з

 

 

.

T

д

f

r

і

1

 

 

 

 

^ j - ^ 2 . f j ^ 1 ^ ^ 3 J +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z.,

J

j г)3 гі2 rfS fl0c8

І

+

7 i j

г » 3 - с ; 1 ^ з Н - / 2

| г»8г»2 dx3 ) -f-

 

 

 

\

2

-Is

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z.jZ-2

 

 

Z-2 j

v^Oi dx-.

 

dx2

 

Z-i j"

 

v3v2dxz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

5

f

2

.

,

,2

U**

 

 

 

 

 

 

^

 

(15.10)

 

 

 

 

•77

- ^ T J ^ ^

8

+

- J - - - s f - = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

(52

Здесь / 2 определяется

аналогично /і,

т. е. h—- L 2

 

dxz'

 

Пусть

U — с р е д н е е по глубине Я

значение вектора

скорости по­

тока, a Ui и £У2 проекции вектора U на направления

x i и х2 . Оче­

видно,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

•»I = C L / C O S ( U , xl)=cU,

 

 

v2=cUcos(U,

 

x2)=cU2,

где

с — функция координат такая, что

 

 

 

 

 

 

 

 

 

 

 

 

Г*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

с

 

dx3=H.

 

 

 

 

 

 

Следовательно,

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

х$

 

 

 

 

 

 

 

 

 

 

•f*

 

 

 

 

 

 

 

 

 

 

 

 

 

|

<ovdx3—HUх,

 

[

•o\dxz=a.

HU\,

 

 

 

 

 

 

z

 

 

 

 

г

 

 

 

 

 

 

 

 

 

 

J

v22dx3=a'HUl

 

 

)

v^dx^a/HUiUi,

 

 

(15.11)

 

 

 

 

 

 

 

a ' =

- ^ J

C 2 rf,V 3 .

 

 

 

 

7*

99

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