УДК 662.9
ИССЛЕДОВАНИЕ ТЕЧЕНИЯ И ЭНЕРГООБМЕНА ЗАКРУЧЕННЫХ ПОТОКОВ ЖИДКОСТИ В ГИДРАВЛИЧЕСКОМ ГЕНЕРАТОРЕ ТЕПЛА
© 2011 В. Н. Белозерцев
Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет)
Рассматриваются вопросы исследования и энергообмена закрученных потоков жидкости в вихревом гидравлическом генераторе тепла. Получены выражения для расчёта распределения скоростей и давления по его длине и радиусу. Подтверждена картина формирования свободного и вынужденного вихрей. Получены уравнения для расчёта подогрева жидкости.
Гидравлический генератор тепла, закрученный поток жидкости, уравнение движения, уравнение неразрывности, распределение, давление, скорость, температура.
Традиционными теплогенераторами, нашедшими широкое применение в технике и быту, являются преобразователи энергии, использующие теплоту сгорания различного рода топлив (жидкого, твёрдого и газообразного). Их недостатками являются загрязнение окружающей среды продуктами сгорания, значительные финансовые затраты при транспортировке и хранении энергоносителей.
Создание, разработка и исследование теплогенераторов, реализующих в своей работе иные принципы преобразования различных форм энергии в тепловую, и экологически чистых представляет значительный интерес. Такими теплогенераторами являются гидравлические вихревые теплогенераторы (рис.1.) В них электрическая энергия (привода насоса) преобразуется в кинетическую энергию поступательного и вращательного движения потока жидкости. Жидкость, приводимая в движение насосом, поступает через входное устройство 1 в завихрителе 2. Благодаря тангенциальному вводу 1 и профилю завихрителя (спираль Архимеда) жидкость закручивается, приобретая вихревой, спиралевидный характер движения. Выйдя из корпуса 3, через развихритель 4 и выходное отверстие 5 она поступает в
насос 6 и из него вновь подается во входное устройство гидравлического теплогенератора. Таким образом, в процессе работы жидкость движется по закольцованному контуру: теплогенера-тор-насос-теплогенератор. Характер течения жидкости, распределение скоростей и давлений внутри устройства будут определять интенсивность её подогрева. Поэтому необходимо знать особенности распределения скоростей и давлений в рабочей зоне теплогенератора, уметь их рассчитывать при заданных условиях на входном устройстве.
Вопросам течения закрученных потоков жидкости посвящено много работ. Как правило, в них рассматривается изотермическое течение жидкости. В вихревом гидравлическом генераторе тепла жидкость подогревается, изменяются её теплофизические характеристики, что влияет и на распределение её параметров внутри устройства. Но доля подогрева жидкости за один проход через теплогенератор будет составлять сотые доли градуса, поэтому зависимости распределения скоростей и давлений в рабочей зоне энергообмена, полученные для изотермической модели, будут здесь также справедливы.
Рис. 1. Общий вид и гидравлическая схема внутреннего контура теплогенератора: 1- входное устройство; 2-завихритель; 3-корпус; 4-развихритель; 5-выходное устройство; 6-насос; 7-электродвигатель; 8-бак
Подогрев жидкости может быть определён, если воспользоваться уравнением Бернулли для вязкой несжимаемой жидкости, вращающейся относительно некоторой оси [1]:
2 2 V -V
------2 + К
2 с
(1)
Откуда удельная энергия, затрачиваемая на преодоление сил вязкого трения, будет равна:
2 2 2 2 Ь _ Р1 ~ Р2 + и1 ~ и2 + и2 ~ и1 (2)
с р 2 2 ’
здесь р1, р2, и1, и2, и1 , и2 - среднеинтегральные значения статического давления, окружных и осевых скоростей соответственно во входном и выходном сечениях теплогенератора. Полагая, что вся энергия, затрачиваемая на преодоление сил вязкостного трения, превращается в теплоту и идёт на подогрев жидкости, можно найти температуру воды на выходе из теплогенератора t2 из выражения
К = С ( ¿2 - ¿1 ) •
Откуда следует: К
і — і +——
2 91^ с
(3)
(4)
где С = С(^)-теплоёмкость воды. Восполь-
зоваться уравнениями (2), (3), (4) можно,
если известны составляющие правой час-
ти уравнения (2) в выходном сечении. Их
можно наити из уравнения неразрывности и движения вязкой несжимаемой жидкости, в основу которого положен закон внутреннего трения Ньютона в применении к жидкой частице в цилиндрической системе координат % в, z [1]:
V—
Эи
1 Э Р р Э г
Эи +1 Эи и
ЭГ г Эг г
(5)
(6)
Эи 1 ЭР Г Э2и 1 Эи
V —-------+УТ —2 +---------
Эи р Э <Т г Эг
Э ( ) Э ( )
—(щ+—( ги)—0,
Эг Э
где, согласно рис.2, осреднённые компоненты скорости: w (по радиусу), и (по углу в), и (по оси I). Соотношения между ними для закрученного потока задаются в рамках одной из полуэмпирических теорий турбулентности [2].
Рис. 2 Составляющие скорости и характерные радиусы в поперечном сечении закрученного потока в теплогенераторе
2
и
г
В работе [3] получены распределения окружной и (г, г), осевой V (г, г)
скоростей, циркуляции
Гн (г)
г и статиче-
ского давления по радиусу и длине трубы Р(гг) :
ч(г, г) —
Гн (г)
1 - ехр
..2 Л
- Яет
4цт г
(7)
1 - ехр
- Яе
(Я - г )
4ц г
2
Г — ги — Гн ()
1 - ехр
/ 2 Л
г
4ц Т г
(8)
_ г
V (г, г)—V+Яо^-
ҐЯ л2
г
V У
ехр
''та 2 Ч
Яе , г
4 Яг
1 - ехр
2
Яе т г
4 Яг
+ Яо Яе т Гн — х
Яе т г 4 Яг
2
г
Р(Г, 2) — Ря + р |
и(г, г)2
Яет г 2Яг
Я
г
2
(9)
7-І
— Р я +
+р
х
х
х
1 - ехр
2
- Яет
1 - ехр
Яе
4ц т2
(я - г) 4ц г
(10)
х
/-І
2
-2
р(г)=рн=соті; Яо —
Гн
V Я2
мых функцией ГН (г), можно получать и различные функции изменения окружных скоростей по длине теплогенератора. Например: при потенциальном свободном вихре на входе имеем ГН(г)=соп$Х=ГН. Тогда
Г
и —
1 - ехр
При г <<
4ц Т г 4ц тг
(11)
в приосевой зоне
где Цт — Я Яе г , ц — Я Яе - турбулентная и кинематическая вязкости; Я - радиус трубы; Яе - критерий Рейнольдса; Гн(г)— гин(г) - циркуляция на входе в трубу; ря -давление в полом вихревом жгуте закрученного потока; г0 - радиус цилиндрического разрыва сплошности вблизи оси вращения, при г<г0 давление
жидкость вращается как твёрдое тело с
Г /
угловой скоростью /4ц г. При
г >> Л14ц тг - ближе к стенкам теплогенератора движение безвихревое, каким оно было в начальном сечении. Уравнение (7) теоретически обосновывает составной «вихрь Ранка». Анализ решения уравнений (7) и (8) показал, что при любом промежуточном состоянии жидкости на входе, какими бы законами ни задавался входящий в теплогенератор закрученный поток, через определённое расстояние от начального створа вследствие вязкой диффузии и диссипации энергии формируется течение, внутренне присущее закрученному потоку.
В имитационной математической модели было проведено исследование течения жидкости при двух различных условиях на входе в теплогенератор. В первом случае величина циркуляции на входе задавалась как для потенциального течения ГН=у вхтвх=соп$Х и как твёрдое тело
Гн —-
Я
-г•
число Росс-
би; Е1(...) - интегрально-показательная функция.
Анализ (7) и (8) показывает, что при различных начальных условиях, задавае-
Для решения этой задачи для двух различных условий была разработана программа расчёта в среде МаШсаё. В этой программе были проведены расчёты для гидравлического теплогенератора с заданными: площадью сопла на входе в улитку 26-16-10-6 м2, радиусом трубы камеры 26-16-10-6 м2, длиной трубы 0,5 м и радиусом на входе гвх = 0,058 (расстояние от оси входного устройства до оси
корпуса), скоростью на входе ьвх — 2,4
м
2
х
г
Г
х
х
х
х
г
3
г
-2
г
вх
Разработаны дополнительные программные функции, позволяющие автоматически переопределять табличные данные р = р ^), С = С ^), V = V ^) (плот-
ности, теплоёмкости и кинематической вязкости) и реализовывать циклически процедуру расчёта осреднённых по площади сечения значений \(г, г); и(г, г); р(г, г) (и изменение температуры за несколько проходов жидкости через теплогенератор. Дополнительные процедуры помещены в отдельный файл МаШсаё, расчётные документы используют вспомогательные функции по ссылке. Результаты расчёта представлены на рис. 3. Исследуемая имитационная модель теплогенератора показала, что задание циркуляции на входе существенно влияет на физическую картину течения и практические расчёты. Так, потенциальное движение на входе, задаваемое ГН=сопв1 оказывается практически неприменимо для описания структуры потока в двух зонах: вблизи оси вращения и твёрдых стенок трубы, в которых существенно воздействие вязкости. Более близкие результаты даёт имитационная математическая мо-
дель при задании циркуляции
Г н ( г)
генератор - насос - теплогенератор. Уточнить подогрев для каждого конкретного циркуляционного контура теплогенератора можно с помощью выражения
А кс —
(12)
і—1
здесь
А Нс
доля удельной энергии, по-
шедшей на нагрев всей жидкости; ^ Ц -
і—1
суммарная удельная энергия, пошедшая на нагрев воды в контуре труб и потери в окружающую среду.
Предположим, что длина присоединительных труб равна длине теплогенератора (в действительности она в несколько раз больше). Величина относительного теплового потока, отводимого от внутренней поверхности теплогенератора и присоединительных труб за счёт теплопередачи, будет составлять п (4 - ^ )
Ч, =-
1
1п
d0
1
(13)
а/ • d1
г по
второму случаю, когда рассматривается вращение жидкости во входном сечении как вращение твёрдого тела с постоянной угловой скоростью. Результаты, представленные на рис. 3, полностью подтверждают наличие свободного вихря у периферии трубы и вынужденного вихря, вращающегося как твёрдое тело. Рассчитанные значения температуры (рис.3, г) получены в математической имитационной модели для адиабатно-изолированной жидкости на внешней поверхности теплогенератора, а также без учёта теплоты, пошедшей на подогрев корпуса теплогенератора, труб, обеспечивающих круговое движение по рабочему гидравлическому контуру теплогенератора, подогрева всей жидкости, находящейся в системе тепло-
В этом выражении ¿1, - темпера-
тура воды на входе в теплогенератор и температура окружающей среды на момент времени %; d1, d2 - внутренний, внешний диаметры; 1 - коэффициент теплопроводности стенки; аг - коэффициент теплоотдачи внешнего теплоносителя.
Тогда тепловой поток О — 2 • Ч, • I, (14)
а доля энергии, пошедшей на нагрев воды,
(15)
АКс — Кс
т.Т
Подогрев воды в теплогенераторе с учётом отвода теплоты в стенку д / = ДС. (16)
с
— г=ь ;
Относительный радиус гЖ, %
х о ; 2=112 ; — / 1. /4
— 2=0
м/с ■ 2=1/21,;
б
г=ь
1-10 ^ 8-Ю5
е" 6-ю5
и <■
§ 4-10 й _
« 2-Ю5
20 40 60 80 100
Относительный радиус г®, %
— г=о
г 1- 2Т.; 7-7.
Рис. 3 Распределение по радиусу и в различных сечениях теплогенератора: а - окружной скорости; б - осевой скорости; в - статического давления; г- температуры жидкости в выходном сечении
а
в
г
Масса жидкости в теплогенераторе и трубопроводах будет приблизительно равна:
тт = 2 р Ут , (17)
здесь р - плотность воды; Ут - внутренний
объём теплогенератора. Тогда при смешении подогретой жидкости из теплогенератора и жидкости, находящейся в трубах,
температура смешения
+ т(г~1
t с
тт 12 + т212
тт + т2
(18)
здесь тт , т2 -массы жидкостей в теплогенераторе и трубопроводах; ^, ^-1 - температуры жидкостей на выходе из теплогенератора и температура жидкости в трубопроводах, г - номер прохода жидкости через теплогенератор. Тогда подогрев
жидкости с учётом теплообмена и присоединённой массы жидкости
д^=(Км -1;). (19)
Из (18) видно, что температура жидкости после смешения будет существенно зависеть от суммарной присоединённой массы жидкости, находящейся в закольцованном контуре теплогенератор-насос-теплогенератор. В табл.1 представлены результаты расчёта температуры смеси и температурного напора Дt воды для различных присоединённых масс при первом проходе жидкости. Как видно из расчётов, величина будет существенно ниже при учёте теплообмена с окружающей средой. При наличии данных о присоединённой массе всех трубопроводов и насоса это значение снизится еще более.
Таблица 1. Результаты расчёта подогрева жидкости
t1 4 At t 4i Ahc At t °C CM > At”
°c °c °C с Вт м Дж кг °C m2 m2
mT 2mT 3mT mT 2mT 3mT
20 20.12 0.12 0.21 274 440.35 0.105 20.052 20.035 20.026 0.052 0.035 0.026
В результате проведённых исследований можно сделать следующие выводы:
1. Разработанная имитационная математическая модель подтверждает физическую картину течения жидкости внутри теплогенератора: наличие двух вихревых потоков, периферийного - свободного и приосевого - вынужденного. Оба потока вращаются в одну сторону, но линейно движутся в разных направлениях.
2. Математическая модель позволяет рассчитывать эпюры окружных, осевых скоростей и статического давления по радиусу и длине теплогенератора, а также подогрев жидкости для адиабатной модели теплогенератора с учётом теплообмена
в окружающую среду, данных по присоединённой массе воды в трубопроводах.
Библиографический список
1. Абрамович, Г. Н. Прикладная газовая динамика [Текст] / Г. Н. Абрамович.- М.: Гос. теор.-техн. изд-во, 1953.600 с.
2. Аксенов, Ю. Н. Экспериментальное исследование движения винтового потока в водосбросных сооружениях [Текст]: автореф. дис.. канд. техн. наук./ Ю. Н. Аксенов. - Харьков, 1969.-16 с.
3. Волшаник, В. В. Закрученные потоки в гидротехнических сооружениях. [Текст] / В. В. Волшаник, А. Л. Зуйкова, А. П. Мордасов. - М.: Энергоатомиздат, 1990.-280 с.
RESEARCH OF CURRENT AND POWER INTERCHANGE OF SWIRLED FLUID STREAMS IN A HYDRAULIC HEAT GENTRATOR
© 2011 V. N. Belozertcev
Samara State Aerospace University named after academician S. P. Korolyov
(National Research University)
The paper deals with the analyses of power interchange of swirled fluid flows in a vertical hydraulic heat generator. Expressions are derived for the calculation of velocity and pressure distribution over its length and radius. The pattern of free and forced vortex formation is confirmed. Equations are obtained to calculate fluid heating.
Hydraulic heat generator, swirled flow of fluid, equation of motion, continuity equation, distribution, pressure, velocity, temperature.
Информация об авторе Белозерцев Виктор Николаевич, кандидат технических наук, доцент кафедры теплотехники и тепловых двигателей, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). E-mail: teplotex ssau@bk.ru. Область научных интересов: рабочие процессы тепловых двигателей и энергетических установок.
Belozertsev Viktor Nicolayevich, Ph.D., associate professor, department of heat engineering and heat engines, Samara State Aerospace University named after academician S.P. Korolyov (National Research University). E-mail: teplotex ssau@bk.ru. Area of research: working processes of heat engines and power plants.