Математическое моделирование процесса конденсации серного ангидрида в атмосфере, прилегающей к тепловой электрической станции

3. Предложенный алгоритм при нулевой длине трассы распространения и отсутствии турбулентных или любых других возмущений обеспечивает удовлетворительную точность регистрации суммарного топологического заряда группы особых точек, только если их количество не слишком велико (<60).

4. При использовании алгоритма для регистрации параметров пучка в турбулентной

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

СПИСОК ЛИТЕРАТУРЫ

1. Soskin M.S., Vasnetsov M.V. Singular optics // Prog. Opt. -2001. - V. 42. - P.219-276.

2. Siegman A.E. Lasers. - Mill Valley, Calif.: University Science, 1986. - 1304 p.

3. Young C.Y., Gilchrest Y.V., Macon B.R. Turbulence induced beam spreading of higher order mode optical waves // Optical Engineering. - 2002. - V. 41. - № 5. - P. 1097-1103.

4. Palacios D.M., Rozas D., Swartzlander G.A., Jr. Observed scattering into a dark optical vortex core // Phys. Rev. Lett. - 2002. -V. 88. - P. 103902.

Поступила 27.11.2012 г.

УДК 536.253

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССА КОНДЕНСАЦИИ СЕРНОГО АНГИДРИДА В АТМОСФЕРЕ, ПРИЛЕГАЮЩЕЙ К ТЕПЛОВОЙ ЭЛЕКТРИЧЕСКОЙ СТАНЦИИ

Д.В. Гвоздяков, В.Е. Губин

Томский политехнический университет E-mail: dim2003@tpu.ru; gubin@tpu.ru

Представлены результаты численной оценки интенсивности процесса конденсации серного ангидрида в воздушном пространстве, прилегающем к тепловой электрической станции. Определены размеры ядер капель кислоты, которые в процессе седиментации могут выпадать на поверхность Земли.

Ключевые слова:

Математическое моделирование, тепловая электрическая станция, дымовые газы, загрязнение, серный ангидрид, конденсация, частица, серная кислота.

Key words:

Mathematical modeling, thermal power plant, smoke gases, pollution, sulphurous anhydride, condensation, particle, sulfuric acid.

Введение

Проблема антропогенного загрязнения атмосферы Земли многие годы является актуальной [1, 2]. Одним из основных источников загрязняющих веществ являются объекты теплоэнергетического комплекса, такие как тепловые электрические станции (ТЭС) [1, 2]. С каждым годом объемы эмиссии загрязняющих веществ в воздушный бассейн планеты увеличиваются [3-5], хотя, следует отметить, прогнозы мировой антропогенной эмиссии оксидов серы и азота неоднозначны (табл. 1). Наличие в воздушном пространстве Земли таких загрязняющих веществ как 802, 8О3, N0^ при их взаимодействии с компонентами атмосферного воздуха, может приводить к образованию соответствующих кислот [2]. Взаимодействие кислотных дождей с биосферой, сооружениями, техникой, как правило, заканчиваются для них губительными последствиями [1, 5].

Таблица 1. Фактические выбросы оксидов серы и оксидов азота в 2000 г. и прогнозы мировой антропогенной эмиссии [3]

Источник Оксиды серы, млн т Оксиды азота, млн т

2000 202Б 20Б0 2100 2000 202Б 20Б0 2100

1РСС (Межправительственная комиссия по изменению климата) 1992 [5], сценарий 1Б92а 178 240 308 296 122 174 210 273

WEC (Всемирный энергетический совет) 1993 [6], сценарий В - 132 - - - 89 - -

Клименко В.В., Клименко А.В. [3] 14Б 13Б 100 Б0 114 126 114 87

1B5

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

• 5Я<0,1-106 м, ядра Айткена;

• 0,1-10-6<5Я<1-106 м, большие;

• 5Я >140-6 м, «гигантские».

При достижении размеров 0,540-3...2,540-3 м возможно выпадение кислотных осадков на поверхность Земли в процессе седиментации [1, 2] как на территориях, прилегающих к ТЭС, так и на отдаленных территориях.

В работе рассматривается серный ангидрид 803, который образуется в газоходах котлов ТЭС путем частичного окисления 802 (до 5 % от общей доли в02) при сжигании высокосернистых топлив и относится к классу умеренно опасных продуктов сгорания натуральных топлив ТЭС. Несмотря на высокую опасность антропогенной эмиссии кислотообразующих веществ в воздушный бассейн и достаточно большое внимание, уделяемое проблеме кислотных дождей [1, 2], до настоящего времени не изучены процессы образования капель кислот в атмосфере Земли.

