Научная статья на тему 'Модель расчета параметров потока на входе в расширение'

Модель расчета параметров потока на входе в расширение Текст научной статьи по специальности «Механика и машиностроение»

CC BY
70
23
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
открытый водный поток / параметры потока / бурный равномерный поток / инженерный метод / двухмерный открытый водный поток / плоскость годографа скорости / интеграл Бернулли / линия сопряжения / крайняя линия тока / эквипотенциаль / произвольная линия тока

Аннотация научной статьи по механике и машиностроению, автор научной работы — Коханенко Виктор Николаевич, Папченко Игорь Владимирович, Папченко Наталья Геннадиевна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — Коханенко Виктор Николаевич, Папченко Игорь Владимирович, Папченко Наталья Геннадиевна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

The subject of our research work is reaching of the analytical decision one of the actual problems in Hydraulics for open stream waters, i.e. the calculation problems of stream parameter in channel expansion. For the decision of the given problem, we used the engineering calculation method with the moving binary measured equations for open water streams on the plane of graph year speed.

Текст научной работы на тему «Модель расчета параметров потока на входе в расширение»

УДК 532.543

МОДЕЛЬ РАСЧЕТА ПАРАМЕТРОВ ПОТОКА НА ВХОДЕ В РАСШИРЕНИЕ © 2008 г. В.Н. Коханенко, И.В. Папченко, Н.Г. Папченко

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

The subject of our research work is reaching of the analytical decision one of the actual problems in Hydraulics for open stream waters, i.e. the calculation problems of stream parameter in channel expansion. For the decision of the given problem, we used the engineering calculation method with the moving binary measured equations for open water streams on the plane of graph year speed.

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

Целью настоящей работы является получение аналитического решения одной из актуальных в гидравлике задач для открытых потоков, т.е. задачи расчета параметров потока в русле с расширяющимися стенками (рисунок).

Пусть бурный равномерный поток (F0 > 1) с параметрами У0, h0 и шириной - «Ь» сопрягается с потоком в расширяющемся призматическом русле, имеющем горизонтальное дно и прямолинейные стенки, наклоненные к продольной оси симметрии потока под углом - 8 ст.

Необходимо рассчитать геометрические и кинематические параметры потока внутри и на границе области G.

Для решения задачи воспользуемся инженерным методом расчета с использованием системы уравнений движения двухмерных в плане открытых водных потоков в плоскости годографа скорости в виде [1]:

Зф Зг

h0

1-3г Зу

2H

о г(1 -г)2

Зф _ 2h 0 г Зу

39 = H0 1 -г Зг ;

39

(1)

где ф = ф (т, 8) - потенциальная функция; у = у (т, 8) -функция тока; 8 - угол, характеризующий наклон местной скорости частиц потока к положительному направлению оси симметрии; т - квадрат скоростного

коэффициента - X; г = X2, X = -

V

V - модуль

^/2gH;

скорости; g - ускорение силы тяжести; Н 0 - посто-

V2

янная в интеграле Бернулли, Н0 = —— + h0^, V0,h0 -

2 g

параметры трехмерного потока на входе в расширение.

к-я пиния тока

т-я эквипотенциаль

План течения потока с расширением прямолинейными стенками: V0 - величина скорости потока на входе в расширение; h0 - глубина потока на в ходе в расширение; F0 = V02 /gh0 - критерий Фруда

Система уравнений (1) согласно теории [1,2] имеет ряд регулярных решений. Из этого ряда для решения задачи в настоящей работе воспользуемся решением в виде:

гь т_т-| + Ci0. т 2h L 1

Ф = -

о

2H

ln-

1 _ X 1 _ X

(2)

+ C 20.

Именно решение (2) может удовлетворять модели потока в русле с расширяющимися прямолинейными стенками.

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

Считаем, что сопряжение пространственного и двухмерного в плане потоков происходит по дуге окружности АВС радиуса R. Определим параметры потока в точках А и В при условии, что угол расширения стенки 6 ст не должен превышать угла поворота потока на его выходе в расширение при свободном растекании 6 тах , т.е. должно выполняться условие [3]:

