Определение основных характеристик форсунки камеры смешения кондиционера с использованием CFD-моделирования

Полный текст

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

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

В работе [6] представлено экспериментальное исследование ряда центробежных форсунок оросительных камер промышленных кондиционеров (рис. 1).

 

Рис. 1. Центробежная форсунка: 1 ‒ входной патрубок; 2 ‒ шнековый завихритель; 3 ‒ камера закручивания; 4 ‒ выходное сопло

Fig. 1. Centrifugal nozzle: 1 ‒ inlet nozzle; 2 ‒ screw swirler; 3 ‒ twist chamber; 4 ‒ outlet nozzle

 

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

• радиус сопла – 2,25 мм;
• средний радиус захода шнека – 12,00 мм;
• количество заходов шнека – 3;
• ширина винтового канала – 3,80 мм;
• глубина винтового канала – 4,00 мм;
• угол наклона шнека – 75°.

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

 

Рис. 2. Геометрия расчетной модели

Fig. 2. Simulation Model Geometry

 

На рис. 3 показана сеточная модель расчетной области, объем сетки составляет 1,5 млн. ячеек. Параметры сетки модели следующие:

• размерность модели – трехмерная;
• тип сетки ‒ гексаэдральная с многогранными призматическими слоями, соединенными через полиэдральный буферный слой;
• количество ячеек ‒ 1,5 млн. шт.;
• высота первого пристеночного призматического слоя – 0,05 мм;
• минимальное количество пристеночных призматических слоев – 5 шт.;
• максимальный коэффициент роста высоты пристеночных слоев – 1,2;
• минимальный размер ячейки – 0,25 мм;
• максимальный размер ячейки – 2 мм;
• ортогональное качество ячеек:
  
  • минимальное – 0,15;
  • максимальное – 1,00.

 

Рис. 3. Сеточная модель расчетной области

Fig. 3. Grid model of the calculation area

 

Для определения основных характеристик центробежной форсунки были проведены численные исследования в программном продукте ANSYS Fluent [7].

Для моделирования движения жидкости в камере смешения форсунки применялся метод объема жидкости (метод Эйлера) Volume of Fluid (VOF) [8]. В качестве параметра, отражающего картину двухфазного течения, здесь применяется параметр объемной доли жидкости в ячейке С (если С = 1, то ячейка заполнена жидкостью, если С = 0, то ячейка заполнена газом). Граница раздела ‒ ячейка с параметром С = 0,5.

При расчете использована модель турбулентности «k-omega SST». Расчет проводился с использованием пространственной дискретизации Compressive, которая обеспечивает наибольшую устойчивость расчета и сходимость по балансам масс и локальным значениям параметров потока. Сходимость расчета определялась по уровням невязок, балансам масс на входе и выходе модели, локальным значениям параметров потока.

Результаты представлены на рис. 4‒8. На рис. 4 показано распределение полного давления в миделевом сечении расчетной модели; на рис. 5‒7 ‒ распределения абсолютной, осевой и тангенциальной скоростей в миделевом сечении расчетной модели.

 

Рис. 4. Распределение полного давления

Fig. 4. Total pressure distribution

 

Рис. 5. Распределение абсолютной скорости потока

Fig. 5. Absolute flow rate distribution

 

Рис. 6. Распределение осевой составляющей скорости потока

Fig. 6. Distribution of axial component of flow velocity

 

Рис. 7. Распределение тангенциальной составляющей скорости потока

Fig. 7. Distribution of tangential component of flow velocity

 

Рис. 8. Распределение относительной концентрации жидкости

Fig. 8. Distribution of relative fluid concentration

 

На рис. 8 показано распределение относительной концентрации жидкости в миделевом сечении расчетной модели. Угол раскрытия жидкостного факела расчетной модели составляет 76,6° для всех режимов. В случае натурного эксперимента изменение угла раскрытия жидкостного факела составляет 75‒78° в рассмотренном диапазоне режимов работы центробежной форсунки. Из сравнения результатов расчетного и натурного экспериментов можно заключить, что CFD-моделирование с высокой степенью достоверности отражает течение в исследуемой форсунке.

