Прогнозирование длины активной зоны конверсионных плазмохимических реакторов

УДК 621.314.211

ПРОГНОЗИРОВАНИЕ ДЛИНЫ АКТИВНОЙ ЗОНЫ КОНВЕРСИОННЫХ ПЛАЗМОХИМИЧЕСКИХ РЕАКТОРОВ

СЕМЕНОВ В. К., д-р техн. наук

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

Ключевые слова: конверсионные плазмохимические реакторы, активная зона, стримерная корона, фронтальные волны ионизации.

CONVERSION PLASMA-CHEMICAL REACTOR CORE REGION LENGTH FORECASTING

V.K. SEMEONOV, Ph.D.

The work represents the method of determining conversion plasma-chemical reactor core region length on the base of streamer crown, ionization frontal waves and barrier discharges.

Key words: conversion plasma-chemical reactor, core region, streamer crown, ionization frontal waves.

Защита окружающей среды от токсичных газообразных примесей во всем мире представляет собой важную экологическую проблему. Наиболее распространенными являются технологии очистки, основанные на использовании сорбционных, каталитических и термических процессов. Общими недостатками названных методов являются их неуниверсальность, высокая энергоемкость, крупные габариты газоочистного оборудования и значительные потери в окружающую среду. Все это приводит к тому, что защита атмосферы от оксидов серы и азота, содержащихся в дымовых газах ТЭС, в настоящее время в основном осуществляется рассеиванием их дымовыми трубами высотой 180, 250 и даже 320 м. Такой подход ни в коей мере не является радикальным решением проблемы. В последние годы продолжаются интенсивные исследования по поиску и разработке новых более эффективных технологий. Большие возможности здесь открываются на основе использования методов химии высоких энергий, в частности плаз-мохимии [1, 2].

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

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

Для оценки длины активной зоны реактора примем предположение о том, что стримеры, а также химически активные частицы равномерно распределены по объему зоны. Кроме того, реальный многоканальный процесс химических превращений компоненты, подлежащей конверсии, аппроксимируем некоторым модельным уравнением чистой гибели

сИп\/ - йУп) = -к

где п - концентрация молекул; й - коэффициент диффузии; к - константа скорости реакции; т - формальный порядок химической реакции (редко бывает выше третьего).

Ради определенности далее анализ будем вести для плоскопараллельного движения

д(пи) 5(пу) ^,д2п д2пч , т

+ + = -кпт.

дх ду "дх2 ду2

Для реакторов с большим расходом газа (реакторы импульсной короны, фронтальных волн ионизации) диффузией частиц вдоль га-

зового потока можно пренебречь, по сравнению с направленным движением:

д(пи) + д(пу) _ 0 5_п =_кпт дх ду ду2

Граничными условиями задачи задаются концентрация компонента на входе в реактор п = п0 и отсутствие потока частиц через боковые стенки реактора. Так как нас интересует зависимость концентрации частиц от продольной координаты реактора, то найдем приближенное решение исходного уравнения методом интегральных соотношений. Для этого проинтегрируем уравнение (3) поперек разрядного канала. С учетом граничных условий задачи имеем

— Гпис1у = _ьГптс1у,