6 <6

= Ci +(V3 _ 1)2,

C1 = arctg

3X 0 1 _V3arctg

1 _ X n

3x о _ 1 _L 1_Xо V3

где x 0 =

V2

2 gH о

Из условия постоянства расхода потока на линии сопряжения потоков АВС, получим равенство:

Vohob = 2hBVBR0 ст

Используя геометрическое соотношение

R = ,

2sin 6 ст

равенство (5) запишем в виде:

6

УЖ« =

(3)

0"0 - 'lBr B ■ р.

Sin 0 с

(4)

Рассматривая равенство (4) как уравнение для определения параметров потока hB, Ув с учетом интеграла Бернулли:

^ + ^т=Н 0,

приходим к условию: hв = h0, Ув = У0,

справедливому при малых углах поворота стенки 6 т.е. при выполнении условий: 6ст < 6о;

0

: sin 0 с

В случае углов 6 ст < 6 для определения параметров потока hв, Ув необходимо решить систему уравнений:

6„„

hoVo = hBVB

sin 6,

H о = hB +

2g

(5)

и в качестве решения выбрать такие корни системы, чтобы Ув было больше У0 и ближайшее к нему.

Система уравнений (5) сводится к приведенному кубическому уравнению:

V3 -2ёИоУ + 2Я = 0.

6 ст

Определив, таким образом, параметры потока на линии сопряжения АВС, решим далее задачу определения параметров потока в произвольной точке D, пересечения Ъй линии тока и т-й эквипотенциали. Заметим, что крайняя линия тока должна удовлетворять условию отсечения от оси симметрии 50 % расхода.

Согласно условию (2) и условию [4] =

& ^

для крайней линии тока, омывающей боковую стенку, получим уравнение:

У кр = С10 ст +

С 2 H о 2h п

Vb

[ln X_x] = ^f.

2

Так как параметр т изменяется вдоль стенки, необходимо потребовать выполнения условия: С 2 = 0.

Следовательно, у кр = С16 ст = , где

С1 = УЬ.

1 26 ст

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

У k = кС16 ст = С16 k ,

где 6к = k6ст, 0 < k < 1 .

k = 1 соответствует крайней линии тока, проходящей вдоль стенки;

k = 0 соответствует оси симметрии потока. В таком случае произвольная точка D в области течения потока определяется ^й линией тока и т-й эквипотенциалью.

В плоскости годографа скорости по параметрам 6 k и т м фиксируется линия тока и эквипотенциаль:

2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

х

0

у k = ci0 k;

Ф M =

C1h0

2Hn

ln-

1 -Т.

1 -Т.

Для перехода в физическую плоскость течения потока используем комплексную дифференциальную связь между физической плоскостью и плоскостью годографа скорости [5]:

d ( x + iy) = ^ d ф + i-^0- d yj 1 e

(6)

Полагая в формуле (6) d у = 0, т.е., идя вдоль линии тока, получим систему дифференциальных уравнений вдоль произвольной линии тока:

dx = -

cos 8

V2gH0

dy = -

sin 8

V2gH0

-d ф;

-d ф,

(7)

где т = -

V2

2 ёН о

Для фиксации произвольной линии тока в системе уравнений (7) полагаем 6 = 6 к. В результате из системы (7) следует система дифференциальных уравнений:

dx = -

cos 6,

■у d ф;

dy = -

sin 6,

(8)

V2gH0

-yT dф.

Учитывая, что

dф = -

C1h0 2H

ln-

1 — Т 1 — Т

C1h0 1 — 3т 'щ0 т(1 — т) 2

d т ,

и, интегрируя систему (8), получим в квадратурах:

x n x Ь —

yD — Ук =

C1h0 Т M cos 6 k J Т В 3т —1

2H 0 Т 32 (1 — Т)

C1h0 Т M sin6 k J Т В 3т —1

2H 0 Т 32 (1 — т)

-d т;

(9)

-d т,

