Научная статья на тему 'ИССЛЕДОВАНИЕ ВЛИЯНИЯ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ НАПРАВЛЯЮЩЕГО АППАРАТА НА ЭНЕРГЕТИЧЕСКУЮ ЭФФЕКТИВНОСТЬ ПНЕВМАТИЧЕСКОГО ЗАТВОРА ТЕПЛОМАССООБМЕННЫХ РЕАКТОРОВ'

ИССЛЕДОВАНИЕ ВЛИЯНИЯ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ НАПРАВЛЯЮЩЕГО АППАРАТА НА ЭНЕРГЕТИЧЕСКУЮ ЭФФЕКТИВНОСТЬ ПНЕВМАТИЧЕСКОГО ЗАТВОРА ТЕПЛОМАССООБМЕННЫХ РЕАКТОРОВ Текст научной статьи по специальности «Механика и машиностроение»

CC BY
0
0
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
агрессивные газы / реактор / пневматический затвор / численное моделирование / запирающая способность / энергетическая эффективность / hazardous gases / reactor / pneumatic gate valve / numerical simulation / sealing capability / energy efficiency

Аннотация научной статьи по механике и машиностроению, автор научной работы — Хабарова Дарья Федоровна, Исмагилов Александр Рашидович, Исмагилов Денис Рашидович

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

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

Похожие темы научных работ по механике и машиностроению , автор научной работы — Хабарова Дарья Федоровна, Исмагилов Александр Рашидович, Исмагилов Денис Рашидович

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

THE INFLUENCE OF THE GEOMETRIC PARAMETERS OF THE GUIDE VANES ON THE ENERGY EFFICIENCY OF A PNEUMATIC GATE VALVE FOR HEAT MASS TRANSFER REACTORS

This article presents the results of a study of a vortex-type pneumatic gate valve. The valve is used to seal the raw material loading channels of heat and mass transfer reactors. The article gives a method for calculating a pneumatic vortex-type gate valve, developed to determine the optimal operating parameters of the device for given geometric parameters and to profile the optimal geometry of the valve for given operating parameters. The article provides an analysis of the mathematical model of the operation of a pneumatic vortex-type gate valve, from which the method for calculating the device was derived. It is shown that the mathematical description of the gas dynamics in the valve does not take into account the gas velocity components at the nozzle exit, which makes it inapplicable for the analysis of the influence of the guide vane configuration on the characteristics of the pneumatic gate valve. Numerical modeling of the gas flows in the valve shows that an increase in the thickness and number of blades in a certain range leads to a deterioration in the sealing capability of the valve and, as a consequence, to an increase in the mass flow rate and gas pressure required for sealing the reactor cavity and reducing the energy efficiency of the device. It is shown that a decrease in the blade angle within the range under study also leads to an increase in the required parameters of the gas for sealing the reactor cavity. The results show that the current method for calculating the operating and geometric parameters of pneumatic vortex-type gate valves does not allow the main geometric parameters of the valve – installation angle, thickness and number of guide vanes – to be calculated.

Текст научной работы на тему «ИССЛЕДОВАНИЕ ВЛИЯНИЯ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ НАПРАВЛЯЮЩЕГО АППАРАТА НА ЭНЕРГЕТИЧЕСКУЮ ЭФФЕКТИВНОСТЬ ПНЕВМАТИЧЕСКОГО ЗАТВОРА ТЕПЛОМАССООБМЕННЫХ РЕАКТОРОВ»

Научная статья УДК 621.542

DOI: 10.14529/power230410

ИССЛЕДОВАНИЕ ВЛИЯНИЯ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ НАПРАВЛЯЮЩЕГО АППАРАТА НА ЭНЕРГЕТИЧЕСКУЮ ЭФФЕКТИВНОСТЬ ПНЕВМАТИЧЕСКОГО ЗАТВОРА ТЕПЛОМАССООБМЕННЫХ РЕАКТОРОВ

Д.Ф. Хабарова, khabarovadf@susu.ru, https://orcid.org/0000-0003-0251-0138 А.Р. Исмагилов, ismagilovar@susu.ru, https://orcid.org/0000-0001-6538-0727 Д.Р. Исмагилов, ismagilovdr@susu.ru, https://orcid.org/0000-0002-9312-295X Южно-Уральский государственный университет, Челябинск, Россия

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

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

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

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

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

