Научная статья на тему 'Water Flow Parameters on the Symmetry Axis and Extreme Line of Current'

Water Flow Parameters on the Symmetry Axis and Extreme Line of Current Текст научной статьи по специальности «Физика»

CC BY
13
4
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Вестник МГСУ
ВАК
RSCI
Область наук
Ключевые слова
mathematical model / hydrodynamics analysis / two-dimensional water flow / open-channel hydraulics / analytical solution / математическая модель / гидродинамический анализ / двумерный поток воды / гидравлика открытого канала / аналитическое решение

Аннотация научной статьи по физике, автор научной работы — Ольга Александровна Бурцева, Мария Сергеевна Александрова

Introduction. An analysis of mathematical models of two-dimensional planned flows was carried out. Such flows are characterized by local depth-averaged velocities and local depths at each point of the flow. The mathematical model formation of the water flow is based on its division into several sections. There is a section where the flow parameters (velocity, depth, width) are kept constant at the stage of flow exit from the pipe — the inertial front. The purpose of the article and its relevance are defined. Materials and methods. By introducing dimensionless complexes on the basis of π-theorem, the formula for the length of inertial front of the water flow at its spreading from a rectangular pipe into a wide diverting channel is derived. An analogy from gas dynamics is used, namely, the transition to the plane of the velocity hodograph. Using the velocity hodograph, the distribution of depths and velocities of the flow along its longitudinal axis of symmetry and along the extreme line of current was obtained. The main computation tasks for the flow parameters have been formulated. Results. Numerical calculations of the formulated main tasks for determining flow parameters are described. Comparison with experimental data is given and the adequacy of the refined mathematical model of a two-dimensional planned flow is confirmed. Conclusions. The resulting formula for the length of the inertial front makes it possible to achieve the desired error in calculating the parameters of the water flow. With flow expansions up to 5, the relative error of the ordinates and flow velocities does not exceed 7–10 %. Calculation formulas and implemented programs will allow HTS designers to quickly and accurately determine the boundaries, speed and depth of free flow on the culvert.

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

Параметры водного потока на оси симметрии и крайней линии тока

Введение. Проведен анализ математических моделей двумерных плановых потоков. Такие потоки характеризуются местными осредненными скоростями по глубине и местными глубинами в каждой точке потока. Формирование математической модели водного потока основано на его разделении на несколько участков. На этапе выхода потока из трубы имеется участок, где сохраняются постоянными параметры потока (скорость, глубина, ширина) — инерционный фронт. Определена цель статьи и ее актуальность. Материалы и методы. Введением безразмерных комплексов и на основе π-теоремы выведена формула длины инерционного фронта водного потока при его растекании из прямоугольной трубы в широкое отводящее русло. Использована аналогия из газовой динамики, а именно переход в плоскость годографа скорости потока. С использованием годографа скорости получены законы распределения глубин и скоростей потока вдоль его продольной оси симметрии и вдоль крайней линии тока. Сформулированы основные задачи расчета параметров потока. Результаты. Описаны числовые расчеты сформулированных основных задач определения параметров потока. Приводится сравнение с экспериментальными данными и подтверждается адекватность уточненной математической модели двумерного планового потока. Выводы. Полученная формула длины инерционного фронта позволяет добиться желаемой погрешности расчета параметров водного потока. При расширениях потока до 5 относительная погрешность ординат и скоростей потока не превышает 7–10 %. Расчетные формулы и реализованные программы позволят проектировщикам ГТС быстрее и точнее определить границы, скорость и глубину безнапорного потока над водопропускной трубой.

Текст научной работы на тему «Water Flow Parameters on the Symmetry Axis and Extreme Line of Current»

ГИДРАВЛИКА. ГЕОТЕХНИКА. ГИДРОТЕХНИЧЕСКОЕ СТРОИТЕЛЬСТВО

RESEARCH PAPER / НАУЧНАЯ СТАТЬЯ UDC 626.31

DOI: 10.22227/1997-0935.2023.8.1262-1271

Water Flow Parameters on the Symmetry Axis and Extreme Line of Current

Olga A. Burtseva, Maria S. Aleksandrova

Platov South-Russian State Polytechnic University (NPI); Novocherkassk, Russian Federation ABSTRACT

Introduction. An analysis of mathematical models of two-dimensional planned flows was carried out. Such flows are characterized by local depth-averaged velocities and local depths at each point of the flow. The mathematical model formation of the water flow is based on its division into several sections. There is a section where the flow parameters (velocity, depth, width) are kept constant at the stage of flow exit from the pipe — the inertial front. The purpose of the article and its relevance are defined.

Materials and methods. By introducing dimensionless complexes on the basis of n-theorem, the formula for the length of inertial front of the water flow at its spreading from a rectangular pipe into a wide diverting channel is derived. An analogy from gas dynamics is used, namely, the transition to the plane of the velocity hodograph. Using the velocity hodograph, the distribution of depths and velocities of the flow along its longitudinal axis of symmetry and along the extreme line of current was obtained. The main computation tasks for the flow parameters have been formulated.

Results. Numerical calculations of the formulated main tasks for determining flow parameters are described. Comparison with experimental data is given and the adequacy of the refined mathematical model of a two-dimensional planned flow is confirmed. Conclusions. The resulting formula for the length of the inertial front makes it possible to achieve the desired error in calculating the parameters of the water flow. With flow expansions up to 5, the relative error of the ordinates and flow velocities does not exceed 7-10 %. Calculation formulas and implemented programs will allow HTS designers to quickly and accu-^ g rately determine the boundaries, speed and depth of free flow on the culvert.

U 3

> 10 KEYWORDS: mathematical model, hydrodynamics analysis, two-dimensional water flow, open-channel hydraulics,

W о

с и to «о

analytical solution

Acknowledgements. The authors are grateful to Doctor of Technical Sciences, Professor of the South Russian State Polytechnic University (NPI) named after M.I. Platov Viktor Nikolaevich Kokhanenko for a number of remarks, which were subsequently eliminated. The authors will be grateful to the reviewers for comments and advice that will help improve the paper,

