Preview

Вестник Концерна ВКО «Алмаз – Антей»

Расширенный поиск

Теоретические основы рабочего процесса вихревой трубы в режиме вакуум-насоса (вихревого эжектора)

https://doi.org/10.38013/2542-0542-2022-2-63-73

Полный текст:

Аннотация

Составлена уточненная физическая модель работы вихревой трубы в режиме вакуум-насоса (вихревого эжектора). Учитывает количество энергии, передаваемой от эжектирующего к эжектируемому потоку газа. Составлена замкнутая система уравнений, описывающая рабочий процесс вихревого эжектора. На базе этой системы уравнений составлены две основные методики расчета вихревого эжектора: методика расчета основных геометрических размеров и методика расчета характеристик при известных геометрических размерах. Показано влияние вязкости на энергообмен в вихревом эжекторе. Совпадение расчетных и экспериментальных данных удовлетворительное. 

Для цитирования:


Кузнецов В.И., Макаров В.В., Шандер А.Ю. Теоретические основы рабочего процесса вихревой трубы в режиме вакуум-насоса (вихревого эжектора). Вестник Концерна ВКО «Алмаз – Антей». 2022;(2):63-73. https://doi.org/10.38013/2542-0542-2022-2-63-73

For citation:


Kuznetsov V.I., Makarov V.V., Shander A.Yu. Fundamentals of the working process of a vortex tube operating in the vacuum pump (vortex ejector) mode. Journal of «Almaz – Antey» Air and Space Defence Corporation. 2022;(2):63-73. https://doi.org/10.38013/2542-0542-2022-2-63-73

ВВЕДЕНИЕ

Газовый эжектор или вихревой компрессор (вихревой эжектор) – простейшие и распространенное газодинамическое устройство, применяемое в разнообразных отраслях промышленности, в частности в авиа- и ракетостроении, вакуумной технике и различных экспериментальных аэродинамических установках [1][2][8].

Основным достоинством эжектора как компрессора является отсутствие движущихся и трущихся деталей, что существенно при работе с горячими либо агрессивными средами.

ПОСТАНОВКА ЗАДАЧИ

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

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

МАТЕРИАЛЫ И МЕТОДЫ ИССЛЕДОВАНИЯ

Существуют несколько физико-математических моделей, объясняющих работу эжектора [1][2][9]. Основной недостаток данных физических моделей состоит в том, что по ним нельзя составить замкнутую математическую модель. Для получения единственного решения вводятся эмпирические формулы на основании экспериментальных исследований. Как правило, эти формулы справедливы в той области, в которой были проведены исследования.

В ряде работ эжекторы называют компрессорами без движущихся частей, но ни в одной работе нет уравнений, показывающих механизм передачи энергии от эжектирующего газа к эжектируемому [5][4][10].

На основании теоретических и экспериментальных исследований других авторов и собственных представлена физическая модель вихревого эжектора, в которой учитывается обмен работой и теплотой между эжектирующим и эжектируемым газами. Введение в систему уравнений, описывающих рабочий процесс вихревого эжектора, обмена работой и теплотой позволяет замкнуть систему уравнений [8][11][17].

На основании решения замкнутой системы уравнений составлены две методики:

  • методика расчета оптимальных геометрических размеров вихревого эжектора (приложение 1);
  • методика расчета характеристик вихревого эжектора при известных геометрических размерах (приложение 2).

ФИЗИЧЕСКАЯ МОДЕЛЬ РАБОТЫ ВИХРЕВОГО ЭЖЕКТОРА

Вихревой эжектор – устройство, в котором полное давление эжектируемого газа увеличивается за счет энергообмена с эжектирующим газом. Эжектирующий газ имеет более высокое полное давление и совершает работу над эжектируемым газом силами вязкости, что приводит к повышению давления эжектируемого газа. В результате взаимодействия потоков на выходе из вихревого эжектора (ВЭ) образуется смесь, имеющая полное давление выше, чем давление эжектируемого газа, но ниже полного давления эжектирующего газа [1].

Предлагается следующая физическая модель рабочего процесса вихревого эжектора.

Эжектирующий (высоконапорный) газ истекает из сопла в камеру смешения (рис. 1).

Рис. 1. Схема вихревого эжектора: 1 – канал для входа эжектирующего газа; 2 – канал для входа эжектируемого газа; 3 – камера энергообмена; 4 – канал для выхода смеси газа