Целью настоящего исследования является численное моделирование процесса конденсации серного ангидрида в атмосфере Земли на поверхности ядер конденсации размером от 0,110-6 до 5,010-6 м.

Постановка задачи и метод решения

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

Механизм (1) образования серной кислоты [9-10, 2] основан на взаимодействии, при котором серный ангидрид 803 поглощается водяным паром.

Я03+Н20^Н2Я04.

На рис. 1 представлена область решения задачи в начальный момент.

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

Рис. 1. Область решения задачи в начальный момент: 1) воздух плюс пары серного ангидрида; 2) ядро конденсации

Процесс переноса энергии, количества движения и массы в рассматриваемых условиях описывается системой нестационарных дифференциальных уравнений в частных производных [11, 12]. Считается, что температура и концентрация на источнике (дымовая труба ТЭС) не зависят от времени. Для моделирования изменения основных искомых функций использованы нестационарные двумерные уравнения смешанной конвекции в приближении Буссинеска [13].

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

60 тт 60 60

—+ и -------+ V --=

6т дХ 6У

1 (620 620) Ог

- + —н+-

Яе КдХ2 6У2

Яе2

60

6Х

62Т 62Т _

2 +--------2 = О,

6Х2 6У

601 + и ,601 + V .601 =

6т 6Х 6У

1 (6201 + 60Л + е^Л

Рг- Яе ' І 6Х2

6У ) рСрМ'

6Сз + и .6С1+у .6Сз

6т

6Х

6У

1

( 62С3 'І6Х2

6У2

р-и0

1 60 620 620,

Бо 6т 6Х2 6У2 ’

где X, Y - координаты декартовой системы координат; т - время; 0 - температура; С3 - концентрация серного ангидрида; и, V - скорости в направлении оси X и Y, соответственно; Т - функция тока; О - вихрь скорости; 2 - тепловой эффект реакции, Дж/кг; WK - массовая скорость конденсации

серного ангидрида, кг/(м3-с); Ьх -геометрический размер частицы (ядра конденсации), м; р - плотность, кг/м3; СР -удельная теплоемкость, Дж/(кг-К); и0 - скорость внешней среды, м/с; Re -число Рейнольдса; Рг - число Прандтля; 8е - число Шмидта; Fo - число Фурье; 1 - воздух плюс пары серного ангидрида; 2 - ядро конденсации; 3 - серный ангидрид.

Для расчета скорости конденсации использовалась формула [13, 14]:

=р

р - р0ец> [- —

\2kRT

I М

где РА - атмосферное давление, Па; Р0 - предэкспо-нента, Па; Е - энергия активации процесса конденсации, Дж/кг; Я - универсальная газовая постоянная, Дж/(моль-К); М - молярная масса, кг/моль; в - коэффициент конденсации.

Выражение (7) устанавливает, что скорость конденсации пропорциональна разности двух односторонних максвеловских скоростей конденсации [13]. Таким образом, вводятся два допущения:

1) конденсирующаяся среда у поверхности ядра конденсации макроскопически неподвижна; 2) состояние конденсирующейся среды может быть описано локальной максвеловской функцией [14].

Начальные условия для системы уравнений (2)-(6) в момент времени т=0 имеют вид: 0(ХД,0)=0, С (ХД,0)=0,

¥(ХД,0)=0, и (ХД,0)=0,

0(ХД,0)=0, V (ХД 0)=0.

Граничные условия для (2)-(6) уравнений:

X = 0,0 < У < Н1: — = 0, С3 = 0, Т = 0,0 = 0;

дХ

X = 0, Н1 < У < Н: © = ©1,

ду=и4, -дХ= -^ ^ = С0;

X = А,0<У<н : = ^-О:-^,

1 1 дХ Я1 дХ АТ •Я1

дС — • I 01 =02, —1 , у/ = 0 = 0;

1 2 дХ Б :р г

дС д^ дО

Х = Ь,0 < У < Н: 0 = 01, —3 = 0,-----------= 0,-= 0;

1 дХ дХ дХ

У = 0,0<Х<1Х: — = 0, С3 = 0,Т = 0,0 = 0;

дУ

У = 0, Ц < Х < Ь: 0 = 01, С3 = С0,

дЖ. = и ^ = -удУ 4, дХ 4;

У = Н1,0 < Х < А: ^ =^ ^ - О-кА,