где хм - значение параметра в точке D на пересечении к-й линии тока и ьй эквипотенциали; 6 к - значение параметра 6 в точке D; хЕ,, уЕ1 - координаты точки D; хк, ук - координаты точки К на начальной эквипотенциали.

Из геометрических соотношений можно определить координаты хк, ук , если известен радиус эквипотенциали R и угол линии тока 6к по формулам:

(xk = R cos 8 k - R cos 8 ст = R (cos 8 к - cos 8 ст j;

Ut = R sin 8 k.

Беря интегралы в системе (9) и пользуясь формулой Ньютона-Лейбница, окончательно получим выражения для определения координат точки D:

x п = хь +

C1h0 2H

cos 8

C\h0 ■ p. yD = yk +T1rsm6к

2Hn

Значения параметра х м определяем из формулы h

хм = 1 —м, где hм - глубина потока в точке М на

Н о

оси симметрии потока. Параметры 6 к и hм задаем произвольно в долях угла 6 ст и глубины hB.

Из системы (10) в частном случае 6 к = 0, т.е. на оси симметрии потока следует:

у D = 0;

C1h с

2H

2

2

LM

(1 —ТM ) ТВ2 (1 — ТB ).

(11)

Уравнение (11) представляет собой закон распределения глубин и скоростей потока вдоль его оси симметрии, так как, задаваясь х м , определяются

глубина hм и скорость Vм по формулам:

Vм =х ;

hм = Н0 (!"хм)•

Из системы уравнений (10) следует, что:

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

параметр х ^ 1, h ^ 0, V ^ Vmax = ^2ёН0 ;

- при малых углах расширения 6 ст < 6 0 параметры потока на линии сопряжения потоков АВС равны V0, h0, т.е. поток деформируется без изменения скоростей и глубин; таким образом, образуется мостооб-разный участок АВСО, ограниченный прямой линией и дугой окружности;

- степень падения глубин вдоль оси симметрии потока зависит от параметра х 0, равно как и от числа

V 2

Фруда Fe = -О- = gh с

2т О 1 —ТГ

- при х м ^ 1, хм, Ум . Вдоль граничной линии тока, как для частного случая, из системы (10) следует:

2

т

х

х

x

Т

х

2

х

d

2

2

x г, _ x и +

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

C1h о 2Hn

cos 8,

(!"тм) тв (1 "тв).

yD = Ук +

C1h0 2Hn

sin 8,

(! "тM) тв2 (1 "тв).

ГДе Хк = 0 Ук = R Sin 8 ст = "2.

Результаты экспериментов, проведенных на экспериментальной установке в НГМА, показывают достаточную для практики расчета гидротехнических сооружений адекватность модели реальному течению потока.

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

Метод, используемый в работе, позволяет вместо применяемых численных методов и полученных ре-

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

В последующих работах покажем, как трансформируется решение задачи в случае учета сил трения потоку и уклона дна русла с расширением.

Литература

1. Коханенко В.Н., Волосухин Я.В., Ширяев В.В., Коханенко Н.В. Моделирование одномерных и двухмерных открытых водных потоков / Под общ. ред. д-ра техн. наук, профессора В.Н. Коханенко. Ростов н/Д., 2007.

2. Ширяев В.В., Мицик М.Ф., Дуванская Е.В. Развитие теории двухмерных открытых водных потоков / Под общей ред. В.В. Ширяева. Шахты, 2007.

3. Емцев Б.Т. Двухмерные бурные потоки. М., 1967.

4. Высоцкий Л.И. Управление бурными потоками на водо-

сбросах. М., 1977.

6 июня 2008 г.

Коханенко Виктор Николаевич - д-р техн. наук, профессор, зав. кафедрой Донского государственного аграрного университета, пос. Персияновский. Тел. 3-55-21.

Папченко Игорь Владимирович

пос. Персияновский

ассистент Донского государственного аграрного университета,

Папченко Наталья Геннадиевна - ассистент Донского государственного аграрного университета, пос. Персияновский

2

2

т

2

2

т

i Надоели баннеры? Вы всегда можете отключить рекламу.