Ключевые слова: агрессивные газы, реактор, пневматический затвор, численное моделирование, запирающая способность, энергетическая эффективность

Благодарности. Исследование выполнено за счет гранта Российского научного фонда № 22-29-20084, https://rscf.ru/project/22-29-20084/

Для цитирования: Хабарова Д.Ф., Исмагилов А.Р., Исмагилов Д.Р. Исследование влияния геометрических параметров направляющего аппарата на энергетическую эффективность пневматического затвора тепломассообменных реакторов // Вестник ЮУрГУ. Серия «Энергетика». 2023. Т. 23, № 4. С. 86-96. DOI: 10.14529/power230410

Original article

DOI: 10.14529/power230410

THE INFLUENCE OF THE GEOMETRIC PARAMETERS OF THE GUIDE VANES ON THE ENERGY EFFICIENCY OF A PNEUMATIC GATE VALVE FOR HEAT MASS TRANSFER REACTORS

D.F. Khabarova, khabarovadf@susu.ru, https://orcid.org/0000-0003-0251-0138 A.R. Ismagilov, ismagilovar@susu.ru, https://orcid.org/0000-0001-6538-0727 D.R. Ismagilov, ismagilovdr@susu.ru, https://orcid.org/0000-0002-9312-295X South Ural State University, Chelyabinsk, Russia

Abstract. This article presents the results of a study of a vortex-type pneumatic gate valve. The valve is used to seal the raw material loading channels of heat and mass transfer reactors. The article gives a method for calculating a pneumatic vortex-type gate valve, developed to determine the optimal operating parameters of the device for given geometric parameters and to profile the optimal geometry of the valve for given operating parameters.

© Хабарова Д.Ф., Исмагилов А.Р., Исмагилов Д.Р., 2023

The article provides an analysis of the mathematical model of the operation of a pneumatic vortex-type gate valve, from which the method for calculating the device was derived. It is shown that the mathematical description of the gas dynamics in the valve does not take into account the gas velocity components at the nozzle exit, which makes it inapplicable for the analysis of the influence of the guide vane configuration on the characteristics of the pneumatic gate valve.

Numerical modeling of the gas flows in the valve shows that an increase in the thickness and number of blades in a certain range leads to a deterioration in the sealing capability of the valve and, as a consequence, to an increase in the mass flow rate and gas pressure required for sealing the reactor cavity and reducing the energy efficiency of the device. It is shown that a decrease in the blade angle within the range under study also leads to an increase in the required parameters of the gas for sealing the reactor cavity.

The results show that the current method for calculating the operating and geometric parameters of pneumatic vortex-type gate valves does not allow the main geometric parameters of the valve - installation angle, thickness and number of guide vanes - to be calculated.

Keywords: hazardous gases, reactor, pneumatic gate valve, numerical simulation, sealing capability, energy efficiency

Acknowledgments. The study was supported by the Russian Science Foundation grant No. 22-29-20084, https://rscf.ru/project/22-29-20084/

For citation: Khabarova D.F., Ismagilov A.R., Ismagilov D.R. The influence of the geometric parameters of the guide vanes on the energy efficiency of a pneumatic gate valve for heat mass transfer reactors. Bulletin of the South Ural State University. Ser. Power Engineering. 2023;23(4):86-96. (In Russ.) DOI: 10.14529/power230410

Введение

Многие технологические процессы на промышленных предприятиях металлургической, химической, нефтяной, газовой и смежных отраслях сопровождаются выхлопами агрессивных газов, ухудшающих экологическую обстановку окружающей среды [1-5]. Общеизвестные промышленные технологии получения современных материалов включают химические реакции под давлением при высоких температурах, сопровождающиеся фазовыми переходами [6-15]. Такие процессы протекают в тепломассообменных аппаратах - реакторах - в составе барботажных колонн, аппаратов гетерогенного катализа и химического синтеза, фьюминг-печей, конвертеров, печей для выплавки стали погружным факелом и др. [1-3, 6, 7, 9].

Рассматриваемые технологические процессы сопровождаются загрузкой сыпучего материала во внутреннюю полость реактора, находящуюся под избыточным давлением, через специальный канал загрузки, снабженный, как правило, механическим затвором. При загрузке запорное механическое устройство открыто, а при протекании химических реакций - закрыто. Применение механического запорного устройства сопряжено с нежелательными выбросами газов из полости реактора в процессе загрузки сырья, а также через зазоры затвора в его закрытом состоянии [2, 3, 8].

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

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

