Численное моделирование потока в осесимметричном конфузоре на входе в первую ступень центробежного компрессора
Завацкая А. В., Данилишин А. М., Зуев А. В.
Выбор конструктивных параметров неподвижных элементов компрессора оказывает существенное влияние на его основные технико – экономические показатели – КПД, напор, массу и габариты. Основные требования, предъявляемые к проточной части входного устройства, сводятся к следующим: минимальные потери давления; равномерное поле скоростей и давлений в выходном сечении; минимальные осевые и радиальные размеры; технологичность конструкции[1].
В настоящей работе проведено исследование течения газа в осесимметричном криволинейном канале конфузора всасывающей камеры первой секции центробежного компрессора от сечения k—k до сечения 0-0 с целью определения оптимальной формы канала, обеспечивающей снижение потерь давления и равномерность потока на выходе.
Для оценки возможности совершенства указанного канала проведены расчеты течения газа при различных значениях конфузорности. Схема осесимметричного конфузора приведена на рис.1.
Геометрические параметры конфузора: , , , , , , коэффициент конфузорности , .Численное моделирование потока проводилось с помощью программного пакета ANSYS CFX 14.5[3]. В ходе работы была построена гексаэдрическая блочно-структурированная сетка, полученная после расчетов на сеточную независимость от числа элементов. Количество элементов сетки – 800000.
Для расчета была принята модель турбулентности SST – двухпараметрическая модель для сдвиговых напряжений. Она объединяет в себе две модели турбулентности: в пограничном слое используется модель , а в ядре потока — . Это обеспечивает получение качественных результатов как в ядре потока, так и вблизи стенок[4].
Исследовались конфузоры с различной степенью конфузорности. Расчеты проводились для совершенного газа – воздуха и для реального попутного нефтяного газа.
Полное давление в сечении k—k на входе в конфузор определялись двумя способами (формулы 1, 2).
(первый способ) | (1) |
(второй способ) | (2) |
В результате были получены значения коэффициента потерь во входном конфузоре и построены графики зависимости коэффициента потерь от коэффициента конфузорности.
Для анализа полученных результатов рассматривались коэффициенты сопротивления в круглой трубе в зависимости от числа Рейнольдса по данным Никурадзе и опытные данные о коэффициентах потерь в конфузоре компрессора для воздуха[5].
При определении числа Рейнольдса коэффициент кинематической вязкости для воздуха принимался равным , а для газа . В таблице 1 приведены результаты расчета для воздуха при определении потерь давления в конфузоре по параметрам в сечении 0-0 (первый способ). Значения коэффициента потерь давления в конфузоре, соответствующие опытным данным, в зависимости от параметров и [5], представлены в таблице 2.
Таблица 1. Результаты расчета для воздуха при определении потерь давления в конфузоре первым способом
nконф | Pk-k* | Po* | ρ | Co | Ϛ | Comax | Comin | ΔСo | Ck-k | Rek-k | Reо-о | λтр |
— | кПа | кПа | кг/м3 | м/c | — | м/c | м/c | — | м/c | -,
105 |
-,
105 |
— |
1,5 | 97,267 | 97,147 | 1,077 | 98,11 | 0,023 | 122 | 40,74 | 0,83 | 65,94 | 1,442 | 2,207 | 0,015308 |
1,83 | 97,255 | 97,148 | 1,077 | 97,94 | 0,021 | 122 | 40,74 | 0,83 | 55,22 | 1,484 | 2,204 | 0,015269 |
2,0 | 97,249 | 97,148 | 1,077 | 98,09 | 0,02 | 122 | 40,74 | 0,83 | 52,84 | 1,519 | 2,207 | 0,015229 |
Таблица 2. Значения коэффициентов потерь давления, соответствующие коэффициентам конфузорности, определяемые по данным [5]
nконф | Ϛ |
— | — |
1,5 | 0,048 |
1,83 | 0,047 |
2,0 | 0,045 |
Таблица 3. Результаты расчета для воздуха при определении потерь давления в конфузоре вторым способом
nконф | Pk-k* | Po* | ρ | Ск-к | Ϛ | Comax | Comin | Co | ΔСo | Rek-л | Reо-о | λтр |
— | кПа | кПа | кг/м3 | м/c | — | м/c | м/c | м/c | — | —,
105 |
-,
105 |
— |
1,5 | 97,266 | 97,146 | 1,077 | 65,94 | 0,051 | 122,2 | 40,74 | 98,108 | 0,83 | 1,442 | 2,2074 | 0,015308 |
1,83 | 97,249 | 97,142 | 1,077 | 55,22 | 0,065 | 122,2 | 40,74 | 98,098 | 0,83 | 1,484 | 2,2072 | 0,015265 |
2,0 | 97,2369 | 97,136 | 1,077 | 52,84 | 0,067 | 122,2 | 40,74 | 98,101 | 0,83 | 1,519 | 2,2073 | 0,015229 |
Как видно, в результате расчета в ANSYS коэффициенты потерь давления по первому способу имеют ту же тенденцию изменения в зависимости от конфузорности канала компрессора, что и в опытных данных [5], но меньше опытных значений коэффициентов потерь.
В таблице 3 приведены результаты расчета по программе ANSYS при определении коэффициента потерь по второму способу.
Полученные значения коэффициента потерь давления отличаются от полученных по первому способу не только по величине, но и по характеру изменения при увеличении конфузорности канала, но более близки по значениям к коэффициентам потерь давления, получаемым по данным [5].
Расчеты течения сжимаемого в компрессоре газа в конфузоре представлены ниже в таблицах 4, 5.
Таблица 4. Результаты расчета для газа при определении потерь давления в конфузоре первым способом
nконф | Pk-k* | Po* | ρ | Со-о | Ϛ | Comax | Comin | ΔСo | Ск-к | Rek-k | Reо-о | λтр |
— | МПа | МПа | кг/м3 | МПа | МПа | кг/м3 | м/c | — | м/c | —,
105 |
—,
105 |
— |
1,5 | 2,5846 | 2,5821 | 16,047 | 109,29 | 0,026 | 132 | 44,01 | 0,81 | 72,32 | 3,999 | 6,215 | 0,006656 |
1,83 | 2,5842 | 2,5818 | 16,049 | 109,03 | 0,025 | 132 | 44,01 | 0,81 | 60,36 | 4,101 | 6,201 | 0,006642 |
2 | 2,5841 | 2,5818 | 16,048 | 109,12 | 0,024 | 132 | 44,01 | 0,81 | 57,72 | 4,195 | 6,206 | 0,006626 |
Таблица 5. Результаты расчета для газа при определении потерь давления в конфузоре вторым способом
nконф | Pk-k* | Po* | ρ | Ск-к | Ϛ | Comax | Comin | Co | ΔСo | Rek-k | Reо-о | λтр |
— | МПа | МПа | кг/м3 | м/c | — | м/c | м/c | м/c | — | —,
105 |
—,
105 |
— |
1,5 | 2,5846 | 2,5821 | 16,345 | 72,32 | 0,059 | 132 | 44,01 | 109,29 | 0,81 | 3,999 | 6,215 | 0,006656 |
1,83 | 2,5841 | 2,5817 | 16,404 | 60,37 | 0,080 | 132 | 44,01 | 109,04 | 0,81 | 4,101 | 6,201 | 0,006641 |
2,0 | 2,5839 | 2,5816 | 16,413 | 57,73 | 0,084 | 132 | 44,01 | 109,13 | 0,81 | 4,195 | 6,206 | 0,006626 |
Как видно, характер изменения коэффициентов потерь давления при различной конфузорности канала для воздуха и для сжимаемого в компрессоре газа аналогичен, однако, значения коэффициента потерь давления в конфузоре для газа получаются выше, чем для воздуха. Можно заметить, что согласно кривым Никурадзе коэффициент трения на стенках канала для воздуха выше, чем для газа.
На следующем этапе будут проведены расчеты с использованием других моделей турбулентности. Это позволит сравнить полученные значения коэффициента потерь и оценить применимость каждой модели для решения данной задачи.
ЛИТЕРАТУРА:
- Селезнев К. П., Подобуев Ю. С., Анисимов С. А. Теория и расчет турбокомпрессоров. – Л.: Машиностроение. Ленинград, 1968. – 408 с.
- Завацкая А. В. Расчет и конструирование центробежного компрессора. Курсовой проект. СПбПУ КВиХТ. 2016. – 60 с.
- ANSYS CFX 14.5: Users Manual [Электронный ресурс] / ANSYS Inc., 2014.
- Любимов А. К., Шабарова Л. В. Методы построения расчетных сток в пакете ANSYS ICEM CFD [Электронный ресурс]: Электронное методическое пособие. – Нижний Новгород: Нижегородский госуниверситет, 2011. – 25 с.
- Идельчик И. Е. Справочник по гидравлическим сопротивлениям. М.: Машиностроение, 1975. – 104 с.