CC BY
12
УДК 502/504:532.543
Е. В. ДУВАНСКАЯ, М. Ф. МИЦИК
Государственное образовательное учреждение высшего профессионального образования «Южно-Российский государственный университет экономики и сервиса»
В. Н. КОХАНЕНКО
Федеральное государственное образовательное учреждение высшего профессионального образования «Донской государственный аграрный университет»
РЕШЕНИЕ ЗАДАЧИ СВОБОДНОГО РАСТЕКАНИЯ В ПЛОСКОСТИ ГОДОГРАФА СКОРОСТИ
Предлагается новый подход к описанию крайней линии тока в плоскости годографа скорости для бурного свободнорастекающегося стационарного потока воды. Предлагаемое решение удовлетворяет условиям непрерывности параметров течения на выходе из трубы в нижний бъеф. Показана адекватность нового решения предыдущему почти на всем плане течения.
Крайняя линия тока, прямоугольная труба, горизонтальное отводящее русло.
The article describes a new approach to the description of the extreme current line in the plane of the velocity hodograph for torrential freely spreading stationary water flow. The proposed solution meets the requirements of continuous flow parameters at the pipe outlet downstream. There is shown an adequacy of the new solution to the previous one on the almost whole flow plan.
Extreme current line, rectangular pipe, horizontal discharge channel.
В работах [1, 2] было получено уравнение математической физики в частных производных второго порядка для функции тока бурного свободнорастекающегося потока:
у = —= const при 9 = 0; т = т0,
(3)
д [ 2т ду! 1 - 3т
дт [1-х дт J
2
д у
2т (1 -т)2 де2
= о,
(1)
где Ь - ширина водопропускной трубы; у0 -модуль вектора скорости на выходе из трубы в широкое отводящее русло.
В связи с этим на выходе потока из трубы нарушается непрерывность течения по параметрам т, 9.
где у = у (т; 9) - функция тока; т - аргумент,
изменяющийся в пределах — <Х0<Х< 1, являющийся квадратом модуля вектора скорости жидкой частицы; 9 - угол наклона вектора скорости жидкой частицы к продольной оси симметрии потока (рис. 1).
В вышеуказанных работах было найдено приближенное аналитическое решение уравнения (1) для верхней крайней линии тока (которую в дальнейшем будем называть просто крайней линией тока)
yi = А
sin е i ■
т 2
(2)
где А - постоянная, определяемая из граничных условий задачи.
Однако решение (2) не удовлетворяет уравнению крайней линии тока:
Рис. 1. План свободного растекания бурного потока
Целью проведенных авторами исследований был поиск приближенных решений уравнения (1) в окрестности точки К(9 = 0; т = т0), удовлетворяющих условию непрерывности по параметрам (т, 9) для крайней линии тока в плане растекания потока.
Уравнение (1) будем решать так:
с \
у 2 = sin 0
A+в 0 1 1
(х-хс )2
(4)
где а > 0.
Кривую вида (4) будем рассматривать, учитывая, что для крайней линии тока (4) выполняется условие sin 0
Г = 1-
lim-
(с-^с )2
(5)
Из (5) следует, что в окрестности точки т = т0 имеет место равенство
sin 0 = (т-То )2 + o | (т - То )1 |,
sin 0-5(1 -х)а ( Ла + lim-v , 7 = В (1 — х с) .
(7)
В = ■
Ve b
2(1-ТсУ ' (8)
Рассмотрим теперь поведение функции (4) при т ^ 1. На крайней линии тока получим:
с Л
lim у 2 = limsin 0
х—1 х—1
A+В (1 -х)
1 + 1
(х - Хс )2
A sin 0m 1
(9)
-+ С = A sin 0„
где 9 - угол между вектором скорости и осью симметрии потока Ох на бесконечности.
Учитывая условие (3), из (9) имеем:
A = ■
Vo b
2вШ 0^ (10)
Такому же значению постоянная А равна и в выражении (2) для функции
С учетом полученных значений постоянных А и В уравнение для крайней
линии тока будет иметь следующий вид:
( \
1 . (1-ОГ
у 2 =
= Veb sin 0
2
. (11)
л2 sin 0_ (1 -т0 Т(т-т0)2
Исследуем свойства кривой при изменении аргумента т. На крайней
линии тока из (11) получим:
( \
veb = \Ь 22
отсюда
-sin I
1
(1 -х)а
Vх 2sin 0max
(1 -хс )а(х-хс )
sin 0 = ■
Х1 (1 - Х С )а(Х-ХС )2sin 0m
где о - бесконечно малая функция относительно своего аргумента,
отсюда
lim sin 0 = Си lim 0 = С. (6)
Х—Хс Х—Хс ^ '
Найдем значения постоянных A и B. Вычислим предел функции (4) при
т ^ V
Asin 0 lim у = lim-1--+
Х— Хп Х—Хп -
С с т-2
_' (12)
(1 -Х0 )"(Х-Х0 )2 +(1 -х)"х28т бшах
Выясним поведение угла 9 в зависимости от т. Для этого сначала определим поведение функции 1/в1п 9 по т. Из (12) получим:
1
1
(1 -Х)а
sin 0
Х 2 sin
0 max (1 -Х0 НХ-Х0 )
1 '
2
(13)
Найдем производную от выражения (13) по т:
(х-хо )2
Поскольку на крайней линии тока выполняется условие (3), то из (7), что
d
1
sin 0
1
2х 2 sin 0
max
(1 -Х)а 1 (а(Х -хс )+1 -х )
(14)
(1 -хс)а (Х-Х0 )
Поскольку выполняются соотношения
1 <хс< х< 1,
2 с
(15)
то каждое из слагаемых правой части равенства (14) при условиях (15)
отрицательно, следовательно, функция 1/sin 9 убывает, соответственно sin 9 -монотонно возрастающая функция. Таким образом, на крайней линии тока возрастанию аргумента т соответствует монотонное возрастание угла 9.
Сравним взаимное расположение кривых у1 и у2 в плоскости годографа скорости. Для этого кривую (2) представим так:
- f (х)-^
sin 9
(16)
Х' sin 9max
Аналогично кривую (11) представим в виде
1 /• / \ 1
- f2 С0 = --+
sin 9
х2 sin 9
(1 -х)а
(17)
(1 -х0 )а(х-х0 )
Исследуем теперь, в каких границах кривые /1(т) и /2(т) мало отличаются между собой. Такая ситуация будет иметь место в том случае, когда мало будут отличаться углы при одинаковых значениях т. Нетрудно видеть, что
f2 f1 f1
< f2 - f = (1 "*)"
> 0.
(18)
(1 -х0 ) (х-х0 )2
Определим интервал по т, на котором выполняется оценка
(1 -х)"
_Ь_¿_ < £
(1 -Х0)"(Х-Х0>1 ' (19)
где е - погрешность аппроксимации, которая в практических расчетах обычно принимается в пределах от 3 до 5 % .
Пусть сначала переменная т удовлетворяет условию
1
1 >х>1 — т (1 — х0)ва. (20)
Тогда выполнятся следущие оценки:
1 -х< m (1 -х0 )e";
(х-Х0 )2 >\¡1-Х0 - m (1-Х0 )е" -
1 I-Т (21)
-(1 -х 0 )2 V1 - me" .
Оценим теперь разность f1 - f :
' (1 -Х0)'
f2 - f <
m e i
(1 -X0T(1 -X0)2 V1 -me e
- m
(1 -X0)
1
- me
m
Потребуем выполнения оценки:
e
)
< e.
(1 — х0) 2 "\/1 - твс
При этом должно быть выполнено следующее условие:
( - Л
т2а < (1 — х0) 1 — тва . (22)
V )
Условие (22) будет выполнено, если т 2а<(1 — х0 )(1 — т ). (23)
Решая неравенство (23) относительно а, получим: 1п (1 — х0 )+1п (1 — т )
21п т
При проектировании дорожных водопропускных сооружений число Фру-да на выходе из трубы Кг0 используется в диапазоне от 1,5 до 5, что соответствует изменению параметра т0 в пределах от 3/7 до 5/7. Поскольку
а >
(24)
3 5 2 4
-<х0<-, то -< 1 - х0 <-.
7 0 7 7 0 7
(25)
Учитывая неравенства (25) и полагая, что т = 6/7, из (24) получим:
1п (2 ) — 21п (7 )
а > —^-= 10,3751,
2 (1п 6 — 1п7)
тогда
а = 1,1. (26)
При неограниченном увеличении значения а функция /2(т) будет как угодно близка к функции /1(т), соответственно функция у2 будет близка к функции у1.
Таким образом, функция выражением (11) определяется теперь однозначно. Выясним, в каком диапазоне изменения переменной х в физической плоскости находится крайняя линия тока, если
1
х0 <х< 1 -т(1 -х0)ва' (27)
На отрезке изменения переменной
х е
х0; 1 - т (1 )е£
кривые и
существенно отличаются, соответственно они отличаются и в плоскости годографа скорости на отрезке [0; хм], который соответствует условиям (27). Для определения границ изменения аргумента х на крайней линии тока, соответствующих условиям (27), воспользуемся формулой из [1], определяющей абсциссу точки М по известному значению т :
м
|Х„
Лкп
2H
где J1 -
Ah0 /
i (J, - kJ2 - kJ3
42gHoy 1 2 3
1 + х
х(1 -х)
х
■ + ln-
1 -х
; к = sin2 0ma;
J -
1 -х
■ + ln
(28)
1 -х
33 - 1п
1 -х
где g - ускорение силы тяжести.
Следует отметить, что формула (28) справедлива для кривой (2), а диапазон изменения аргумента т у функции ^ отличается:
х <х<х (29)
где значение тк > т0 и определяется как решение системы [1]:
Sin 0K = ХК^п 0max '
cos 0
K
(30)
хК (1 -хк) хo2 (1 -х0)
при этом значение 0
^ m
формуле [3]:
0max = (V3 - 1)^2 + arCt^^3!
P (1 -х0 )
вычисляется по
3х0-1
- хп
Значение тм определяем из (27):
хм -1 -т(1 -х0)sa. (32)
Полагая, что в (31) s = 0,03, вычисляем:
'м
-1" 14 (0).
(33)
Упростим выражение коэффициента в формуле (28), используя известные [1] выражения для глубины к0 и скорости у0:
К = н о (1 -хо);
1
v0 -х 0
Следовательно,
Ah
v,
bbH0 (1-TQ )
2 H^^lgH^ 2 sin 0mxJlgHo2Ho
1
- b (1 -х0 )х2 4sin 0 .
(34)
Таким образом, абсцисса точки M представляет собой выражение вида
хм = ^СО-
Выражение для функции F(t0) получается, если вычислить тк из системы
(30), а величины 0max и тм - по формулам
(31), (33). Такие преобразования были реализованы в среде прикладных математических программ Maple 9.5. График функции F(t0) представлен на рис. 2.
Отрицательные значения абсциссы xM при т0 < 0,58 объясняются следующим образом: формула (28) для нахождения границ по x определялась для кривой у1, задаваемой уравнением (2), значения переменной т удовлетворяют условию (29), т.е. даже при
b
0,7 0,5 0,3 0,1
Х- 0,1 b
0\---0,5 0,6 0,7 0,8 0,9 1 х0
Рис. 2. График функции F(t0)
хм -
к
т0 = 1/3 значение тк = 0,58.
При проектировании дорожных водопропускных сооружений значения параметра т0 используются в пределах от 3/7 до 5/7, или т0 е [0,42; 0,72]. (35) Из рис. 2 нетрудно видеть, при таких значениях параметра т0 выполняется условие х/Ь < 0,1.
Таким образом, кривые у1 и у2 при условиях (35) практически совпадают в плане растекания потока, за исключением малого участка потока вблизи выходной трубы, ширина которого не превосходит 0,1Ь (рис. 3).
для крайней линии тока А =
(/ - /)
/1
100%
не превосходит 3 %.
Решение у2 удовлетворяет условию непрерывности течения на выходе из трубы (при т = т0, 9 = 0). Угол 9 на кривой у2 монотонно возрастает с ростом т.
Уравнение крайней линии тока удовлетворяет граничным условиям свободного растекания бурного потока как на выходе из трубы, так и на бесконечности.
Рис. 3. Сравнение крайних линий тока уг и у2 в плане течения
При значениях т0 от 3/7 до 5/7 (т.е. в тех диапазонах значений, которые используются при проектировании дорожных водопропускных сооружений) и значениях х > 0,1Ь относительное рассогласование выражений
1. Мицик М. Ф., Косиченко Н. В., Лемешко М. А. Метод с использованием годографа скорости применительно к расчету параметров бурного двухмерного потока // Математическое и компьютерное моделирование естественнонаучных и социальных проблем: сб. статей IV Междунар. научн.-техн. конф. молодых специалистов, аспирантов и студентов. - Пенза: Приволжский Дом знаний, 2010. - С. 130-141.
2. Моделирование одномерных и двухмерных открытых водных потоков: монография / В. Н. Коханенко [и др.]. - Ростов на Доду: Изд-во ЮФУ, 2007. - 168 с.
3. Емцев Б. Т. Двухмерные бурные потоки. - М.: Энергия, 1967. - 212 с.
Материал поступил в редакцию 17.09.10 Дуванская Елена Викторовна, кандидат технических наук, доцент E-mail: delvik2004@list.ru
Мицик Михаил Федорович, кандидат технических наук, доцент E-mail: m_mits@mail.ru
Коханенко Виктор Николаевич, доктор технических наук, профессор, зав. кафедрой «Техническая механика и гидравлика» E-mail: kosmi74@yandex.ru
