УДК 644.1
И.И. Полосин - доктор технических наук, профессор, заведующий кафедрой отопления и вентиляции
С.Н. Кузнецов - кандидат технических наук, доцент
А.В. Портянников - инженер
А.В. Дерепасов - инженер
Воронежский государственный архитектурно-строительный университет (ВГАСУ)
ПОСТРОЕНИЕ ЭФФЕКТИВНОГО ВОЗДУХООБМЕНА ДЛЯ ПОМЕЩЕНИЙ ПРОИЗВОДСТВ С МЕСТНОЙ ВЫТЯЖНОЙ ВЕНТИЛЯЦИЕЙ С ИСПОЛЬЗОВАНИЕМ МАТЕМАТИЧЕСКОЙ МОДЕЛИ
АННОТАЦИЯ
Рассмотрена двухмерная модель распространения стационарных воздушных потоков в помещении, основанная на системе стационарных уравнений неразрывности и Навье-Стокса с использованием стандартной к- e модели турбулентности. Математическая модель и алгоритмы ее решения реализованы на ЭВМ в виде программы. Проанализированы результаты расчетов воздухообменов в помещении с местной вытяжной вентиляцией и выбрана наиболее эффективная схема воздухообмена.
КЛЮЧЕВЫЕ СЛОВА: Вентиляция, вредные вещества, распределение воздуха, математическая модель.
I.I. Polosin - doctor of technical sciences, professor, head of Heating and Ventilation department
S.N. Kuznetsov - candidate of technical sciences, associate professor
A.V. Portjannikov - engineer
A.V. Derepasov - engineer
Voronezh State University of Architecture and Civil Engineering (VSUACE)
CONSTRUCTION OF EFFECTIVE AIR EXCHANGE FOR FACTORY PREMISES WITH LOCAL EXHAUST VENTILATION USING MATHEMATICAL MODEL
ABSTRACT
A two-dimensional model of distribution of stationary air streams in premise based on a system of steady-state equations of continuity and Navier-Stokes using a standard к- e turbulence model has been considered. A mathematical model and algorithms of its solution have been implemented as a computer program. Results of calculations of air exchange in premise with local exhaust ventilation have been analyzed and the most effective air exchange scheme has been selected.
KEYWORDS: Ventilation, harmful substances, distribution of air, mathematical model.
Разработка энергоэффективных схем организации Уравнение неразрывности:
воздухообмена в гальванических цехах требует учета
э
характера движения воздуха в помещении, который ____(ри ) = 0, (1)
зависит от выбранного типа воздухораспределения. Эхг-
Стоимость и трудоемкость исследований
воздухообмена в натурных условиях или на воздушно- где р — пл°тн°стъ кг/м3; X. — /-я
тепловых моделях очень высоки, поэтому в настоящее пространственная координата, м; и. — компоненты
время для решения этой задачи все чаще используется скорости течения воздуха, м/с.
метод математического моделирования [5]. Уравнения Навье-Стокса'
Математическая модель вентиляционного процесса
распространения воздуха основана на уравнении э э Э
неразрывности, системе осредненных по Рейнольдсу -------(ри ¡и .) = —Э— +--
стационарных уравнений Навье-Стокса и уравнениях Эх. Эх1 Эх.
к-е модели турбулентности [1]. Рассмотрим основные уравнения модели.
э (r__ )_ dp + э (Эщ_ +эЦ£ _ 2 ^ дщ_
Эх. dxt 3 у Эх,.
V jj
№
Эх,.
р uuj _m,
Эи, + duj Эх, Эх,.
у J , у
3
pk
Эх,
(2)
5„, (3)
где р - давление, Па; т - динамическая вязкость, кг/мс; т - турбулентная динамическая вязкость кг/мс; g - ускорение свободного падения, м/с2; к -кинетическая энергия турбулентности, м2/с2; ] = 1, 2, 3;
8у - символ Кронекера.
Перенос кинетической энергии турбулентности в стандартной к-е модели определяется уравнениями:
Эх (рЧ )_ЭТ
Эх Эх
Э
VV
Эк
Эх
у J 0
+ Gk + Gb _pe, (4)
Э ( ) Э
—(pe_j )_ — Эх, Эх,
VV
m+— ае
Л Эе Л
Эх,
+ рСДе_
у j у
е2 е
:_pC2-^r= + Cie f C3eGb, к + VM-e к
(5)
где е - скорость диссипации энергии турбулентности, м2/с3; С - константы модели турбулентности; т, - турбулентная динамическая вязкость, кг/мс; ок и ое - турбулентные числа Прантля для к и е.
Скорость образования турбулентности:
Gk _ 2m,
I
2
Эх
\ Ґ +m
I
i> j
V
Эи. Эи,
—L + —J Эх. Эх.
J 1
2
Gb _mg1 Эт
p Эх3
Величина C определяется выражением:
(6)
(7)
C1 _ даах
0,43,
h
h + 5
<?к
Л = Р^-.
е
Турбулентная динамическая вязкость определяется уравнением:
C к2 m_ pQ—.
Константа Ст определяется выражением:
Cm_-
(8)
А + A
kU*
где
U _,/SS +WjW,
Wj _ W ij 2ej ,
Wj _ 1
iJ 2
Эи, Эи
г г
Эх, Эх,
Константы А0 и As определяются как: Ао _ 4,04,
As _ л/б cos j,
где
j_ 3 cos 1 (f6W),
W-
SjSjkSki
S _ 1 V 2
Эи, Эи.
____j_ +____i_
Эх, Эх,
J ,
Для констант к-е модели турбулентности принимались следующие значения [1]:
Си = 1,44; Сз£ = 1,40; С2 = 1,9; ак = 1,0; ае = 1,2.
Рассматриваются два типа границ: твердые стенки и свободные границы. Скорость и кинетическая энергия турбулентности на границе равны нулю, для давления выполняются условия Неймана. Коэффициент диссипации турбулентной энергии определяется непосредственно функциями пограничного слоя.
На свободной границе давление равно атмосферному, нормальные к свободной поверхности градиенты компонент скорости равны нулю. Граничными условиями для входящего потока будут условия Дирихле со значениями переменных в окружающей среде, а для выходящего потока - условия Неймана.
Для решения системы уравнений (1-8) использован метод конечных разностей с предварительным расщеплением их по физическим процессам [2] и использованием метода фиктивных областей. Подобные методы решения подробно рассмотрены в [2, 3, 4]. Решение строится с помощью двухциклической схемы покомпонентного расщепления с использованием неявной монотонной разностной схемы второго порядка.
1
е
Полученные алгоритмы решения уравнений (1-8) реализованы на ЭВМ в виде программы на языке С++.
В настоящее время уже известны подобные программные продукты, разработанные в СПбГАСУ для численного моделирования аэрации корпуса электролиза алюминия, позволяющие установить характер движения воздуха в сечении модуля корпуса. Разработанная нами математическая модель и программа, в отличие от известных, позволяют получить устойчивое решение в зонах больших градиентов скоростей течений воздуха вблизи бортовых отсосов. При помощи разработанной программы изучалось распределение двухмерных стационарных воздушных потоков в поперечном сечении цеха с местными отсосами с целью определения влияния способа подачи приточного воздуха на циркуляцию воздушных потоков.
В качестве объекта исследования рассмотрено двухмерное сечение цеха бесконечной длины с промышленными ваннами наружных размеров: ширина - 1100 мм и высота - 1000 мм. Ванны размещались в два ряда в поперечном сечении цеха и оборудованы двухсторонним бортовым отсосом. Помещение цеха имеет ширину 18,0 м и высоту 12,0 м. Расчетная кратность воздухообмена составляла 10 ч-1. Температура ограждающих конструкций помещения и приточного воздуха принята равной 18 °С.
Исследовалось три варианта подачи приточного воздуха на высоте 0,8к (9,6 м) от пола:
- плоскими струями с подачей воздуха вертикально вниз;
- плоскими струями с использованием перфорированного воздуховода круглого сечения, с углом раскрытия струи 45° и подачей воздуха вниз;
- плоскими струями с подачей воздуха вертикально вниз по центру помещения и под углом 20° к вертикали по сторонам помещения.
Схемы подачи приточного воздуха приведены на рис. 1.
Рис. 1. Схемы подачи приточного воздуха:
1— плоская струя; 2, 3 — плоская струя с подачей воздуха под углом 20° к вертикали;
4 — подача через перфорированный воздуховод
Удаление загрязненного воздуха из помещения цеха осуществлялось бортовыми отсосами от ванн в объеме 90 % от общего количества приточного воздуха и 10 % из верхней зоны через отверстие в перекрытии цеха (см. рис. 2).
Рис. 2. Схема помещения. Стрелками показаны места подачи и удаления воздуха: 1-1, 2-2, 3-3 — сечения, в которых исследовалась вертикальная составляющая скорости воздуха
Анализ полученных результатов, приведенных на рис. 3-4, позволяет сделать ряд выводов. Поток воздуха от воздухораспределителей движется вниз, взаимодействует с зонами циркуляции и ограждающими конструкциями, отражается от пола, разворачивается и, отдавая часть воздуха местным отсосам, поднимается вверх. Во всех вариантах по обе стороны от нисходящих приточных струй образуются зоны циркуляции. Размер зон циркуляции по горизонтали примерно равен половине расстояния между воздухораспределителями, а по вертикали — высоте установки воздухораспределителей.
Вдоль стен помещения во всех случаях наблюдается низкоскоростной восходящий воздушный поток, направленный к местам удаления воздуха в верхней зоне помещения.
Сравнение двух первых вариантов подачи приточного воздуха показало, что использование плоских струй с подачей воздуха вертикально вниз даёт схожую картину распределения воздушных потоков и скорости движения воздуха, как в объеме помещения, так и в рабочей зоне.
В реальном помещении струи приточного воздуха от воздухораспределителей находятся в условиях стеснения строительными конструкциями и технологическим оборудованием. Воздух вокруг струй находится не в состоянии «покоя», а движется навстречу, в результате чего профиль скоростей струи сильно деформируется. Скорость воздуха на оси струй зависит от типа воздухораспределителя на расстояниях до 4,6 м от воздухораспределителя. При больших расстояниях от воздухораспределителя скорость воздуха на оси струй слабо зависит от типа воздухораспределителя.
Местные отсосы сильно изменяют направление движения струи. Так, при расстоянии от воздухораспределителя 5,6 м отклонение оси струи может достигать 0,7 м, а в рабочей зоне — 2 м. В
Рис. 3. Подача воздуха плоскими струями вертикально вниз: а) векторное поле скоростей воздушных потоков;
б) вертикальная составляющая скорости воздуха:__________- на высоте 8 м от уровня пола (сечение 1-1);
_ _ - на высоте 6 м от уровня пола (сечение 2-2);.......- на высоте 4 м от уровня пола (сечение 3-3)
Рис. 5. Подача приточного воздуха плоскими струями вертикально вниз по центру помещения и под углом 20° к вертикали по сторонам помещения: обозначения те же, что и на рис. 3
результате чистый воздух поступает не в рабочую зону, а к местным отсосам и удаляется ими.
Кроме 1 и 2 вариантов, исследовались также другие способы подачи приточного воздуха. Наиболее эффективным оказался вариант с подачей приточного воздуха плоскими струями вертикально вниз по центру помещения и под углом 20° к вертикали по сторонам помещения (см. рис. 5).
Чистый воздух подается на рабочие места и только затем поступает в зону действия местных отсосов. Восходящие потоки находятся над ваннами и выносят вредные вещества из рабочей зоны в верхнюю зону помещения.
Таким образом, установлено, что при взаимодействии со строительными конструкциями и местными отсосами происходит сильная деформация струй. Видимо, корректнее говорить не о струях, а о воздушных потоках.
Для исследования концентрации вредных веществ в рабочей зоне важно знать, как в рабочую зону поступает чистый воздух и как удаляется из нее загрязненный. Первые два варианта создают неэффективный воздухообмен, т.к. чистый воздух удаляется местными отсосами и не поступает в рабочую зону. В третьем случае чистый воздух поступает в рабочую зону и только затем удаляется местными отсосами. Восходящие потоки находятся над ваннами и выносят вредные вещества из рабочей зоны в верхнюю зону помещения.
Предложенный подход позволяет рассчитать и выбрать систему воздушных потоков в помещениях, наиболее эффективно использующую приточный воздух и повысить эффективность действиявентиляции.
Литература
1. Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. —М.: Энергоатомиздат, 1984. — 152 с.
2. Марчук Г.И. Методы вычислительной математики.
— М.: Наука, 1977. — 456 с.
3. Иванов М.Я., Крупа В.Г., Нигматулин Р.З. Неявная схема С. К. Годунова повышенной точности для интегрирования уравнений Навье-Стокса // Журнал вычислит. математики и математич. физики, 1989, т. 29, № 6. — С. 888-901.
4. Мелькумов В.Н., Кузнецов И.С. Расчет нестационарного поля концентраций двухкомпонентной газовой смеси // Вестник ВГТУ. Сер. Энергетика, 2006, т. 2, № 6. — С. 125-127.
5. Посохин В.Н., Зиганшин А.М. Численное моделирование конвективных течений над протяженными теплоисточниками // Известия КазГАСУ 2005, № 1 (3). — С. 136-140.