CC BY
10
ÄUNIVERSUM:
№ 2 (95)_ТЕХНИЧЕСКИЕ НАУКИ_Февраль. 2022 г.
DOI - 10.32743/UniTech.2022.95.2.13126
МОДЕЛИРОВАНИЯ МАКСИМАЛЬНОЙ СКОРОСТИ ПОТОКА ВЫЗЫВАЮЩЕЙ КАВИТАЦИЮ И РЕЗКОЙ ПЕРЕСТРОЙКИ ПОТОКА
Худайкулов Савет Ишанкулович
д-р техн. наук, профессор, Научно-исследовательский институт ирригации и водных проблем,
Республика Узбекистан, г. Ташкент E-mail: S. I.Xudaykulov@mail. ru
Муминов Ойбек Алишер угли
докторант,
Ферганский политехнический институт, Узбекистан, г. Фергана
SIMULATION OF THE MAXIMUM FLOW RATE CAUSING CAVITATION AND ABRUPT RESTRUCTURING FLOW
Savet Khudaykulov
Doctor of Technical Sciences, Professor, Research Institute of Irrigation and Water Problems,
Uzbekistan, Tashkent
Muminov Oybek Alisher ugli
Doctoral student, Fergana Polytechnic Institute, Uzbekistan, Fergana
АННОТАЦИЯ
В статье исследуется число кавитации, которое становится равным коэффициенту давления, в обтекающем потоке в том месте, где достигается максимальная скорость, возникает кавитация, которая может привести к резкой перестройке всего течения жидкости. Если число кавитации меньше коэффициента давления, то безразмерное число кавитации приобретает существенное значение как определяющий безразмерный параметр. Приводится величина коэффициента давления в точках поверхности тела.
ABSTRACT
The article examines the number of cavitation that becomes equal to the pressure coefficient, in the flowing flow at the point where the maximum velocity is reached, cavitation occurs, which can lead to a sharp restructuring of the entire fluid flow. If the cavitation number is less than the pressure coefficient, then the dimensionless cavitation number becomes essential as a defining dimensionless parameter. The value of the pressure coefficient at the points of the body surface is given.
Ключевые слова: вибрация, резонанс, поперечные волны, волны возмущения, зона кавитации, начальное сечение поворота трубы, ось быстротока, радиус кривизны при повороте.
Keywords: vibration, resonance, transverse waves, disturbance waves, cavitation zone, initial section of pipe rotation, velocity axis, radius of curvature during rotation.
Введение. Известно, что жидкости, встречающиеся в природе и применяемые в технике, содержат взвешенные твердые частицы и растворенные газы. В тех местах потока, где давление падает, происходит нарушение сплошности течения и образуется область, заполненная пузырьками, внутри которых находятся пары жидкости или газ, выделившийся из раствора. Это явление называется кавитацией. Начальную стадию кавитации можно трактовать как
явление закипания жидкости при понижении давления. При дальнейшем понижении давления мелкие пузырьки объединяются и в потоке возникают большие полости - каверны, заполненные выделившимися из жидкости газами и парами жидкости. (рис. 1)
Постановка задачи. Величину давления
можно рассматривать как физическую характеристику, которая не влияет на движение жидкости при
р ^ ра. При р = ра в жидкости может возникать
Библиографическое описание: Худайкулов С.И., Муминов О.А. МОДЕЛИРОВАНИЯ МАКСИМАЛЬНОЙ СКОРОСТИ ПОТОКА ВЫЗЫВАЮЩЕЙ КАВИТАЦИЮ И РЕЗКОЙ ПЕРЕСТРОЙКИ ПОТОКА // Universum: технические науки : электрон. научн. журн. 2022. 2(95). URL: https://7universum. com/ru/tech/archive/item/13126
№ 2 (95)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
февраль, 2022 г.
кавитация, оказывающая существенное влияние на законы движения жидкости. Кавитация может возникнуть, например, вблизи минимального сечения в трубке с пережатием, в поршневом насосе, когда давление за поднимающимся поршнем стремится к нулю, а также при обтекании различных тел потоком жидкости.
Особенно такие движения потока соответствуют движению дисперсной смеси в всасывающих насосах для поднятия воды на опеределенную высоту (рис.1) в наливные Каркидонского водохранилища. Затем исследуются выше приведенными условиями с помощью уравнения движения для смеси жидкости, который имеет вид [3]:
K =
8/24
1
(1)
R2
Ol
fi Ol
Здесь /о - вязкость несущей жидкости; Я - радиус трубы. Если исследуемая смесь жидкости состоит из жидких и твердых фаз, совершающих взаимопроникающие движения, то вязкость смеси отличается от вязкости жидких (несущей) фаз и зависит от объёмного содержания твердых фаз и физических свойств взвеси смеси и т.д.
Расчетные формулы для коэффициента взаимодействия, приводимые в работе [2,4], полученные на основе формулы Эйнштейна, выглядят следующим образом:
K =
&fi№o
1
(2)
R
Qcm
Ol
fi (1 + 2,5 f2)
Для барометрического давления В на отметке установки рабочего колеса турбины получим
1
^ = * - К -у
з2
см 2
см 5
2g 2g
2 2 2 2 _ h + wcmx _ WCM 2 + UCM2 _ исмх _ h
w,2 _5 ^ cmw, x_2
2 g
(3)
где Их - теоретическая высота отсасывания;
g2 Q2
зсм 2 зсм 5_h
2g 2g
*'сми> 2-5 — динамическое разрежение,
обусловленное формой и размерами отсасывающей трубы;
2 2 2 2 w _ Wo + и о _ и
смх см 2 см 2 смх
2 g
1CMW,x_2 - динами-
ческое разрежение, обусловленное размерами и формой лопастей рабочего колеса турбины. Нужно отметить, что многочлены зависят, кроме указанного, еще и от режима работы турбины. Уравнение (3) можно записать иначе:
Px Pn
У У
H ~ H
В _ hx _
v у.
H
з
см 2
Lg 2g
з2
-см5 _ h +
hw 2_5 +
2 2 2 2 wcmx _ wcm 2 + ucm 2 _ ucmx _ fo
(4)
2 g
l о
cmw x_2
Из уравнения (4) находим условие возникновения кавитации: Кавитация в турбине появляется,
если абсолютное давление Е^ в точке х будет равно
7
давлению парообразования, т.е. когда правая часть уравнения (4) обратится в нуль:
1 H
1 H
Г ^2
з
пт О
Pn
В _ hx
V У j
з
2
ПЛ1 'S
2g 2g
hw 2_5 +
2 2 2 2 W _ W о + и о _ и
rr пт x rr nm 9
2 g
h
'cmw x_2
^ 0
В первой скобке заключены величины, характерные для данной установки: местное барометрическое давление В ; высота расположения рабочего колеса над уровнем нижнего бьефа,
называемая высотой отсасывания Ьх; давление парообразования — ; напор Н в данной установке.
?
У
№ 2 (95)
A, UNI
те)
UNIVERSUM:
технические науки
февраль, 2022 г.
Рисунок 1. Влияние кинетической составляющей давления на распределение избыточного давления на повороте
а - прямой водовод, давление в сечение гидростатическое; б - изогнутый водовод, за счёт кинетической составляющей, избыточное давление увеличивается на внешной стенке и уменьшается на внутренней.
Выражение в правой части уравнения - кавита-ционный коэффициент установки:
а
в - h - En.
H
(5)
Во второй скобке уравнения (4) содержатся квадраты скоростей в абсолютном, переносном и относительном движении и потеря энергии между точками х и 5, отнесенная к напору.
Потеря напора Н^, может быть выражена через
2
квадрат скорости смеси жидкости Зсм в любой точке
потока и соответствующий этой скорости коэффициент сопротивления
hwx-5 =
2 g
(6)
Скорости воды в потоке смеси пропорциональны корню квадратному из величины используемого напора ]ГH , где ]г - гидравлический КПД турбины. Поэтому все скорости могут быть представлены в виде скоростного коэффициента, соответствующего данной скорости в рассматриваемой точке. Скоростной напор, соответствующий в той же точке, имеет вид
З2
= K зЛг H.
2 g
(7)
В условиях кинематических подобных (изогональных) режимах работы для турбин данной
серии все скоростные коэффициенты К9 и коэффициенты сопротивления д имеют постоянные значения, С* для турбин данной серии и изогональных режимов остается неизменным. Таким образом, кавитационный коэффициент турбины С представляет собой отношение наибольшего динамического разрежения на лопасти рабочего колеса к используемому напору ]гH.
Если при каком-то режиме рх = рп, то динамическое разрежение имеет максимальное значение и возникает кавитация. При этом кавитационный коэффициент турбины для этого режима будет иметь критическую величину
а* = B - - h
а
У Л гH Л г
(8)
Из этого выражения можно получить условие для максимально допустимой высоты отсасывания Нх,
при которой обеспечивается бескавитационная работа турбины
H = К < в - ^ - а* л гH
У
(9)
Нетрудно также расчетным путем установить точку х, в которой можно ожидать наибольший вакуум. Поэтому в практике гидромашиностроения кавитационный коэффициент С определяется путем испытания моделей гидротурбин на кавитационном стенде.
У
№ 2 (95)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
февраль, 2022 г.
Приведем таблицу сопоставления значений кинетической энергии в точках 2 и 5 со значениями кинетической энергии, вычисленной по формуле (3), где с учетом взаимодействий фаз многофазного дисперсного потока при /х = 0,4 ; /2 = 0,6 таких значениях [1,2,3,5].
Установлено, при наличии отсасывающей трубы под рабочим колесом создается дополнительное понижение давления, пропорциональное удельной
кинетической энергии
з2
см 2
2 gH
оставшейся в потоке
на выходе из рабочего колеса. Величина этой энергии зависит от состава смеси, системы турбины, ее коэффициента быстроходности и пропускной способности, характеризуемой значением приведенного расхода Q'l, из рабочего колеса, выраженные в
процентах от рабочего напора для турбин различной быстроходности. [5,6,8] Ориентировочные значения удельной кинетической энергии на выходе приведены в табл. 1.
Таблица 1.
Сопоставление вычисленных значений удельной кинетической энергии на выходе водовыпуского сооружения по формуле (3) с данными водохранилища
Система турбин H, м Ol, л / с n1, об / мин ns Q2 2 100% 2 gH По формуле (4) 2 100% 2 gH Q2 з 100% 2 gH По формуле (7) 32 см 5 100% 2 gH
Поворотно- 5.0 2250 180 930 53.0 52.7 4.0 3.8
лопастные 10.0 2010 170 830 43.0 42.8 2.8 2.5
15.0 1850 160 755 36.0 34.8 2.3 2.2
20.0 1750 150 695 32.0 2.1 1.9 1.8
30.0 1430 140 580 24.0 20.8 1.4 1.3
40.0 1240 130 500 20.0 19.2 1.1 0.9
50.0 1110 120 435 18.0 17.6 0.9 0.8
60.0 1040 112 395 17.0 15.6 0.8 0.75
70.0 940 107 355 16.0 15.6 0.7 0.7
80.0 830 105 350 14.0 14.0 0.65 0.65
Радиально- 30.0 1400 78 320 11.0 10.9 1.5 1.4
осевые 45.0 1370 77 315 10.5 10.0 1.4 1.2
75.0 1250 75 290 8.5 8.0 1.2 1.1
115.0 1030 72 255 6.0 5.8 0.8 0.7
170.0 650 69 195 2.5 2.3 0.4 0.3
230.0 420 67 150 2.0 1.8 0.2 0.14
310.0 280 65 120 1.8 1.6 0.1 0.8
00.0 200 62 95 1.6 1.4 0.05 0.04
500.0 150 60 80 1.5 1.2 0.02 0.01
Сопоставление значений кинетической энергии в точках 2 и 5, которые приведены в табл. 1.2, показывает, если в отсасывающих трубах тихоходных турбин удельная кинетическая энергия в потоке за рабочим колесом составляет 1.5-3% от рабочего напора, то по мере увеличения быстроходности турбины эта цифра увеличивается, доходя у быстроходных поворотно лопастных турбин до 50% и более от всей энергии. Но из-за наличия внутренних взаимодействий фаз возрастают внутренние сопротивления, которые отрицательно влияют на кинетическую энергию потока, увеличивая тепловую энергию. Взаимодействие фаз дисперсной смеси уменьшает высоту отсасывания, и часть выделяемой энергии преобразовывается в тепловую энергию [7, 9].
п 3
В ряде случаев отношение — определяется кия,
нематическими условиями задачи, в частности, ниже увидим, что так обстоит дело при непрерывном потенциальном обтекании тел неограниченным потоком идеальной несжимаемой жидкости. В этом случае максимальная скорость Змах, достигается на границе текущей жидкости, т. е. на поверхности тела и
отношение
3
3
зависит только от геометрических
свойств поверхности тела и его ориентации относительно скорости набегающего потока.
При обтекании тел поступательным потоком безразмерные характеристики поля скоростей в иде-
№ 2 (95)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
февраль, 2022 г.
альной несжимаемой жидкости определяется системой безразмерных параметров х, у, £ , а, р, где й
й й й
-характерный размер тела, а, р-углы, задающие ориентацию тела относительно скорости набегающего потока. Безразмерное отношение 9 - не зависит от
9.
скорости, плотности и давления в набегающем потоке и получается постоянным при фиксированных
X У £
безразмерных координатах — ,£_,_, а, р.
й й й
Максимальное значение 9мах , соответствует
вообще одной вполне определенной точке на поверхности тела. При учете сжимаемости в случае адиабатических движений совершенного газа получится:
К=f
к
a, ß, X, У, £,
d d d a„
Максимальной скорости 9 частиц жидкости в потоке соответствует минимальное давление р .
Величину 2(рГСТ - р)
рК
в точках поверхности тела
называют коэффициентом давления и обозначают через ср .
На основания формулы для коэффициента давления, соответствующего точке минимального давления,
можно написать:
CPmn
_2( p
ГСТ pmm
рК
К
(10)
Наступление кавитации определяется условием
2(PГСТ - pd )
рК2
X
(11)
Безразмерное число
X = 2( p ГСТ
Pd
рК
называется числом кавитации. Число кавитации определяется заданными условиями обтекания. Значения % зависят от давления в бесконечности через
ргсг, - которое зависит от глубины погружения тела в жидкость. При фиксированной разности ргсг - рй число кавитации % резко падает с увеличением скорости набегающего потока 9..
Выводы:
1. В тот момент, когда число кавитации становится равным коэффициенту давления, в обтекающем потоке в том месте, где достигается максимальная скорость, возникает кавитация, которая может привести к резкой перестройке всего течения жидкости. Если число кавитации меньше коэффициента давления, то безразмерное число кавитации приобретает существенное значение как определяющий безразмерный параметр.
2. В этом случае число кавитации необходимо вводить наряду с числом Рейнольдса и числом Фруда в качестве основного параметра, характеризующего гидродинамический поток, и основного критерия подобия при моделировании.
3. При движении в жидкости любого профиля при увеличении его скорости неизбежно наступление кавитации. Кавитация наступает тем позже, чем ближе к единице отношение максимальной скости к скорости потока в бесконечности, т.е. чем меньше профиль возмущает поток.
4. Кавитация может возникнуть не только при увеличении скорости данного профиля, но и при уменьшении давления в бесконечности, которое зависит от глубины погружения тела в жидкость. Очевидно, что с погружением на глубину, когда р растет, наступление кавитации затрудняется.
1
1
c
Список литературы:
1. Бегимов У.И. Худайкулов С.И. «Дисперс аралашмалар окимидаги кавитация ва пульсация жараёнларнинг шаклланиш конуниятини моделлаштириш ва сув омборларига куллаш». Монография. 146 б. Бухоро 2019й.
2. Худайкулов С.И., Нишонов Ф.Х. «Математические модели гидравлического удара в гидросооружениях и производственных комплексах» Ташкент - 2017. 146 с.
3. Хамидов А.А., Худайкулов С.И. Теория струй многофазной вязкой жидкости. "ФАН" 2003. 138 с.
4. Khudaykulov S.I, Yakhshibaev D.S, Usmonov A. Nishonov F.X. НЖ: Change in concentration of collector waters along the flow length taking intoaccount the difference in densities Asian Journal of Research, Japan, Osaka. № 3, (3) 2019, pp 39-43 Scientific Journal Impact factor (№24), IF=5.1
5. Khudaykulov S.I., Jovliev U.T. Nishonov F.X. Modeling of disperse mixture flow with the formation of a bubble zone bringing hydro system to vibration. Ж: Asian Journal of Research, Japan, Osaka. № 3, (3) 2019, pp 159 -162 Scientific Journal Impact factor (№24), IF=5.1.
6. УсмоноваН.А., ХудайкуловС.И., УсмановА.Х. Мухандислик коммуникация кувурларида кавитация ва пульсация хосилбулишини динамикасини моделлаштириш. Научно-технический журнал ФерПИ. 2019 . Том 23 . спец. вып. № 3.С. 74-79.
№ 2 (95)
UNIVERSUM:
ТЕХНИЧЕСКИЕ НАУКИ
февраль, 2022 г.
7. Худайкулов С.И., Н.А. Усмонова, А.А. Усманов Установление начальных интервалов и прогнозов возникновения пульсации и кавитации пульсирующего многофазного потока в трубопроводах каркидонского водохранилища Научно-технический журнал ФерПИ. 2019 . Том 23 . спец. вып. № 3. С. 183-187.
8. Худайкулов С.И., Негматуллоев З.Т., Бегимов У.И. Течение дисперсной смеси в трубе с наличием магнитного поля. Узбекистан ^ишлок; ва сув хужалиги журнали илмий иловаси. "AGRO-ILM» 2020 № 1. 86-89.
9. Худайкулов С.И., Усманов А.А., Усмонова Н.А. «Характеристика защитных оболочек водосброса Каркидонского водохранилища» O'zbekiston respublikasi madaniyat vazirligi o'zbekiston davlat san'at va madaniyat institutining farg'ona mintaqaviy filiali "oriental art and culture" ilmiy-metodik jurnali "o'zbekistonda ilm-fan va ta'lim"№2-sonli konferensiya to'plami (8-qism) 24 Aprel 2020 yil Qo'qon shahri 341-346с.
10. Ishankulovich, Khudaykulov Savet, Usmonova Nodirakhon Akramovna, and O'tbosarov Shuhratjon Rustamjon ugli. 2021. "Simulation of the Lift of Two Sequential Gate Valves of the Karkidon Reservoir". Middle European Scientific Bulletin 18.
11. Jovliev Oktam Temirovich, Mamlakat Khasanovna Khujakulova, Nodirakhon Akramovna Usmonova and Shuhratjon Rustamjonugli O'tbosarov. 2021. "Modeling the Theory of Liquid Motion Variable on the Way Flow". Middle European Scientific Bulletin 18 (November):455-61.
12. Abdulkhaev Z.E. (2021). Protection of Fergana City from Groundwater. Euro Afro Studies International Journal, 6,
13. Erkinjonovich A.Z., & Mamadaliyevich M.M. (2021, May). Water Consumption Control Calculation In Hydraulic Ram Device. In E-Conference Globe (pp. 119-122).
14. ugli Mo'minov O.A., Maqsudov R.I., & qizi Abdukhalilova S.B. (2021). Analysis of Convective Finns to Increase the Efficiency of Radiators used in Heating Systems. Middle European Scientific Bulletin, 18, 84-89.
15. Мадрахимов М.М., & Абдулхдев З.Э. (2019). Насос агрегатини ишга туширишда босимли сув узатгичлардаги утиш жараёнларини х,исоблаш усуллари. Фаргона Политехника Институти Илмий-Техника Журнали, 23(3),
16. Mo'minov O.A., & O'tbosarov Sh, R. Type of heating radiators, principles of operation and theoretical analysis of their technical and economic characteristics.
17. Mamadalievich M.M., & Erkinjonovich A.Z. Principles of Operation and Account of Hydraulic Taran. JournalNX, 1-4.
18. Абдукаримов Б.А., Муминов О.А., & Утбосаров Ш.Р. (2020). Оптимизация рабочих параметров плоского солнечного воздушного обогревателя. In Приоритетные направления инновационной деятельности в промышленно сти (pp. 8-11).
19. Madraximov M.M., Abdulkhaev Z.E., & Orzimatov J.T. (2021). Gidravlik taran qurilmasining gidravlik hisobi. Scientific progress, 2(7), 377-383.
20. Usarov M., Ayubov G., Usarov D., & Mamatisaev G. (2022). Spatial Vibrations of High-Rise Buildings Using a Plate Model. In Proceedings of MPCPE 2021 (pp. 403-418). Springer, Cham.
21. Abdulkhaev Z.E., Abdurazaqov A.M., & Sattorov A.M. (2021). Calculation of the Transition Processes in the Pressurized Water Pipes at the Start of the Pump Unit. JournalNX, 7(05), 285-291.
22. Abobakirovich A.B., Mo'Minov Oybek Alisher, O.G., & O'G'Li, S.M. A. (2019). Calculation of the thermal performance of a flat solar air heater. Достижения науки и образования, (12 (53)), 9-11.
23. Абдулхаев З.Э., & Сатторов А.М. (2020). Central pump case adjustment by changing the rotation frequency. Актуальные научные исследования в современном мире, (6-1), 20-25.
24. Mo'minov O.A., Abdukarimov, B.A., & O'tbosarov, S.R. (2021). Improving support for the process of the thermal convection process by installing reflective panels in existing radiators in places and theoretical analysis. In Наука и инновации в строительстве (pp. 47-50).
25. Mirsaidov M., Usarov M., & Mamatisaev G. (2021). Calculation methods for plate and beam elements of box-type structure of building. In E3S Web of Conferences (Vol. 264). EDP Sciences.
26. Koraboevich U.M., & Ilhomidinovich M.G. (2021, June). Calculation of the free vibrations of the boxed structure of large-panel buildings. In " online-conferences" platform (pp. 170-173).
70-81.
56-60.
