экология воздушной среды и космического пространства
ecology of air atmosphere and space
Статья поступила в редакцию 10.06.13. Ред. рег. № 1677
The article has entered in publishing office 10.06.13. Ed. reg. No. 1677
УДК 621.1:502:519.6
ГИДРОДИНАМИКА ВРЕДНЫХ ВЫБРОСОВ В АТМОСФЕРУ Р.А. Амерханов, К.А. Гарькавый, А.С. Кириченко
Кубанский государственный аграрный университет 350044 Краснодар, ул. Калинина, д. 13 Тел.: (861) 221-58-54, e-mail: energyksau@mail.ru
Заключение совета рецензентов: 15.06.13 Заключение совета экспертов: 18.06.13 Принято к публикации: 20.06.13
Приведены математические модели конвективной диффузии вредных выбросов в атмосферу. Рассматриваются стационарное и ламинарное движение среды; турбулентная диффузия, приводятся нелинейные уравнения диффузии. Изложены методы решения сформулированных математических моделей.
Ключевые слова: источники вредных выбросов, ламинарное стационарное движение, турбулентная диффузия, нелинейные уравнения, методы решения.
HYDRODYNAMICS OF HARMFUL WASTES INTO ATMOSPHERE R.A. Amerkhanov, K.A. Garkavij, A.S. Kirichenko
Kuban State Agrarian University 13 Kalinina str., Krasnodar, 350044, Russia Tel.: (861) 221-58-54, e-mail: energyksau@mail.ru
Referred: 15.06.13 Expertise: 18.06.13 Accepted: 20.06.13
There were considered the theoretical methods of convective diffusion of harmful wastes into atmosphere. There were presented the stationary and laminar medium movement, non-linear equations of diffusion are cited. There were given an account of solving methods of formulated non-linear differential equations of diffusion.
Keywords: sources of harmful wastes, laminar stationary movement, convective diffusion of a substance, non-linear equations, methods of solving.
Роберт Александрович Амерханов
^ f
Константин Алексеевич Гарькавый
Сведения об авторе: профессор кафедры электротехники, теплотехники и возобновляемых источников энергии Кубанского гос. аграрного университета, д-р техн. наук. Заслуженный работник высшей школы РФ, почетный работник высшего профессионального образования РФ, почетный работник по науке и технике РФ, награжден орденом Почета.
Основной круг научных интересов: энергосбережение естественных ресурсов при использовании нетрадиционных и возобновляемых источников энергии в агропромышленном комплексе.
Публикации: около 340.
Сведения об авторе: профессор кафедры электротехники, теплотехники и возобновляемых источников энергии Кубанского гос. аграрного университета, канд. техн. наук. Премия - Олимп науки Кубани.
Основной круг научных интересов: энергосбережение естественных ресурсов при использовании нетрадиционных и возобновляемых источников энергии в агропромышленном комплексе.
Публикации: 63.
i -4
Анна Сергеевна Кириченко
Сведения об авторе: аспирант Кубанского гос. аграрного университета.
Основной круг научных интересов: энергосбережение естественных ресурсов при использовании нетрадиционных и возобновляемых источников энергии в агропромышленном комплексе.
Публикации: 5.
Международный научный журнал «Альтернативная энергетика и экология» № 06/1 (127) 2013 © Научно-технический центр «TATA», 2013
Широкий класс технологических процессов, при которых имеют место химические реакции и физические превращения, сопровождается большим количеством вредных для окружающего мира выбросов. Значительная часть вредных ингредиентов возникает и попадает в атмосферу в результате горения. Кроме того, следует иметь в виду, что вредные вещества создаются в значительных количествах не только антропогенным путем, но также в результате природных явлений: лесных пожаров, нарушений в работе электростанций и т.д.
Надо признать, что в настоящее время происходит непрерывное нарастание загрязняющего эффекта. Это объясняется тем, что промышленность страны, обеспечивающая основной вклад в ВВП, по-прежнему ориентирована на ресурсо- и энергозатратное производство. Использование возобновляемых источников энергии и вторичных энергетических ресурсов недопустимо мало. Другой причиной служит то, что значительная доля выбросов происходит от коммунального хозяйства.
На сегодняшний день одной из основных задач является оценка вредных выбросов в атмосферу на основе анализа теоретических методов. Это позволит определить конкретные пути решения экологической проблемы [1, 2].
Будем считать, что нарушение равновесия в среде связано с изменением концентрации вредных веществ от точки к точке и с наличием в ней макроскопического конвективного движения. Предположим, что температура смеси является постоянной, а градиент давления достаточно малым.
Таким образом, перенесение вещества в подвижной среде будет обусловлено лишь двумя разными механизмами. Во-первых, взвешенные газоаэрозольные частицы увлекаются средой и переносятся вместе с ней, во-вторых, при наличии разницы концентраций в среде возникает молекулярная и турбулентная диффузия. Совокупность обоих процессов называют конвективной диффузией вещества в воздушной или жидкой среде. Запишем дифференциальное уравнение, которое описывает диффузное перенесение в подвижной несжимаемой среде. Предположим, что процесс происходит в стационарном и ламинарном режиме [3].
В данном пространстве можно провести некоторые мнимые поверхности, в каждой точке которых в момент времени т концентрация постоянна. Такие поверхности называют «изоповерхность». Основная гипотеза математической теории диффузии заключается в том, что диффузионный поток вещества через единицу площади изоповерхности в единицу времени изнутри наружу равняется
J = - D = - D grad C = - D dC dn dx.
(1)
где Б - коэффициент молекулярной диффузии, которая зависит в общем случае от концентрации С и тем-
пературы Т; п - внешняя нормаль к поверхности. Знак минус указывает на то, что поток вещества направлен в сторону уменьшения концентрации вещества.
Если взвешенное вещество находится в подвижной со скоростью V среде, то последняя захватывает его в своем движении. При этом рядом с диффузионным движением возникает конвективный поток (масса вещества, которая переносится средой через единицу площади в единицу времени): 3 к = CV. Полный поток вещества представляет собой сумму конвективного и диффузионного потоков:
J = Jk + Jd = VC - D grad C .
(2)
Выделим в рассмотренной сплошной среде произвольный объем Ж, ограниченный поверхностью & Пусть в рассмотренном объеме находится источник вещества, мощность которого равняется Q (мг/с).
Из закона сохранения следует, что изменение массы вещества в объеме Ж в единицу времени равняется сумме массы вещества, которое проходит через поверхность, ограничивающую данный объем, и притока вещества от источника. Принимая во внимание это положение и используя поверхностное превращение интеграла по формуле Остроградского-Гаусса [4], будем иметь:
дС дт
+ div VC = div( D grad С) + G.
div( VC) = V grad С + С div V .
или
- с -
VAiC Ж
(3)
В прямоугольной декартовой системе координат уравнение конвективной диффузии в подвижной ламинарной среде запишется следующим образом:
дС диС ЗуС дwC -+-+-+-=
дт дх ду дх
д ( дСЛ д ( дС Л д ( дС Л „
= — \ Б-1^ — 1 Бу-\ + —\ Я-1 + G, (4)
дх { х дх ) ду { у ду ) дх У 2 дх )
где и, V, V! - соответствующие компоненты скорости ветра подвижной сплошной среды, которые в общем случае являются функциями координат и должны быть определены из гидродинамической задачи; Бх, Бу, Б2 - коэффициенты молекулярной диффузии в соответствующих направлениях.
Представим
(5)
В силу несжимаемости среды, V = 0, уравнение концентрации принимает вид
— + V ягаа С = ШУ( б ^ С) + G (6)
дт
dC dC dC dC ——+ ——+ v——+ w—— = дт dx dy dz
= ±Г dCY^r Щ + ±(вя д£) + G. (7)
dx v dx ) dy V y dy ) dz v dz
Э 0 I
3 CL
International Scientific Journal for Alternative Energy and Ecology № 06/1 (127) 2013
© Scientific Technical Centre «TATA», 2013
Р.А. Амерханов, К.А. Гарькавый, А.С.Кириченко. Гидродинамика вредных выбросов в атмосферу
Обобщив уравнение для случая турбулентной подвижной среды, представим концентрацию вещества и скорость ветра в виде суммы среднего значения и малых пульсационных составляющих:
С = С + С'; и = и + и'; V = V + у'; V = V + . (8)
Причем
1 х+Дх i х+Дх
С = — I Cdт; u = — I Дт J Дт J
1 т + Дт
"Дт
с + дт
J vdt;
т+Дт
I udт; Дт т
Дт
1 т+Дт
w = -— J wdт
(9)
т+Дт т+Дт
С = — f Cdт = 0; м7 = — f uf dт = 0;
At J At J
1 т+Дт i т+Дт
— f V dt = 0; w = — f w'dт = 0. (10)
Дт Дт
В большинстве случаев молекулярным переносом возможно пренебречь, и в этом случае будем иметь:
дС -дС -дС -дС
--+ и--+ V--+ V-=
дт дх ду дх
дх
De
дС дх
Л
(
D„
дС
Л
дУ { У дУ
(
дх
De
дС дх
Л
+ G. (11)
С(х, y, z, т) I т=0 = f (х, y, z):
(12)
Предельное условие второго рода заключается в задании плотности потоку вещества на предельную поверхность как функции времени:
- D.
дС
дп
= (XS , yS , ZS , т
(14)
Самый простой случай однородного предельного условия второго рода состоит в постоянстве потока вещества на поверхность:
- D.
дС
дп
= ц = const.
(15)
При исследовании распространения газоаэрозольных примесей в атмосфере считают, что на поверхности земли выполняется условие отображения (поток вещества на поверхности равняется нулю):
- D„
дС дп
= 0.
(16)
Следовательно, конвективная диффузия описывается уравнением в частных производных второго порядка с переменными коэффициентами. Для его решения необходимо знать распределение вещества в данной области в начальный момент времени (начальное условие), геометрию рассматриваемой области (безграничной, полуграничной или ограниченной среды) и закон взаимодействия исследуемого взвешенного вещества с предельной поверхностью (предельное условие). Начальное условие записывается в следующем виде:
где fx, y, z) - известная функция координат.
В качестве предельного условия при бесконечности задают, что концентрация при удалении от источника стремится к фоновому значению С0 или к нулю: С = С0, при l ^ го, где l - расстояние от источника.
На предельной поверхности условия выглядят сложнее. При задании концентрации вещества на предельной поверхности S в любой мо
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.