В работе [15] приводится метод расчета пневматического запорного устройства вихревого типа. Согласно работе [15], предложенный метод расчета позволяет решить как задачу определения оптимальных режимных параметров устройства при заданной геометрии, так и обратную задачу - определения оптимальных геометрических параметров при заданных режимных. Однако исследований, посвященных верификации и апробации научных основ данного метода, в научно-технической литературе обнаружить не удалось.

Постановка цели и задач

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

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

Научная новизна

В результате исследования установлено, что угол установки, толщина и число лопаток направляющего аппарата существенно влияют на требуемые для запирания давления в реакторе пара-

метры активного газа и, как следствие, энергетическую эффективность затвора.

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

Теоретическая часть

Принципиальная схема пневматического затвора вихревого типа показана на рис. 1 [15].

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

Пневмозатвор является работоспособным при условии полного запирания загрузочного канала при отсутствии истечения выхлопных газов. Рабо-

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

Вихрь формируется за счёт отклонения активного потока лопастями в количестве п толщиной 5, высотой Н, установленными под углом а к плоскости среза сопла. В результате истечения активного потока из кольцевого сопла с лопастями на выходе из него поток имеет осевые их3 и тангенциальные иф3 составляющие скорости с общим массовым расходом т. Формирование потока происходит при воздействии давления на входе в сопло р\ и температуры среды Ть а также давления на выходе из сопла рр (давление в реакторе) и температуры газов в реакторе Тр. Интенсификация активного потока происходит при увеличении

Рис. 1. Расчетная схема вихревого пневмозатвора: n - число лопаток; 5 - толщина лопатки; R3 - радиус запираемого отверстия; H - высота лопатки Fig. 1. Schematic diagram of a vortex-type pneumatic gate valve: n - number of blades; 5 - blade thickness; R3 - radius of the sealed hole; H - blade height

массового расхода активного газа m, тангенциальной составляющей скорости иф3, давления питания рь температуры Т1. Вышеперечисленные параметры связаны между собой уравнениями математической модели, на основе которой разработан метод расчета пневматического затвора, приведенный в работе [15]:

ф3

2 ( pp - f2 )

RPTP

_ p2 R2T2

(1)

'лх3

(R3 + H )2

, Р2Мф3

(pp - p2)--— +

%(R3 + H)

2 Р2Мф3 in ррмф3

2

2ßpp

Р2«ф3