1 1 дХ А1 дХ АТ •Я1

дС — • I 01 =02, —1 , у/ = 0 = 0;

1 2 дХ Б :р г

У = Н,0<Х<Ь: 0 = 01,= 0,— = 0, — = 0.

1 дУ дУ дУ

где и4, V4 - значения компонент скорости на границе, соответствующей источнику ввода серного ан-

гидрида; С0 - безразмерная концентрация серного ангидрида на границе области.

Разностные аналоги дифференциальных уравнений решены локально-одномерным методом [15-16]. При решении задачи использовался алгоритм, разработанный [17, 18] для решения задач сопряженного теплопереноса в областях с локальными источниками энергии и массы.

Начальная концентрация серного ангидрида в воздушном пространстве принималась соответствующей значениям для летнего времени года [19]. Начальная концентрация серного ангидрида составляет 3,0 мг/м3. При постановке задачи были приняты следующие допущения:

1) не учитывалось испарение капли в процессе движения;

2) не принималось во внимание влияние солнечной радиации;

3) ядро конденсации имеет форму куба.

Результаты и обсуждения

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

Таблица 2. Результаты численных исследований

t, с 50-106, м

0 0,03 0,1 0,5 0,8 1,0 5,0

5Ч-103, м

600 0,012 0,045 0,110 0,178 0,230 0,920

1800 0,025 0,089 0,270 0,445 0,440 1,380

3600 0,036 0,118 0,400 0,587 0,580 1,790

5400 0,045 0,138 0,471 0,652 0,680 2,080

7200 0,052 0,155 0,520 0,719 0,750 2,310

9000 0,059 0,169 0,562 0,776 0,820 2,510

10800 0,064 0,181 0,588 0,826 0,880 2,670

Примечание: 80 - размер ядра конденсации в начальный момент времени, м; (5Ч - размер сформировавшейся частицы в момент времени ^ с.

Установлено, что при начальном размере ядра конденсации 1,010-6 м уже через 1800 с образуется частица размером более 1,010-3 м. Известно, что частицы соответствующего размера относятся к типу «грубодисперсных» [9] (являются дождевыми каплями). При начальном размере ядра

5.010-6 м формируется частица размером более

1.010-3 м за 600 с с момента начала процесса конденсации. Дальнейшее увеличение начального размера ядра может приводить к образованию частиц с размером более 2,510-3 м. Такие частицы размером 1,010-3 м и более способны выводится из

0 2000 4000 6000 8000 10000 12000

5

6

Рис. 2. Изменение размеров частиц (м): 1) 8=0,03-10-6; 2) 8=0,1-10-6; 3) 8=0,5-10-£; 4) 8=0,8-10-6; 5) 8=1,0-10

6) 8=5,0-10-6

1

2

3

4

атмосферы за счет процесса седиментации [1, 2, 9]. В результате тепломассопереноса в области дымового факела тепловой электрической станции образовавшиеся частицы размером менее 1,010-3 м с начальным размером ядра конденсации 0,0310Л..0,810-7 м могут осаждаться на поверхности более крупных частиц [1, 2]. Таким образом, ускорение процесса формирования атмосферных осадков может очень быстро привести к выпадению дождевых капель на территориях, прилегающих к тепловой электрической станции.

На рис. 2 графически показан рост частицы в зависимости от времени при различных начальных размерах ядер конденсации.

Известно, что на фиксированном расстоянии от тепловой электрической станции при определенном классе устойчивости атмосферы можно выделить территории с максимальными значениями концентраций веществ антропогенного происхождения [20]. Таким образом, целесообразным является дальнейший анализ процесса образования ча-

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

Заключение

Проведен численный анализ процесса конденсации серного ангидрида на поверхности ядра конденсации в воздушном пространстве, прилегающем к ТЭС. Установлено, что частицы, способные выводиться из атмосферы Земли за счет процесса седиментации, образуются при начальном размере ядра конденсации 1,010-6 м и более за время 0,5.3 ч. Частицы с меньшим размером способны переноситься на большие расстояния от источника антропогенных компонент за счет конвекции и диффузии.

СПИСОК ЛИТЕРАТУРЫ

1. Заиков Г.Е., Маслов С.А., Рубайло В.Л. Кислотные дожди и окружающая среда. - М.: Химия, 1991. - 144 с.

2. Израэль Ю.А., Назаров И.М., Прессман А.Я., Ровин-ский Ф.Я., Рябошапко А.Г. Кислотные дожди. - Л.: Гидроме-теоиздат, 1989. - 270 с.

