ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 07.07.2024
Просмотров: 146
Скачиваний: 1
Исходные предпосылки работы Дж. Смита [93] близки к пред посылкам Энгелунда. Так же как Энгелунд, он рассматривал пло ское осреднениое турбулентное движение жидкости с постоянным коэффициентом турбулентной вязкости. Расход наносов опреде лялся по формуле М. Ялина. Записав выражение для касательного напряжения по (3.75), Смит без достаточных оснований отбросил
на линии возмущенного дна производную dw и в конечном ре
зультате пришел к заключению о невозможности установить до минирующую длину волн. Наряду с этим, работа Смита содержит обстоятельное решение задачи о сопротивлении донных волн.
В своей совокупности исследования, о которых говорилось в § 3.1 и 3.2, позволяют высоко оценить вклад, внесенный методом малых возмущений в решение проблемы о происхождении волно образного рельефа подвижного дна. Будучи применен после многих лет исканий п догадок, метод малых возмущений поставил, нако нец, решение этой проблемы на научную основу.
§ 3.3. В О З М У Щ Е Н И Я с Б О Л Ь Ш И М И Д Л И Н А М И волн
Формирование плёсов и перекатов, меандрирование и ветвление русла есть процессы, начало которым кладет неустойчивость русла по отношению к возмущениям с большими длинами волн. Длины этих волн имеют порядок ширины потока. Современная динамика русловых потоков не располагает, однако, средствами математи ческого описания плановых деформаций русел. Поэтому вопрос о начальной устойчивости русла по отношению к длинноволновым возмущениям удается пока что поставить лишь применительно к следующей сильно схематизированной модели. Пусть дан прямо линейный равномерный поток в русле прямоугольного сечения весьма большой ширины В с подвижным дном и недеформируе мыми берегами. По дну могут перемещаться гряды и рифели. Ус тойчивость течения будем исследовать в приближении плановой задачи движения руслового потока, т. е. определяя поле скоростей с точностью до средних скоростей на вертикалях. Вариациями ско ростей, обусловленными грядовым рельефом дна, будем пренебре гать.
Таким образом, гряды и рифели не будут играть какой-либо другой роли, кроме роли элементов макрошероховатости. Пусть координатная плоскость (х, у) расположена горизонтально, ось х направлена вдоль оси потока, ось г вверх, средний уклон дна мал. Наложим на плановые линии тока малое синусоидальное воз мущение. Следствием искривления линий тока будут отклонения свободной поверхности от плоскости (рис. 3.10). Вдоль береговых линий будут чередоваться повышения и понижения свободной по верхности. На участках, где свободная поверхность понижается, течение будет ускоряться и произойдут размывы дна. Там, где сво бодная поверхность повышается, течение будет замедляться и про
изойдет отложение наносов. В результате возникает шахматная цепочка углублений дна (на рис. 3.10 они показаны штриховкой) и шахматная цепочка отложений. Этот возмущенный рельеф дна уси лит искривление струй, возмущения будут нарастать, пока не сфор мируется устойчивая извилистая линия наибольших глубин. В ра ботах по морфологии рек в этих случаях говорят о русловом меандрировании или о меандрировании меженного потока, понимая под этим извилистость течения внутри меженного русла. Если бе рега легко поддаются размыву, то следующим шагом будет ис кривление берегов, т. е. переход к пойменному меандрированию.
Можно себе представить, что в русле очень большой ширины параллельно развиваются два или больше возмущений планового движения. Поток, а с ним и русло разделятся по ширине на не сколько фрагментов, каждый из которых будет повторять картину,
Рис. 3.10. Возмущения свободной поверхности при искривлении плановых линий тока.
изображенную на рис. ЗЛО. Этот процесс положит начало делению русла на рукава.
Постановка и решение задачи об устойчивости русла в описан ном виде принадлежит Р. Калландеру [59]. Перейдем к математи ческой формулировке задачи. Исходной будет система уравнений планового движения потока в деформируемом русле (2.27) — (2.30). Перепишем эту систему, исключив высоту дна с помощью равен ства 2 s = z,„ — /г и представив продольный градиент потерь энергии,
в соответствии с первым из равенств (2.25), в виде |
|
|
/ / = - |
Jfo |
(ЗЛ00) |
h ’ |
||
где то— проекция на ось х вектора |
касательного |
напряжения на |
дне (индекс х у величины то, а также у проекции вектора удель ного расхода наносов qs будем опускать). Получим следующие уравнения:
dz^j |
|
dU \ |
(ЗЛ01) |
|
дх |
~~ |
dt J’ |
||
|
||||
dz^ ___ |
|
(ЗЛ02) |
||
~ду |
g |
|
||
|
|
63
|
|
|
h U ) |
d (h V ) |
elk |
|
|
(3.103) |
|
|
|
д (дх |
дѵ |
dt = о , |
|
) =0 . |
|
|
т М - |
dqs |
д |
д (zw |
h |
(3.104) |
||
|
дх |
ду М + |
dt |
|
||||
Невозмущенное и возмущенное движения должны удовлетво |
||||||||
рять системе |
(3.101) — (3.104). В невозмущенном движении имеем: |
|||||||
U = Uü, |
У= 0, |
h = h , то=(т0)о, <7s=?so, |
zw= z w{0 ) — І0х, |
причем |
||||
|
|
|
|
Ы о= ^ рѴ 0. |
|
|
|
(3.105) |
Возмущенное движение определяется равенствами: U=Uo+U', |
||||||||
V = V', |
h = h o + h', |
То= (т0)о+т', qs= |
qs0 + q', |
zlo= z w(0) — I0x + |
+z'w. Перенеся эти равенства в систему (3.101) — (3.104), исполь
зовав (3.105) и отбросив члены, нелинейные по переменным за дачи, получим
1 |
( тоLo |
-г |
дСГ |
dU’ |
(3.106) |
дх |
Рho |
— g h ho |
U0 дх |
dt |
|
ду |
|
|
g |
, |
( |
Uo |
\ |
4)/U |
fd i/o |
ax |
H |
dt |
||
|
dU’ |
. n |
^h’ |
?h0 |
|||||||||||
|
"Q |
. |
дх |
|
U ° |
dx |
|
hr, |
|
V ’ . |
dh’ |
— О |
|||
|
|
|
|
dy |
‘ |
dt |
|||||||||
|
|
|
|
* |
|
|
|
-4 - '*0 |
|
|
|
||||
1 |
1 |
|
\ d4s |
|
|
[ |
gsO |
|
d V ’ |
\ |
|
dzw |
|
dh’ |
|
|
1 |
|
|
|
|
|
|
1 |
|
|
|
||||
|
дх |
|
|
1 |
u 0 |
|
ду |
) |
dt |
|
dt |
||||
|
|
|
|
|
|
|
, |
|
|||||||
|
— E |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
)■, (3.107)
(3.108)
0.(3.109)
Чтобы замкнуть систему (3.106) — (3.109), свяжем возмущения касательного напряжения и расхода наносов с возмущениями про дольной скорости. Малость амплитуд возмущений позволяет вос пользоваться линейными соотношениями:
|
|
|
|
|
Ч>- -m \U , |
|
|
|
|
|
(3.110) |
|||
|
|
|
|
|
Ь |
- = т М ' . |
|
|
|
|
(3.111) |
|||
При учете этих |
соотношений |
и |
равенства |
(3.105) |
система |
|||||||||
(3.106) — (3.109) переписывается в следующем виде: |
|
|
||||||||||||
dz.. |
■ |
Ч |
пцЦ’ |
|
|
|
|
|
|
|
|
|
|
|
дх |
(то)о |
|
£ |
- ] + |
|
т |
№ |
+ 3 |
4 - |
(ЗЛ12) |
||||
|
дг„ |
|
=/о |
V |
|
|
|
|
d V ’ |
d V ’ |
|
(3.113) |
||
|
ду |
|
и. |
|
|
|
дх |
|
у |
|
||||
|
|
|
|
|
о |
|
|
дѴ> |
|
|
|
|
|
|
, |
dU’ |
■U, |
|
dh’ |
|
-Ао |
|
dh’ |
=o, |
|
(3.114) |
|||
|
|
0 |
|
|
=0. |
|||||||||
Л°~дЗГ |
Qso____ |
|
|
ду |
|
, |
dt |
dh’ |
(3.115) |
|||||
|
|
|
0 |
дх |
|
|
|
|
|
|
||||
dU ’ |
|
|
|
d V ’ |
|
dz., |
|
|
|
|||||
|
(1 _ « ) f/( |
|
|
|
|
|
|
|
||||||
т* -д Т |
|
|
ду |
|
|
|
dt |
dt |
|
|
||||
|
|
|
о |
|
|
|
|
|
|
64
Возмущения элементов движения представим выражениями:
|
і А /„ (x—cUat) |
(3.116) |
|
U'—S (y)U0e |
п° |
||
V' = |
і А-/„ (х—Ы/оО |
(3.117) |
|
T(y)U 0e |
ll° |
|
|
|
1 |
(x-ci/at) |
(3.118) |
h’— R (у) ІіФ |
й° |
||
. |
i |
!„ (x—cUd) |
(3.119) |
Zv—Z (y) hQe |
0 |
Здесь S, Т, R, Z — безразмерные комплексные функции. Волно вое число k и комплексная скорость с также безразмерны. Вслед ствие соотношений (3.110) и (3.111) возмущения касательного на пряжения и расхода наносов находятся в фазе с возмущениями продольной скорости. Возмущения V', h' и z'w могут быть смещены
по фазе относительно U'.
Подставив выражения для возмущений (3.116) — (3.119) в урав нения (3.112) — (3.115), после дифференцирования и сокращения экспоненты получим:
Гm\Uo |
I |
)] s - - R - |
|
L Hob |
1■ i k Fr0(l —c |
||
|
|
|
^ |
[ \ - \ - i k |
Fr0(l — с)] T - . = — |
©IoJ>| |
|
|
|||
s + (1 — c ) R = i |
Ao |
|
|
A/Q |
|
QsO
S + c { R - - Z ) - 1 ( l - - e) Uо/?о
\ - i |
k Z — |
■ 0 , |
|
(3.120) |
|
d Z |
’ |
|
(3.121) |
|
dy |
|
||
|
|
|
||
dT |
|
|
|
(3.122) |
dy |
У |
|
|
|
|
ho |
dT |
|
(3.123) |
|
kJ0 |
dy |
• |
|
|
|
Система уравнений (3.120) — (3.123) дает возможность решить задачу об устойчивости планового движения потока в эродируемом русле. Опуская громоздкие преобразования, приведем общий ход решения.
Исключая из уравнений (3.120), (3.122) и (3.123) величины 5
dT
и R, получаем уравнение, связывающее Z с производной
dy
Решая это уравнение совместно с (3.121), исключаем Z и получаем обыкновенное дифференциальное уравнение для безразмерной амп литуды поперечной скорости
сРТ |
А Г =0, |
(3.124) |
|
дуг |
|||
|
|
где Л — параметр, сложным образом зависящий от элементов не возмущенного движения, волнового числа k и комплексной скоро сти с. Интерес представляют собственные значения параметра А,
5 Зак. № 550 |
65 |
т. е. те его значения, при которых уравнение (3.124) имеет нетри виальные (не нулевые) решения. Чтобы найти собственные зна чения Л, используем граничное условие обращения амплитуды Т в нуль на линиях берегов
7’ ( - | - ) = 7 ( ---- |
§ - )= 0 . |
(3.125) |
Из условий (3.125) следует, что собственные значения Л дол жны быть вещественными и положительными
R eA > 0 , Іш А =0.
Это сразу приводит к решению
Т—А cos даг-g -, /і= 1 , 3, 5, . . . . |
(3.126) |
где А — постоянная. Собственные значения Л определяются выра жением
(3.127)
Приравняв правую часть (3.127) выражению Л через элементы невозмущенного движения, волновое число k и скорость с и от бросив некоторые малые члены, получаем уравнение, в котором ве щественная и мнимая части могут быть разделены. Это дает два алгебраических уравнения для определения сти с,; тем самым мы решаем задачу об устойчивости.
Уравнения для сг и Сі ставят эти величины в зависимость от сле дующих шести безразмерных параметров, из которых первые че тыре характеризуют невозмущенный поток, а два последние — воз мущения:
Fr0, |
Т1Ъ\ Up |
m2 |
4 s0 . |
|
(то)о |
AQ ’ |
7/QAQ |
|
k, tva |
AQ |
|
|
|
IQB ■ |
|
Число n, входящее в последний параметр, представляет собою число полос с шахматными дорожками экстремумов zw и zs, на ко торые разбивается возмущенный поток. Значение п = 1 отвечает меандрированию, значения /г > 1 — делению на рукава.
Если касательное напряжение на дне определять по формуле
Шези, то |
ml =2gp |
halo |
Если расход наносов выра |
|
и г |
||||
|
|
|
||
жается степенной зависимостью вида |
(2.21), то mz = m-yr—. |
|||
|
|
|
и о |
Таким образом, зная величины £Л>, /го, /о и qso и построив связь между qSQи ІІ0, мы можем найти все четыре параметра, характе ризующие невозмущенный поток.
66