В начальный момент работы ВЭ эжектирующий газ начинает течение из сопла, а эжектируемый газ находится в состоянии покоя. Под действием центробежных сил эжектирующий газ движется по периферии камеры энергообмена 3, создавая разряжение на оси камеры. Под действием разности атмосферного давления и давление разряжения на оси камеры 3 воздух из атмосферы поступает через сопло 2 в камеру 3. В камере 3 силами вязкости происходит энергообмен между эжектирующим и эжектируемым газами. После завершения энергообмена смесь газов поступает в сопло 4 для дальнейшего использования по назначению. Трение и теплообмен на стенках корпуса не учитываются, так как они учитываются КПД сжатия ηc и расширения ηр.

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

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ РАБОЧЕГО ПРОЦЕССА ВИХРЕВОГО ЭЖЕКТОРА

Сумма энергий эжектирующего и эжектируемого газов равна энергии газа, выходящего из вихревого эжектора

(1)

(2)

где N1 – энергия эжектирующего газа, Вт;
G1 – расход эжектирующего газа, Кг/с;
L – удельная энергия эжектирующего газа, Дж/кг.

(3)

где ηр – суммарный КПД расширения (отношение реальной работы расширения к изоэнтропной).

Энергия эжектируемого газа после подведения к нему энергии от эжектирующего газа силами вязкости за счет разности линейных скоростей

(4)

где G2 – расход эжектируемого газа, кг/с;
L2 – удельная энергия эжектируемого газа после совершения над ним работы эжектирующего газа, Дж/кг.

(5)

где η*c – суммарный КПД процесса сжатия (отношение изоэнтропной работы сжатия к реальной).

Суммарная энергия газа на выходе из вихревого эжектора

(6)

где G3 – расход газа на выходе из вихревого эжектора, кг/с;
L3 – удельная энергия газа на выходе из вихревого эжектора, Дж/кг.

(7)

Cтепень понижения полного давления эжектирующего газа

(8)

где Р01 – полное давление эжектирующего газа на входе в эжектор, Па;
Р03 – полное давление газа на выходе из эжектора, Па.

Степень повышения полного давления эжектируемого газа за счет энергообмена силами вязкости с эжектируемым газом

(9)

где РН – давление среды, откуда поступает эжектируемый газ, Па.

Степень понижения полного давления газа на выходе из вихревого эжектора

(10)

Уравнения (2)–(10) подставляются в уравнение (1):

(11)

Уравнение (11) является уравнением энергии газа в механической форме. В тепловой форме это уравнение имеет вид [11][13][15]

(12)

где i0j = Cpj Т0j – полная энтальпия газа, Дж/кг;
Cpj – теплоемкость при постоянном давлении j-го газа, Дж/(кг•К);
Т0j – полная температура j-го газа, К;
J = 1–3 – параметры газа на входе эжектирующего газа, эжектируемого газа и на выходе из вихревого эжектора соответственно.

При Cpj = Cp = const и G2/G1 = n, уравнение (12) преобразуется к виду

(13)

Полные температуры эжектирующего и эжектируемого газов на входе в эжектор могут иметь одну и ту же температуру, т.е. Т02 = Т01. Тогда из уравнения энергии в тепловой форме (13) следует, что

(14)

C учетом уравнения (14) уравнение энергии в механической форме (11) можно преобразовать следующим образом:

(15)

В уравнении (15) одна неизвестная величина. Это полное давление газа на выходе из вихревого эжектора Р03.

Для нахождения полного давления Р03 раскрываются скобки уравнения (15):

(16)

Все члены уравнения (16) умножаются на величину :

(17)

Уравнение (17) преобразовывается к виду:

(18)

где

;

Полное давление газа на выходе из вихревого эжектора находится путем решения квадратного уравнения (18):

(19)

В уравнении (19) знак плюс перед корнем квадратным взят потому, что полное давление газа на выходе из вихревого эжектора не может иметь отрицательное значение.

На входе в вихревой эжектор скорости эжектирующего и эжектируемого газов имеют различные значения. За счет разности скоростей возникают касательные напряжения, приводящие к снижению скорости эжектирующего газа и повышению скорости эжектируемого газа. Следовательно, силами вязкости кинетическая энергия передается от эжектирующего газа к эжектируемому [1]. Таким образом, эжектирующий газ совершает работу над эжектируемым газом, в результате чего полное давление эжектирующего газа падает, а эжектируемого – растет. Обмен работой приводит также и к изменению полной температуры газа. Полная температура эжектирующего газа падает, эжектируемого – растет [1][12][14].

Изменение полных температур эжектирующего и эжектируемого газов можно определить с помощью известных уравнений термодинамики [6][7]:

(20)

(21)

где T01p – полная температура эжектирующего газа после завершения энергообмена с эжектируемым газом, К;
T02p – полная температура эжектируемого газа после завершения энергообмена с эжектирующим газом, К.

Энергообмен заканчивается тогда, когда полное давление эжектирующего и эжектируемого газов будут равны полному давлению их смеси на выходе из вихревого эжектора, т.е. P03.

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

(22)

Тепловой поток и температуру T03 на выходе из эжектора можно определить по формулам теории теплопередач или по формуле (13).

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

Касательные напряжения можно определять эмпирической зависимостью Ж. Буссинеска, по гипотезам Прандтля, Тейлора, А. Ферри, Колмогорова [16][3].

Эмпирическая зависимость Ж. Буссинеска

Для вычисления турбулентного касательного напряжения пользуются соотношением, аналогичным закону трения Ньютона

(23.1)

где K – коэффициент турбулентной вязкости, который определяется экспериментально.

Гипотеза Прандтля

Исходя из выражения «кажущегося» напряжения трения

(23.2)

Прандтль показал, что

(23.3)

где l – длина пути смешения.

Кроме того, Прандтлем была найдена полуэмпирическая зависимость, которую назвали «новая формула Прандтля».

В этом случае напряжение трения

(23.4)

где H – экспериментальная величина;
– скорости потока на границах зоны смешения.

Гипотеза А. Ферри

Зависимость для коэффициента турбулентной вязкости

(24)

где Vmax – максимальная скорость в струйном слое, b – ширина струйного слоя, H1 – постоянная.

Между коэффициентом турбулентной вязкости и параметрами течений была установлена связь [16]

(25)

где H2 – эмпирическая постоянная; ρ, ρm – соответственно плотности в струе и на внешней границе потока;

– скорости потоков.

Гипотеза Колмогорова

Турбулентная вязкость характеризуется кинетической энергией пульсаций и некоторым масштабом турбулентности

(26)

где
– величина, пропорциональная кинетической энергии турбулентности;
L – масштаб турбулентности, соответствующий среднему размеру турбулентных вихрей;
k = 0,2 – эмпирическая постоянная. Кроме того, величину касательных напряжений можно определить и по формуле [17]

(27)

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

Расчетные и экспериментальные данные приведены на рисунке 2 [15][18].

ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ

Вышеприведенная математическая модель рабочего процесса вихревого эжектора замкнута. По ней можно составить две методики расчета:

  • методика расчета геометрических характеристик вихревого эжектора для заданных термодинамических характеристик;
  • методика расчета термодинамических характеристик вихревого эжектора при известных геометрических параметрах.

Расчеты по данным методикам дают удовлетворительное совпадение с экспериментальными данными собственными и других авторов [2][8][15].

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

Обмен работой и теплотой с эжектируемым газом приводит к падению полного давления и полной температуры эжектирующего газа. На выходе из эжектора полное давление и полная температура эжектирующего и эжектируемого газов выравниваются.

ЗАКЛЮЧЕНИЕ

Представлены результаты теоретических исследований авторов работы и сравнение их с результатами экспериментов других авторов [2][8][15].

Показано, что в вихревом эжекторе в результате взаимодействия эжектируемый газ повышает свое полное давление и полную температуру за счет совершения над ним работы со стороны эжектирующего газа. В результате обмена работой полное давление и полная температура эжектирующего газа падает. В результате энергообмена возникает разность термодинамических температур, что приводит к возникновению теплового потока от высокотемпературного потока к низкотемпературному. Передача энергии от эжектирующего газа к эжектируемому идет силами вязкости за счет разности угловых скоростей [1].

Приложение 1.

Методика расчета оптимальных геометрических размеров вихревого эжектора для получения заданной степени эжекции

Основные данные берутся в соответствии со схемой, показанной на рисунке 1.

Исходные данные:

Определить: F1 – площадь сопла для входа эжектирующего газа (м2); F2 – площадь сопла для ввода эжектируемого газа (м2); F3 – площадь сопла для выхода смеси из эжектора (м2); l – длина эжектора (м), P03 – полное давление газа на выходе из эжектора.

Расчет

1. Коэффициент a

2. Коэффициент b

3. Коэффициент с

4. Полное давление газа на выходе из вихревого эжектора

5. Полная температура газа на выходе из эжектора

6. Площадь среза сопла эжектирующего газа

7. Эжектирующее сопло может быть круглым, квадратным или прямоугольным. Все определяется конструктивными особенностями эжектора.