O rr thereby raising its level.

>

FOR CITATION: Burtseva O.A., Aleksandrova M.S. Water Flow Parameters on the Symmetry Axis and Extreme Line of Crrent. Vestnik MGSU [Monthly Journal on Construction and Architecture]. 2023; 18(8):1262-1271. DOI: 10.22227/1997£ o 0935.2023.8.1262-1271 (rus.).

O <u v '

<u <u

§ о Corresponding author: Olga A. Burtseva, [email protected]

<o > S

S = ™ §

от " от Е

.Е о

Параметры водного потока на оси симметрии и крайней линии тока

£ с Ольга Александровна Бурцева, Мария Сергеевна Александрова

й ° Южно-Российский государственный политехнический университет (НПИ) имени М.И. Платова

о Е (ЮРГПУ (НПИ) имени М.И. Платова); г. Новочеркасск, Россия

§ ° -

^ ^ АННОТАЦИЯ

ОТ "£= Введение. Проведен анализ математических моделей двумерных плановых потоков. Такие потоки характеризуются

— 2 местными осредненными скоростями по глубине и местными глубинами в каждой точке потока. Формирование

^ • математической модели водного потока основано на его разделении на несколько участков. На этапе выхода

О ^ потока из трубы имеется участок, где сохраняются постоянными параметры потока (скорость, глубина, ширина) —

О инерционный фронт. Определена цель статьи и ее актуальность.

^ Ш Материалы и методы. Введением безразмерных комплексов и на основе п-теоремы выведена формула длины инер-

X ционного фронта водного потока при его растекании из прямоугольной трубы в широкое отводящее русло. Использована

X

с

_ аналогия из газовой динамики, а именно переход в плоскость годографа скорости потока. С использованием годографа

¡^ ¡^ скорости получены законы распределения глубин и скоростей потока вдоль его продольной оси симметрии и вдоль

Ш > крайней линии тока. Сформулированы основные задачи расчета параметров потока.

Результаты. Описаны числовые расчеты сформулированных основных задач определения параметров потока.

© О.А. Бурцева, М.С. Александрова, 2023 Распространяется на основании Creative Commons Attribution Non-Commercial (CC BY-NC)

Приводится сравнение с экспериментальными данными и подтверждается адекватность уточненной математической модели двумерного планового потока.

Выводы. Полученная формула длины инерционного фронта позволяет добиться желаемой погрешности расчета параметров водного потока. При расширениях потока до 5 относительная погрешность ординат и скоростей потока не превышает 7-10 %. Расчетные формулы и реализованные программы позволят проектировщикам ГТС быстрее и точнее определить границы, скорость и глубину безнапорного потока над водопропускной трубой.

КЛЮЧЕВЫЕ СЛОВА: математическая модель, гидродинамический анализ, двумерный поток воды, гидравлика открытого канала, аналитическое решение

Благодарности. Авторы выражают благодарность д.т.н., профессору Южно-Российского государственного политехнического университета (НПИ) имени М.И. Платова Виктору Николаевичу Коханенко за ряд замечаний, которые в дальнейшем были устранены. Авторы будут признательны рецензентам за замечания и советы, которые помогут улучшить статью, тем самым повысить ее уровень.

ДЛЯ ЦИТИРОВАНИЯ: Бурцева О.А., Александрова М.С. Water flow parameters on the symmetry axis and extreme current line // Вестник МГСУ. 2023. Т. 18. Вып. 8. С. 1262-1271. DOI: 10.22227/1997-0935.2023.8.1262-1271

Автор, ответственный за переписку: Ольга Александровна Бурцева, [email protected].

INTRODUCTION

The subject of research is open stationary water flows behind an unpressurised opening of rectangular cross-section, with a small vertical velocity component, namely, two-dimensional planar flows. The theory and methods for solving problems of planar hydraulics are most fully described in monographs by G.I. Sukhomel, I.I. Levi and I.A. Sherenkov [1], B.T. Emtsev [2] and VN. Kokhanenko [3].

Often, when designing hydraulic structures, it is necessary to know the characteristics of open water flows, for which the Froude criterion is greater than one. These streams are: with free spreading of the flow from open channels or free-flow pipes into a wide channel; during the flow of flows in the connections of channels of different widths; in open spillways behind hydroelectric power plants; in case of river floods behind culvert road structures; behind small bridges. According to the nature of the flow of a turbulent water flow through a non-pressure pipe, non-pressure and semi-pressure flow regimes are considered.

Open water flow in a suddenly opening channel can have velocities substantially exceeding the permissible velocities for the unstrengthened part of the channel. As shown by field surveys conducted by the authors [1, 4], the main cause of destruction of hydraulic structures is dangerous erosion of the downstream. Therefore, to calculate the anchorage of the outlet channel, hydraulic engineering structures designers need information on the flow parameters in the vicinity of the outlet edge of the pipe. Therefore, the relevance of this work is confirmed.

A water flow model behind an unconfined outlet has several sections [3-8]. In works [3-21] the parameters and shape of free-diffusion region of turbulent flow at its outflow from unpressurized pipe into wide diversionary channel is determined experimentally and analytically. In [22, 23] a formula for inertial front length of flow at the tube outlet is given. However, there is no sufficient justification for this formula.

The purpose of this paper is to derive the formula for the inertial front length on the basis of dimension

theory and n-theorem, as well as testing its appropriateness in calculating the parameters of water flows.

MATERIALS AND METHODS

There are basic units of measurement in any system of units. They are introduced from experience with the help of standards. In the SI, for example, the basic units are the meter, the second, and the kilogramme. The expression of an arbitrary unit of measurement in the basic units is called the dimensionality. For each basic unit, a notation is introduced: L is length, T is time, M is mass, etc. The theory of dimensionality proves [24, 25], that the dimension of any quantity is a degree monomial of the form = Ll • T • Mm ... and is called the dimension formula. The statement sought follows from the fact that the ratio of two numerical values of a physical quantity should not depend on the choice of scales for the basic units of measurement.