^Р3 H (2Я3 + H— 2Р3

m = Мх3Р3 A3 = Мх3Р3H (2R3 + H

Р3жН (2R3 + H)

-- v^

tga = _ii;

иф3

tg2a = -Р2(1 + y)

2Р3(2 + y)

Рр Р2

1 + ln

2 Л Рриф3

2ßpp

y

2(1 + y)2

(2)

(3)

(4)

(5)

где

р2, р3, Рр - плотность газа в окружающей среде, на срезе кольцевого сопла и в реакторе соответственно; р2 - давление окружающей среды;

А3 = пИ (2R3 + И) - пЪИ = пИ (2R3 + И - эффективная площадь среза сопла; Я3 - радиус канала загрузки;

в - коэффициент расхождения давления на срезе сопла и давления в реакторе; + И) - и5 , п5

V = -

n(2R3 + H)

= 1 —

n(2R3 + H)

- коэффициент стеснения лопатками;

Y =

R3 H

относительный радиус канала загрузки.

Количество выхлопных газов в загрузочном канале зависит от давления запираемого потока выхлопных газов из реактора в загрузочный канал рр, температуры Тр и радиуса загрузочного канала Я3. Увеличение этих параметров увеличивает количество выхлопных газов в загрузочном канале.

Анализ уравнений (1)-(5) показывает, что максимальное запираемое давление в реакторе определяется величиной тангенциальной составляющей скорости активного газа на срезе кольцевого сопла: при ее увеличении максимальное запираемое давление растет. Данная зависимость, рассчитанная при абсолютном давлении в загрузочном канале р2 = 100 кПа, абсолютном давлении на входе в сопло р1 = 375 кПа, абсолютных температурах в загрузочном канале и на входе в сопло Т1 = Т2 = 293 К, абсолютной температуре в реакторе Тр = 673 К, радиусе канала загрузки Я2 = 146 мм показана на рис. 2. Следует отметить, что монотонно возрастающая зависимость рр от иф3 не является ожидаемой с точки зрения функционального назначения, создаваемого тангенциальной составляющей скорости вихря, при увеличении диаметра которого возможна газодинамическая картина потоков, не способствующая запиранию канала, - вихри локализуются на удаленной от канала загрузки периферии.

Рис. 2. Взаимосвязь величины тангенциальной составляющей скорости активного газа и максимального запираемого противодавления в реакторе Fig. 2. Relationship between the tangential component of the active gas velocity and the maximum sealed back pressure in the reactor

p

4

+

Heat engineering

Рис. 3. Влияние режимных и геометрических параметров затвора на величину тангенциальной составляющей скорости активного газа: а - угла установки лопатки направляющего аппарата;

b - массового расхода активного газа Fig. 3. The influence of the operating and geometric parameters of the valve on the value of the tangential component of the active gas velocity: a - of the installation angle of the guide vane; b - of mass flow rate

of active gas

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

Зависимости тангенциальной составляющей скорости активного газа иф3 от угла установки лопатки а и массового расхода активного газа т показаны на рис. 3. Данные зависимости рассчитаны при тех же параметрах, которые описаны к рис. 2.

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

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

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

Практическая часть

В качестве исходных данных для численного моделирования приняты следующие геометриче-

ские размеры вихревого пневмозатвора, обозначение которых приведено на рис. 1:

- ширина сопла Н = 0,88 мм;

- радиус канала загрузки сырья R2 = 150 мм;

- угол установки лопасти а = 80° (случаи а = {40°; 60°; 80°} рассмотрены для конструкции с 8 лопатками и толщиной лопаток 2 мм);

- количество лопаток 8, 12 и 16;

- толщина лопатки 5 = 2 мм (случаи 5 = {1 мм, 2 мм, 3 мм} рассмотрены для конструкции с 8 лопатками и углом установки лопасти а = 80°).

Также приняты следующие физические условия работы вихревого пневмозатвора:

- абсолютное давление воздуха перед соплом р1 = 400 кПа;

- абсолютное давление в реакторе pp = 100 кПа и 102 кПа;

- барометрическое давление в цехе p2 = 100 кПа;

- температуры T = T2 = 293 К, Тр = 673 К.

- массовый расход активного газа m = 0,577 кг/с.

При проведении моделирования вводились

следующие допущения и упрощения:

- моделирование проводилось в стационарной постановке (установившееся течение);

- стенки конструкции являются гладкими, не-разрушаемыми и непроницаемыми;

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

В качестве математической модели, описывающей движение вязкого турбулентного газа, принимается система уравнений Навье - Стокса, осредненная по Рейнольдсу (RANS) [16-19], замкнутая полуэмпирической моделью турбулентности Realizable k-e [20-24].

Семейства двухпараметрических диссипатив-ных k-e моделей турбулентности хорошо описывают свойства свободных сдвиговых течений моделью турбулентности. Realizabl k-e модель по

сравнению со стандартной k-e моделью имеет два существенно важных отличия [24]:

- Realizable k-e модель содержит альтернативную формулировку для турбулентной вязкости;

- модифицированное уравнение переноса для скорости диссипации e было получено из точного уравнения для переноса среднеквадратичных колебаний завихренности.

Данная модель удовлетворяет точным математическим ограничениям по напряжениям Рей-нольдса, вытекающим из физики турбулентной жидкости.

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

Твердотельная геометрическая модель конструкции с основными геометрическими размерами показана на рис. 4.

Расчетная схема модели вихревого пневмоза-твора для CFD-анализа приведена на рис. 5.

Сформированная, на основе геометрической модели (рис. 5а), расчетная сеточная модель (рис. 5b)

состоит из тетраэдрических элементов (размер элемента 10 мм на поверхности каналов, 1 мм на поверхности каналов в области лопаток, 30 мм на поверхности реактора, 0,5 мм на поверхностях лопаток) и призматических элементов вблизи стенок (высота первого слоя - 0,5-0,02 мм в зависимости от зазора вблизи поверхности стенок, количество слоев - 2-6). Общее количество элементов составляет 13,97 млн.

Сеточная модель имеет следующие критерии качества расчетных элементов:

- минимальное значение коэффициента ортогональности 3 • 10-05;

- минимальное значение коэффициента пропорциональности 1390;

- минимальный объем расчетного элемента 1,71-10-13 м3, максимальный объем расчетного элемента 2,83-Ш-6 м3.

На рис. 6 показаны распределения параметров потока в продольном сечении при противодавлении в реакторе рр = 0 кПа и 2 кПа (по избыточной шкале) для конструкции с 12 лопатками толщиной 2 мм.

На рис. 6Ь по векторному полю скоростей видно, что при противодавлении в реакторе 2 кПа

Рис. 4. Твердотельная геометрическая модель конструкции с основными

геометрическими размерами Fig. 4. Solid geometric model of the structure with basic geometric dimensions

поток активного газа выходит через сечение канала загрузки сырья. Это свидетельствует о срыве рабочего режима устройства - запирание полости реактора при противодавлении 2 кПа не происходит.

Результаты трёхмерного моделирования при различных количествах лопастей п, их толщин 5 и углов установки а сведены в таблице.

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

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

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

Практическая значимость

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

Рис. 5. К CFD-анализу параметров затвора: а - расчетная схема конструкции; b - конечно-элементная расчётная модель Fig. 5. To CFD analysis of valve parameters: a - calculating diagram of the structure; b - finite element calculation model

Исмагилов Д.Р. направляющего аппарата на энергетическую эффективность...

Рис. 6. Распределение параметров потока в продольном сечении (сверху вниз) результирующей скорости (локальный диапазон изменения, м/с), избыточного давления (ати), температуры (К), векторное поле скоростей (м/с): а - при противодавлении в реакторе 0 кПа; b - при противодавлении в реакторе 2 кПа Fig. 6. Distribution of flow parameters in the longitudinal section (from top to bottom) resulting velocity (local range of change, m/s), excess pressure (atm), temperature (K), vector velocity field (m/s): a - at a back pressure in the reactor of 0 kPa; b - at a back pressure in the reactor of 2 kPa

Результаты моделирования Simulation results

Количество лопаток n Толщина лопаток 5, мм Угол установки лопаток а, град. Избыточное давление на входе в канал подвода активного газа, кПа Массовый расход газа через канал загрузки (выхлоп), кг/с

Противодавление в реакторе (давление сопротивления на границах области реактора, кроме верхней поверхности), кПа

0 2 0 2

8 2 40 - 530,3 - -1,019

8 2 60 - 328,3 - -0,863

8 1 80 259 270,3 1,49 -0,692

2 262,8 273,3 1,50 -0,698

3 255,6 267,3 1,52 -0,752

12 2 80 261,7 271,3 1,66 -0,771

16 2 80 267,8 280,5 1,68 -0,811

Заключение

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

Численным моделированием газодинамических процессов в пневматическом затворе вихревого типа показано, что уменьшение угла установ-

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

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

1. Основы металлургического производства / В.А. Бигеев, К.Н. Вдовин, В.М. Колокольцев и др. СПб.: Лань, 2017. 615 с.

2. Спиридонов Е.К., Гришина Е.А., Подзерко А.В. Характеристики и расчет эжекционного пневмоза-твора // Инновация, экология и ресурсосберегающие технологии на предприятиях машиностроения, авиастроения, транспорта, сельского хозяйства: тр. IX Междунар. науч.-техн. конф. Ростов н/Д: ИЦ ДГТУ, 2010. С. 784-790.

3. Гришина Е.А., Спиридонов Е.К., Подзерко А.В. Рабочий процесс и конструкции эжекционного пневмозатвора // Динамика машин и рабочих процессов: сб. докл. всерос. науч.-техн. конф. 2009. С. 41-46.

4. Guerrero A.P., Paredes G.E. Linear and Non-linear Stability Analysis in Boiling Water Reactors. The Design of Real-Time Stability Monitors. Woodhead Publishing Series in Energy, 2019. 465 p. DOI: 10.1016/C2017-0-01640-3

5. Potential environmental pollution from copper metallurgy and methods of management / G. Izydorczyk, K. Mikula, D. Skrzypczak et al. // Environmental Research. 2021. Vol. 197. P. 111050. DOI: 10.1016/j.envres.2021.111050

6. Roderick I.L.G. Chapter 4.1 - Process Kinetics, Fluid Flow, and Heat and Mass Transfer, in Process Metallurgy // Treatise on Process Metallurgy. Publisher: Elsevier, Netherlands, 2014. P. 428-443. DOI: 10.1016/B978-0-08-096984-8.00007-0

7. Dynamic chaos of imaging measurements for characterizing gas-liquid nonlinear flow behaviour in a metallurgical reactor stirred by top-blown air / K. Yang, X. Zhang, G. Yang et al. // The Canadian Journal of Chemical Engineering. 2023. Vol. 1. P. 1-17. DOI: 10.1002/cjce.25077

8. Pearson M., Martinez Fonte D. Use of Remotely Operated Shut-Off Valves for Emergency Isolation of Pressure Vessels and Hazardous Fluids // Proceedings of the 62nd Conference of Metallurgists. 2023. P. 107-119. DOI: 10.1007/978-3-031-38141-6_12

9. Yamazaki Y. Chapter 3 - Gasification Reactions of Metallurgical Coke and Its Application - Improvement of Carbon Use Efficiency in Blast Furnace // Gasification for Practical Applications. Publisher: InTech, Croatia, 2012. P. 51-84. DOI: 10.5772/51680

10. Blanchard J.N., Brunet Y., Merlen A. Influence of a counter rotating vortex pair on the stability of a jet in a cross flow: an experimental study by flow visualizations // Experiments in Fluids. 1999. Vol. 26. P. 63-74. DOI: 10.1007/s003480050265

11. Gopalan S., Abraham B.M., Katz J. The structure of a jet in cross flow at low velocity ratios // Physics of Fluids. 2004. Vol. 16. P. 2067-2087. DOI: 10.1063/1.1697397

12. Ben M.M., Malcangio D., Mossa M. Vertical dense jet in flowing current // Environmental Fluid Mechanics. 2018. Vol. 18. P. 75-96. DOI: 10.1007/s10652-017-9515-2

13. Farokhi S., Taghavi R., Rice E.J. Effect of initial swirl distribution on the evolution of a turbulent jet // AIAA Journal. 2012. Vol. 27 (6). P. 700-706. DOI: 10.2514/3.10168

14. Illyas S.M., MuthuManokar A., Kabeel A.E. Experimental and Computational Study on Effect of Vanes on Heat Transfer and Flow Structure of Swirling Impinging Jet // Journal of Applied Fluid Mechanics. 2023. Vol. 16 (2). P. 205-221. DOI: 10.47176/JAFM.16.02.1296

15. Спиридонов Е.К., Гришина Е.А., Подзерко А.В. К расчету пневмозатворов // Вестник ЮУрГУ. Серия «Машиностроение». 2011. № 11 (228). С. 4-11.

16. Cottet G.-H., Koumoutsakos P.D. Chapter 2 - Vortex Methods for Two-Dimensional Flows // Vortex Methods: Theory and Practice. Cambridge University Press, USA, 2000. P. 10-54.

17. Anderson J.D. Computational Fluid Dynamics. The basics with applications. 3rd ed. Springer, Germany, 2009. 332 p.

18. Girault V., Raviart P.A. Finite Element Methods for Navier-Stokes Equations: Theory and Algorithms. Springer Berlin, Heidelberg, 1986. 376 p. DOI: 10.1007/978-3-642-61623-5

19. Lee C., Jiang X. Flow structures in transitional and turbulent boundary layers // Physics of Fluids. 2019. Vol. 31 (11). P. 11301. DOI: 10.1063/1.5121810

20. Shaheed R., Mohammadian A., Gildeh H.K. A comparison of standard k-e and realizable k-e turbulence models in curved and con-fluent channels // Environmental Fluid Mechanics. 2019. Vol. 19 (3). P. 543-568. DOI: 10.1007/s10652-018-9637-1

21. Sun Y., Xu G., Shi Y. Numerical investigation on noise reduction of rotor blade-vortex interaction using blade surface jet blowing // Aerospace Science and Technology. 2021. Vol. 116. P. 106868. DOI: 10.1016/j.ast.2021.106868

22. Versteeg H.K., Malalasekera W. An introduction to computational fluid dynamics the finite volume method. 2nd ed. Bell and Bain Ltd, Scotland, 2007. 353 p.

23. Simulation of gas-liquid two-phase flow in metallurgical process / W. Bo, S. Shiyi, R. Yanwei et al. // ActaMetall Sin. 2020. Vol. 56. P. 619-632. DOI: 10.11900/0412.1961.2019.00385

24. A new k-e eddy viscosity model for high Reynolds number turbulent flows / T. Shih, W. Liou, A. Shabbir et al. // Computers & Fluids. 1995. Vol. 24. P. 227-238.

References

1. Bigeev V.A., Vdovin K.N., Kolokol'tsev V.M. et al. Osnovy metallurgicheskogo proizvodstva [Basics of metallurgical production]. St. Petersburg: Lan' Publ.; 2017. 615 p. (In Russ.)

2. Spiridonov E.K., Grishina E.A., Podzerko A.V. [Characteristics and calculation of ejection pneumatic valve]. In: Innovatsiya, ekologiya i resursosberegayushchie tekhnologii na predpriyatiyakh mashinostroeniya, aviastroeniya, transporta, sel'skogo khozyaystva: tr. IX Mezhdunar. nauch.-tekhn. konf. [Innovation, ecology and resource-saving technologies at enterprises of mechanical engineering, aircraft manufacturing, transport, agriculture: proceedings IX International scientific-technical conference]. Rostov-on-Don: DSTU Publ.; 2010. P. 784-790. (In Russ.)

3. Grishina E.A., Spiridonov E.K., Podzerko A.V. [Working process and designs of ejection pneumatic valve]. In: Dinamika mashin i rabochikh protsessov: sb. dokl. vseros. nauch.-tekhn. konf. [Dynamics of machines and work processes: collection of reports of the All-Russian scientific and technical conference]. 2009. P. 41-46. (In Russ.)

4. Guerrero A.P., Paredes G.E. Linear and Non-linear Stability Analysis in Boiling Water Reactors. The Design of Real-Time Stability Monitors. Woodhead Publishing Series in Energy; 2019. 465 p. DOI: 10.1016/C2017-0-01640-3

5. Izydorczyk G., Mikula K., Skrzypczak D., Moustakas K., Witek-Krowiak A., Chojnacka K. Potential environmental pollution from copper metallurgy and methods of management. Environmental Research. 2021;197:111050. DOI: 10.1016/j.envres.2021.111050

6. Roderick I.L.G. Chapter 4.1 - Process Kinetics, Fluid Flow, and Heat and Mass Transfer, in Process Metallurgy. In: Treatise on Process Metallurgy. Publisher: Elsevier, Netherlands; 2014. P. 428-443. DOI: 10.1016/B978-0-08-096984-8.00007-0

7. Yang K., Zhang X., Yang G., Li M.,Wang H., Xiao Q. Dynamic chaos of imaging measurements for characterizing gas-liquid nonlinear flow behaviour in a metallurgical reactor stirred by top-blown air. The Canadian Journal of Chemical Engineering. 2023;1:1-17. DOI: 10.1002/cjce.25077

8. Pearson M., Martinez Fonte D. Use of Remotely Operated Shut-Off Valves for Emergency Isolation of Pressure Vessels and Hazardous Fluids. In: Proceedings of the 62nd Conference of Metallurgists. 2023. P. 107-119. DOI: 10.1007/978-3-031-38141-6_12

9. Yamazaki Y. Chapter 3 - Gasification Reactions of Metallurgical Coke and Its Application - Improvement of Carbon Use Efficiency in Blast Furnace. In: Gasification for Practical Applications. Publisher: InTech, Croatia; 2012. P. 51-84. DOI: 10.5772/51680

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

10. Blanchard J.N., Brunet Y., Merlen A. Influence of a counter rotating vortex pair on the stability of a jet in a cross flow: an experimental study by flow visualizations. Experiments in Fluids. 1999;26:63-74. DOI: 10.1007/s003480050265

11. Gopalan S., Abraham B.M., KatzJ. The structure of a jet in cross flow at low velocity ratios. Physics of Fluids. 2004;16:2067-2087. DOI: 10.1063/1.1697397

12. Ben M.M., Malcangio D., Mossa M. Vertical dense jet in flowing current. Environmental Fluid Mechanics. 2018;18:75-96. DOI: 10.1007/s10652-017-9515-2

13. Farokhi S., Taghavi R., Rice E.J. Effect of initial swirl distribution on the evolution of a turbulent jet. AIAA Journal. 2012;27(6):700-706. DOI: 10.2514/3.10168

14. Illyas S.M., MuthuManokar A., Kabeel A.E. Experimental and Computational Study on Effect of Vanes on Heat Transfer and Flow Structure of Swirling Impinging Jet. Journal of Applied Fluid Mechanics. 2023; 16(2):205-221. DOI: 10.47176/JAFM.16.02.1296

15. Spiridonov E.K., Grishina E.A., Podzerko A.V. To calculation of pneumoshutters. Bulletin of the South Ural State University. Ser. Mechanical Engineering. 2011;11(228):4—11. (In Russ.)

16. Cottet G.-H., Koumoutsakos P.D. Chapter 2 — Vortex Methods for Two-Dimensional Flows. In: Vortex Methods: Theory and Practice. Cambridge University Press, USA; 2000. P. 10-54.

17. Anderson J.D. Computational Fluid Dynamics. The basics with applications. 3rd ed. Springer, Germany; 2009. 332 p.

18. Girault V., Raviart P.A. Finite Element Methods for Navier-Stokes Equations: Theory and Algorithms. Springer Berlin, Heidelberg; 1986. 376 p. DOI: 10.1007/978-3-642-61623-5

19. Lee C., Jiang X. Flow structures in transitional and turbulent boundary layers. Physics of Fluids. 2019;31(11):11301. DOI: 10.1063/1.5121810

20. Shaheed R., Mohammadian A., Gildeh H.K. A comparison of standard k-e and realizable k-e turbulence models in curved and con-fluent channels. Environmental Fluid Mechanics. 2019;19(3):543-568. DOI: 10.1007/s10652-018-9637-1

21. Sun Y., XuG., Shi Y. Numerical investigation on noise reduction of rotor blade-vortex interaction using blade surface jet blowing. Aerospace Science and Technology. 2021;116:106868. DOI: 10.1016/j.ast.2021.106868

22. Versteeg H.K., Malalasekera W. An introduction to computational fluid dynamics the finite volume method. 2nd ed. Bell and Bain Ltd, Scotland; 2007. 353 p.

23. Bo W., Shiyi S., Yanwei R., Shuyong C., Wangjun P., Jieyu Z. Simulation of gas-liquid two-phase flow in metallurgical process. ActaMetall Sin. 2020;56:619-632. DOI: 10.11900/0412.1961.2019.00385

24. Shih T., Liou W., Shabbir A., Yang Z., Zhu J. A new k-e eddy viscosity model for high Reynolds number turbulent flows. Computers & Fluids. 1995;24:227-238.

Информация об авторах

Хабарова Дарья Федоровна, канд. техн. наук, доц., кафедра гидравлики и гидропневмосистем, Южно-Уральский государственный университет, Челябинск, Россия; khabarovadf@susu.ru.

Исмагилов Александр Рашидович, канд. техн. наук, доц., кафедра гидравлики и гидропневмосистем, Южно-Уральский государственный университет, Челябинск, Россия; ismagilovar@susu.ru.

Исмагилов Денис Рашидович, старший преподаватель, кафедра летательных аппаратов, ЮжноУральский государственный университет, Челябинск, Россия; ismagilovdr@susu.ru.

Information about the authors

Darya F. Khabarova, Cand. Sci. (Eng.), Ass. Prof., Department of Hydraulics, Hydraulic and Pneumatic Systems, South Ural State University, Chelyabinsk, Russia; khabarovadf@susu.ru.

Alexandr R. Ismagilov, Cand. Sci. (Eng.), Ass. Prof., Department of Hydraulics, Hydraulic and Pneumatic Systems, South Ural State University, Chelyabinsk, Russia; ismagilovar@susu.ru.

Denis R. Ismagilov, Senior Lecturer, Department of Aircraft, South Ural State University, Chelyabinsk, Russia; ismagilovdr@susu.ru.

Статья поступила в редакцию 16.09.2023; одобрена после рецензирования 10.11.2023; принята к публикации 10.11.2023.

The article was submitted 16.09.2023; approved after review 10.11.2023; accepted for publication 10.11.2023.

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