CC BY
20ТЕХНИЧЕСКИЕ НАУКИ
Qurbonalizoda S.S.
Dr., Deputy Director of cascade Varzob 's HPPs OSHC «Barqi Tojik», Tajikistan
Gulakhmadov A.A. Dr. Chief specialist of energy department Ministry of energy and water resources of the Republic of Tajikistan Kobuliev Z. V. Dr., professor, ccorresponding-member, Director Institute of water problems, hydropower and ecology of Academy of sciences of the Republic of Tajikistan Kholmuhamadzoda S.A. PhD student, Head of energy department Ministry of energy and water resources of the Republic of Tajikistan
ANALYSIS OF THE SPATIAL STRUCTURE OF THE FLOW IN HIGH PRESSURE HYDRAULIC TURBINES
Summary. The article discusses the study of the velocity field and pressure conducting on the hydraulic turbine type R0500 (Francis turbine radial-axial 500 pressure) and the basic hydrodynamic equations used in the process of mathematical modeling of turbulent flow are based on the fundamental laws of conservation: mass, momentum, energy. This paper presents the distribution of flow rates before and after the impeller of high-pressure R0500, as well as the loss of flow energy in the spiral chamber and guide vanes.
Key words: HPP, hydraulic turbine, spiral chamber, impeller, Francis turbine, energy.
Introduction
In the 1950s and 1960s of the last century, in connection with the widespread development of Siberian rivers, as well as the construction of a cascade of hydro power stations in Ukraine, Central Asia, the Caucasus and the Baltic states, the development of experimental and computational-theoretical works in the development of flowing parts and impellers received a great impetus hydraulic turbines of various speeds.
During these years, the domestic nomenclature (industry standard) was developed. The efficiency level of the developed model hydraulic turbines reached 8789% at the optimum.
With the development of two-dimensional and three-dimensional methods for calculating blade systems in the years 70-80, it was possible to increase the efficiency to 91-92%.
But still, as in the 50-60 years, for example, when developing the impellers of the Krasnoyarsk or Sayano-Shushenskaya hydro power plants, up to 20 impellers and their modifications were tested on a model installation, the experimental refinement of the flow parts and impellers was the main and most laborious part of the job [9-11].
With the widespread introduction of mathematical modeling based on computer technologies into the process, which made it possible to increase the level of optimal efficiency to 93 ^ 94% (with Dlmod = 0.5 m), doubts arose about the importance and need for further experimental studies of the flow structure [2,7].
The paper substantiates the feasibility of experimental studies and their importance for improving the mathematical models used in computer technology.
The relevance of the study
The basic equations of hydrodynamics used in the process of mathematical modeling are based on the fundamental laws of conservation: mass, momentum, angular momentum, energy.
In the process of developing methods for mathematical modeling of turbulent flows, simplifications are introduced into the process which is under consideration.
First, diffusion and convective transport is carried out as the simultaneous transfer of mass, energy, momentum and angular momentum.
Naturally, the desire to take into account in the process of mathematical modeling this complex transfer process. However, it is known that the Navier-Stokes equation is a differential analog of the momentum conservation law, and the continuity equation is a differential analog of the mass conservation law. The law of conservation and transfer of angular momentum are not directly applied in mathematical modeling.
Secondly, a system of time-averaged Reynolds and continuity equations is used to describe the behavior of a turbulent flow.
Averaging the flow over a selected cross section or over a time interval can be carried out in various ways:
Wschodnioeuropejskie Czasopismo Naukowe (East European Scientific Journal) #3(55), 2020
22
For example, when determining the average speed, you can set the task a) average speed should provide the same flow rate that occurs in a real stream, or b) same amount of movement (momentum), or c) same moment of momentum (angular momentum), or d) the same kinetic energy.
Of course, these will be different values of the averaged speeds.
Averaging the Navier-Stokes equations (equations of momenta), we obtain an adequate model of force interaction on an elementary liquid volume with a violation of the energy balance and the balance of the angular momentum as applied to this liquid volume.
An attempt to clarify the problem is to use the mathematical modeling of the equations of balance and dissipation of turbulent kinetic energy (the k-s model of turbulence).
As known, writing the Navier-Stokes equations in a conservative (energy) form by multiplying term by term (the tensor form of each term) is averaged over energy.
Furthermore, the equation of balance of the total kinetic energy is obtained, and the equation of balance of the averaged kinetic energy is subtracted from this equation.
The influence of the angular momentum on the process, as a rule, is not taken into account.
Subsequent, in the process of implementing the calculation on a computer, the simplifications of the mathematical model are not amenable to complete analysis from the point of view of the introduced errors and the accuracy of the results obtained.
Thirdly, the influence of the boundary conditions is so significant on the results obtained, both in terms of the accuracy of the averaged flow parameters specified in the input section and in terms of such turbulence characteristics as scale, intensity, and other parameters [1].
As a result of the above mentioned, the experimental verification of the flow parameters obtained as a result of mathematical modeling, as well as the verification of the accuracy of setting the boundary conditions, is the main criterion for the accuracy of the calculated results.
Flow structure analysis
In the case of solving the direct problem of the flow around the impeller blades, the actual flow form at the entrance to the blade system, and for some (channel) methods and at the exit from it, it is necessary to know in order to set the boundary conditions. With the known of the distribution of velocities and pressures in the characteristic sections of the turbine cavity, one can find the corresponding averaged characteristics, which is primarily necessary for determining the energy balance of the hydraulic turbine.
The analysis of the energy balance thus compiled allows us to study a number of important issues in the theory of the working process: the degree of influence of individual elements of the flow part and the establishment of the most effective ways to increase speed, energy and cavitation indicators.
The study of the velocity and pressure fields was carried out on a hydraulic turbine bench in the flow part of the radial-axial turbine model with impellers Francis turbine R0500 I-2b and RO500/683 (Russian standard P0500 I-26 h P0500/683) [3-6]. The main elements of a hydraulic turbine model are characterized by the following parameters:
1) Circular spiral camera with a coverage angle in plan 9=345° and the radius of the input section
p=133mm (arated t= = 0.6), right rotation calculated
Vwp
by law vur = const;
2) Stator with 13 profiled columns (including spiral tooth);
3) Profile of guide vanes normalized, asymmetric, positive curvature, ^ = 1.2; ^ = 0.08; z0 = 24;
4) Impellers (Di=400 mm): R0500 I-2b (Francis turbine), z = 17, ^ = 0.686; R0500/683, z = 19,
^ = 0.686;
D±
5) Suction pipe (draft tube) - straight conical with a taper angle 5°, height h = 3.18DX.
The velocity and pressure field in the turbine cavity was measured in two sections: before of the impeller (section I-I) (fig. 1) and after it (section II-II).
From the input section of the spiral
Figure 1. The location of the dimensional sections in the turbine cavity and the decomposition scheme of the absolute velocity vector
In the plan, the cylindrical section I-I is located to the output section of the spiral (in the direction of rotation of the impeller) at an angle (p= 142°30'^160°30' = 18°) at the impeller RO500 I-2b and at an angle (=142°30'^162°(40 = 19°30) at the impeller R0500/683. Section II-II in plan is at an angle (=276° to the output section of the spiral.
When assessing the influence of the shape of the flow in a spiral on its structure in front of the impeller, it became necessary to additionally study the flow in section III-III, which is located at an angle (= 112°30' to the input section of the spiral and passes through the center of the circle of the section. The flow was studied using probes in the cavity of a hydraulic turbine with an impeller РО500I-2b in 22 heads of pressure, and with the impeller R0500/683 (Francis turbine) - in six heads of pressure H = 4.5 ^ 5.5m.
The following can be cited as the main regularities of the structure of the averaged flow (Fig. 2 - 6):
1. The energy loss of the flow in the spiral chamber and the guiding vanes for the studied zone of regimes does not significantly depend on the regime of operation of the turbine and amount to about 2 ^ 3.5%. At the upper and lower rings of the guide vane, the total energy is 1 ^ 2% less than in the central part.
2. The nature of the unevenness in height of the guide vane averaged in the circumferential direction of the flow velocity is determined primarily by the value of a0 and does not depend on the number of rotations of the impeller. This is indicated by dependency matching
%»). Vrm.
, ; , ; , (given to 1 m head and averaged in Qi Qi Qi
the circumferential direction) at a0 =17 mm and differences nj (n\ = 80; 75; 65; 55 and 46,5 min-1).
By comparing these graphs with similar ones in the case of a cavity without an impeller or with other impellers, it was established that the pattern of change in speed before of the impeller along the height of the
Wschodnioeuropejskie Czasopismo Naukowe (East European Scientific Journal) #3(55), 2020
24
guide vanes is specific for hydraulic turbines of this speed. With increase ao, uneven values
QÍ
QÍ
QÍ
the height of the guide vane is leveled.
3. When n' changes over a wide range (from values nj = 80min-1 untill n' = 46.5min-1) the angle change graph a^ ( between the direction of absolute speed and its peripheral component) the height of the guide vane practically does not change and only at changing a0 angle a^ is changing, increase with increasing a0.
If we compare the nature of the change in the angle a with the same openings a0 for the case with an impeller of RO500 nj = 65min-1 and without an impeller, we find that in the absence of an impeller, the angle a is less in 2°r4°.
Thus, the influence of the lattice of the impeller blades on the flow behind the guide vanes is established. This influence consists in some equalization of the flow unevenness in the circumferential direction and in the height of the guide vanes and in an increase (by 2°r4°) of the angle.
It was found that uneven flow patterns
(W°';
, , , is a function of the angle a flow Q/ Q/ 6
averaged over the entire cylindrical surface and the
dependence is obtained a from opening the guide vane.
Thus, in advance, when profiling a new impeller,
it is possible to take into account the uneven flow along
the height of the guide vane.
4. As a result of the analysis of graphs of changes in the angle /?(0) flow (between relative and peripheral speed) height of the guide vane for various a0 and n' we can conclude that when changing nj ranging from n' = 80 min-1 until n' = 46.5 min-1 angle ^ varies accordingly from the values 22° at the bottom ring of the guide vane, 20° in the central part and 27° at the upper ring when n'j = 80 min-1 until 112° at the lower ring, 133° in the central part and 122° at the upper ring whennj = 46.5 min-1.
Thus, at turbine operating conditions other than optimal, the flow can flow onto the impeller blade system at angles of attack varying from -50° r -60° go 50° - 60°. In all likelihood, this is one of the main reasons for the sharp change in efficiency when changing nj at a0=const. Such unfavorable conditions of flow inflow limit the possibility of installing high-pressure hydraulic turbines at hydro power plants with significant pressure fluctuations, as well as their operation at starting (low) pressures.
At n' = c onst (nj = 65 min-1) when it changes a0 (from a0=9 mm until a0=29 mm) angle ft from 59° r 78° until 43° r 47°. Smaller values /? occur at the lower ring of the guide vane, and larger v alues at the upper.
5. Vortex flow structure determined by /2r(0), ^z(0), does not change significantly when
changing the opening of the guide vane (nj = const) and almost does not change at a0=const and variable
nJ.
Figure 2. The distribution offlow rates in front of the impeller Francis turbine RO500I-2b (D1 = 400mm, n't= 65min-1)
Figure 3. The distribution offlow rates in front of the impeller Francis turbine (Russian standard RO5QQI-2b) (D1 = 400mm, a0= 17mm)
As a result of the analysis of the parameters of the stream averaged over the cylindrical surface in front of the impeller, the following regularities were established.
1) Lines a = const coincide with the lines ao=const and the difference between d flow and aHa (angle of the midline of the profile of the blade on the output edge), almost equal to zero (0° -r- 0.5°).
2) In the studied area of the turbine, the angle P varies within 30° ^ 120°.
3) The energy loss in the spiral chamber and guide vane varies slightly (from 2% to 3.5%) when changing the turbine operating regime from ao=9 mm to ao=29 mm and from n' = 46.5 min-1 to n' = 80min-1 and are less important at large values nt.
4) For small values n\ (n\ = 46.5 + 55min-1) EK«En«50%, for large values n'and Q', ek«30% and En«70% from full energy.
+40 +30 +20 +10 0 -10 -20 -30 -40 -50 -60 -70 -80 -90
+40 +30 +20 +10 0 -10 -20 -30 -40 -50 -60 -70 -80 -90
Figure 4. Flow characteristic in section I-I before the impeller Francis turbine (Russian standard R0500I-2b) (D1 = 400mm)
5) Quantities (vur)j decrease with increasing study area of the turbine (9mm < a0 < 29mm; nj and $ from (^U^)J = 1.6 to (^U^)j = 1.05 in the 46.5 min-1 < nj < 80 min-1).
Figure 5. Average flow parameters in front of the impeller Francis turbine (Russian standard R0500I-2b)
Ilj= 65 O' r
1 2 3 4 5 6 7 ti 1 iu II 12 13 14
X
1 2 3 4 5 6 7 I 9 10 11 12 13 14
X
X
X
x
X
x
X
a0= 17 mm (D 1=400mm)
1 2 3 4 5 6 7 8 9 10 11 12 13 14
1 2 3 4 5 6 7 8 9 10 11 12 13 14
Figure 6. Flow parameters behind the impeller Francis turbine (Russian standard R0500I-2b)
In the optimal regime, there is a minimum value of the kinetic energy of the flow in the suction pipe of the order of 3-4%. At the a0=9 mm and a0=25 mm, nj = 65 min-1, ek«7%, and at the a0=29 mm, nj = 65 min-1, the kinetic energy of the flow reaches 10-12% of the energy of the input streamflow.
The potential energy of the flow in the suction pipe changes insignificantly when the turbine operating regime changes (n'andQj), as a rule, it has a negative
value (pressure is less than atmospheric) and in the optimal regime En ■ is about 1-2%.
Thus, the energy loss of the flow in the suction pipe at the optimum regime is equal to 1^2% and increase to 4% at ao=25 mm and nj = 65 min-1, and to « 8^9% at the ao=29 mm, nj = 65min-1; also to « 9^10% at ao=17 mm and nj = 80min-1. Basically, these losses are associated with the energy of the
Wschodnioeuropejskie Czasopismo Naukowe (East European Scientific Journal) #3(55), 2020
circulation component of the velocity
(Vuj)2
2g
which at
the optimum regime is close to zero.
The optimum universal characteristics of the best impellers are achieved with a positive, close to zero, value of the averaged circulation at the exit of the impeller.
The presented on figures 2-4 graphs of the distribution of flow parameters in a cylindrical section R=212 mm (for Di=400 mm) impeller Francis turbine R0500I-2b (modification P0400I-2b - stock number
29
3515) can be used in the development of new blade systems for pressure 400^500 m, as well as for mathematical modeling of the flow as boundary conditions.
It should be noted that the use of modern three-dimensional methods of mathematical modeling of the flow, based on rather complex and rigorous turbulence models, requires the same high degree of accuracy of the task, as the averaged flow parameters, and the characteristics of the flowing turbulence as boundary conditions. in the input section.
1 - spiral chamber; 2 - guiding vanes; 3 - impeller; 4 - shaft; 5 - draft tube Figure 7. General view of a section along a radial-axial hydraulic turbine
Conclusions
This article was concluded by the process of developing methods for mathematical modeling of turbulent flows. Firstly, diffusion and convective transport as the transfer of mass, energy, momentum and angular momentum. In the process of mathematical modeling is this complex transfer process. However, it is known that the Navier-Stokes equation is a differential analog of the momentum conservation law, and the continuity equation is a differential analog of the mass conservation law. The law of conservation and transfer of angular momentum are not directly applied in mathematical modeling. Also in the optimal regime, the minimum kinetic energy of the flow in the draft tube and the flow behind the impeller in the spiral chamber and the guide vanes changes slightly (from 2% to 3.5%) with a change in the turbine operating regime.
References
1. Barlit V.V. // Gidravlicheskiye turbiny (Hydraulic turbines). - Kiev, Vishcha school, 1977.
2. Barlit V.V., Borodayevskiy V.V., Yakovleva L.K. // Metodicheskiye ukazaniya k kursovoy rabote, k kursovomu i diplomnomu proyektam po gidroturbinam i obratimym gidromashinam «Vybor tipa gidroturbiny, nasos-turbiny pri pomoshchi universal'nykh kharakteristik, opredeleniye osnovnykh parametrov i postroyeniye yeye ekspluatatsionnoy kharakteristiki s ispol'zovaniyem EVM», (Methodical instructions for the term paper, for the term and graduation projects on hydraulic turbines and reversible hydraulic machines "Choosing the type of hydraulic turbine, pump-turbine using universal characteristics, determining the main parameters and constructing its operational characteristics using computers"). - Khar'kov, KHGPU, 1994.
3. Barlit V.V., Borodayevskiy V.V., Yakovleva L.K. // Metodicheskiye ukazaniya k kursovoy rabote, k kursovomu i diplomnomu proyektam po gidroturbinam i obratimym gidromashinam «Raschet i postroyeniye ravnoskorostnogo meridional'nogo potoka v polosti lopastnoy (RO, DPL, PL) gidromashiny pri pomoshchi EVM», (Methodical instructions for the term paper, for
the term and graduation projects on hydraulic turbines and reversible hydraulic machines "Calculation and construction of an equal-speed meridional flow in the cavity of the blade (RO, DPL, PL) hydraulic machines using a computer"). - Khar'kov, KHGPU, 1994.
4. Kovaleva. L. // Spravochnik po gidroturbinam (Handbook of hydraulic turbines) Mashinostroyeniye, 1984.
5. Borodayevskiy V.V., Yakovleva L.K. // Metodicheskiye ukazaniya k kursovoy rabote, k kursovomu i diplomnomu proyektam po gidroturbinam i obratimym gidromashinam «Raschet na prochnost' lopatki napravlyayushchego apparata», (Methodical instructions to the term paper, to the term and graduation projects on hydraulic turbines and reversible hydraulic machines "Strength calculation of the blades of the guide vanes"). - Khar'kov, KHGPU, 1994.
6. Mironov K.A. Yakovleva L.K., Gulakhmadov A.A. // Vybor tipa gidroturbiny pri proyektirovanii vysokonapornoy GES (The choice of the type of hydraulic turbine when designing a high-pressure hydro power plant). - Promislova gidravlika i pnevmatika. - Vinnitsya: VNAU. - №1 (43). - 2014. -p. 51-54.
7. Mironov K.A. Yakovleva L.K., Gulakhmadov A.A. // Uskoreniye protsedury proyektirovaniya i rascheta rabochikh koles radial'no-osevykh gidroturbin (Speeding up the design and calculation of impellers for radial-axis hydraulic turbines). Sbornik nauchnykh trudov SWorld. - Odessa - 2013. - Tom 6. - p. 23-27.
УДК 005.8:004.73
8. Potetenko O.V., Drankovskiy V.E. Krupa Ye.S., Vakhrusheva O.S., Gulakhmadov A.A. // Vikhrevaya struktura potoka i analiza razlichnykh matematicheskikh modeley potoka v kanalakh vysokonapornykh radial'no-osevykh gidroturbin RO 400, RO 500 i RO 600 (Vortex structure of the flow and analysis of various mathematical models of flow in the channels of high-pressure radial-axis hydraulic turbines RO 400, RO 500 and RO 600). Skhidno-Evropeys'kiy zhurnal peredovikh tekhnologiy. - Kharkiv, 2012. - № 3/8 (57). - p. 50-57.
9. Bryzgalov V.I. // Iz opyta sozdaniya i osvoyeniya Krasnoyarskoy i Sayano-Shushenskoy gidroelektrostantsiy (From the experience of creating and developing the Krasnoyarsk and Sayano-Shushenskaya hydro power plants). Proizvodstvennoye izdaniye. — Krasnoyarsk: Sibirskiy ID «Surikov», 1999. — 560 p.
10. Yefimenko A.I., Rubinshteyn G.L. // Vodosbrosnyye sooruzheniya Sayano-Shushenskoy GES (Spillway facilities of the Sayano-Shushenskaya hydro power plant). SPb.: Izdatel'stvo OAO «VNIIG im. B. Ye. Vedeneyeva», 2008. — 511 p.
11. Dvoretskaya M.I., Zhdanova A.P., Lushnikov O.G., Sliva I.V. // Vozobnovlyayemaya energiya. Gidroelektrostantsii Rossii (Renewable energy. Hydro Power Plants of Russia SPb.: Izdatel'stvo Sankt-Peterburgskogo politekhnicheskogo universiteta Petra Velikogo, 2018. — 224 p.
Marshak O. I.
Head of Department of Information Technologies of Shipbuilders College of the AdmiralMakarov National University of Shipbuilding
SYSTEMATIC RESEARCH FOR ENHANCEMENT OF COMPETITIVENESS OF SHIPBUILDERS COLLEGE IN THE SPACE OF PROFESSIONAL PRE-HIGHER EDUCATION.
Маршак Олена 1ллгвна
3aeidyea4 вiддiлення тформацшних технологш
Коледжу корабелiв Нацюнального yuieepcumemy кораблебудування iMeui adMipana Макарова
СИСТЕМНЕ ДОСЛ1ДЖЕННЯ ЩДВИЩЕННЯ КОНКУРЕНТОСПРОМОЖНОСТ1 КОЛЕДЖУ КОРАБЕЛ1В У ПРОСТОР1 ФАХОВО1 ПЕРЕДВИЩО1 ОСВ1ТИ
Summary. This work is dedicated to solving an important scientific and practical task of analytical forecasting of processes of professional pre-higher education development. For this matter a system-analytical research has been made. This research includes formulation of criteria for improving efficiency of innovational and research activities of Information technology department of Shipbuilders Collage of National university of shipbuilding and building cause and effect diagrams for the matter of situation forecast.
Анотащя. Робота присвячена виршенню важливо! науково-прикладно! 3aAa4i аналиичного прогнозування процеав розвитку фахово! передвищо! освгги, для чого виконано системно -аналтгачне дослвдження, яке включатиме у себе формулювання критерив тдвищення ефективносп шновацшно! та дослвдницько! дiяльностi вщдшення тформацшних технологш Коледжу корaбелiв Нацюнального ушверситету кораблебудування iменi aдмiрaлa Макарова (м. Микола!в) з побудовою для прогнозу ситуаци причинно-нaслiдкових дiaгрaм.