n-theorem. A relation, independent of the choice of units, between n dimensional quantities, k of which have independent dimensions, can be represented as a relation between (n-k) quantities which are dimen-sionless combinations of n dimensional quantities.

Formula for the length of the inertial front as itflows from a rectangular pipe into a wide diversion channel. Let us consider the problem of deriving a formula for the length of the inertial front X'D. The need to derive the formula is justified by the fact that along the inertial front all flow parameters retain their values. Hence, the formulas [3, 5] for calculating the flow parameters will change.

The process of free flowing of supercritical uncon-fined potential flow into a wide diversionary channel is determined by four parameters with:

• culvert width b, cm;

• depth of flow in the pipe h0, cm;

• initial velocity of the flow V0, cm/s;

• free fall acceleration g = 981, cm/s2.

We introduce dimensionless complexes characterizing the spreading process:

< n

8 8 TIS

kK

G Г

S 2

0 ся

t CO

1 ф

У 1

J to

u-

^ I

Ф 3 О

oi

о 3

b H

V2

F = — 0 gK

CO CO

3 ю

Ф 6 > §6

ф )

ii

® 8

. DO

■ T

s □

s У с о <D * »00

О О 10 10 U W

where

V 2

H o = -L + h 2g

(1)

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

the constant in D. Bernoulli's integral, 8 (radian) is an angle characterizing the slope of the velocity vector to the longitudinal axis of flow symmetry; F0 = 2t0 / (1 - t0 ). Froude's criterion and the square of the velocity coefficient t0 = V02/lgH0, 0max is the maximum flow angle.

Thus, the number of independent units is k = 2, the dimensional quantities n = 4. The number of independent dimensionless complexes is n - k = 2, and there is a single equation linking these complexes.

Consider the fall time of the water and the time of its inertial motion. Let h0 = 10 cm, V0 = 100 cm/s, XT = 4 cm, then the free fall time of a liquid particle from a height of h0 will take the form:

*f =

J± =0.145 s, 1,000 \ 50

and the time ofinertia travel tt = XD/V0 = 4/100 = 0.04 s. The inertia time of a water particle is 4 times shorter than its fall time. Thus, the water, before it has time to fly away, comes out of the pipe as a jet. Therefore, the greater the depth of the flow, the greater the length

To derive the structure of the formula X'D, we use of the inertial front, i.e. XT « h0/b , and X'T « hjH0

formula (1) in the form:

From the triangle ABC can be written (see Figure):

