
книги из ГПНТБ / Картвелишвили Н.А. Потоки в недеформируемых руслах
.pdfРассмотрим единственный член
г
в (15.10), содержащий компоненты тензора напряжений. Положим, аналогично тому, как это было сделано при одномерной идеализа ции,
^ , = g - t f ( / + Л , |
( 1 5 Л 2 ) |
где величина /* учитывает гидравлические сопротивления, а /** —
количество движений пульсаций. Можно предположить, опять-таки как и в случае одномерной идеализации, что влияние слагаемого /**
(и /** для уравнения по направлению xz) выражается в некотором увеличении корректива количества движения а', т. е. что величину /** можно отбросить, если вместо а ' ввести полный корректив коли
чества движения а. К гидравлическим |
же сопротивлениям |
/* мы |
вернемся ниже. |
|
|
Положим V3 = fVj., где / — некоторая |
функция координат, |
анало |
гичная функции с. На основании (15.6) можно написать |
|
|
W I F + Z - ^ - . |
• - & Ч " £ - - g - |
|
и далее |
|
|
гx2
x*
j |
v^dx^HU^Vy^L |
|
|
+ ^ w ) , |
|
/ = 1 , 2, |
||
|
h=-RT$ |
§ |
f(xi> |
x*> 6. |
i)a\dx3, |
|
||
|
|
г |
хг |
|
|
|
|
|
Рп = |
- Ж . [ |
I с (*і> |
x*> 5. |
0 / C * i . *2. |
E, |
t)d%dx3, |
||
|
.'z |
x3 |
|
|
|
|
|
|
P I I I = - 7 7 2 " I |
J с 2 ( л " і > |
^2. |
S. |
x2, |
Є, |
t)d$,dx3, |
z хг
1 |
г |
|
|
|
|
Pv=77-J |
c(xlt |
x2, xa, |
t)f(xu |
x2> xa, |
t)dxz, |
z |
|
|
|
|
|
$Vl = ~ff\ |
£ 2 ( - * l , |
x2> x3> |
fyf(Xl> |
X2> X3, |
t)dxz. |
z |
|
|
|
|
|
Обозначим через Єї, Єг, єз направляющие косинусы вертикали, направленной вниз. Тогда Fi=g&i, Fz^g&z, F$=g&3. Наконец, обра тим внимание на то, что поперечная скорость оз по крайней мере на порядок меньше скорости U, поэтому в (15.10) можно пренебречь членами, содержащими vz3. Учитывая все сказанное и полагая, как
и в одномерной гидравлике, что корректив |
количества |
движения а |
|||||||
и коэффициенты |
Pi, ... , pvi постоянны, |
можно привести уравнение |
|||||||
(15.10) к виду |
|
|
|
|
|
|
|
|
|
- ^ [ . , + в а ( / , + ^ Г - - ^ - ) ] + в - Я / + |
dHU, |
|
|||||||
dt |
|
||||||||
dLxL2HU\ |
|
|
|
dLxHU\U2 |
|
дЦ |
них— |
||
L\L2 |
dx. |
|
L,L |
дхо. |
|
дх. |
|||
|
|
L,LL |
|
dL2 |
HUl+A^O, |
|
(15.13) |
||
|
|
2 |
дхі |
|
|
|
|
|
|
|
|
\ |
|
|
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
d* |
|
|
|
|
|
|
|
dx7 |
Ьк '-2 |
• ~kdxj |
|
(Pi. ^ г + Р ш |
w)+ |
|||
|
1 |
d |
|
L,H\U2 |
^„ |
dH |
P . u |
^ |
|
|
LiL2 |
dx2 |
dt |
|
|||||
|
|
|
|
|
1 |
d |
|
|
|
dt |
|
|
|
|
|
dxx |
|
|
|
dH |
|
|
|
h_ |
|
|
dH |
+ |
|
J vi W |
|
L\L2 |
dx. |
|
|
|
|||
|
|
|
|
|
|||||
|
|
|
|
|
dH |
|
|
|
|
|
dx2 |
|
|
dt |
|
|
|
|
Если в уравнении (15.13) заменить индекс 1 индексом 2, а ин декс 2 индексом 1, то получится второе динамическое уравнение двумерного потока. Выписывать это уравнение в общем виде нет необходимости.
Интегрируя уравнение неразрывности (15.4) по Н и используя (15.6), получаем соответствующее уравнение для двумерного по тока
dL2HUx |
дЦНЦ2 |
(15.14) |
|
dx. |
dx2 |
||
|
101
Для избыточного давления на дно р+— р% из (15.7) нетрудно по лучить следующее выражение:
|
н |
+ |
+ |
1 |
дН |
|
|
L\L<i \ |
дхх |
|
|
|
дН |
|
(15.15) |
дхо |
|
|
|
|
|
|
|
Член N+ в этой формуле есть изменение давления |
на дно потока |
||
за счет турбулентности |
и вязкости, т. е. очень |
малая |
по сравнению |
с остальными членами величина, которой можно пренебречь. Если считать, что разрывная прочность жидкости практически равна нулю, как это всегда и бывает в обычных, не лабораторных, усло виях, то для отсутствия отрывов потока от дна необходимо и доста точно, чтобы правая часть выражения (15.15) была положительной при всех Л'І и А'2. В дальнейшем всюду предполагается, что это ус ловие выполнено.
Полученные уравнения двумерного потока занимают примерно такое же положение по отношению к уравнениям гидродинамики и одномерной гидравлики, какое занимают уравнения теории оболо чек по отношению к уравнениям теории упругости и сопротивления материалов.
Рассмотрим теперь гидравлические сопротивления. Уравнение (15.13) с точностью до постоянного множителя есть проекция урав нения равновесия действующих в потоке сил на координатное на правление Xi. Следовательно, выражение (15.12) есть проекция не которой силы (опять-таки с точностью до постоянного множителя)
на направление х\, а значит и член /* в (15.13) должен |
рассматри |
|||||
ваться как проекция на направление хі некоторого |
вектора j , по |
|||||
рождаемого гидравлическими сопротивлениями. |
Направление по |
|||||
следнего совпадает с направлением вектора скорости |
U, т. е. имеют |
|||||
место формулы |
|
|
|
|
|
|
У u\ |
+ u\ + ul ' |
л * = ; |
и > |
+ ul |
(15.16) |
|
|
у u\ + ul |
|
||||
где |
-і'* |
|
|
|
|
|
|
|
дН |
|
|
|
|
|
|
|
|
|
(15.17) |
|
|
|
|
|
|
|
|
Если бы линии |
тока и траектории |
вектора U |
были |
прямоли |
нейны, то величина / была бы вполне аналогична уклону гидравли ческого трения в одномерном случае, т. е. можно было бы написать, аналогично формуле (13.21), j = U2/C2H, где С —коэффициент Шези, определяемый через глубину потока Н и коэффициент шеро
ховатости дна п, например по формуле Маннинга |
С = Н'Іс/п. |
В этом |
|||
случае касательные |
напряжения на |
плоскостях, |
ортогональных |
||
ко дну и проходящих |
через траектории |
вектора |
U, |
равны |
нулю. |
Но при искривлении этих траекторий касательные напряжения на линейчатых поверхностях, ортогональных дну, для которых траек тории вектора U служат направляющими, не могут оставаться рав ными нулю. Если бы эти напряжения оставались равными нулю, в потоке не могли бы существовать водоворотные зоны: все эти зоны не получали бы энергию от основного потока и, затрачивая свою кинетическую энергию на преодоление сопротивлений, порождае мых дном, должны были бы в конце концов превратиться в застой ные зоны, тогда как в действительности наблюдаются устойчивые водоворотные зоны. К выводу о том, что искривление траекто рий вектора усредненной скорости вызывает возникновение допол
нительных касательных напряжений, |
|
||||
приводит |
также |
факт |
существова |
|
|
ния местных потерь напора, не за |
|
||||
висящих |
от шероховатости |
стенок |
|
||
потока: эти потери возникают толь |
|
||||
ко там, где искривляются траектории |
|
||||
вектора скорости. |
|
|
|
||
Таким |
образом, необходимо допу |
|
|||
стить, что j имеет еще и второе сла |
|
||||
гаемое, |
прямо |
пропорциональное |
|
||
(так же как и первое) |
U2 и обратно |
|
|||
пропорциональное радиусу |
кривизны |
|
|||
R линий тока вектора U, т. е. |
|
||||
|
1 |
т |
|
(15.18) |
Рис. 16. К формуле Вейсбаха для |
|
|
|
|
||
|
g\R\ |
|
|
потери напора на закруглении |
|
|
|
|
|
|
|
Параметр % в этом выражении |
трубы. |
||||
можно |
определить |
из |
формулы |
|
Вейсбаха для коэффициента сопротивления в закруглениях трубы (рис. 16). Эта формула, приводимая в любом гидравлическом спра
вочнике, имеет вид |
|
С = [ 0 , 1 3 Ц - 0 , 1 6 5 ( - | - ) 3 ' 5 ] - ^ . |
(15.19) |
Рассматривая бесконечно тонкую трубку, охватывающую тра |
|
екторию вектора скорости, следует положить d/R-^-0, |
тогда фор |
мула (15.19) приобретает вид |
|
=0,262—, |
|
Потеря напора на пути RQ равна |
|
и, следовательно, дополнительный гидравлический уклон, вызывае мый искривлением линии тока, определяется по формуле
h |
0,262 |
1 ^ = 0 , 0 4 1 7 - ^ , |
||
RQ |
|
|||
|
gR |
' |
gR ' |
т. е. Л = 0,0417.
По отношению к выражению (15.18) необходимо сделать два за мечания.
1. Мы предполагали, что траектория вектора скорости есть пло ская кривая. В общем случае это — пространственная кривая (когда дно не плоское). Следует ожидать, что выражение (15.18) сохра няет силу и для пространственной кривой с радиусом кривизны R .
2. Выражение (15.18) было бы связано с действительными ка сательными напряжениями, если бы траектории вектора v локаль ной усредненной скорости лежали на ортогональной ко дну линей чатой поверхности, проходящей через траекторию вектора U сред ней скорости. В общем случае это условие не выполняется. Поэтому выражение (15.18) отражает только эффект некоторых фиктивных напряжений, который должен быть таким же, как и для действи тельных напряжений. Поэтому может оказаться, что эксперимен тальное исследование двумерной идеализации приведет к некото рому изменению численного значения коэффициента X.
С учетом выражения (15.18) имеем
Til |
111 |
|
|
fi=uui(-w+-gTm)- |
|
|
|
( 1 5 - 2 0 ) |
|||
Поскольку |
|
|
|
|
f/з |
|
|
|
|
|
|
|
R |
|
|
|
|
|
|
|
|
|
= |
— |
, |
|
|
(15.21) |
|
где w — модуль векторного произведения |
вектора |
скорости |
на век |
||||||
тор ускорения |
|
|
|
|
|
|
|
|
|
w=y (и2и3- |
и3и2)2+{и3их- |
|
ихи3)2+(ихи2-и2их)\ |
||||||
|
|
|
|
|
|
|
|
|
(15.22) |
то |
|
|
|
|
|
|
|
|
|
,«_ |
UUi |
, I |
UjW |
' |
-*_UU2 |
і |
Ujw |
/іцоо\ |
|
J) |
D |
' g ' U2 |
J 2 |
D ~i~ g ' |
U2 ' |
U J - z , - 7 |
|||
Если в качестве Ui, U2, U3 взять компоненты полного ускорения, |
|||||||||
то R будет радиусом кривизны траекторий вектора |
U, а если опери |
||||||||
ровать с субстанциональными |
компонентами ускорения, то R будет |
радиусом кривизны линий тока. Напряжения, вызывающие гидрав лический уклон /, можно увязать с кривизной траекторий. Поэтому, учитывая что
I = 1 |
d L \ _ |
д І 2 _ |
ди^ |
_ д и 2 _ д и г |
п |
3 |
' дх3 |
дхй |
дх3 |
дх3 |
дх3 |
напишем
dUx |
1 U l |
dUx |
, u2 |
dUx |
|
|
m 2 |
dL2 |
|
, uxu2 |
dLx |
' |
||||||
dt |
1 |
Ц |
dxx |
1 |
L 2 |
dx |
2 |
|
|
dx |
x |
1 |
LxL2 |
dx |
2 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(15.24) |
||
dU-> і |
u2 |
dU2 |
|
ux |
dU2 |
|
|
и] |
dlx |
|
|
. UXU2 |
dL2 |
' |
||||
dt |
1 |
І2 ' |
dx2 |
1 |
Ц |
dx |
x |
|
|
LXL2 |
dx |
2 |
|
' |
L\L2 |
dx |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(15.25) |
||
|
|
Of |
dU3 |
, |
Ui |
dUz |
, |
Uo |
dx |
|
• |
|
(15.26) |
|||||
|
|
dt |
|
L 2 |
|
|
||||||||||||
|
|
|
1 |
Ц |
dx |
x |
1 |
2 |
|
|
|
|
||||||
Существование |
дополнительных |
|
сопротивлений, |
вызываемых |
||||||||||||||
искривлением |
|
линий |
тока, |
приводит |
|
к |
ряду |
новых |
вопросов |
даже |
||||||||
в одномерном случае. На этих |
вопросах мы частично |
остановимся |
||||||||||||||||
в следующем параграфе, а более подробно — в главах V I и V I I . |
|
|||||||||||||||||
16. О замыкании системы уравнений |
Рейнольдса |
|
|
|
Как уже указывалось в § 12, конструирование гипотез, позволяющих замк нуть уравнения Рейнольдса, оправдывается только для тех задач, которые прин
ципиально |
требуют |
трехмерной трактовки и притом |
не могут быть рассмотрены |
в рамках |
теории идеальной жидкости. Таких задач современная гидрологическая |
||
и гидротехническая |
практика выдвигает множество. |
Сюда относятся, например, |
все задачи, связанные с течениями в глубоких водохранилищах, где нельзя абстрагироваться ни от распределения скоростей в плане, ни от их распреде ления по вертикалям, ибо на разных глубинах течения могут быть направлены даже в противоположные стороны.
Гипотезам о замыкании уравнений Рейнольдса посвящена большая литера тура, в частности, капитальная монография [45]. Попробуем построить гипо тезу, которая давала бы более полное замыкание для частного случая руслового потока, включая сюда и течения в водохранилищах. Это значит, что из всех при
чин возникновения турбулентности выделяются только две: влияние |
стенок по |
|||
тока и влияние |
кривизны линий тока (т. е. в конечном |
счете |
местных |
сопротив |
лений). Случаи, |
когда турбулентность возникает вне |
этих |
причин |
(например, |
в океанических течениях), а также случаи отрывного обтекания стенок предла гаемой гипотезой не охватываются. Мы поставим себе целью не только замк нуть уравнения Рейнольдса, но и удовлетворить краевым условиям на твердых границах и свободных поверхностях. Необходимо специально подчеркнуть, что
речь |
идет только о гипотезе, а не о чем-либо |
большем, |
и во |
всяком |
случае |
|
ке о |
теории. Поэтому |
излагаемую ниже конструкцию ни в коем |
случае |
нельзя |
||
считать чем-то бесспорным п окончательным. |
|
|
|
|
||
|
Ко всякой новой |
гипотезе предъявляются |
следующие |
требования: |
1) она |
не должна противоречить уже установленным и проверенным фактам и теориям; 2) она должна объяснять факты, которые не могли объяснить прежние гипо тезы и теории. Но, построив новую гипотезу, удовлетворяющую этим требова
ниям, мы всегда должны быть готовы к тому, что могут |
появиться |
новые факты, |
|||
не укладывающиеся |
в эту гипотезу, |
и тогда ее придется |
обобщать |
или заменять |
|
другой, |
более совершенной гипотезой. |
|
|
||
Для |
построения |
интересующей |
нас гипотезы введем понятие |
призматиче |
ского потока, понимая под этим поток, в котором линии тока усредненных мест ных скоростей есть параллельные прямые. Очевидно, что необходимым условием призматичности потока является призматичность русла и однородность шерохо ватости дна и стенок вдоль образующих. Из § 11 ясно, что это условие не явля ется одновременно достаточным.
Концепции классической гидравлики в вопросе о сопротивлениях безусловно верны для двух видов призматических потоков с однородной по периметру
шероховатостью: 1) в круглой трубе, 2) в бесконечно широком прямоугольном русле. Эти случаи будут служить в качестве эталонных.
Как указывалось в предыдущем параграфе, искривление линий тока и траек торий вектора локальной усредненной скорости должно вызывать дополнительные касательные напряжения. Поэтому можно предположить, что в непризматическом потоке существуют две системы напряжений: первая, связанная с шероховато стью границ потока, аналогичная существующей в призматическом потоке, и вто рая, зависящая от кривизны линий тока и не связанная с шероховатостью не посредственно.
Мы рассмотрим сначала призматические потоки, что приведет к некоторым отправным положениям, а затем обратимся к общему случаю непризматнческого потока. При этом мы будем всюду иметь в виду поток с развитой турбу лентностью и пренебрегать вязкими напряжениями.
Основной принцип учения о гидравлическом турбулентном |
трении выра |
жается известной формулой |
|
тп |
|
в которой U — средняя по сечению скорость потока; у — удельный |
вес жидкости; |
т 0 — среднее касательное напряжение на смоченной поверхности; |
С — коэффици |
ент Шезн. Классическая гидравлика коэффициент С ставит в зависимость от
гидравлического потока |
R и коэффициента |
шероховатости |
смоченной поверхно |
||
сти п. Так как в силу |
бесспорных |
законов |
механики xc=yRI, |
где / — гидравли |
|
ческий уклон |
потока, |
то из (16.1) |
вытекает известная формула для расхода |
||
Q = UF=K У7, |
в которой |
|
|
|
K-=CFYR |
(16.2) |
— пропускная способность русла. |
|
Нетрудно заметить, что формула (16.2) при изложенной |
трактовке коэффи |
циента Шези в некоторых случаях приводит к невозможным для призматиче ских потоков результатам. Действительно, пусть К а и Кб — пропускные способ-
|
6) |
|
V |
|
V77Y |
/У777777777777, |
77777777777/V |
Рис. 17. К основному парадоксу класси ческой гидравлики.
ности труб, имеющих сечения, представленные соответственно на рис. 17 а и б. Если определять С по формуле Маннпнга, то элементарные выкладки приводят к тому, что
|
К6=*чКа, |
« = |
/ _ 0 ± ^ > £ |
|
|||
|
|
b |
|
|
. |
Л |
|
|
|
— |
|
• |
Х = ^ - - |
(16.3) |
|
При малых |
а, т. е. при малых |
Ь, ср<1, т. е. Кб<Ка, в то время как физи |
|||||
чески очевидно, |
что КвЖа, |
если |
в |
потоке |
нет поперечных |
(циркуляционных) |
|
течений. Такое |
противоречие |
не дает |
себя |
знать на практике |
по причинам, кото- |
рые будут ясны из дальнейшего. Но тем не менее оно свидетельствует по край ней мере о неполноте классических концепций.
Попробуем рассмотреть более детально касательные напряжения в призма тическом потоке. Примем за координатные линии л:, линии тока вектора усред
ненной |
скорости, |
а за |
координату х.\ |
расстояние вдоль этих линий. |
Тогда |
|||||||
д/дхі = 0 |
для |
всех |
переменных, |
кроме |
давления, |
Li = l, |
v2=v3=0. |
Для |
потока |
|||
в поле силы |
тяжести |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
p ^ i — ? |
xz- |
|
= — Р £ 7 |
= — i f 7 |
|
|
|
|
и из (12.17) |
имеем |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
(/-3^12) |
+ • |
|
з)] = |
77 - |
|
(16.4) |
|||
|
|
|
І2Ц |
дхл |
|
|||||||
|
|
|
і дх2 |
|
|
|
|
|
|
|
||
Уравнения, получающиеся из |
(12.17) |
при £=2 и £=3, содержат |
р и |
опреде |
ляют отличия закона распределения давления по сечению потока от гидростати
ческого закона. Эти уравнения |
пока не потребуются. |
|
|
|
|
|
||||||||||||
Координатные |
линии |
х2 |
и |
х3 |
лежат в |
плоскости |
сечения потока. Примем |
|||||||||||
за линии |
-Ї2 изотахи, |
а за |
линии |
х3 |
семейство |
линий |
ортогональных |
изотахам. |
||||||||||
Это есть траектории вектора градиента мо |
|
|
|
|
|
|
|
|||||||||||
дуля скорости ||v|;|=Ui, мы их будем назы |
|
|
|
|
|
|
|
|||||||||||
вать просто линиями |
градиента |
(рис. 18). |
|
|
|
|
|
|
|
|||||||||
Рассмотрим напряжение Т12, т. е. напряже |
|
|
|
|
|
|
|
|||||||||||
ние, действующее |
|
в направлении |
координат |
|
|
|
|
|
|
|
||||||||
ных линии х\ по |
площадкам, |
нормальным |
|
|
|
|
|
|
|
|||||||||
к координатным |
линиям |
х% Иначе |
говоря, |
|
|
|
|
|
|
|
||||||||
Ті2 — касательное |
|
напряжение |
на цилиндри |
|
|
|
|
|
|
|
||||||||
ческих поверхностях, |
для которых |
линии |
|
|
|
|
|
|
|
|||||||||
градиента есть направляющие. Для потока |
|
|
|
|
|
|
|
|||||||||||
в круглой трубе или в |
бесконечно |
широ |
|
|
|
|
|
|
|
|||||||||
ком прямоугольном русле вследствие сим |
|
|
|
|
|
|
|
|||||||||||
метрии |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
( ^ 1 2 ) = |
0, |
|
^,2 = |
0. |
(16.5) |
Рис. 18. |
Изотахи и линии |
гради |
|||||||||
дх2 |
|
|
|
|
|
|
|
|
|
|
ента модуля скорости в призмати |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Естественно |
предположить, |
что COOT |
|
|
ческом потоке. |
|
|
|||||||||||
ношения (16.5) имеют место н |
в |
общем |
|
|
|
|
|
|
|
|||||||||
случае, так как иначе существовали |
бы силы, приводящие |
к тому, что на линиях |
||||||||||||||||
градиента |
было бы dvddx2=^=0, что противоречит физическому смыслу линии гра |
|||||||||||||||||
диента. Точно такие же соображения |
приводят к тому, |
|
что Тіз=0 при |
х3—0, |
||||||||||||||
т. е. в точке максимума |
скорости. Тогда для касательных |
напряжений |
Тіз в на |
|||||||||||||||
правлении |
линий |
|
д.'і по цилиндрическим |
поверхностям, |
для которых образующие |
|||||||||||||
есть изотахи, получим из (16.4) |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
*із (*и. *з) = |
и |
|
J 9 > |
Х з ) |
J h |
<*2. 6) Ц (*а , |
£) d€. |
|
(16.6) |
||||||||
Здесь интегрирование выполняется вдоль координатной линии х3. Для каса |
||||||||||||||||||
тельного напряжения на смоченной |
поверхности |
(х3=х%) |
эта формула |
имеет вид |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
Х* |
|
|
|
|
|
|
|
|
|
"із ( * 2 . •**) = |
^ |
,J2 |
х |
^ |
j |
Ц {х2, |
х3) Ц |
(х2, |
|
xs) dx5. |
|
(16.7) |
Выделим двумя бесконечно близкими линиями градиента элемент площади сечения потока (рис. 18). В данном случае — это некоторый криволинейный треугольник ОАВ. Площадь элемента равна
-г*
А/ р = j L 2 (х2, х3) L 3 (х2, х3) |
dx2dx3, |
а твердая часть его периметра составляет |
-V*) dx2 • |
|
||||
|
|
А В = |
L 2 (х2, |
|
|
|
Введем длину |
|
|
|
|
|
|
|
AF |
1 |
-*•'* |
|
||
|
Г |
|
|
|||
1=1 |
{ ^ - Ш |
- = Ц { Х 2 |
, . Ч ) |
J |
М *2 . *зНз ( . *2 . |
X3)dxz |
и перепишем |
формулу |
(16.7) так: |
|
|
|
|
|
|
1,3 (х2, .г*) = |
f Л (А-2). |
(16.8) |
В случае равномерного потока в бесконечно широком прямоугольном русле глубиной Я и в случае потока в круглой трубе радиусом г длина I равна соот ветственно Н и г/2, т. е. в обоих случаях она совпадает с гидравлическим радиусом R. Средняя же скорость V* в элементе совпадает со средней скоростью U по сечению. В потоках с иной формой сечения длина I зависит от х2 и ее можно рассматривать как обобщенное понятие гидравлического радиуса. Это подсказывает следующую гипотезу, по аналогии с (16.1),
4 3 (*2 . •*«) |
( 1 6 < 9 ) |
ТС2 •
Здесь |
С — коэффициент |
Шези, |
рассчитанный |
по |
обобщенному |
гидравличе |
||||||||
скому радиусу /, т. е. зависящий от х«. Очевидно, что |
|
|
|
|
|
|||||||||
|
|
|
|
г* |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
] |
«і (*2. х3) |
L 2 (JC2 , л-3) I 3 (х2, |
л-3) rfx3 |
|
|
|
||||
|
c7* = £ / * ( * 2 ) = - |
f |
|
|
|
|
|
|
• |
|
(16.Ю) |
|||
|
|
|
|
|
r* |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
J-2 (-^2. |
л-3) I 3 (x2, |
x3) |
dx3 |
|
|
|
||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
Средняя скорость |
потока |
равна |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
-і"» |
|
|
|
|
|
|
и = |
-р- [ І'І dF = |
-p- ф [f* (л-,) rfx2 |
J /,, (x2, |
x3) L 3 |
(x2, x3) |
dx3 |
= |
|||||||
|
|
F |
|
|
|
|
О |
|
|
|
|
|
|
|
|
|
|
= l^-j)CY7ldx2L2{x2, |
|
|
|
|
|
|
(16.11) |
||||
В этих |
формулах |
ф —символ |
интегрирования |
по замкнутым координатным |
||||||||||
линиям х2. |
Так как L2(x2, |
xse)dx2=d%, |
где d% — элемент |
смоченного |
периметра |
|||||||||
потока, то из (16.11) |
получается следующая |
формула |
для пропускной |
способно |
||||||||||
сти: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
К=§Сіу1 |
|
|
|
|
|
|
|
(16.12) |
|
В этой формуле, в отличие от (16.2), пропускная способность зависит от |
||||||||||||||
распределения |
скоростей |
в потоке |
и, следовательно, |
от формы русла. |
|
|
||||||||
К формуле |
(16.12) мы вернемся в § 26 в связи с исследованием устойчиво |
|||||||||||||
сти призматического потока. Там же рассмотрим потоки в призматических |
руслах |
с различным коэффициентом шероховатости на различных частях смоченного пе риметра.
Попробуем обобщить высказанную гипотезу на непризматические потоки, имея в виду при этом не только касательные напряжения, а все компоненты тензора турбулентных напряжений (тензора пульсации), вызванные как шерохо ватостью дна и стенок потока (первая система), так и кривизной линий тока (вторая система). В призматическом потоке действует только первая система напряжений, связанная, согласно предположениям предыдущего параграфа, с ли
ниями тока |
вектора |
v |
и ортогональными |
к ним траекториями градиента нормы |
||v|| этого |
вектора. |
В |
общем случае |
непризматического потока траектории |
grad |Jv|| не ортогональны к линиям тока v. Это ясно из того, что в общем слу
чае не существует |
семейства поверхностей, ортогональных |
к линиям тока. |
Чтобы |
||||
убедиться |
в последнем, заметим, что если |
существует поверхность |
ф(Хі, х2, |
х3) = |
|||
=0, |
ортогональная |
к линиям тока v, то |
нормаль к этой |
поверхности должна |
|||
совпадать |
по направлению с вектором v, т. е. должно быть |
|
|
|
|||
|
|
|
v = Xgrada. |
|
(16.13) |
||
|
Отсюда |
|
|
|
|
|
|
|
|
rot v = |
rot (X grad <р) = X rot grad <p -4- [grad X, grad <?] |
= |
|
||
|
|
|
— [grad X, grad a] = |
[grad X, v). |
|
|
|
|
А так как векторное произведение перпендикулярно |
своим |
сомножителям, |
||||
то |
отсюда |
следует |
(v, r o t v ) = 0 . |
|
(16.14) |
||
|
|
|
|
||||
|
Это и есть условие существования поверхностей, ортогональных |
линиям |
тока. |
||||
|
Имея |
в виду общий случай, в котором условие (16.14) не выполняется, по |
пробуем найти обобщение понятия линий градиента, сохраняя требование орто гональности этих новых линий к линиям тока v.
Выделив определенную линию тока v, зафиксируем на ней некоторую точку 0. Примем эту точку 0 за начало координат некоторой декартовой прямоуголь ной системы О-ЇІ-^З и проведем через эту точку плоскость а, ортогональную линии тока. Уравнение линии тока есть
dx, |
dx, |
dx, |
( 1 6 Л 5 ) |
v, |
v2 |
v3 |
|
Уравнение же плоскости а, проходящей через точку 0 (рассматривается бесконечно малая окрестность начала координат), имеет вид
A, dx\ 4- А2 dx2 + A3 dx3 = 0. Эта плоскость ортогональна линии тока, если
|
Щ _ |
Щ _ |
Уз |
|
|
|
|
||
|
А, |
А2 |
|
А3 |
|
|
|
||
Таким образом, уравнение плоскости а есть |
|
|
|
||||||
|
v, dx, + v2 dx2 |
+ v3 dx3 — 0. |
|
|
|
||||
Рассмотрим теперь вектор |
grad |
||v||. Очевидно, что он совпадает |
по напра |
||||||
влению с вектором grad v2. Составляющие |
последнего по осям |
х,, х2, |
х3 |
обозна |
|||||
чим С О О Т В е Т С Т В е Н Н О Через gt, g2, ёз- |
|
|
|
векторы v и grad v2 |
|
|
|||
Уравнение плоскости р, в которой |
лежат |
и |
которая |
||||||
проходит через начало координат, есть |
|
|
|
|
|
|
|||
|
a, dx, + а2 dx2 |
+ я 3 |
dx3 = 0, |
|
|
|
|||
где |
|
|
|
|
|
|
|
|
|
|
а, = v2g3 |
— v3g2, |
а2 |
= |
v3g, — v,g3, |
|
|
|
|
|
|
«3 = |
^1^2 — |
v2g3. |
|
|
|
||
Плоскости а и fi пересекаются по прямой |
|
|
|
|
|||||
|
dx, |
|
dx2 |
|
|
dx3 |
|
|
(16.16) |
a3v2 — a2v3 |
a,v3— |
a3v, |
a2v,—a,v2 |
|
|
||||
|
|
|
|||||||
Если переходить |
от одной линии тока к другой и от одной точки 0 к другой, |
||||||||
то уравнение (16.16) |
будет определять кривые, ортогональные |
к линиям |
тока v. |