сх I {

0 0

Ь2кпт_1

где Ь =-0--безразмерный параметр.

V

Уравнение (4) записано в безразмерном виде. В целях упрощения обозначений безразмерные величины обозначены теми же символами, что и размерные. В качестве масштабов приняты: ширина (при симметричном профиле скоростей полуширина) канала И, исходная концентрация примеси п0 и гидродинамический масштаб скорости V/И .

В цилиндрических координатах приведенное уравнение имеет следующий вид:

С 1 1

— ГпигСг = _Ь ГптгСг.

сИ {

0 0

Задавая профиль концентрации в виде п

— = _ВГ, (5) с1х

где

В = •

Ь| фт (у)Су

_0_

1 .

|иф(у)с1у

Вид решения уравнения (5) зависит от формального порядка реакции т. При реакции первого порядка зависимость п(х) носит экспоненциальный характер:

п

При втором и третьем порядке реакции решения соответственно имеют вид

п

1 + В(п0х'

п(х) =

+ 2В2п2х

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

п = п

пр

для длины в размерном виде: 1

и > 1п—|у(у)ф(у)с1у ппр 0

I >•

к| Ф(у)с1у

Здесь функция у

С учетом этого формула (10) принимает

вид

п„

I >•

< и > 1п— п

пр

(12) к

Полагая для примера к = 10 1/с, п0/ппр = 100, <и> = 2 м/с, получим Ь > 0,92 м.

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

' кИ2 п0 ДЛ

< Ре >;-; —-; К; А

V пп

^ = Г И

'пр

Здесь <Ре> - среднее по сечению канала реактора число Рейнольдса; V - кинематическая вязкость газа; к - константа скорости реакции

„ Сг первого порядка; К =•

< Ре >

- конвективный

0

0

п

0

1 -X

Рг

параметр; А =-— - параметр, характери-

< Ре >

зующий эффективность отвода тепла через боковые стенки реактора; Сг, Рг - число Грасс-гофа и число Прандтля, соответственно.

В реакторах барьерного разряда, где скорость протока газа мала, при больших градиентах концентрации может оказаться, что продольной диффузией пренебречь уже нельзя. В этом случае уравнение для ^-функции будет иметь вид

с^

a— -pfm =0,

dx2 dx где

i

V < Re > j<p(y)y(y)dy о

а =

i

Dj9(y)dy

P = -

kh2nm-1 |9m(y)dy

о_

i '

Dj9(y)dy

Решение нелинейного уравнения (14) будем искать приближенным методом Галеркина. Суть метода заключается в том, что решение уравнения ищется в виде отрезка ряда, представляющего собой линейную суперпозицию некоторых базисных функций, удовлетворяющую граничным условиям. Коэффициенты разложения находятся из интегральных соотношений, выражающих ортогональность невязки к каждой базисной функции. При удачном выборе функций достаточно точный результат получается при отрезке ряда, содержащем несколько функций. Так как нас интересует оценка длины реактора, то можно ограничиться одной базисной функцией. Выбор ее очевиден. Уравнение (14) при т = 1 имеет точное решение: концентрация обезвреживаемого компонента падает вдоль реактора по экспоненциальному закону. При т ф 1закон падения концентрации будет близок к экспоненциальному, поэтому решение уравнения (14) будем искать в виде

f(x) = exp(-yx)'

После элементарных вычислений получим

2р

а а Y =--+ J— +

2 V 4 m +1

Таким образом, длина активной зоны реактора должна удовлетворять неравенству пп

hln-

L >•

n

пр

Y

Полагая для примера к = 10 1/с, п0/ппр = 100, Ре = 30, И = 2 мм, получим I > 15 см. В безразмерном виде длина активной зоны реактора должна определяться следующими критериями подобия:

( ьи2„т-1 „ Л

L = f h

< Re >; -; D

kh2nm-1

; K; a

Ппр

В качестве примера на рис. 1 приведены результаты численного и аналитического расчета для реакции второго порядка.

1

0,8 _ 0,6

0.4.

0,2

0

0 5 10 15 20 х

Сравнение результатов численного и аналитического решения

Сравнение результатов расчета по формуле (18) с результатами численного решения уравнения (14) для конкретных данных показало, что для реакций первого порядка относительная погрешность не превосходит 0,01, тогда как для реакций второго и третьего порядка она составляет и 0,1 и 0,2, соответственно.

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

1. Высоковольтные электротехнологии / О.А. Ано-шин, А.А. Белогловский, И.П. Верещагин и др. - М.: МЭИ, 2000.

2. Семенов В.К., Мешалкин Д.В. Экспериментальное исследование плазмохимических реакторов фронтальных волн ионизации // Изв. Вузов: Ядерная энергетика. -2000. - № 2. - С. 50-53.

3. Семенов В.К. К вопросу теплообмена в вертикальных каналах газоразрядных аппаратов // ТВТ. - 1998. -№ 3. - Т. 36. -С. 503-507.

Семенов Владимир Константинович,

ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор кафедры атомных электростанций, телефон (4932) 26-99-15, e-mail: npp@aes.ispu.ru

о