CC BY
42
ным 1/2, для элементов, смежных с ним, - по 1/4, для шести элементов, находящихся с подветренной стороны, - по 1/6.
Принимая 0 = 0,1 и учитывая вышеизложенные правила и расчетные соотношения (1), (3), (4), получим ячеечно-графическую схему клеточного автомата в условиях изменившихся метеопараметров через минимальный интервал времени А т.
Новое значение концентраций в базовой ячейке составит С0(г+ 1)= 1,098.
Предложенный клеточный автомат реализован в виде программного приложения на языке С#. Пример работы программы продемонстрирован на рис. 3, 4. Показаны исходные данные (рис. 3) и результаты моделирования процесса распространения примеси загрязняющего вещества (рис. 4) в результате гипотетического аварийного выброса на территории Новомосковской акционерной компании «Азот».
Развитие предложенного математического аппарата и программного обеспечения будет связано с внесением изменений в правила клеточного автомата для учета неоднородности подстилающей поверхности, определенного сценария изменения метеопараметров, настройкой параметров клеточного автомата для различных условий внешней среды, а также с возможностью моделирования загрязнения воздуха источниками различных типов.
УДК 519.68:544.65
В.А. Василенко, А.С. Князева, Э.М. Кольцова
Российский химико-технологический университет им. Д.И. Менделеева, Москва, Россия
РАЗРАБОТКА ПРОГРАММНОГО МОДУЛЯ ДЛЯ
МОДЕЛИРОВАНИЯ ПРОЦЕССОВ, ПРОТЕКАЮЩИХ В ТОПЛИВНЫХ ЭЛЕМЕНТАХ
A mathematical model of DMFC is developed. It taking into account electrochemical and convection processes in polymer electrolyte membrane fuel cell (PEMFC). Electrodes voltage-current curves, expressions for the profiles of proton concentration across membrane thickness and time are obtained. Values of proton diffusion coefficients, convection rate and electrochemical reactions constants are derived thanks to model.
Представлена математическая модель, учитывающая электрохимические и конвективные процессы, протекающие в ионно-обменной мембране топливного элемента с прямой подачей метанола. Получены зависимости для определения концентрации протонов по толщине мембраны и во времени, вольт-амперных характеристик электродов и поляризационных кривых. Получены значения коэффициентов диффузии протонов, скорости конвективного переноса и константы электрохимических реакций, протекающих на электродах.
В настоящее время технология на основе топливных элементов является перспективным направлением малой энергетики. Экологически чистое и эффективное производство электроэнергии привлекает все больше
разработчиков и производителей в эту область. Топливные элементы, как правило, представляют интерес, как источники тока для автотранспорта, портативных устройств и промышленных энергоустановок. Ключевой частью топливного элемента является мембранно-электродный блок. К настоящему моменту в качестве „электролита" широкое применение получила полимерная электролитная мембрана (ПЭМ). Наиболее распространённая ПЭМ - плёнка политетрафторэтилен]-перфторсульфоновой поликислоты, известная под торговым именем №йоп® фирмы БиРогй™ (США). На сегодняшний день, процессы, протекающие в ионообменных мембранах, достаточно хорошо изучены [1-2].
В данной работе исследовался транспорт протонов через ПЭМ при различных внешних условиях на примере двух типов мембран №йоп® N115 толщиной 120 мкм и N117 толщиной 183 мкм фирмы БиРот1М (США). Исследования проводились на базе Технического университета г. Карлсруэ (Германия). Часть экспериментальных исследований, проводимых в рамках данной работы, ранее были уже опубликованы [3].
Электрохимические реакции, протекающие на электродах, могут быть представлены следующим образом:
Реакция на аноде:
СН3ОН(ж) + Н20(ж) >С02(Г) + 6Н+ + бе" (1)
Реакция на катоде:
_1,502(г) +6Н+ +6е-—^ЗН20(ж)_(2)
Суммарная реакция:
СН30Н(ж)+1,502(г) —» СО
2(Г)+2Н2°(Ж)' исе11-1,21¥ (3)
Под действием платино-рутениевого катализатора на аноде метанол реагирует с водой до образования углекислого газа и протонов, при этом освобождаются электроны. Протоны диффундируют через мембрану к катоду и там, под действием платинового катализатора с кислородом воздуха реагируют до образования воды.
Массоперенос ионов водорода в мембране осуществляется посредством диффузии и конвекции. С учетом фрактальности структуры мембраны и идей фрактальной геометрии процесс может быть описан по аналогии с [4] уравнением с дробной производной по времени:
дуС . дС. . д2С. .
-н±_ + УУ!}-У = в -
д^ дх 1 Эх2 ' ( )
где С + — концентрация протонов, молъ/м , у - «порозность» фрактальной
н
структуры, v - скорость конвективного потока протонов, а/ с; х - координата по толщине мембраны в ячейке, м; - модифицированный коэффициент диффузии протонов, лГ/с, рассчитывается по уравнению, предложенному в [4]:
(5)
S - площадь кластера, открытого для протекания, лГ ; Si - площадь сечения межкластерного канала, а/2; £ — относительная концентрация "зубьев" на канале длиной /, % = S\/S.
Движущей силой конвективного массопереноса протона через мембрану является разность потенциалов на электродах (градиент электростатических напряжений dcp/dz и градиент давлений dp/dz ).
Согласно уравнению Шлёгла [5] скорость конвективного потока:
v = -k^c'" -F .^--h-.ÈL (6)
ц н+ dz |х dz' где к - электрокинетическая проницаемость мембраны, лГ; кр -
гидравлическая проницаемость мембраны, лГ; ц - вязкость жидкости в порах мембраны, кг/м/с; F - постоянная Фарадея, Кл/молъ\ z - зарядовое число реакции; R - универсальная газовая постоянная, Дж/К-молъ; Т -температура ячейки, К.
Для решения дифференциальных уравнений необходимы дополнительные условия. Искомая функция, обозначенная как концентрация протонов Сн+, является функцией времени и координат, поэтому запишем
начальное условие, характеризующее её распределение в пространстве в начальный момент времени. Для уравнения (4) оно имеет вид:
CH+(t = 0-x) = Cl+; (7)
Граничные условия для уравнения (4) учитывают образование протонов в ячейке на аноде, а также уход заряженных частиц за счёт реакции на катоде.
Левое граничное условие:
ÔC .
~ D U +v!L^CH+ |х=0 = U (8)
у их п
W, = ехр^^^ССДзОД , (9)
где w\ - скорость образования протонов по реакции на аноде, люлъ/лГ-с; к" -стандартная константа скорости электрохимической реакции на аноде, а/ с; аА - коэффициент переноса на аноде; щ - перенапряжение на аноде, В. Правое граничное условие:
" D U +v'L^CH+ U = w2 \X=L (10)
у ^ П
<П>
где м>2 - скорость образования воды по реакции на катоде, люлъ/лГ-с, к" -константа скорости электрохимической реакции на катоде, а/ с; ак - коэффициент переноса на катоде; г]К - перенапряжение на катоде, В.
Для решения уравнения с дробной производной по времени (4) использовался численный метод, предложенный в [6].
Табл. 1. Значения кинетических констант и коэффициентов диффузии
Параметры Подготовленные мембраны Неподготовленные мембраны
Nafion 115 Nafion 117 Nafion 115 Nafion 117
Т-60°С Т-70°С Т-60°С Т-75°С Т-60°С Т-60°С
Dg, M2/cg 8,93* Ю-9 1,0* Ю-9 1,55* 10~9 1,62* Ю-9 2,65* Ю"10 2,24* Ю"10
К1, м/с 2,7* 10"4 2,8* 10"4 2,7* 10"4 2,7* 10"4 2,7* 10"4 2,7* 10"4
К2, м/с 3,0* 10"6 3,5* 10"6 3,0* 10"6 3,7* 10"6 3,0* 10"6 3,0* 10"6
Vg, м/с'3 7,8* 10"6 7,3*10"6 5,81* 10~6 5,81* 10"6 1,18*10 5 1,22*10"5
На основе математической модели был разработан генетический алгоритм поиска кинетических параметров процесса: констант скорости электрохимических реакций на электродах, модифицированного коэффициента диффузии протонов водорода, "порозность" мембраны. Алгоритм написан на языке программирования С++. Включает в себя параллельные и последовательные вычисления. Программа работает под операционными системами семейства Unix.
Табл. 2. Значения перенапряжения и силы тока
Параметры Подготовленные мембраны Неподготовленные мембраны
Nafion 115 Nafion 117 Nafion 115 Nafion 117
Т-60°С Т-70°С Т-60°С Т-75°С Т-60°С Т-60°С
>7.4, В 0.096 0.098 0.096 0,100 0.095 0,095
т, в 0,067 0,068 0,067 0.069 0.065 0,064
/', мА/см2 339 342 349 349 316 344
/, А 7,81 7,86 8,03 8.03 7,26 7,91
В результате моделирования были получены значения кинетических констант и коэффициентов диффузии для различных типов мембран (см. табл. 1), зависимости изменения концентрации протонов во времени и по толщине мембраны.
Из таблицы 1 видно, что коэффициенты диффузии у подготовленных мембран выше, чем у неподготовленных, причем у мембраны Nafion N117 больше, чем у Nafion N115. Таким образом, использование подготовленной мембраны Nafion N117 является предпочтительным. В таблице 2 приведены данные для рассчитанных значений перенапряжения и силы тока для различных мембран.
4000 1
3500 —-3000 2500 —-
Ъ
5 2000
1500 1000
500 —-
0--
0,00
перенапряжение, В
Рис. 1. Катодные и анодные поляризационные кривые неподготовленной мембраны N115
Из таблицы 2 видно, что в случае использования подготовленных мембран N 117 сила тока, вырабатываемого ТЭ тока выше, чем при использовании неподготовленных мембран.
На рис. 1 приведены поляризационные анодные и катодные кривые для неподготовленной мембраны N 115. Из графика видно, что анодный ток в 2 раза меньше катодного, что говорит о падении протонной проводимости с уменьшением влажности материала. Эта зависимость определяется тем фактом, что протонная проводимость мембраны осуществляется благодаря диссоциации сульфогрупп в присутствии воды, входящей в ее структуру, с последующим переносом образовавшихся гидратированных форм протона по протонгидратной оболочке. Проводимость такой мембраны в значительной степени определяется количеством адсорбированной воды [2].
Было рассчитано омическое сопротивление мембраны, даны рекомендации по поиску протонообменных мембран с наибольшим выходом по току.
Работа выполнена в рамках государственного контракта № 16.513.11.3039 от 12 апреля 2011 г.
Библиографические ссылки
1. Ярославцев, А.Б. / А.Б. Ярославцев, В.В. Никоненко // Российские нано-технологии, 2009. Том 4. № 3-4. С. 44.
/
анод /
катод
/
У
2. Иванчёв, С.С. / С.С. Иванчёв, С.В. Мякин// Успехи химии, 2010..№ 79. Том 2. С. 118.
3. Philip Scharfer /Philip Scharfer, Wilhelm Schabel, Matthias Kind.// Journal of Membrane Science. V. 303 (2007). P. 37.
4. Nigmatullin R R. Phys. stat. sol. (b), 1986. V. 133. P. 425.
5. Schlögl, R./ R.Schlögl // Berichte der Bunsengesellschaft für Physikalishe Chimie, 1966. № 70/ S. 400-414.
6. Кольцова, Э.М./ Э.М. Кольцова, В.А. Василенко, B.B. Тарасов // Журн. физ. Химии, 2000. Т. 74. № 5. С. 954.
УДК 502.55.001.18
А.П. Виноградов, Т.В. Савицкая, A.B. Горанский
Российский химико-технологический университет им. Д.И. Менделеева, Москва, Россия
БАЗА ДАННЫХ ПО ПОЖАРО- ВЗРЫВОБЕЗОПАСНОСТИ ОПАСНЫХ ПРОИЗВОДСТВЕННЫХ ОБЬЕКТОВ
The database's structure on the dangerous industrial objects's fire and explosion hazard is considered. The database functional purpose is given. Examples of information filling are introduced.
Рассмотрена структура базы данных по пожаро- взрывобезопасности опасных производственных объектов. Дано функциональное назначение базы данных. Приведены примеры информационного наполнения
На предприятиях химической, нефтехимической и нефтеперерабатывающей промышленности постоянно осуществляется интенсификация технологических процессов и производств, связанная с укрупнением единичных мощностей установок и аппаратов, а в результате наличием в них больших запасов взрыво-, пожароопасных веществ, возрастанием температур и давлений протекания процессов. Отсюда, как следствие, возможен рост числа аварий с все более тяжелыми социальными, экономическими и экологическими последствиями.
Пожарная безопасность и взрывобезопасность опасных производственных объектов (ОПО) регламентируются различными нормативными документами, но применительно к свойствам опасности технологических процессов химически-опасных объектов (ХОО) они должны рассматриваться комплексно, исходя из свойств веществ (технологических сред), параметров технологических процессов, совместного действия поражающих факторов пожаров и взрывов.
Для систематизации и хранения данных о пожаро- и взрывобезопасности технологических процессов, установок и оборудования химически
