CC BY
69
УДК 613.6:621.43
МАТЕМАТИЧЕСКОЕ И ЭКСПЕРИМЕНТАЛЬНОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ РАСПРОСТРАНЕНИЯ ОКСИДОВ УГЛЕРОДА И ИЗБЫТКОВ ТЕПЛОТЫ В ГАЗОВОЗДУШНОЙ СРЕДЕ ПОМЕЩЕНИЯ
Б.Ч. МЕСХИ, Е.И. МАСЛОВ, А.Н. СОЛОВЬЁВ, Ю.И. БУЛЫГИН, Д.А. КОРОНЧИК
(Донской государственный технический университет)
Представлены результаты математического и экспериментального моделирования процессов распространения угарного газа и избытков теплоты в воздушной среде гаражного помещения. На основе теоретикоэкспериментального метода получены адекватные экспериментальным данным поля подвижности воздуха, концентрации оксидов углерода и температуры внутри исследуемого помещения.
Ключевые слова: загазованность, избытки теплоты, концентрация, температура, вентиляция.
Введение. Снижение эмиссии вредных веществ (ВВ) в рабочую зону операторов участков обкатки машиностроительных предприятий является актуальной задачей [1]. С другой стороны, выбросы при горячей обкатке машин «тёплые», что приводит к образованию в загазованной рабочей зоне областей с неудовлетворительными параметрами микроклимата. В результате формируется целый комплекс опасных и вредных производственных факторов (ОВПФ), таких как загазованность, повышенная подвижность воздуха и температура. Каждый из ОВПФ необходимо учитывать, так как они оказывают определённое влияние друг на друга, что в конечном итоге сказывается на точности моделирования.
Постановка задачи. Математическое моделирование процессов распространения вредных веществ в воздушной среде вентилируемых помещений связано с рядом трудностей.
Во-первых, результаты физического моделирования исследуемых процессов позволяют уточнить как сами модельные уравнения, так и их начальные и граничные условия.
Во-вторых, доминирующее влияние на формирование полей загазованности, подвижности воздуха и температуры оказывают коэффициенты турбулентной диффузии (обмена) и теплопроводности газовой смеси в помещении, входящие в систему дифференциальных уравнений. В системе используются дифференциальные уравнения, описывающие процессы конвективно-диффузионного тепломассопереноса, введенные в работах М.Я. Берлянда [2], Г.И. Марчука [3], С.В. Патанкара [4].
Поэтому вполне обоснованной выглядит поставленная в работе задача - исследовать поведение математической модели, описывающей распространение угарного газа и избытков теплоты в воздушной среде заранее выбранного модельного объекта с простой геометрической формой и известным внутри стационарным источником загрязнений.
В качестве такого объекта - газовой камеры - был выбран типовой металлический гараж. Источником загрязнений являлся автомобиль ВАЗ 2104, отработанные газы от которого были отведены в центр гаражного помещения.
Математическая модель тепломассопереноса вредных веществ и избытков явной теплоты.
Система уравнений в общем виде
С
а
-+
аС-с, ),а( V. с, = *ла^+д )^|+
дх
+
+
ду
ду
дС,Л
ду
дР
дz
+д'(( ап-+Д )^1+£в' ) •*(" - г )■
дс
(1)
дР дР
и = —, V = —, ^ = — дх ду дz
д2Р д2 Р д2 Р п
+----г- +---Г- = 0,
дх2
С д +— ду
ду дz
дТ 1 д(и • Т) 1 д(у• Т) | д((* + ws)• Т)^__^
дх
(2)
+
дх
+
ду
+
дz
(А + А^ + А р)
дТ_
дх
+
[^(а + а + а р^+д-.^(а + а, + а р)~ ^+^ () ^( г- Г]),
где
А
га лп
(3)
(4)
Ci - концентрация /-го ВВ в помещении, мг/м3; t - текущее время, с; и, V, ш - компоненты вектора скорости воздушной среды, м/с; - скорость оседания/подъёма примеси/газов; АП = (АХ, АУ, А2) - коэффициенты турбулентной диффузии в помещении, м2/с;
D = фХ, DУ, DZ) - коэффициенты диффузии /-го газа, м2/с; 2/ (I) - интенсивность выброса /-го вредного вещества от у-го источника в помещении, кг/с; 5(г - г) - дельта-функция Дирака; г. = (ху, у., zj) - координаты источника выброса, м; р - плотность газовой смеси, кг/м3; Ср - удельная изобарная теплоемкость, Дж/(кг ■ К); А - коэффициент теплопроводности /-го газа, Вт/(м^К); А( - коэффициент турбулентной теплопроводности, Вт/(м^К); А - коэффициент радиационной теплопроводности, Вт/(м^К); д. - интенсивность внут-
ренних источников тепла, Вт; Рга - диффузионное число Прандтля, Рн - турбулентное число Прандтля.
Граничные условия для уравнения (1):
- на боковых стенках:
(5)
в последнем равенстве указывается, что стенки помещения, в котором проводятся исследования, являются непроницаемыми для концентрации загрязняющего вещества;
- на входной и выходной границе:
дС=-л. • с
дп Ап
д
[и ,и >0
\ п ’ п ■■ _
где q = ^ , это означает, что поток примеси зависит от ее концентрации, скорости ип и
t0, Un ^ 0*
направления воздушного потока, индуцированного работой вентиляции.
Граничные условия для уравнения (2):
- на твердых стенках
8P
= 0, (7)
сп
где n - единичный вектор внешней нормали. Это означает, что стенки исследуемого помещения являются непроницаемыми для воздушных потоков;
- на выходной границе
CP
— Ч,, (8)
сп
где ип - известное значение скорости;
- на входной границе (границе втекания воздушного потока в помещение)
P = const (9)
(условие Дирихле).
Граничные условия для уравнения (3): начальные условия заключаются в задании распределения поля значений температуры в начальный момент времени (t = to), т. е. предшествующий расчетному. Они должны быть заданы в виде функции
Т = f (^ У, z) - (10)
для пространственной задачи;
- на боковых стенках
f=T'K-T-)- (11)
данное уравнение обозначает, что со стенками помещения ведется теплообмен. Здесь должны быть заданы коэффициент теплоотдачи а и температура окружающей тело среды Тс;
- на входной границе (границе втекания воздушного потока в помещение)
Т = const - (12)
известное и постоянное значение температуры втекаемого в помещение потока воздуха;
- на выходной границе — = -— -(Тп), что означает, что поток температуры примеси за-
dn Xt п
висит от ее температуры, скорости un и направления воздушного потока, индуцированного работой вентиляции.
В ходе экспериментальных работ граничные условия уточняются.
Коэффициенты турбулентного обмена Ап. При турбулентном течении потока и турбулентном переносе коэффициенты Ап пропорциональны диссипируемой энергии и определяющему размеру рассматриваемого источника. Для вентилируемых помещений, зная величину энергии, поступающей в помещение, и определяющий размер источника, можно определить коэффициент турбулентного обмена в плоскости (x, y), м2/с:
AX = Ay = 0,25-в^ • 1П4з, (13)
где в - кинетическая энергия воздушных потоков в помещении; 1П - определяющий размер помещения, равный для прямоугольных отверстий и сечений неправильной формы
П =^, м, (14)
где F- площадь сечения помещения, перпендикулярная направлению движения воздушных потоков.
Коэффициент турбулентного обмена по высоте помещения находится из зависимости, предложенной М. Я. Берляндом для наибольших классов устойчивости воздушной среды:
Аш = А=1 ' (2І'2і) , (15)
что допустимо для исследуемого помещения. Здесь Аг=1 - коэффициент турбулентного обмена на
высоте 1 м, г1 = 1 м, г- текущая высота помещения, м.
Кинетическая энергия воздушных потоков в помещении в определяется из формулы
8 = вПС + ВТС + вДП ~ вВС , (16)
куда входят соответственно энергии приточных струй вПС, тепловых струй вТС, движущихся
п и вытяжных струй вВС . Величины вПС ,
предметов в ния ВВ и находятся из формул:
вТС, вДП учитываются в модели рассеива-
2
g - Q - г 1 + п
С -Т-р.
4
в ДП =
1 <Яч а < а Р 1 X
1_ 2 - Си ] _3600_
(17)
(18) (19)
где Ь/УП - кратность воздухообмена в помещении свободным объёмом УП; V - средняя скорость выхода воздуха из приточных отверстий, м/с; £, - коэффициент местного сопротивления на выходе из воздухораспределителей; Q - теплонапряжённость объема, Вт/м3; д = 9,8 м/с2; z - расстояние от полюса до рассматриваемой точки, м; п = 0,8 - экспериментальная константа; Ср - теплоемкость воздуха на удалении от источника, ; Т0 - температура
воздуха в помещении на удалении от источника; рв - плотность воздуха на удалении от источника, £, - коэффициент аэродинамического сопротивления движущегося предмета; Fn -площадь поперечного сечения движущегося предмета; \п - скорость движущегося предмета; х - средняя продолжительность движения предмета с данной скоростью в течение 1 ч; Gn - масса воздуха в объеме помещения; так как Gn = У -р, где V- объем помещения, то:
в ДП =
1 <Яч а < а 1 X
2 - V _3600_
(20)
Характеристика источника выделения вредностей
В уравнение (1) входит Q/ ^) - интенсивность выброса /-го вредного вещества от /-го источника. В нашем случае источник загрязнения - это двигатель автомобиля ВАЗ 2104, который работает на холостом ходу. При этом с отработанными газами через шланговый отсос в воздух гаража поступают токсичные компоненты. Как показано в работе [5], одним из основных токсич-
вТС
ных компонентов, который учитывается при расчёте потребных воздухообменов на вентиляцию, является оксид углерода.
Для определения его массового выброса необходимо экспериментально определить концентрацию СО в отработавших газах (ОГ) источника и измерить расход отработавших газов автомобиля. Для одиночного стационарного источника функция Q/(t) может быть определена по формуле
О.(О = СД ' °ОГсрД , г/с, (21)
^/У> 3600
где С д - концентрация ВВ в ОГ от /-го источника загрязнения, г/м3; ООГсрд - расход отработавших газов двигателя автомобиля на холостом ходу, м3/с.
Явное теплопоступление в исследуемое помещение учитывается выражением
О. = срОг - ООГ - (ТОГ - То) - Ад, кВт, (22)
ХХ ХХ
где СрОГ _ изобарная теплоемкость ОГ на режиме холостого хода, Дж/кг ■ К; Тог - температура ОГ на холостом ходу, К; Т0 - температура газовоздушной среды в помещении, К;
Ад - количество источников.
Экспериментальные исследования. С целью обоснования адекватности разработанной математической модели распространения угарного газа, избытков теплоты и подвижности воздуха внутри заранее выбранного модельного объекта с простой геометрической формой и известным стационарным источником загрязнения были проведены экспериментальные исследования.
Одновременно в заданные моменты времени замерялись концентрации угарного газа, температура газовоздушной смеси и скорость движения воздуха непосредственно в разных точках объёма исследуемой газовой камеры. Для оценки пространственного распределения вредных веществ по помещению применялся стационарный метод отбора проб (функции устройств для отбора и оперативного анализа проб выполняли электрохимические датчики газового анализа, которые с прибором размещались в фиксированных точках помещения).
Натурные исследования процессов проводились в гаражном помещении стандартных размеров (длина - 6,40 м, ширина - 3,30 м, высота - 2,30 м), схема которого (вид сверху) представлена на рис. 1.
Рис. 1. Схема размещения измерительных точек в гаражном помещении
Измерения концентраций и температур газовоздушной среды за разные промежутки времени проводились в пяти заданных точках камеры, где указаны их геометрические координаты: «а» (1.10,1.10,1.5), «Ь» (1.10,2.20,1.5), «с» (3.20,1.65,1.5), «й» (5.20,1.10,1.5), «е» (5.20,2.20,1.5). Геометрические параметры источника загрязнения соответствуют координатам точки «с» - центр помещения (см. рис. 1).
Согласно разработанной «Программе испытаний», замеры концентраций угарного газа в воздухе газовой камеры осуществлялись на высоте 1,5 м в характерных точках, представленных на схеме исследуемого объекта (см. рис. 1).
Измерения концентраций угарного газа в воздухе газовой камеры проводились в холодный период года (ноябрь) в течение часа по следующему графику (табл. 1).
Таблица 1
Последовательность взятия проб в течение часа
№ п/п Номер пробы на газовый анализ Время, прошедшее с начала эксперимента, с
1 I 0
2 II 175
3 III 325
4 VI 475
5 V 625
6 VI 1260
7 VII 2100
8 VIII 2700
9 IX 3300
10 X 3600
Для полноты исследования физической картины образования и распространения температурных полей и полей концентраций ВВ в камере испытания проводились для трёх характерных случаев: 1) закрытые ворота (вентиляционные отверстия заклеены); 2) открыта одна створка ворот; 3) ворота приоткрыты (щель размером 0,47 м).
Приборное оснащение эксперимента. Результаты измерений
Для измерения концентраций ВВ применялся газоанализатор Drager X-am 5000 (Германия): метод измерения - электрохимический; диапазон показаний - CO: 0-2000 ppm, NO2: 0-50 ppm, SO2: 0-50 ppm (ррт - мг/м3); предел допускаемой основной погрешности - 5%.
Микроклиматические условия в помещении (температура и подвижность воздуха) измерялись согласно схеме, представленной на рис. 1, на высоте 1,5 м. Для измерений скорости движения воздушных потоков и температуры газовоздушной среды использовался высокоточный дифференциальный манометр Testo 521. Он имеет внутренний сенсор с измерительным диапазоном от 0,1 до 100 ГПа и является прибором для измерений скорости потока от 0,1 до 100 м/с. Применяется совместно с трубкой Пито. Модель Testo 521-1 со встроенным сенсором давления имеет погрешность 0,2% от полной шкалы, а модель Testo 521-2 со встроенным сенсором давления имеет погрешность 0,1% от полной шкалы (рис. 2).
О О
=Г (ТЗ о_ 1— X си / // // Точка а Точка b — — Точка с
X о /у »/ I/ I/ #/ Точка d • Точка е
15 175 325 475 625 1260 2100 2700 ЗЗОО 3600 Время, с
Рис. 2. Динамика изменения концентраций СО в газовой камере (точки а-е см. рис. 1)
Измерения температуры производились в холодный период года внутри помещения на высоте 1,5 м. Использовались ртутные термометры с точностью измерений до 0,1 °С. Зависимости температуры газовоздушной смеси в камере в заданных измерительных точках представлены на рис. 3.
Точка Ь Точка с - — Точка сі ------ Точка е
15 175 325 475 2100 2700 3300 3600
Время, с
Рис. 3. Динамика изменения температуры газовоздушной среды в камере (точки а-е см. рис. 1)
Экспериментальное определение характеристик источника загрязнений
Массовый выброс СО от источника загрязнений QCO () определялся экспериментально.
Для проведения измерений использовался 4-компонентный газоанализатор СО-СН-СО2-О2 "АВТОТЕСТ-01.03-М" с вычислением коэффициента избытка воздуха Л. Краткие технические характеристики анализатора газа АВТОТЕСТ 01.03 М: диапазон измерения содержания оксида углерода СО 0-7 об. %, абсолютная погрешность ± 0,2%, относительная погрешность ± 6%.
Экспериментальное определение коэффициентов турбулентного обмена
Как следует из выражения (13), определяющее влияние на коэффициенты турбулентного обмена кроме геометрических размеров помещения оказывает кинетическая энергия газовоздушных потоков в помещении. Точное определение данных коэффициентов невозможно, так как в формулы (17) - (19) входят трудноопределимые экспериментально величины: кратность вентиляции помещения; коэффициент аэродинамического сопротивления и ряд других. Поэтому в задачи экспериментальных исследований входила возможность оценки вышеприведённых величин.
Для корректного определения энергии приточных струй (17) в газовой камере были проведены экспериментальные работы по определению кратности вентиляции помещения гаража для двух исследуемых случаев: 1) открыта одна створка ворот; 2) ворота приоткрыты (щель размером 0,47 м). С этой целью вначале создавалась максимальная загазованность в помещении, а затем контролировалось время "спада" концентрации СО до фоновых значений. Параллельно измерялась скорость вытекания газовоздушной смеси из камеры. Зная площади поперечных сечений отверстий вытекания газов 5 (2,56 м2 и 1 м2), определяли потребный воздухообмен и кратность вентиляции (табл. 2).
/
/
/
""ЛГ
Таблица 2
Потребный воздухообмен и кратность вентиляции
Исследуемый случай Время проветривания, Тпр , ч Скорость вытекания газов из камеры, , м/с (м/ч) Потребный воздухообмен камеры, Цт =\р • S •и, м3/ч Кратность воздухообмена, к - ^°тр Vп ’
Одна створка открыта 0,37 1 (3600) 3410 75
Щель размером 0,47 м 0,52 1,5 (5400) 2700 58
Затем по формулам (14) - (19) были рассчитаны энергии приточных и тепловых струй с учётом определяющих размеров помещения для каждого из случаев (табл. 3).
Таблица 3
Энергии приточных и тепловых струй
Исследуемый случай Энергия приточных струй, Є пс , м2/с3 Энергия тепловых струй, ЄГС, м2/с3 Суммарная энергия в камере Є , м2/с3 Определяющий размер помещения, 1 -\[Ё, м
Одна створка открыта 0,75 0,000227 0,750227 1,6
Щель размером 0,47 м 0,58 0,000227 0,580227 1,0
Наконец, в соответствии с формулами (4), (13) и (15) определили значения коэффициентов турбулентной диффузии и турбулентной теплопроводности в помещении для исследуемого угарного газа для трёх изучаемых случаев. Результаты расчётов сведены в табл. 4.
Таблица 4
Коэффициенты турбулентной диффузии и теплопроводности
Исследуемый случай (вариант расчёта) Коэффициенты турбулентной диффузии помещения камеры, м2/с Коэффициенты турбулентной теплопроводности помещения камеры, Вт/м^К Коэффициент турбулентной диффузии газа, м2/с Коэффициент турбулентной теплопроводности газа, Вт/м^К
А - А лпх - лпу АП2 V - КУ Оса \со
Закрытые ворота 0 0 0 0 15,240-ь 26,2 •Ю"3
Открыта одна створка ворот 0,42 Апг-1,5 - 0,3 504 К-1,5 - 360 15,240-ь 26,2 •Ю-3
Щель размером 0,47 м 0,2085 АШ-1,5 - 0,3 240 Ка-1,5 - 360 15,2^ 10-6 26,2 •Ю-3
Коэффициенты турбулентной диффузии угарного газа Dco рассчитаны по формуле [6]:
^ 4,3• 10“7 • Т3/2 П 1
~ Р-(уА + ^В3)2 умА + М~В' (23)
где МА, МВ - молярные массы СО и воздуха соответственно, г/моль; VА, VВ - мольные объёмы
СО и воздуха, см3; p - давление газовоздушной смеси, Па; T - абсолютная температура газовоздушной смеси, К.
Результаты модельных расчётов и сравнение их с экспериментальными данными.
В связи с тем, что на практике в различных помещениях всегда работает общеобменная или активная вентиляция, то случай 1 (закрытые ворота - вентиляционные отверстия заклеены) являет-
ся наиболее подходящим для теоритических исследований процессов диффузии газа, без учета вентиляционных потоков. Данный случай помогает выяснить влияние газов на диффузионную составляющую уравнений системы, описывающей процесс тепломассопереноса, поэтому было принято решение подробнее рассмотреть данный вариант исследования.
Результаты модельных расчётов концентраций и температур газовоздушной среды внутри гаражного помещения для первого случая - закрытые ворота (вентиляционные отверстия заклеены) - представлены на рис. 4-7.
Точка е
Точка й
Рис. 4. Динамика изменения концентраций СО в газовой камере (точки а-е см. рис. 1)
т Точка <1
Рис. 5. Динамика изменения температуры газовоздушной среды в камере (точки а-е см. рис. 1)
Как следует из графиков зависимостей (см. рис. 4, 5), процесс роста температуры в замкнутом объёме в значительной степени инерционный по сравнению с быстрым ростом концентрации угарного газа.
Также следует отметить, что концентрации ВВ быстро выравниваются по величине в разных частях исследуемого помещения (см. рис. 4), в отличие от температур газовоздушной смеси,
870
где имеет место значительное расхождение между параметрами источника (точка «с» - максимальная температура) и другими расчётными точками (см. рис. 5). Очевидно, что последнее обстоятельство связано со сложными процессами теплообмена газовой камеры с окружающей средой.
Картина сформированных полей концентраций ВВ и температур газовоздушной среды исследуемого помещения в вертикальном «срезе» на момент времени 2700 с представлена на рис. 6, 7. Как видно из результатов модельных расчётов, по высоте помещения имеет место чёткое расслоение концентраций СО, причём их максимум приходится в центре помещения у потолка. Тёплый угарный газ поднимается вверх, скапливаясь у потолка. В результате концентрации отличаются у пола и потолка в 2 раза. Несколько иное расслоение имеет место по температурам в помещении. Картина имеет симметричную форму, влияние на которую оказывает теплообмен с окружающей камеру средой. Необходимо отметить, что исследования проводились в «холодный» период года и разница температур ОГ и атмосферного воздуха составляла около 40-50 °С.
Рис. 6. Поле концентраций СО в газовой камере на высоте 1,5 м от пола (точки а-е см. рис. 1)
Рис. 7. Поле температур газовоздушной среды в камере на высоте 1,5 м от пола
(точки а-е см. рис. 1)
Наконец, на рис. 8, 9 даны графики изменения концентраций и температур внутри закрытого помещения камеры во времени, полученные как экспериментально, так и в результате расчёта по разработанной математической модели.
7000
£ бооо
Сн
. 5000
рГ
§■ 4000 Л
В"3000
X
о
а 2000
X ,о
юоо о
15 175 325 475 625
Время, с
Рис. 8. Сравнение теоретически рассчитанных концентраций СО в газовоздушной среде камеры по модели с экспериментальными данными (точки а, с см. рис. 1). Случай - ворота закрыты
ЗЬ,0
0,0
15 175 325 475 625 126021002700 33003600 Время, с
Рис. 9. Сравнение теоретически рассчитанных температур газовоздушной среды камеры по модели с экспериментальными данными (точки а, с см. рис. 1). Случай - ворота закрыты
Как видно из рис. 8, отклонение экспериментальных данных от расчётных концентраций СО в воздушной среде исследуемой камеры составляет не более 2-3%, что соответствует пределу допускаемой основной погрешности измерительной аппаратуры - 5% в интервале времени 0-175 с. Дальнейшее расхождение расчётных и экспериментальных зависимостей в интервале 175-625 с связано с тем, что газоанализатор Drager X-am 5000 имеет ограниченный предел измерений по СО до 2000 ppm и не позволяет определить дальнейшее возрастание концентраций.
На рис. 9 представлены графики динамики изменения температуры газовоздушной среды исследуемой камеры в течение часа, полученные как в результате численного моделирования, так и при проведении натурного эксперимента. Сравнение произведено для двух характерных точек помещения.
Как следует из рис. 9, отклонение экспериментальных данных от расчётных температур газовоздушной среды камеры составляет не более 2-3 °С, что является вполне допустимым для моделирования таких достаточно сложных явлений тепломассопереноса веществ в замкнутом пространстве.
Выводы.
1. Разработана математическая модель процессов распространения газообразного вещества и избытков теплоты в воздушной среде гаражного помещения.
2. Доказана адекватность представленной математической модели путём сравнения результатов численных расчётов искомых параметров внутри исследуемого помещения с параметрами, определёнными экспериментально на основе построенной физической модели.
3. Показана возможность применения разработанной модели при исследовании загазованности производственных помещений и для определения параметров микроклимата в рабочих зонах.
Библиографический список
1. Математическая модель процессов распространения вредных веществ и избытков явной теплоты в производственных помещениях / Б.Ч. Месхи [и др.] // Вестн. Донск. гос. техн. ун-та.
- 2010. - Т. 10. - С. 1031-1038.
2. Берлянд М.Я. Прогноз и регулирование загрязнения атмосферы / М.Я. Берлянд. - Л.: Гидрометеоиздат, 1985. - 272 с.
3. Марчук Г.И. Математическое моделирование в проблеме окружающей среды / Г.И. Марчук. - М.: Наука, 1982. - 320 с.
4. Патанкар С.В. Численные методы решения задач теплообмена и динамики жидкости / С.В. Патанкар. - М.: Энергоатомиздат, 1984. - 149 с.
5. Моделирование процессов переноса и ассимиляции вредных веществ в загазованном помещении участка обкатки / Л.Н. Алексеенко [и др.] // Вестн. Донск. гос. техн. ун-та. - 2009.
- Спецвыпуск. Технические науки. Ч. I. - С. 56-69.
6. Конечно-элементное моделирование активной вентиляции загрязненных помещений / Л.Н. Алексеенко [и др.] // V Всероссийская школа-семинар. Математическое моделирование и биомеханика в современном университете, 1-5 июня 2009 г. - С. 9.
Материал поступил в редакцию 10.05.11.
References
1. Matematicheskaya model' processov rasprostraneniya vredny'x veshhestv i izby'tkov yavnoj teploty' v proizvodstvenny'x pomeshheniyax / B.Ch. Mesxi [i dr.] // Vestn. Donsk. gos. texn. un-ta.
- 2010. - T. 10. - S. 1031-1038. - In Russian.
2. Berlyand M.Ya. Prognoz i regulirovanie zagryazneniya atmosfery' / M.Ya. Berlyand. - L.: Gi-drometeoizdat, 1985. - 272 s. - In Russian.
3. Marchuk G.I. Matematicheskoe modelirovanie v probleme okruzhayushhej sredy' / G.I. Marchuk. - M.: Nauka, 1982. - 320 s. - In Russian.
4. Patankar S.V. Chislenny'e metody' resheniya zadach teploobmena i dinamiki zhidkosti / S.V. Patankar. - M.: E'nergoatomizdat, 1984. - 149 s. - In Russian.
5. Modelirovanie processov perenosa i assimilyacii vredny'x veshhestv v zagazovannom po-meshhenii uchastka obkatki / L.N. Alekseenko [i dr.] // Vestn. Donsk. gos. texn. un-ta. - 2009.
- Speczvy'pusk. Texnicheskie nauki. Ch. I. - S. 56-69. - In Russian.
6. Konechno-e'lementnoe modelirovanie aktivnoj ventilyacii zagryaznyonny'x pomeshhenij / L.N. Alekseenko [i dr.] // V Vserossijskaya shkola-seminar. Matematicheskoe modelirovanie i biomex-anika v sovremennom universitete, 1-5 iyunya 2009 g. - S. 9. - In Russian.
MATHEMATICAL AND EXPERIMENTAL SIMULATION OF CARBON OXIDES AND EXCESS HEAT PROPAGATION IN GAS-AIR INDOOR ENVIRONMENT
B.C. MESKHI, E.I. MASLOV, A.N. SOLOVYEV, Y.I. BULYGIN, D.A. KORONCHIK
(Don State Technical University)
Results of the mathematical and experimental simulation of the distribution of carbon monoxide and excess heat in the air of the garage premises are presented. Fields of air movement, carbon oxide concentration and temperature indoor comparable to the experimental data have been obtained on the ground of the theoretical and experimental method.
Keywords: air pollution, excess heat, concentration, temperature, ventilation.
