Научно-инжиниринговая группа

Компрессорная, вакуумная, холодильная техника и системы транспорта и переработки газа

Создана специалистами кафедры "Компрессорная, вакуумная и холодильная техника" и Научной школы компрессоростроения Санкт-Петербургского политехнического университета Петра Великого

Численное моделирование потока в осесимметричном конфузоре на входе в первую ступень центробежного компрессора

Завацкая А. В., Данилишин А. М., Зуев А. В.

Выбор конструктивных параметров неподвижных элементов компрессора оказывает существенное влияние на его основные технико – экономические показатели – КПД, напор, массу и габариты. Основные требования, предъявляемые к проточной части входного устройства, сводятся к следующим: минимальные потери давления; равномерное поле скоростей и давлений в выходном сечении; минимальные осевые и радиальные размеры; технологичность конструкции[1].

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

Для оценки возможности совершенства указанного канала проведены расчеты течения газа при различных значениях конфузорности. Схема осесимметричного конфузора приведена на рис.1.

Рис.1. Схема осесимметричного конфузора[2]

Рис.1. Схема осесимметричного конфузора[2]

Геометрические параметры конфузора: , , , , , , коэффициент конфузорности  , .

Численное моделирование потока проводилось с помощью программного пакета ANSYS CFX 14.5[3]. В ходе работы была построена гексаэдрическая блочно-структурированная сетка, полученная после расчетов на сеточную независимость от числа элементов. Количество элементов сетки – 800000.

Для расчета была принята модель турбулентности SST – двухпараметрическая модель для сдвиговых напряжений. Она объединяет в себе две модели турбулентности: в пограничном слое используется модель , а в ядре потока — . Это обеспечивает получение качественных результатов как в ядре потока, так и вблизи стенок[4].

Исследовались конфузоры с различной степенью конфузорности. Расчеты проводились для совершенного газа – воздуха и для реального попутного нефтяного газа.

Полное давление в сечении kk на входе в конфузор определялись двумя способами (формулы 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

 

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

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

 

ЛИТЕРАТУРА:

  1. Селезнев К. П., Подобуев Ю. С., Анисимов С. А. Теория и расчет турбокомпрессоров. – Л.: Машиностроение. Ленинград, 1968. – 408 с.
  2. Завацкая А. В. Расчет и конструирование центробежного компрессора. Курсовой проект. СПбПУ КВиХТ. 2016. – 60 с.
  3. ANSYS CFX 14.5: Users Manual [Электронный ресурс] / ANSYS Inc., 2014.
  4. Любимов А. К., Шабарова Л. В. Методы построения расчетных сток в пакете ANSYS ICEM CFD [Электронный ресурс]: Электронное методическое пособие. – Нижний Новгород: Нижегородский госуниверситет, 2011. – 25 с.
  5. Идельчик И. Е. Справочник по гидравлическим сопротивлениям. М.: Машиностроение, 1975. – 104 с.