Диаметр сопла:

8. Расход эжектируемого газа

9. Площадь среза сопла эжектируемого газа

где q2) = 1, так как на оси эжектора вакуум и перепад давлений на срезе сопла сверхкритический.

10. Диаметр среза сопла эжектируемого газа

11. Температура эжектирующего газа после совершения работы над эжектируемым газом

12. Температура эжектируемого газа после обмена работой с эжектирующим газом

13. Газодинамическая функция давления газа на выходе из вихревого эжектора

где Р3 = РН – при истечении в атмосферу с докритическим перепадом давления или при регулируемой площади сопла.

14. Коэффициент скорости на выходе из эжектора

15. Расходная газодинамическая функция на выходе из эжектора

16. Площадь среза сопла на выходе из эжектора

17. Диаметр среза сопла на выходе из эжектора

18. Величина касательных напряжений между эжектирующим и эжектируемым газом

19. Статическая температура газа на выходе из эжектора

где τ (λ3) = f 3) = 0,8335.

20. Критическая скорость газа на выходе из эжектора

21. Скорость газа на выходе из эжектора

22. Средняя скорость движения эжектирующего газа от входа до выхода из эжектора

23. Число Рейнольдса

24. Коэффициент сопротивления трения (формула Блазиуса)

25. Плотность газа на входе в эжектор

26. Длина пути по вихревой линии

27. Внутренний диаметр смещения потоков

28. Длина одного шага

29. Число шагов движения потока газа

30. Линейная длина камеры смещения вдоль оси эжектора

Приложение 2.

Методика расчета термодинамических характеристик вихревого эжектора при известных геометрических размерах

Исходные данные:

P01 – полное давление эжектирующего газа, Па;

T01 – полная температура эжектирующего газа, К;

P02 – полное давление эжектируемого газа, Па;

T02 – полная температура эжектируемого газа, К;

k = 1,4 – показатель адиабаты воздуха;

CD = 1003,5 (Дж/кг•К) – теплоемкость воздуха при постоянном давлении;

; (кг•К/Дж)0,5 – коэффициент m = 0,0404 – для воздуха;

R = 287 – газовая постоянная воздуха;

= 0,92 – КПД процесса расширения;

= 0,85 – КПД процесса сжатия;

F1 = 4,0•10–4 м2 – площадь входного сопла эжектирующего газа;

F2 = 4,146•10–4 м2 – площадь входного сопла эжектируемого газа;

F3 = 8,81•10–4 м2 – площадь выходного сопла вихревого эжектора.

Определить: G1, G2, G3, n, P03, T03.

Расчет

1. Расход эжектируемого воздуха

2. Расход эжектирующего газа

3.

4. Коэффициент a

5. Коэффициент b

6. Коэффициент c

7. Полное давление газа на выходе из вихревого эжектора

8. Полная температура газа на выходе из эжектора

9. Газодинамическая функция давления газа на выходе из вихревого эжектора

10. Коэффициент скорости на выходе из эжектора

11. Расходная газодинамическая функция на выходе из эжектора

12. Расход газа на выходе из эжектора

Список литературы

1. Кузнецов В. И., Макаров В. В. Эффект Ранка: эксперимент, теория, расчет. М.: Инновационное Машиностроение, 2017. 376 с.

2. Ардов Ю. К. Новые газовые эжекторы и эжекционные процессы. М.: Издательство физ.-мат. литературы. 2001. 356 с.

3. Уколов А. И., Родионов В. П. Моделирование дефекта внутренней поверхности струйного кавитатора // Вестник Госуниверситета морского и речного флота им. адмирала Макарова, 2018.

4. Уколов А. И., Родионов В. П. Верификация результатов численного моделирования и экспериментальных данных влияния кавитации на гидродинамические характеристики струйного потока // Вестник Госуниверситета морского и речного флота им. адмирала Макарова, 2018.

5. Курников А. С., Мизгирев А. С., Герепкова Е. А. Расчет гидродинамического кавитатора с тороидальной камерой смещения // Вестник ВГАВТ. 2016. Вып. 47.

6. Блинков В. Н., Мелихов В. И., ПарфеновЮ.В., Никонов С. М., Емельянов Д. А., Елкин Н. В., Трубкин Е. И. Экспериментальное исследование расходно-напорной характеристики эжектора с тангенциальным подводом пассивной среды для системы аварийного охлаждения активной зоны АЭС с ВВЭР // Фундаментальные исследования. 2013. № 10(5). С. 949–952. URL: http://www.fundamentalresearch.ru/ru/article/view?id=32430