Расходная характеристика также является одной из основных для форсунки. В частности, она позволяет определить зависимость коэффициента расхода форсунки μ_Ф от режимных параметров. Расход воды G_В в натурном эксперименте определен весовым методом [6], в численном эксперименте этот параметр определен в постпроцессоре ANSYS Fluent.

Коэффициент расхода форсунки определяется соотношением

${\mu }_{\mathrm{Ф}}=\frac{G_{\mathrm{В}}}{G_{\mathrm{В}. \mathrm{ИД}.}} ,$

где $G_{\mathrm{В}.\mathrm{ИД}.}=F_C\sqrt{2\cdot {\rho }_{\mathrm{В}}\cdot \mathrm{\Delta }P}$ ‒ идеальный расход воды через форсунку,

здесь F_C – площадь сопла форсунки; ρ_В – плотность воды; ΔP – перепад давления на форсунке.

На рис. 9 показана расходная характеристика исследуемого варианта форсунки, вычисленная весовым методом в натурном эксперименте [6], и расходная характеристика форсунки из CFD-моделирования.

 

Рис. 9. Расходная характеристика форсунки

Fig. 9. Nozzle flow characteristic

 

Результаты определения значений идеального расхода воды и коэффициента расхода форсунки представлены ниже. В табл. 1 приведены результаты определения коэффициента расхода форсунки в натурном эксперименте [6], в табл. 2 ‒ результаты определения коэффициента расхода форсунки в CFD-моделировании.

 

Таблица 1. Значения режимных параметров форсунки в натурном эксперименте

Table 1. Values of operating parameters of the nozzle in the field experiment

ΔP, МПа	GВ , кг/с	GВ.ИД. , кг/с	μФ
0,098	0,092	0,223	0,413
0,196	0,102	0,315	0,324
0,294	0,131	0,386	0,340
0,304	0,128	0,392	0,326

 

Таблица 2. Значения режимных параметров форсунки в численном эксперименте

Table 2. Values of operating parameters of nozzle in numerical experiment

ΔP, МПа	GВ , кг/с	GВ.ИД. , кг/с	μФ
0,051	0,074	0,160	0,465
0,101	0,086	0,226	0,378
0,203	0,104	0,320	0,325
0,304	0,120	0,392	0,306

 

Сравнение результатов натурного и численного экспериментов приведено на рис. 10. Заметна хорошая сходимость результатов натурного и численного экспериментов.

 

Рис. 10. Коэффициент расхода форсунки в натурном и численном экспериментах

Fig. 10. Nozzle flow coefficient in full-scale and numerical experiments

 

Выводы. Проведено исследование характеристик течения жидкости (воды) через центробежную форсунку оросительной камеры промышленных кондиционеров с использованием численного моделирования. По результатам работы можно заключить, что для угла раскрытия жидкостного факела расхождение результатов натурного и численного экспериментов не превышает 2,1 %. Результаты определения коэффициента расхода форсунки натурного и численного экспериментов, как видно из проведенного исследования, соответствуют как качественно, так и количественно. Следовательно, CFD-моделирование обладает хорошей прогностической способностью при проектировании центробежных форсунок оросительной камеры (камеры смешения) промышленных кондиционеров.

Об авторах

Данила Владимирович Зеленцов

Самарский государственный технический университет

Автор, ответственный за переписку.
Email: tgv@samgtu.ru

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

Россия, 443100, г. Самара, ул. Молодогвардейская, 244

Сергей Всеволодович Крашенинников

Самарский государственный технический университет

Email: tgv@samgtu.ru

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

Россия, 443100, г. Самара, ул. Молодогвардейская, 244

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

Дополнительные файлы