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

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

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

уравнения

(37.1)

и (37.2)

к такому

виду:

 

 

 

 

 

 

 

 

1

дН

\

|

аЧ\

 

 

дН ,

a.q{q2

 

дН

Є і +

Є

3 ( / . +

Ц

д Х

[

}

I

gL^H-6

 

 

--дхі

g L i l 2 H i

дХ2

— /* j

1

d9l

\

 

а

 

L ^ 1 | ^ 2 \ •

 

^

 

 

 

d<

г

gLffl2

\

п

V

дх,

дх2 )

^

4 2

дх2

 

 

 

VI

|

92^2 ^ dLj

 

qxq2

 

dl2

 

 

(37.5)

 

 

 

у L2

 

L\ J д х і

 

 

L2

дх2

 

 

 

 

 

 

 

 

 

 

 

1(І2

1

 

д Н

\

 

і

а(&

 

 

дН

.

aq,q2

 

дН

L2

 

2

)

 

2L2H3

 

 

дх2

1

gL^H3

 

 

дхх

 

 

 

 

 

 

 

 

— / *

І

1 .

д?2

І

 

 

а

I7

^?2

і

1 і

 

д?2

J^gLxH

 

д*

^

g L ^ y A 1

дх2

д Х х

)+^~дх7

 

 

 

Я2

,

q\Lx

\

дії

 

 

qxq2

dLj_ ,

 

 

^ g

 

 

^ Lx '

Z 2 j д х а

 

 

L l

д Х і

 

 

 

 

Примем

за координатную

поверхность

(хі, х2)

горизонтальную

плоскость. Тогда

єі = є 2 = 0 ,

є3 = — 1,

и в

уравнении

(37.5) будет

 

 

• l + « 8

( / . + - £ •

 

 

дН

\

,

 

*Я\

 

дН

 

 

 

 

 

 

 

дх,

і

1

, / ?b н з

дхх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gLxL\H

 

 

 

 

 

 

 

 

 

 

7

і I

gL\H6

І

 

 

 

 

 

но

 

 

 

 

 

 

,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2Н3

 

\

с

)

 

 

 

 

 

где c = ^2gH — скорость распространения фронта волны неустано­ вившегося движения на глубине Я. Различного рода хозяйствен­ ные водоемы характеризуются тем, что c~^>U и поэтому величиной (U/c)2, а тем более (Ui/c)2 по сравнению с единицей можно прене­ бречь. Тогда, учитывая, что

_ j

1_ _

дН

=

1_ _ d(z

+ Н) =

1_ _ дх.м

= J

1

Li

дхі

 

L\

дхі

Li

дх,

1

есть уклон свободной поверхности по направлению Xi и что ана­ логичные соотношения имеют место для направления х2, можно, ог­ раничиваясь случаем установившегося движения, придать уравне-

ииям (37.4) —(37.6) такой вид:

 

 

 

 

 

 

 

 

dgi

L i ? L ^ o ,

 

 

(37.7)

 

 

дхх

1

дх2

 

 

 

 

 

 

 

 

 

 

 

 

 

дхп

 

ч\

і

llL2 \ dL2

 

qxq2

dL2

 

(37.8)

У^2^

 

L 2

j

дхх

 

L 2

дхх j •

 

 

 

 

 

r і "-71-72

дН _

* ,

 

a

{

dq2 .

dq2

 

Li " r

L \ І дх2

 

 

Ц ' dx2

 

(37.9)

 

 

 

 

Во многих случаях, как и в одномерной

идеализации,

можно

пренебречь так же силами инерции, тогда уравнения

(37.8) и

(37.9)

будут иметь вид

 

 

 

 

 

 

 

 

 

 

 

Л = У ь

 

/ 2

= / 2 .

 

 

(37.10)

Иногда (например, при исследовании бурных потоков) удобнее совмещать координатную поверхность (xi, Xz) не с горизонтальной плоскостью, а с дном потока. Тогда

где ф — отметка дна. Будем опять рассматривать установившееся движение. Уравнение неразрывности сохраняет форму (37.7), ди­ намические же уравнения будут иметь вид:

 

 

-

 

 

1

^

і

£ І

Я

^

^

і

 

«Я\Ч*

 

 

 

дН

.

 

 

 

 

 

 

 

i i

"d 3 ~1 Г

/

дхі

1" Т

^ , / 1

Я

3

'

<?*2

У 1 " Г

 

 

 

 

 

 

 

 

„ ,2

2 „ 3

 

гіг.

 

.

2, 2 г г3

 

Л , .

 

 

 

 

 

 

j t

y i i

 

 

a g i

 

[ _ £ I

 

і

<&h \

 

dL2

qxq2

_dL2

^ / : 2 2 Я 2

j V 1

^ .

1 4 2

 

dx2

 

 

 

 

L2 ) 0xx

 

-Ц-

S x \ y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(37.11)

 

 

 

 

 

 

1

/

,

 

аЯІ

\

dH

,

 

aqtfo

 

 

 

dH

.* .

 

 

 

 

 

 

L2 [ 4

4 r ^ L f F j

 

dx2

+

 

 

 

 

• - d x 7 — j 2

+

 

 

 

 

 

 

 

dq2

 

.

 

dq2

(

ЧІ

 

,

q%

\

 

д

ц

gig2

 

^/1

2

72

Я

2

j

1 '*2

Л ь " i V l

d x ,

^ А, " г Z 2

j

 

 

 

z,

-53^ j

 

 

 

dx2

-<

 

d j f

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(37.12)

251

Величины j * и j * определяются формулами которые при принятых предположениях дают

 

 

 

 

 

 

І2

 

 

dU2

1

Ш7)

^

aY/2

 

 

dx2

i t

 

UK;

L\L2

1 dr2

 

 

 

LxL2

az.2

 

 

 

 

 

™*

Ьхі

 

 

 

dL,

L^L^H^

 

 

 

dx2

 

 

 

dx

 

dqi

£i_

 

 

 

 

<?x2

L,

 

 

 

 

д і ї

 

 

 

 

 

 

dx2

 

dxi

(15.22) — (15.26),

(37.14)

dxo

u2 dx.

<li . .

Z.2 dx,

dZ,2

(37.15)

Существенное упрощение всех выкладок и расчетов по уравне­ ниям (37.11) и (37.12) получается, если принять єз~-—1. Считая допустимой погрешность в 5%. имеем

е 3

1 Е 1 £ 2

1

s2 + s2

>0,95,

откуда вытекает,

что касательная

к

линии наибольшего падения

рельефа дна должна составлять с горизонтальной плоскостью угол

не более 18°. Принимая

ез = — 1 ,

мы должны одновременно

не де­

лать различия мел-еду:

 

 

 

 

 

 

1) синусом

и тангенсом

угла

наклона к горизонту касательной

к координатной

линии

(т. е. считать, что е.і и е 2 — н е синусы, а тан­

генсы соответствующих

углов);

 

(х\,

х2)

2) длиной дуги кривой

на координатной поверхности

и длиной проекции этой дуги на

горизонтальную плоскость;

 

3) кривизной кривой

на

координатной поверхности {хи хг)

и

кривизной

проекции этой

кривой на горизонтальную плоскость.

 

Правые

части уравнений

(37.11) и (37.12) остаются при этом

без изменений, а левые приобретают соответственно вид

1

, Л

а<їл \

дНdH } ,

Li

dx.

 

dxx

 

 

 

dH

 

dx?

' gL\H*

dx2

aqxq2

dH

(37.16)

gL.L2Hs

dxo

 

аЯ\Яъ

dH

(37.17)

gL]L2H3

dx,

 

Интегрирование уравнений двумерного потока почти во всех случаях, имеющих практическое значение, выполнимо только ко­ нечно-разностными (сеточными) методами. Эта задача громоздка и кропотлива, но принципиальных затруднений для современных вычислительных средств не представляет. Трудности расчета тече­ ний в водохранилищах и озерах возникают не в связи с интегри­ рованием уравнений двумерных потоков, а в связи с постановкой краевых условий. Рассмотрим два случая, в которых эти трудности носят различный характер: 1) водоем с вертикальными ограждаю­ щими стенками, очерченными в плане по достаточно регулярным кривым; 2) естественный водоем с отлогими берегами. В первом случае возможны два варианта краевых условий на ограждающих непроницаемых стенках:

1) нормальная к стенке составляющая Un вектора скорости U равна нулю;

2) вектор скорости U равен нулю.

В первом варианте касательная к стенке составляющая скоро­ сти U, вообще говоря, не равна нулю и, следовательно,'-на криво­ линейных участках стенки нужно допустить существование каса-' тельных напряжений второй системы (§ 16), вызванных кривизной линий тока. Во втором же варианте нужно допустить существова­ ние на стенке касательных напряжений первой системы, связан­ ных с шероховатостью стенки. Но касательные напряжения на ог­ раждающих стенках вообще противоречат идее двумерной идеали­

зации, которая

отводит ограждающим

боковым

стенкам

чисто

геометрическую

роль. Это противоречие

немедленно устраняется

с переходом к идеальной жидкости, когда первый

вариант,

£УП = 0,

полностью согласуется с двумерной идеализацией, или с переходом к трехмерной идеализации, с которой полностью согласуется вто­ рой вариант U = 0 . Но в рамках двумерной идеализации для тече­ ний с гидравлическими сопротивлениями рассматриваемое проти­ воречие неустранимо, и не остается ничего другого, как принять тот или иной из указанных вариантов, закрывая глаза на получаю­ щуюся при этом физическую неувязку. Картины течения, отвечаю­ щие этим двум вариантам краевых условий, отличаются одна от другой только у стенок, на транзитном же потоке выбор варианта краевых условий практически не сказывается/ С вычислительной точки зрения вариант 11 = 0, по-видимому, удобнее.

В замыкающем створе /—1 водохранилища (рис. 66) обычно известна глубина и расход на единицу длины створа. От возможного перекоса водного зеркала в этом створе обычно отвлекаются. До тех пор, пока отсутствует боковое сжатие, как на рис. 66, никаких за­ труднений при этом не возникает. Если же в замыкающем створе имеет место боковое сжатие (рис. 67), то.краевые условия, строго говоря, могут быть поставлены только из гидродинамического ана­ лиза этого сжатия. Однако, учитывая, что локальные особенности потока в замыкающем створе не могут сказываться на общей картине течения в водохранилище, в практических расчетах целесообразнее

видоизменить очертание границ потока так, чтобы получить истече­ ние из водоема без бокового сжатия (пунктир на рис. 67).

Во входном створе 2—2 (рис. 66) должно быть задано только распределение расхода по створу, но, конечно, не глубина, которая должна быть получена в результате расчета.

Рис. 66. К постановке краевых условий для дву­ мерного потока.

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

протекающей через водоем. Но главная трудность заключается

не

 

в этом, а в изрезанное™ берегов есте­

 

ственных водоемов (рис. 68). Рассуж­

 

дая принципиально, мы можем интег­

 

рировать

двумерные

уравнения

для

 

всей области, показанной на рис. 68,.

 

поставив краевое условие: на урезе

 

воды

U = 0

(в данном случае это усло­

 

вие достаточно точно отвечает действи­

 

тельности

и не противоречит двумер­

1

ной

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

Но извилистость

1границы области интегрирования соз­ дает очень серьезные затруднения вы­

Рис. 67. К постановке краевых

числительного

характера,

которые мо­

условий в замыкающем

створе

гут перерасти

в принципиальные

труд­

при наличии

бокового

сжатия

ности. Преодолевать эти трудности, по-

вытекающего

из

водохранили­

ща

потока.

 

видимому, не имеет практического-

А — действительное

положение

смысла,

ибо

в

заливах

образуются

крышки водослива,

В — измененное

зоны с

ничтожными скоростями

(за­

положение.

 

 

 

 

 

стойные

зоны),

которые

не влияют

 

 

 

 

на основной

транзитный

поток.

По­

этому можно заменить фактическую границу потока некоторой

сглаженной

кривой

(пунктир на рис. 68),

на

которой

в каче­

стве краевого условия можно

принять U = 0.

Строгий

анализ

ошибки, получающейся в результате такой

замены,

требующий

отыскания

взаимно

однозначного

отображения

одной

на

другую-

областей, очерченных фактической и сглаженной границами, очень сложен. Практически можно рекомендовать сопоставлять расчеты

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

скорости или по конфигурации линий тока).

Если картина

течения

в какой-либо области между фактическим и

сглаживающим

конту­

ром (штриховка на рис. 68) все же будет представлять практиче­ ский интерес (например, для выяснения времени смены массы воды в такой области), можно произвести дополнительный расчет тече-

Рис. 68. К постановке краевых условий для по­ тока в водохранилище с отлогими естественными берегами.

ния в такой области по упрощенной системе уравнений (37.7) и (37.10) для контура, ограниченного фактической границей и неко­ торой близкой к сглаживающей границе линией тока АВ, получен­ ной в расчете транзитного течения (рис. 68). При этом краевые ус­ ловия должны быть следующими: на линии тока АВ и в сечениях 3—3 и 4—4 скорости равны скоростям, полученным расчетом тран­ зитного потока, на фактической границе U = 0 .

Переход бурного течения в спокойное приводит к необходимости постановки специфических краевых условий или, точнее, условий контакта области, в которой c>U, и области, в которой c<U. По­ ток между двумя близкими линиями тока вектора U можно рас­ сматривать как поток в непризматическом русле, а переход из обла­ сти c<U в область о U — как прыжок в непризматическом русле. При этом, в соответствии со сказанным в § 30, в рамках представле­ ний о медленной изменяемости течения мы должны считать, что прыжок происходит в определенной точке на линии тока, причем глубины по разные стороны от этой точки связаны условием (30.5). Обозначая индексами «с» и «б» величины, относящиеся

к спокойному и бурному течению, можно записать это условие так:

aq2 НІ aq2 НІ

где <7н. = (7б Яб = с/с Я0 . В отличие от одномерного случая, здесь Ц*Ф Ф const, а изменяется вдоль линии тока. Трудность расчетов тече­ ний с переходами из бурного в спокойное состояние связана с тем, что граница, разделяющая области этих состояний, заранее не из­ вестна и ее следует выявить в процессе расчета.

В заключение необходимо заметить, что Ы. М. Вернадский [3, 4]-—первый исследователь плановых задач гидравлики — выделял водоворотные зоны. Тогда, естественно, появляются краевые условия на границе водоворота. Однако необходимость выделения водоворотных зон была вызвана принятой им моделью двумерного (планового) потока, методом решения задачи и ручной техникой вычислений. Описанная здесь модель отражает все свойства потока, способные привести к образованию водоворотных зон, а сеточные методы интегрирования автоматически обнаруживают эти зоны. Поэтому выделения водоворотных зон и постановки специальных краевых условий на их границах для данной модели не требуется.

38. О влиянии рельефа дна на свойства течения. Поворот потока

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

. Направим координатные линии Хі по линиям тока вектора U, вернемся от неизвестных qi и qz к прежним неизвестным Ux и Uz и будем считать выполненными условия, при которых левые части

уравнений

(37.11)

и

(37.12)

имеют вид (37.16) и (37.17). Учиты­

вая, что при выбранной координатной

сетке

Ui = U, Uz = qz=j* =0 ,

получаем

из

(37.8),

(37.11) и (37.12):

 

 

 

 

 

 

 

1

dHU

 

НЦ

dLz

 

 

 

 

 

тхдкт^-цп,дхТ—°-

 

 

(38Л>

 

 

 

Ц I. дхх

" Г V

gH j дхх\ —

 

 

 

\D +

g LXL2

дх2 j

' gL\H

дхх

(d8'2>

 

 

1

 

д(<г +

Н) _

oU2

1

дЦ

 

, „ R ~

 

 

Ц

'

dx2

g

" LXL2

' ~dxT

^Q-6>

Так как

L.] dx^—

ds\,

Z#2 dx2ds2,

 

— —jj—

* -^—-

Ri '

 

 

 

 

 

Z.1Z.2

длг2

 

 

1 _

_

дЦ =

1

 

 

 

 

IjZ.2

'

dxi

R2

 

 

(где dsi, i?i и

ds2,

Rz — элементы дуг

и радиусы

кривизны соот­

ветственно линий тока и ортогональных к ним кривых), то урав­

нения (38.1),

(38.2)

и (38.3)

могут

быть написаны

в

виде

 

 

 

 

 

 

 

 

Ш и

Ии-=0,

 

 

 

 

 

 

 

(38.4)

 

 

 

 

 

dsi

R2

 

 

 

 

 

 

 

 

 

 

 

-ЖГ +

I і

—W)

=

ІТГ +

 

 

 

даті)

 

 

• " S T '

( 3

8

- 5 )

 

 

 

 

 

д (9 4- Я )

 

 

at/2

 

 

 

 

 

 

(38.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть

теперь X i и хо — по-прежнему

произвольные

ортогональ­

ные криволинейные координаты. Уравнение линий тока

(2.2)

в этих

координатах

будет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L\ dx\

 

L 2

dx2

 

 

 

 

 

(38.7)

 

 

 

 

 

Ux

 

 

U2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

производная

по направлению

s

любой

скалярной

функции

F —

=

F(xi, х2)

выразится так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dF

cos

(іГлгі )

dF

 

.

cos

(s, x2)

 

dF

 

 

/ 0

0

o\

 

 

 

 

 

- L T ^ -

- d x

T

+

 

T~2

' ~~дхї '

 

(

3 8 - 8 )

 

В качестве

координат xi

и x2

 

введем

полярные

координаты

в

и р (рис. 69), тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

хх 9,

х2=р,

cfsj == рof©,

ds2=dp,

7-i =

p,

£ 2 = 1 ,

 

 

 

 

 

Ut

= Ue,

Ut=Uf,

i / P / £ / e

= t g a = / ( P >

9),

 

 

 

 

где a — угол между линией тока и координатной линией. Уравнение линий тока тогда будет

«ли -=р/(р, в),

; (38.9)

а уравнение линий, ортогональных к линиям тока, напишется так:

Далее

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

 

 

 

 

1

 

 

.

 

/

 

 

 

 

COSGt =

,

 

sin<x=-

 

 

 

 

dF

 

p '

dQ +

/

dp

 

dF _

p

" d8 +

dp

 

d s l

 

/ 1 + / 2

'

 

^ 2

 

/ 1 " + / 2 "

 

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

формулы

 

 

 

 

 

 

 

 

1

Р2

+ 2р2 — рр

 

 

 

 

 

 

 

 

~~R

( р 2 +

р 2 ) а / 2

 

 

 

 

и уравнений

(38.9) и (38.10)

имеем:

 

 

 

 

 

1 _ ^ 7

9 1

dp

дв

1 _ ^ 1 +

7

+ /

dp

дв

Rl

 

р (1 +

/ 2 ) Э Ь

 

'

А?2

7

 

(1 +

/ 2 ) 3 / *

 

Учитывая

все эти соотношения,

можно

представить

уравнения

(38.4) — (38.6)

в таком виде:

 

 

 

 

 

 

 

 

 

/

р

df

df

 

 

 

 

'

 

т

1

d?

. f , dy?

,і Л(л

at/2U W\ I 1

 

дЯІdH\ ., f

P

d9

' dp

+V

gH ) \ p

'

d6 + 7

j - 0 , '-(38.!.)

,dH ) dp j )

 

Yl+f2

 

 

,'

f

df

, df

£/2 +

 

 

і

X

,

P /

dp+

d9

 

 

,

all <

1

dHU_,fdHU_\

 

m

9 .

 

+177 П

 

 

d e - + 7 - ^ p — >

 

( 3 8 л 2 >

/

d ( y + Я ) . d (у + H)

 

at/ 2 і 1

 

 

 

іd e — + — щ ^ - g v v — и 1 ? 2 ( 3 8 л 3 )

Это и есть искомая форма уравнений движения. Эти уравнения отличаются от уравнений предыдущего параграфа тем, что вместо трех неизвестных Н, Uі, ІІ2 в них фигурируют три неизвестные Н, U, f а также тем, что выбрана конкретная система координат —

полярная.

то к трем неизвестным Н, U, f до­

Если рельеф дна не задан,

бавляется четвертая — ф, и

система трех

уравнений (38.11) —

(38.13) становится не замкнутой. Это значит,

что на поток можно

наложить некоторое дополнительное условие. При проектировании поворотных и рассеивающих виражей в качестве такого условия обычно принимают #=const . При #=const система уравнений

(38.11) — (38.13)

после

преобразования

уравнения (38.12) с по­

мощью уравнения

(38.13) приобретает такой вид:

 

 

dU

 

 

dU

Uf

 

 

р

df

 

 

 

 

 

 

 

 

/

 

 

 

 

(38.14)

д@

 

 

dp

р

 

 

 

1 +

/ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дв

 

 

 

 

 

 

 

 

/ 1

+ / 2

и2-

аЦ

 

 

 

 

 

 

(38.15)

 

 

D

 

g

p

'

dQ

>-!

dp

j

f

ду

.

с5а

 

aU2

I 1

nf

d f

J-

d f

'

(38.16)

 

 

і

+P

 

= 0 .

р

дв

'

др

 

 

 

 

 

 

 

Рассмотрим

поворот

бурного потока. Пусть до поворота, т. е.

до сечения СіСг на

рис. 69,

и после

поворота,

т. е. после

сечения

Рис. 69. Система полярных координат для поворота по­ тока в плане.

ЕіЕг, поток течет в прямоугольных призматических лотках, уклон и ширина которых в общем случае могут быть различны. Границы

потока на повороте, т. е. пунктирные дуги CiEi

и С2Е2, должны

быть заданы. При этом, по условиям безотрывного

обтекания, они,

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

и после поворота. Очевидно, что

эти

дуги будут

линиями тока.

Пусть

уравнения

дуг С\ЕІ

И CZEZ

имеют вид соответственно

 

 

Pi—Pi (в),

Р2=Р 2 (Є) .

(38.17)

Из

уравнения

(38.9) следует краевое условие: на дугах (38.17)

должно быть соответственно

 

 

 

 

 

Pi

(9 )

/ =

Р2 ( В )

(38.18)

 

 

Р і ( 0 )

Р2 ( в )

17*

259

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