3. Клименко В.В., Клименко А.В. Сокращение выбросов малых парниковых газов как альтернатива снижению эмиссии углекислого газа. Ч. I // Теплоэнергетика. - 2000. - № 6. - С. 6-12.

4. Внуков А.К. Защита атмосферы от выбросов энергообъектов: справочник. - М.: Энергоатомиздат, 1992. - 176 с.

5. Climate Change 1992. The Supplementery Report to the IPCC Scientific Assessment / J.T. Houghton, B.A. Callender,

S.C. Varney. Cambridge University Press, 1992.

6. Brendow K. West and East European approaches to sustainable energy // World Energy Concil Journal. Desember. 1993. - P. 41-53.

7. Фелленберг Г. Загрязнение природной среды. Введение в экологическую химию. - М.: Мир, 1997. - 232 с.

8. Грин Х., Лейн В. Аэрозоли - пыли, дымы и туманы. - Л.: Химия, 1969. - 428 с.

9. Терехов В.И., Пахомов М.А. Тепломассоперенос и гидродинамика в газокапельных потоках: монография. - Новосибирск: Изд-во НГТУ, 2008. - 284 с.

10. Андруз Дж., Бримлекумб П., Дишклез Т. и др. Введение в химию окружающей среды. - М.: Мир, 1999. - 271 с.

11. Роуч П. Вычислительная гидродинамика. - М.: Мир, 1980. - 616 с.

12. Берковский Б.М., Полевиков В.К. Вычислительный эксперимент в конвекции. - Минск: Университетское, 1988. - 167 с.

13. Кнаке О., Странский И.Н. Механизм испарения // УФН. -1959. - Т. 68. - № 2. - С. 261-305.

14. Лабунцов А.Д. Физические основы энергетики. Избранные труды по теплообмену, гидродинамике, термодинамике. - М.: Изд-во МЭИ, 2000. - 388 с.

15. Пасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло- и массообмена. - М.: Наука, 1984.- 288 с.

16. Самарский А.А., Попов Ю.П. Разностные схемы газовой динамики. - М.: Наука, 1975. - 424 с.

17. Кузнецов Г.В., Шеремет М.А. Об одном подходе к математическому моделированию тепловых режимов радиоэлектронной аппаратуры и электронной техники // Микроэлектроника. -2008. - Т. 37. - №2. - С. 150-158.

18. Кузнецов Г.В., Шеремет М.А. Сопряженный теплоперенос в замкнутой области с локально сосредоточенным источником тепловыделения // Инженерно-физический журнал. - 2006. -Т. 79. - № 1. - С. 56-63.

19. Гвоздяков Д.В., Губин В.Е., Кузнецов Г.В. Численного оценка формирований условий атмосферного образования серной кислоты в районе расположения тепловой электрической станции // Научно-технические ведомости СПбГПУ. - 2012. - № 2-2. -С. 195-200.

20. Федосов А.А. Математическая модель загрязнения атмосферы выбросами тепловых электрических станций // Альтернативная энергетика и экология. - 2006. - № 9 (41). - С. 98-102.

Поступила 27.05.2013 г.

УДК 532.529.2

СЕДИМЕНТАЦИЯ ЧАСТИЦ В ФЕРРОКОЛЛОИДАХ С УЧЕТОМ ДИНАМИЧЕСКИХ ГРАНИЧНЫХ УСЛОВИЙ

А.Ф. Глухов, В.А. Демин

Пермский государственный национальный исследовательский университет E-mail: demin@psu.ru

Проведено численное исследование процесса осаждения частиц в феррожидкости на основе керосина в тонком вертикальном канале с учетом возможного термодиффузионного перераспределения частиц при неоднородном нагреве. Определены значения термодиффузионного параметра, при которых действие термодиффузии для установленной разности температур на торцах канала становится сопоставимым с седиментационным механизмом расслоения.

Ключевые слова:

Феррожидкость, седиментация частиц, термодиффузия, динамические граничные условия. Key words:

Ferrofluid, sedimentation of particles, thermodiffusion, dynamical boundary conditions.

Большинство задач гидродинамики и тепломассообмена математически сводятся к краевым задачам. Очевидно, что граничные условия определяют вид решений наравне с уравнениями. Однако степень обоснованности у граничных условий на разные неизвестные величины в гидродинамике неодинаковая. В задачах теплообмена краевые условия на температуру выставляются, как правило, физически строго: на границах расчетных областей температура должна быть непрерывной,

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