7. Ершов М. И., Волкова Ю. В., Мунц В. А. Моделирование режимов работы эжектора для рециркуляции анодных газов // Труды второй НТК молодых ученых Уральского энергетического института. Екатеринбург: УрФУ, 2017. С. 62–64.

8. Усс А. Ю., Пугачук А. С., Чернышев А. В., Тухтабуллин Ф. Г. Разработка стенда для визуализации и экспериментального исследования рабочего процесса в вихревом струйном устройстве // Омский научный вестник. Сер. Авиационно-ракетное и энергетическое машиностроение – 2020. – Т. 4, № 3. – С. 45–55.

9. Усс А. Ю., Чернышев А. В., Атамасов Н. В. Разработка метода расчета и создание вихревого струйного устройства для управления потоком газа. 2019.

10. Uss A. Yu., Chernyshev A. V. The Development of the Vortex Gas Presenture Regulator // Procedia Engineering. 2016. Vol. 152.P. 380–388. DOI: 10.1016./j.proeng.2016.07.718

11. Uss A. Yu., Chernyshev A. V. The Development of the Vortex Gas Presenture Regulator // AJP Conference Proceedings. 2017. Vol. 1876. 020025. DOI: 10. 1063/1.4998845

12. Tesat V., Brouckova Z., Kordik J., Travnicek Z., Peszynsku V. Vortex With flow control by syntheticjets // The European Physical Journal Web of Conferences. 2012. Vol. 25. 01092. P. DOI: 10.1051/epjcont/20122501092

13. Усс А. Ю., Чернышев А. В., Атамасов Н. В. Разработка метода расчета и создание вихревого струйного устройства для управления потоками газа // Омский научный вестник. Сер. Авиационно-ракетное и энергетическое машиностроение. 2019. Т. 3. № 2. С. 78–86. DOI: 10.25206/2588-0373-2019-3-2-7886

14. Степин В. А. Вихревые и двухфазные потоки в технологии промывки судовых систем: монография. Северный (Арктический) федеральный университет им. М.В. Ломоносова.

15. Кузнецов В. И., Якумушкин Р. В., Шенбергер А. Ю. Результаты сравнительных испытаний охладителя наддувочного воздуха комбинированного дизеля // Сб. научных статей. Спб: (ВИ ЖДВИ ВОСО) ВА МТО, 2019. С. 131–139.

16. Краснов Н. Ф., Кощевой В. Н., Калугин В. Т. Аэродинамика отрывных течений: учеб. пособие для втузов. М.: Высшая школа, 1988. 351 с.

17. Кутателадзе С. С. Теплопередача и гидродинамическое сопротивление. Справочное пособие. М.: Энергоатомиздат, 1990. 367 с.

18. Халиулин Р. Р., Сыченков В. А., Панченко В. И. Разработка методики расчета эжекторов с малыми коэффициентами эжекции // Известия высших учебных заведений. Авиационная техника. 2018. № 1. С. 65–70.


Об авторах

В. И. Кузнецов
Федеральное государственное автономное образовательное учреждение высшего образования «Омский государственный технический университет»
Россия

доктор технических наук, профессор, профессор кафедры «Авиа- и ракетостроение»,

Омск



В. В. Макаров
Федеральное государственное автономное образовательное учреждение высшего образования «Омский государственный технический университет»
Россия

кандидат технических наук, доцент, доцент кафедры «Авиа- и ракетостроение»,

Омск



А. Ю. Шандер
Федеральное государственное автономное образовательное учреждение высшего образования «Омский государственный технический университет»
Россия

старший преподаватель кафедры «Авиа- и ракетостроение»,

Омск



Рецензия

Для цитирования:


Кузнецов В.И., Макаров В.В., Шандер А.Ю. Теоретические основы рабочего процесса вихревой трубы в режиме вакуум-насоса (вихревого эжектора). Вестник Концерна ВКО «Алмаз – Антей». 2022;(2):63-73. https://doi.org/10.38013/2542-0542-2022-2-63-73

For citation:


Kuznetsov V.I., Makarov V.V., Shander A.Yu. Fundamentals of the working process of a vortex tube operating in the vacuum pump (vortex ejector) mode. Journal of «Almaz – Antey» Air and Space Defence Corporation. 2022;(2):63-73. https://doi.org/10.38013/2542-0542-2022-2-63-73

Просмотров: 69


Creative Commons License
Контент доступен под лицензией Creative Commons Attribution 4.0 License.


ISSN 2542-0542 (Print)