W (O

N N

o o

tv N

00 00

K (V

U 3

> in

E M 2

to 00

. r

« Q

?!

<D <u

O g

o

o g<

o

CO

™ O

o

ro

o

E o

CL ° c

Ln O

s H

o E

fe ° CD ^

Ü (0 №

H„ = h0 r v2 „ H

_2 gh„

^ h„ 2

H „ F„ + 2

H=F+1 ^

Figure shows that:

tga0 =

b/2

x^ '

where we get

X =

b/2

tga0

where a, is the wave angle [3] defined by the expression:

. 1

4F'

X

D A____

-_-

b/2

B X*^ Vo

x'D « x =

2 tga0

then [3]:

Assume that X'D is proportional to X*. From

X =

tga0

2

Let the initial abscissa of the flow X* be directly proportional to the length of the inertial front X'D, and also directly proportional to the flow depth h0 and inversely proportional to sin0 :

J A A max

XD X

X'D x-0, X'D x

X'd D H

T

Sinö„

we get the structure of the formula:

bjF^î i

—---x

2 Sin6max

i

XI

xh A_ № h

b H

„ 2 SinBmax ( + 2y

Whence can be written:

X'd JFTl 1

h

<x -

F0+ 2 sin0n

_iz To determine the initial abscissa of the flow

As a result of regression analysis applied to the experimental data, the authors obtained a formula for the length Next, we will prove that X'D <x h0/ b and X'D « X *1 of the inertial front of the flow when it spreads from a rectangular tube into a wide diversionary channel in the form of:

It's obvious that:

X'D x

SinG „

XD = trunc

VF^i i

F0 + 2 sin0n

+ 1,

1 The sign k denotes the proportionality of the quantities. 1264

1

1

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

which is verified by experiment. Here "1" is taken in length units: for the experiment, it is a centimeter, for nature calculations it is a meter.

Determination the distribution law of the depth and velocity along its longitudinal symmetry axis. The dynamic equations of motion of a two-dimensional open flow are obtained from the equations of L. Euler, supplementing them with terms that take into account the resistance forces. Complications also arise when taking into account the slope of the bottom of the conduit, i.e. the need to take into account the influence of gravity. As a result, the system of equations of motion and continuity of the flow form a system of essentially non-linear partial differential equations, the analytical solution of which has not been found so far.

The study of the flow equations begins with the simplest case.

Where in:

• resistance forces are neglected (it is true for supercritical flows passing through a narrowing, or an expansion, or a bend in the channel, since in these cases, inertial forces have a predominant effect on the flow);

• movement is considered potential;

• the bottom of the flow is assumed to be horizontal.

In this case, the equations of motion of a two-dimensional supercritical flow turn out to be similar to the equations of plane gas motion at supersonic speeds. Therefore, many known results of gas dynamics can be applied to supercritical flows [2]. Let's use the known method of S.A. Chaplygin to move to the velocity hodo-graph plane, that is, we pass from the coordinates in the plane to the variables of the velocity vector and its angle of inclination to the axis of symmetry of the flow. In this case, the nonlinear system of equations is reduced to a linear system of equations in partial derivatives, which greatly facilitates its study.

Using the velocity hodograph plane [2, 3], a differential complex relation between the flow plan and the velocity hodograph plane can be established in the form of a generalizing relationship between the conjugate velocity and the derivative of the complex potential along the coordinate [3]:

dz — dx + idy -1 dq + ih0 dy j 1 eá

(2)

account the internal friction of the flow and its friction against the channel has the form2, 3 [2, 3, 22]:

d<p ~dx

1-3t

¥

2H

t (1 - T)2 59'

d<p _ 2h0 t 5y

59 " H5x '

(3)

One of the solutions to system (3) are the stream functions and the potential function of the form:

¥ = -4 sin6;

T

a h

9 = A

cosO

(4)

H

r1/2

(1 - T)

Assuming along the streamline dy = 0, and taking into account that on the axis of symmetry of the flow 8 = 0 taking into account (2) and (4), we obtain an ordinary differential equation relating dx, dx:

dx =

Ah

(3T -1)

2Hj2gH0 x2a - T)2

(5)

where indicated A = V0b/2sin0max — constant for the entire flow; t is the parameter of the flow kinetics.

Integration of equation (5) taking into account the initial conditions x — XJD+ X , t = t0 allows you to get a dependency of the form:

x = X'D + X +

Ah

2H 0j2gH~0

1 + t , 1 - T

--ln--

1 + T„

T(1 - T)

>( - To)

(6)

The values of the flow velocity and its depth are found by the formulas [6]:

v = t'2j2hg ;

h =hn(\ -t).

(7)

where 0 is the angle characterizing the direction of the velocity vector to the axis Ox; V is the modulus of the liquid particle velocity vector.

The main system of planned equations of motion for stationary, potential water flows without taking into

The following problems [22, 23] can be solved using the formula (6).

Problem 1. By setting the kinetics parameter on the flow symmetry axis, the abscissa of this point on the symmetry axis, counting from the edge of the pipe, can be determined, then the depth at this point and the flow velocity can be found.

Problem 2. By setting the abscissa of a point on the longitudinal axis of flow symmetry, the kinetic

2 Certificate of state registration of computer programs No. 2022618552. Determination of Parameters of a Freely Spreading Flow. Burtseva ОА. etal. 05.12.2022.

3 Certificate of state registration of computer programs No. 2022666655. Determination of Flow Parameters Along the Extreme Current Line. Aleksandrova M.S. 08.29.2022.

< DO

iH

kK

G !

S 2

0 œ t co

1 i

y 1

J to

u-

^ I

3 o

i 3

o i

o 3

co

CO

3

i 6 r §6

i )

® m

» DO ■ £

s □

s y c o <D *

S°S°

2 2 O O 2 2 W W

W (0

N N

o o

tv N

oo oo

K o

u 3

> in

E M 2

U 00

. r

« Q

?!

<D <u

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

O ig

o o

ig< o

CO

™ O

o

ro

o

E o

CL ° c

Ln O

s H

o E

CD ^

■8

ü (0

parameter can be found from equation (5), and further, the flow velocity and depth at that point can be determined.

Problem 3. By setting the flow depth, it is possible to determine the flow velocity, the kinetics parameter, and then from equation (5) the abscissa of this point on the symmetry axis, counting from the pipe edge.

The set tasks are implemented as separate blocks in MathCad environment3.

Specifying the distribution law of the flow parameters along the outermost current line. To determine the coordinates of the outermost current line, the coupling equation between the physical flow plane and the velocity hodograph plane (2) is used.

For the current and potential functions from equation (2), taking into account solution (4), we obtain a system of differential equations that is valid along the outermost flow line [3]:

dx — —

Ah0 cos0 2 H 0T1,2j2gH0

(1 — 3x) cos0 , x sin0 „ dx +--—- d0

2T(1 — T)2 T1/2 1 — T T3/2

dy - —

Ah0 sin0 2 H 0T1,2j2gH0

(8)

(1 — 3x) cos0 , x sin0

-dx +--— d0

2x(1 — t)2 T1'

1 — x x3

A long the outermost flow line the dependence is valid [3, 12]:

sin0

= sin0„

(9)

dT

cos0d0 = 1/2sin0—.

T

(10)

Equation (10), taking into account dependence (9), is reduced to the form:

dx

cos0d0 = 1/2—sin0max.

T

Taking into account the transformations made above, we transform the system of equations (7) to the form:

dx =

Ah

2 H 04lgH0

3t -1 2sin2 0m

dv = Ah

t2(1 - T)2 (1 - T)2 cos0

dT,

H ojlgH,

,t1/2(1 - T).

(11)

Integration of system (11) allows one to obtain the parametric equations of the outermost flow line in the form [3]:

x = XD +

Ahn

2H 0yßgH0

1 + t 2sm20„

- ln

1 + t.

t i1 - t )

2sin2 0„„

1 - t__^ +_=

t Tk i1 - Tk) 1 - Tk cos0

ln

1 - t 1 - t„

b Ah0 sin2 0m y = - + — --m

2 H^V2iH

"k

cos0.

dT, (12)

T1/2 (1 - T) Tk,2(1 - Tk )

Since along the outermost flow line dy = 0, then from the first equation of the system (4) we obtain the connection between the parameters:

System of equations (12) at known parameters 0k and Tk [3] at the corner point of the outlet edge pipe of the outermost flow line allows to determine the coordinates of the very point X, Y of an arbitrary point of the outermost flow line.

Using the system of equations (12) the following problems can be solved.

Problem 4. On the equipotentiality, the kinetic flow parameters and the angle of inclination of the flow velocity vector to the longitudinal axis of symmetry have the same value. By setting the parameter tw at the point W on the flow symmetry axis, it is possible to determine the parameters tc, 0c of the point C at the outermost flow line and then the coordinates xC, yC of this point. And further changing the parameter tw we obtain a set of points outermost flow line.

Problem 5. Assuming in the first equation of the system (12) x = xC, we determine the kinetic parameter tc at the point C of the flow at the outermost flow line and then the angle of the velocity vector to the longitudinal axis of symmetry of the flow 0C at this point. Substituting the found flux kineticity into the second equation of the system (12), we determine the ordinate of the outermost flow line >>C.

Thus, given a flow extension PC = VJb0 it is possible to estimate the distance mismatch for a selected flow cross-section.

The authors of this study are actively pursuing research in this area, improving and refining the proposed algorithms.

RESULTS

Let us consider the tasks. Experimental data on the free spreading of the flow behind a rectangular culvert are borrowed from [3], where the experimental setup, flow measuring instruments, other details of the experiment are described in detail, see Table 1. The article presents only the results of a numerical calculation using the developed programs2, 3.

The flow has the following characteristics:

• initial flow velocity V0 = 147.654 cm/s;

• initial flow depth relative to the bottom h0 = 9.27 cm;

• acceleration of gravity g = 981 cm/s2;

• pipe width b = 16 cm.

T

Table 1. Point characteristics of free spreading of the plan flow

XE, cm 4 24 44 64 71

yf cm 9.5 38 59 76 80

h, cm 8.5 2.66 1.54 1.12 1.09

V, cm/s 151.928 186.461 191.243 192.714 191.49

then find flow velocity. The results of the problem solution are given in Table 4. For the convenience of analysis, the data at the same points as in Table 3 are given.

It should be noted that there is no discrepancy between the calculated data in Table 3 and Table 4. This is due to the fact that the calculation is performed using analytical relationships.

Next we find:

• Froude's number F0 = 2.397;

• hydrodynamic pressure H0 = 20.382 cm;

• initial flow kinetics t0 = 0.545;

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

• inertial front length X'D = 3 cm;

• wave angle at a point M0 a0 = 0.702 r or a0 = 40°23'. 1. Let us set a kinetic parameter on flow symmetry

axis, dividing segment [t0, 1] into equal parts in increments At = 4.548 • 10-3. Then we calculate the abscissa on the symmetry axis corresponding to the value of kinetic flow according to formula (6). Then we determine flow velocities and depths at the corresponding point. The results are given in Table 2. For the sake of convenience, the data for points where abscissa are close to experimental values are given.

Table 2. Determination of the abscissa on the longitudinal axis of the flow according to a given value of the kinetic parameter

Step number t. i X., cm V, cm/s h , cm

70 0.864 24.477 185.831 2.772

83 0.923 43.493 192.087 1.569

89 0.95 66.978 194.907 1.121

90 0.955 73.612 195.373 1.019

Step number t. i X,, cm V, cm/s h , cm

6 0.861 24 185.548 2.834

11 0.924 44 192.179 1.558

16 0.948 64 194.664 1.068

18 0.953 72 195.267 0.948

Table 4. Definitions of the kinetic parameter, point abscissa and flow velocity on the axis of symmetry for a given flow depth

Step number t. i X., cm V, cm/s h , cm

6 0.861 24 185.548 2.834

11 0.924 44 192.179 1.558

16 0.948 64 194.664 1.068

18 0.953 72 195.267 0.948

2. In order to determine the kinetic parameter of the flow, let us set an abscissa on its longitudinal axis with the least deviation with experimental data, for instance 4 cm step. Solving equation (6), we find the only root corresponding to the kinetic parameter at the point on the longitudinal axis of the stream with given abscissa. Then we find the velocity and the depth of the flow at the given point. The results of the problem solution are given in Table 3.

3. The depth of flow h we set, using the calculations in problem 2. The kinetic parameter is determined from formula (7). Then we find abscissa on flow symmetry axis corresponding to this kinetic parameter from equation (6),

Table 3. Determination of the kinetic parameter for a given abscissa

4. Now let us consider the determination of flow parameters along the outermost flow line. From the second equation of system (4) it is seen that an equipotential can be singled out by a specific value of the parameter t on the flow symmetry axis. By setting the parameter t0 for the point M0 (Figure) on the flow symmetry axis, we determine the parameters tc, 0c followed by the coordinates xC, yc of the point C on the outermost flow line corresponding to this equipotentiality. Then by changing the kinetic parameter (e.g. with a constant step) we obtain the set of points of the outermost flow line.

5. Assuming in the first equation of the system (12) x = xC, we determine the kinetic parameter tc at the point C of the flow at the outermost flow line and then the angle of the velocity vector to the longitudinal axis of symmetry of the flow 0c at this point. Substituting the found flux kineticity into the second equation of the system (12), let us determine the ordinate of the outermost flow line yc. The results coincide with the numerical experiment described in the previous paragraph.

The numerical experiment consisted in calculating the values of the kinetic parameter, the angle of inclination of the flow velocity vector to the symmetry axis and the coordinates of the points on the outermost flow line. The results are given in Table 5. In order to determine the adequacy of the parameter calculation algorithm, a relative error of the flow ordinate has been calculated in comparison with experimental data.

CONCLUSIONS

The formation of a mathematical model of water flow is based on its division into several sections. At the stage of flow exit from the pipe there is a section where flow parameters (velocity, depth, width) are kept constant — inertial front. The formula of the inertial front length of the water stream has been derived at its spread from a rectangular pipe into a wide outlet

< n

8 8 is

o

S

c

O w t co

y i

J CD

u-

^ I

DD 3

0 CD

01

o 3

co co

l\J CO

o DD§

r §6

o o

CD

o

C 3

CD )

ii

® 8

» DO

■ T

s □

s y

c o

® S

CO CO

2 2

O O

2 2

W W

Table 5. Definitions of the kinetic parameter, the angle of inclination of the flow velocity vector to the axis of symmetry, and the ordinate of points on the extreme streamline according to the given kinetic parameters on the axis of symmetry of the flow

Step number X., cm t. on the axis of i symmetry TCi on outermost flow line YCi, cm/s \y - y\ 1 C e\ 100% Ya

0 0 0.545 0.545 0.661 8 0

1 4 0.545 0.731 0.791 8.97 5.914

6 24 0.861 0.952 0.946 34.926 9.854

11 44 0.924 0.975 0.963 63.3 6.241

16 64 0.948 0.983 0.969 92.257 17.275

18 72 0.953 0.985 0.97 103.922 22.729

W (0 N N

o o

N N

00 00 K (V U 3 > in C M 2

HQ 00

. r

« q

?!

<D <1J

O £

o

o g<

o

CO

™ O

o

ro

o

E O

CL ° c

Ln O

s H

o E

CD ^

■8

0 (0

channel by introducing dimensionless complexes and on the base of n-theorem.

Obtained formula for the length of the inertial front, allows us to achieve the desired error. At flow extensions up to 5, relative errors of ordinates and flow velocities do not exceed 7-10 %. At steps 16 and 18 the error is more than 10 %. This can be explained by the fact that these points are close to the transition zone. This area should be investigated further.

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

Using analogies from gas dynamics, namely the transition to the flow velocity hodograph plane, it is possible to obtain the depth and velocity distribution laws along the longitudinal axis of symmetry and along the outermost flow line. With the help of system of equations (12) set tasks have been successfully solved. Software has been developed for realization of these algorithms. On the basis of experimental data, the adequacy of the refined mathematical model of two-dimensional planar flow has been proved.

REFERENCES

1. Sherenkov I.A. Experimental studies of the spreading of a turbulent flow behind the outlet heads of culverts. Proceedings of the joint seminar on hydraulic engineering and water management construction. Kharkiv, 1958; 1:39-43. (rus.).

2. Emcev B.T. Two-dimensional turbulent flows. Moscow, Energiya Publ., 1967; 212. (rus.).

3. Kokhanenko V.N., Volosukhin Ya.V., Le-meshko M.A., Papchenko N.G. Modeling of stormy two-dimensional in plane water flows. Rostov-on-Don, Southern Federal University Publ., 2013; 180. (rus.).

4. Konstantinov N.M. Some issues of the hydraulics of the tailwater of small road culverts with free spreading of a turbulent flow. Hydraulics of road culverts. 1969; 255-269. (rus.).

5. Kosichenko N.V. About the petal of the free flow of a turbulent stream into a wide fortified channel. Prirodoobustrojstvo (Environmental Engineering). 2011; 3:58-62. URL: https://cyberleninka.ru/article/n/o-lepestke-svobodnogo-rastekaniya-burnogo-potoka-v-shirokoe-ukreplennoe-ruslo (rus.).

6. Kosichenko N.V. Analysis of the study and refinement of methods for free flow spreading behind non-pressure culverts. Bulletin of the Saratov State Agrarian University. 2011; 9:27-33. (rus.).

7. Kohanenko V.N., Mitsik M.F., Kosichenko N.V. The more accurate equation of the extreme line current in the hodograph plane speed in free flow stormy spreading for unpressurised pipes. Bulletin of Higher

Educational Institutions. North Caucasus region. Technical Sciences. 2013; 1:33-35. URL: https://cyberlenin-ka.ru/article/n/utochnennoe-uravnenie-krayney-linii-to-ka-v-ploskosti-godografa-skorosti-v-zadache-svobodno-go-rastekaniya-burnogo-potoka-za-beznapornymi (rus.).

8. Micik M.F., Kohanenko N.V., Kosichenko N.V. Determination of parameters of a turbulent flow, freely spreading from a pressureless tube in a broad channel without the flooding of the lower canal pond. Privol-zhsky Scientific Journal. 2013; 4:52-56. URL: https:// www.elibrary.ru/item.asp?id=21306122 (rus.).

9. Tang S.L., Antonia R.A., Djenidi L., Abe H., Zhou T., Danaila L., Zhou Y. Transport equation for the mean turbulent energy dissipation rate on the centreline of a fully developed channel flow. Journal of Fluid Mechanics. 2015; 777:151-177. DOI: 10.1017/ jfm.2015.342

10. Anderson W., Barros J.M., Christensen K.T., Awasthi A. Numerical andexperimental study of mechanisms responsible for turbulent secondary flows in boundary layerflows over spanwise heterogeneous roughness. Journal of Fluid Mechanics. 2015; 768:316347. DOI: 10.1017/jfm.2015.91

11. Aranda J.A., Beneyto C., Sanchez-Juny M., Blade E. Efficient design of road drainage systems. Water. 2021; 13:1661. DOI.org/10.3390/w13121661.

12. Anees M.T., Abdullah K., Nordin M.N., Rahman N.N., Syakir M.I., Kadir M.O. One- and two-dimensional hydrological modelling and their uncertain-

ties. Flood Risk Management. 2017; 11:221-244. DOI: 10.5772/intechopen.68924.

13. Nematollahi B., Abedini M.J. Analytical solution of gradually varied flow equation in non-prismatic channels. Iranian Journal of Science and Technology — Transactions of Civil Engineering. 2020; 44(1):251-258. DOI: 10.1007/s40996-019-00316-5

14. Hager W., Castro-Orgaz O. Transcritical flow in open channel hydraulics: from boss to de marchi. Journal of Hydraulic Engineering. 2016; 142(1). DOI: 10.1061/(asce)hy.1943-7900.0001091

15. Hager W. Unconfined expansion of supercritical water flow. Journal of Engineering Mechanics. 1997; 123(5):451-457. DOI: 10.1061/(ASCE)0733-9399(1997)123:5(451)

16. Liu J.L., Wang Z.Z., Zhao Y.F. Explicit equations for critical depth in open channels with complex compound cross sections. Flow Measurement and Instrumentation. 2012; 24:13-18. DOI: 10.1016/j.flow-measinst.2011.12.005

17. Chanson H. Explicit equations for critical depth in open channels with complex compound cross sections: A discussion. Flow Measurement and Instrumentation. 2013; 29:65-66. DOI: 10.1016/j.flowmea-sinst.2012.10.006

18. Castro-Orgaz O., Cantero-Chinchilla F.N. Non-linear shallow water flow modelling over topography with depth-averaged potential equations. Environmental Fluid Mechanics. 2020; 20(2):261-291. DOI: 10.1007/s10652-019-09691-z

19. Li J., Li S.S. Near-bed velocity and shear stress of open-channel flow over surface roughness.

Environmental Fluid Mechanics. 2020; 20(2):293-320. DOI: 10.1007/s10652-019-09728-3

20. Jesusdhas V., Balachandar R., Wang H., Murzyn F. Modelling hydraulic jumps: IDDES versus experiments. Environmental Fluid Mechanics. 2020; 20(2):393-413. DOI: 10.1007/s10652-019-09734-5

21. Leng X., Chanson H. Hybrid modelling of low velocity zones in box culverts to assist upstream fish passage. Environmental Fluid Mechanics. 2020; 20(2):415-432. DOI: 10.1007/s10652-019-09700-1

22. Kokhanenko V., Burtseva O., Kelekhsaev D. Here determination of the free flow maximum expansion into a wide diversion channel of turbulent water flow behind the non-pressure water throughput hole. IOP Conference Series: Materials Science and Engineering. 2020; 753(3):032053. DOI: 10.1088/1757-899X/753/3/032053

23. Kohanenko V.N., Kelekhsaev D.B. Solution of the problem of determining the equation of the extreme streamline and the parameters along it taking into account the section XD. Research Results — 2019: materials of the IV National Conference of the Faculty and Researchers (Novocherkassk, May 14, 2019). Novocherkassk, 113-117. (rus.).

24. Sedov L.I. Methods of similarity and dimension in mechanics. Moscow, Nauka Publ., 1977; 440. (rus.).

25. Romanov A.S., Semikolenov A.V., Taranen-ko S.N., Shahorin A.P. Theory of similarity and dimension. Boundary layer. Moskva, MGTU im. N.E. Baumana, 2011; 48. (rus.).

< П

8 8 iH

kK

G Г

S 2

со со

Received March 30, 2023.

Adopted in revised form on May 25, 2023.

Approved for publication on July 7, 2023.

BioNoiEs: Olga A. Burtseva — Candidate of Technical Sciences, Assistant Professor, Associate Professor of the Department of General Engineering Disciplines; Platov South-Russian State Polytechnic University (NPI); 132 Prosveshcheniya st., Novocherkassk, 346428, Russian Federation; ID RSCI: 161751, Scopus: 650-7800-801, ResearcherlD: ABG-9531-2020, ORCID: 0000-0003-4288-8709; [email protected];

Maria S. Alexandrova — postgraduate student of the Department General Engineering Disciplines; Platov South-Russian State Polytechnic University (NPI); 132 Prosveshcheniya st., Novocherkassk, 346428, Russian Federation; ORCID: 0000-0003-2901-7116; [email protected].

Contribution of the authors:

Olga A. Burtseva — idea, writing an article, translating the text into English, editing the text after computer typing, final conclusions.

Maria S. Alexandrova — compilation of a program for calculating flow parameters, calculation by algorithms, verification of the adequacy of the mathematical model.

All authors have made an equivalent contribution to the publication and declare no conflicts of interest.

СПИСОК ИСТОЧНИКОВ

1. Шеренков И.А. Экспериментальные исследования растекания бурного потока за выходными оголовками водопропускных сооружений // Труды

объединенного семинара по гидротехническому и водохозяйственному строительству. Харьков, 1958. Вып. 1. С. 39-43.

y ->■

J CD I

n

i 3 о

n

со со

n

a g

i §

r §6

i )

m

. DO

■ г

s S

s у с о <D *

S0 S0

2 2 О О 10 10 U W

2. Емцев Б.Т. Двухмерные бурные потоки. М. : Энергия, 1967. 212 с.

3. Коханенко В.Н., Волосухин Я.В., Лемеш-ко М.А., Папченко Н.Г. Моделирование бурных двухмерных в плане водных потоков. Ростов н/Д : Издательство Южного федерального университета, 2013. 180 с.

4. Константинов Н.М. Некоторые вопросы гидравлики нижнего бьефа малых дорожных водопропускных сооружений при свободном растекании бурного потока // Гидравлика дорожных водопропускных сооружений. 1969. С. 255-269.

5. Косиченко Н.В. О лепестке свободного растекания бурного потока в широкое укрепленное русло // Природообустройство. 2011. № 3. С. 58-62. URL: https://cyberleninka.ru/article/n/o-lepestk:e-svobodnogo-rastekaniya-burnogo-potoka-v-shirokoe-ukreplennoe-ruslo

6. Косиченко Н.В. Анализ изучения и уточнения методов свободного растекания потока за безнапорными водопропускными отверстиями // Вестник Саратовского ГАУ. 2011. № 9. С. 27-33.

7. Коханенко В.Н., МицикМ.Ф., Косиченко Н.В. Уточненное уравнение крайней линии тока в плоскости годографа скорости в задаче свободного растека-

й Я ния бурного потока за безнапорными водопропуск-

5 5 ными трубами // Известия вузов. Северо-Кавказский со со регион. Технические науки. 2013. № 1. С. 33-35. URL: о § https://cyberleninka.ru/article/n/utochnennoe-uravnenie-е « krayney-linii-toka-v-ploskosti-godografa-skorosti-v-щ gj zadache-svobodnogo-rastekaniya-burnogo-potoka-za-oj ^ beznapornymi

8. Мицик М.Ф., Коханенко Н.В. Косиченко Н.В. о -Ц Определение параметров свободно растекающегося

• «* бурного потока в окрестности его выхода из безна-£ £ порной прямоугольной трубы в широкое отводя-"£= о щее русло // Приволжский научный журнал. 2013. ~ f № 4. С. 52-56. URL: https://www.elibrary.ru/item. g J asp?id=21306122

cd

4 "g 9. Tang S.L., Antonia R.A., Djenidi L., Abe H.,

po ^ Zhou T., Danaila L., Zhou Y. Transport equation for

z the mean turbulent energy dissipation rate on the centre-

OT E line of a fully developed channel flow // Journal of Fluid

~ § Mechanics. 2015. Vol. 777. Pp. 151-177. DOI: 10.1017/

£ О jfm.2015.342

Й ° 10. Anderson W., Barros J.M., Christensen K.T.,

CO —

о E Awasthi A. Numerical andexperimental study of mech-

6 о anisms responsible for turbulent secondary flows in

O)

boundary layerflows over spanwise heterogeneous от "i= roughness // Journal of Fluid Mechanics. 2015; 768:3167 1 347. DOI: 10.1017/jfm.2015.91 Sj -j 11. Aranda J.A., Beneyto C., Sanchez-Juny M.,

L_ g Blade E. Efficient design of road drainage systems // ® EE Water. 2021. Vol. 13. P. 1661. DOI: 10.3390/

* s£ w13121661

¡3 12. Anees M.T., Abdullah K., Nordin M.N.,

bq £ Rahman N.N., Syakir M.I., Kadir M.O. One- and two-dimensional hydrological modelling and their

uncertainties // Flood Risk Management. 2017. Vol. 11. Pp. 221-244. DOI: 10.5772/intechopen.68924

13. Nematollahi B., Abedini M.J. Analytical solution of gradually varied flow equation in non-prismatic channels // Iranian Journal of Science and Technology — Transactions of Civil Engineering. 2020. Vol. 44. Issue 1. Pp. 251-258. DOI: 10.1007/s40996-019-00316-5

14. Hager W., Castro-Orgaz O. Transcriti-cal flow in open channel hydraulics: from boss to de marchi // Journal ofHydraulic Engineering. 2016. Vol. 142. Issue 1. DOI: 10.1061/(asce)hy.1943-7900.0001091

15. Hager W. Unconfined expansion of supercritical water flow // Journal of Engineering Mechanics. 1997. Vol. 123. Issue 5. Pp. 451-457. DOI: 10.1061/ (ASCE)0733-9399(1997)123:5(451)

16. Liu J.L., Wang Z.Z., Leng C.J., Zhao Y.F. Explicit equations for critical depth in open channels with complex compound cross sections // Flow Measurement and Instrumentation. 2012. Vol. 24. Pp. 13-18. DOI: 10.1016/j.flowmeasinst.2011.12.005

17. Chanson H. Explicit equations for critical depth in open channels with complex compound cross sections: A discussion // Flow Measurement and Instrumentation. 2013. Vol. 29. Pp. 65-66. DOI: 10.1016/j. flowmeasinst.2012.10.006

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

18. Castro-Orgaz O., Cantero-Chinchilla F.N. Non-linear shallow water flow modelling over topography with depth-averaged potential equations // Environmental Fluid Mechanics. 2020. Vol. 20. Issue 2. Pp. 261-291. DOI: 10.1007/s10652-019-09691-z

19. Li J., Li S.S. Near-bed velocity and shear stress of open-channel flow over surface roughness // Environmental Fluid Mechanics. 2020. Vol. 20. Issue 2. Pp. 293-320. DOI: 10.1007/s10652-019-09728-3

20. Jesusdhas V., Balachandar R., Wang H., Murzyn F. Modelling hydraulic jumps: IDDES versus experiments // Environmental Fluid Mechanics. 2020. Vol. 20. Issue 2. Pp. 393-413. DOI: 10.1007/s10652-019-09734-5

21. Leng X., Chanson H. Hybrid modelling of low velocity zones in box culverts to assist upstream fish passage // Environmental Fluid Mechanics. 2020. Vol. 20. Issue 2. Pp. 415-432. DOI: 10.1007/s10652-019-09700-1

22. Kokhanenko V., Burtseva O., Kelekhsaev D. Here determination of the free flow maximum expansion into a wide diversion channel of turbulent water flow behind the non-pressure water throughput hole // IOP Conference Series: Materials Science and Engineering. 2020. Vol. 753. Issue 3. P. 032053. DOI: 10.1088/1757-899X/753/3/032053

23. Коханенко В.Н., Келехсаев Д.Б. Решение задачи определения уравнения крайней линии тока и параметров вдоль нее с учетом участка ХД // Результаты исследований — 2019 : мат. IV Национальной конф. проф.-преп. состава и научных работников

(г. Новочеркасск, 14 мая 2019 г.). Новочеркасск, 25. Романов А.С., СемиколеновА.В., Таранен-

2019. С. 113-117. ко С.Н., Шахорин А.П. Теория подобия и размерно-

24. Седов Л.И. Методы подобия и размерности сти. Пограничный слой. М. : МГТУ им. Н.Э. Баумана,

в механике. М. : Наука, 1977. 440 с. 2011. 48 с.

Поступила в редакцию 30 марта 2023 г. Принята в доработанном виде 25 мая 2023 г. Одобрена для публикации 7 июля 2023 г.

Об авторах: Ольга Александровна Бурцева — кандидат технических наук, доцент кафедры общеинженерных дисциплин; Южно-Российский государственный политехнический университет (НИИ) имени М.И. Платова (ЮРГИУ (НИИ) имени М.И. Платова); 346428, г. Новочеркасск, ул. Просвещения, д. 132; РИНЦ ID: 161751, Scopus: 650-7800-801, ResearcherlD: ABG-9531-2020, ORCID: 0000-0003-4288-8709; [email protected];

Мария Сергеевна Александрова — аспирант кафедры общеинженерных дисциплин; Южно-Российский государственный политехнический университет (НИИ) имени М.И. Платова (ЮРГИУ (НИИ) имени М.И. Илатова); 346428, г. Новочеркасск, ул. Просвещения, д. 132; ORCID: 0000-0003-2901-7116; sergand1957@ gmail.com.

Вклад авторов:

Бурцева О.А. — идея, написание статьи, перевод текста на английский язык, редактирование текста после компьютерного набора, итоговые выводы.

Александрова М.С. — составление программы расчета параметров потока, расчет по алгоритмам, проверка адекватности математической модели.

Все авторы сделали эквивалентный вклад в подготовку публикации и заявляют об отсутствии конфликта

< Ш

интересов. e ф

t О

Ü

со со

n ю a 0

Ф6 >86

ф )

ii

ф 8

Ю DO

■ Т

(Л У

с о

Ф ф

со 00

2 2

О О

2 2

W W

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