Модель термодесорбции водорода из металлоорганической структуры

Тип работы:
Реферат
Предмет:
Физика


Узнать стоимость

Детальная информация о работе

Выдержка из работы

Труды Карельского научного центра РАН № 1. 2013. С. 33−37
УДК 519. 6:539. 2
МОДЕЛЬ ТЕРМОДЕСОРБЦИИ ВОДОРОДА ИЗ МЕТАЛЛООРГАНИЧЕСКОЙ СТРУКТУРЫ
Ю. В. Заика
Институт прикладных математических исследовании Карельского научного центра РАН
Рассматривается дегазация методом термодесорбции (ТДС) металлоорганического образца, предварительно насыщенного водородом. Система мелких пор является каналом переноса, крупные поры играют роль ловушек с обратимым захватом. Представлена математическая модель ТДС-эксперимента.
Ключевые слова: водородопроницаемость, моделирование, краевые задачи.
Yu. V. Zaika. MODEL OF HYDROGEN THERMODESORPTION FROM METAL-ORGANIC FRAMEWORKS
Degassing by the method of thermal desorption (TDS) of an organometallic sample pre-saturated with hydrogen is considered. The system of small pores is the channel for migration, and large pores act as traps with reversible capture. The mathematical model of the TDS-experiment is presented.
Key words: hydrogen permeability, modeling, boundary-value problems.
Введение
Исследования в области водородного материаловедения инициированы прежде всего перспективами водородной энергетики. Рассмотрим прикладную задачу, связанную с проблемой накопления и хранения водорода. Одними из перспективных материалов в этом отношении являются металлоорганические структуры (metal-organic frameworks, MOF)[4, 1, 3], имеющие пористую структуру. В [2] представлена модель переноса водорода в пористом материале. В качестве приложения рассмотрена водоро-допроницаемость пористого вольфрама. Поры считаются настолько большими, что накапливающийся в них водород является
газообразным и приходится оперировать параметром давления Н2 в порах. В MOF характерные размеры пор в диаметре значительно меньше (исчисляются ангстремами — АА), что требует корректировки модели. Приведем подходящую математическую модель для конкретного эксперимента: по-
сле предварительного насыщения водородом материал дегазируется методом термодесорбции. Канал переноса (диффузии) молекул Н2 — система мелких
пор, соизмеримых с размерами молекулы. Большие поры изолированы в том смысле, что молекула Н2 не может попасть сразу из одной большой поры в другую. Поэтому систему больших пор считаем ловушками-накопителями. Вместе с тем, они относительно малы (несколько десятков, А в диаметре).
¦@
Образец является двухслойной пластиной из подложки и исследуемого материала. Образец насыщается водородом при комнатной температуре и медленно охлаждается до криогенной температуры. После установления равновесия происходит в режиме вакуумиро-вания медленный нагрев со стороны подложки, с помощью масс-спектрометра регистрируется выходной десорбционный поток водорода.
Математическая модель
Пренебрегаем торцами и считаем подложку практически непроницаемой для молекул Н2 со стороны слоя МОР. Обозначим через Б площадь одной стороны поверхности. По толщине слоя МОР выделим отрезок [х, х + йх] и рассмотрим элемент объема йУ = Бйх. В дальнейшем по контексту под Б, У понимаем как сами геометрические фигуры, так и численные значения их площади и объема. Введем обозначения, следуя [2]: пр — концентрация больших пор, [пр] = 1/еш3- Ур — объем поры- Бр — ее площадь поверхности. В силу регулярности структуры МОР эти величины являются константами в образце. Тогда прУр — объем пор в единичном кубике, т. е. численно это объемная доля пор в материале. Эта безразмерная величина называется пористостью. Обозначим ее П = прУр. Величина прБр имеет смысл общей поверхности пор в единице объема. Эту объемную плотность поверхности обозначим Б.
Перейдем к уравнениям материального баланса. Объем йУ разбивается на две части: ПйУ — объем пор, (1-ц)йУ — объем, в котором осуществляется диффузионный перенос (система мелких пор). Обозначим через, 1р+ плотность потока молекул Н2 из поры: столько молекул пересекает единичную площадку сферы большой поры изнутри в единицу времени. Аналогично определяем величину — это плотность встречного потока диффундирующего водорода, попавшего в пору-накопитель. Через с (Ь, х) обозначим концентрацию диффундирующего Н2 в системе мелких пор (в канале переноса) в момент времени Ь на глубине х по толщине слоя МОР.
Далее в выкладках опускаем бесконечно малые о (АЬ), о (Ах) в силу последующего предельного перехода АЬ ^ 0, Ах ^ 0. Запишем изменение количества молекул Н2 в объ-
еме (1 — п) йУ в единицу времени (скорость):
(1 — п) йУАс/ АЬ ^ д^(Ь, х)(1 — п) йУ (АЬ ^ 0).
Найдем выражение этой величины через баланс потоков. Площадь плоского сечения равна Б. Саму геометрическую фигуру и ее площадь для простоты обозначаем одним символом. Сечение делится на две части: Бd — через эту площадь идет диффузия, (Б — Ба) — общая площадь кругов, являющихся сечениями крупных пор. Для численного значения Ба имеем Ба = (1 — п) йУ/йх = (1 — п) Б. Через плоское сечение с координатой х (по толщине слоя МОР) входит диффузионный поток Н2 плотности -(Бдс/дх)|х, а через сечение х + Ах выходной поток имеет плотность -(^дс/дх)|х+Ах. Умножив упомянутые плотности на Б а, получим потоки Н2 по каналу мелких пор сквозь сечения х и х + Ах. Разность -Бдс/дх1хБа — (-^дс/дх|ж+дж)Ба определяет изменение количества молекул Н2 в объеме йУ в единицу времени за счет диффузии. Кроме того, молекулы водорода попадают в диффузионный канал переноса из больших пор и захватываются обратно из канала мелких пор. Дисбаланс этих двух встречных потоков равен (,]+ - ,]~)БрпрБ (1х. Здесь прБйх — общее количество больших пор в объеме йУ = Бйх, так что БрпрБйх = БйУ — это общая поверхность больших пор в объеме йУ. Для МОР Б «1.
Составляем баланс молекул Н2 в объеме канала диффузии (1 — п) йУ = (1 — п) Бйх (с учетом АЬ = йЬ ^ 0, Ах = йх ^ 0):
дсБайх = (Ь, х)(1 — п) йУ =
дЬ дЬ
(вд)
V дх) х+Ах V дх) х
Б (IV
+ - '-^р) Брпр ¦ Бйх.
Подставим Ба = (1-п)Б, поделим на йУ = Бйх и устремим йх к нулю: Б = прБр,
Поделив дополнительно на величину (1 — п), окончательно получаем: Т = Т (Ь, х),
§ = дх (°{Т) дхх) + (¦'* - т-п ¦
0
Для определенности зависимость «физикохимических» коэффициентов модели от температуры Т считаем аррениусовской. В частности, в приведенном уравнении В (Т) =
ехр{-Ев/[КТ]}, где Ев — энергия активации, К — универсальная газовая постоянная.
Аналогично составляем материальный баланс для объема пйУ, занятого большими порами (ш — концентрация Н2 в них):
= - (4 — 4) Б * ж = -($ - ъ) Б.
В качестве модели плотностей потоков 3± сквозь поверхность больших пор примем
7+ = к+(т)w (Ь, х) [1 — c (Ь, х) стах,
7- = к-(Т)с (Ь, х) [1 — ш (Ь, х) штЦ ¦
Величина в+ = с/стах имеет смысл степени насыщения канала диффузии. Чем больше концентрация в больших порах и меньше степень насыщения малых пор, тем больше плотность выходного потока молекул Н2 из поры-ловушки в канал диффузии. Соответствующий коэффициент пропорциональности к+ зависит от текущей температуры Т. Аналогично интерпретируется выражение для плотности встречного потока в поры-накопители.
Подчеркнем, что в уравнения переноса входит не только пористость материала п, но и объемная плотность поверхности больших пор Б. При одной и той же пористости величина Б для МОР может варьироваться в широких пределах. Перенос тепла диффундирующими молекулами (эффект Дюфо) и термодиффузию (эффект Соре) считаем второстепенными.
Перейдем к граничным условиям. Подложка практически непроницаема со стороны МОР: дхс (Ь, 0) = 0, Ь ^ 0. Значение х = 0 соответствует началу слоя МОР, Ь = 0 — началу нагрева образца в вакууме. Выходная сторона слоя МОР обращена в вакуум и вследствие пористости материала имеет изрезанный вид. Можно ввести параметр, а шероховатости поверхности [2], равный отношению «истинной» площади поверхности к геометрической. Но численно значение, а оценить трудно. Поступим по-другому. Выделим тонкий приповерхностный слой толщины ?, который в силу изрезанности существенно отличается по физико-химическим свойствам от
регулярных внутренних слоев общей толщины L (? ^ L). Кроме того, в процессе изготовления материала, подготовки и проведения эксперимента неизбежно появление микропримесей и окислов на поверхности. Фактически этот слой является несколько иным материалом по сравнению с объемным. Внешний слой обладает накопительным потенциалом: молекула H2 имеет дополнительную возможность зацепиться за неровности поверхности. Рассмотрим сечение x = L. Площадь этого сечения равна S. Из объема подходит диффузионный поток -SdD (T)дхе1ь. Молекулы H2 могут не сразу десорбироваться, а предварительно «растекаться» в поперечном направлении. Таким образом, десорбция в вакуум идет «сквозь окно» площади S. Выделенный иррегулярный внешний слой настолько мал, что в его пределах объемная концентрация молекул водорода практически не меняется по толщине x? [L, L + ?]. Обозначим ее ce (t). Поскольку имеем дело с физической адсорбцией (абсорбцией в ?-слой), то плотность потока десорбции моделируем реакцией первого порядка: b (T)ci (t). Эффективный коэффициент десорбции b имеет размерность скорости (cm/s). Дисбаланс потоков определяет изменение концентрации в приповерхностном слое:
SicfXt) = -SdD (T)dxc (t, L) — Sb (T)ce (t).
Остается поделить на S и заменить отношение Sd/S на (1 — n). В [2] в аналогичной ситуации (но в больших порах-полостях находится газ) принято ci (t) = c (t, L), что при? ^ 1 соответствует относительно малой поверхностной концентрации q (t) = ?ci (t). Можно расширить возможности модели, полагая ce (t) = h (T)c (t, L) и допуская h «1: к диффундирующему водороду добавляется водород из открытых приповерхностных пор. Обозначая? h = 1/g, формально приходим к динамическому граничному условию, аналогичному условиям баланса потоков
co (t) = g (T)qo (t), ce (t) = g (T)qe (t),
qo (t) = уs (T)po (t) — b (T)q^(t) + D (T)cx (t, 0) ,
qe (t) = ys (T)pe (t) — b (T)q2 (t) — D (T)cx (t,?)
для экспериментального метода проницаемости H сквозь пластину толщины? (у — кинетический коэффициент, s — коэффициент
¦(c)
прилипания, р — давление газообразного водорода). Такая форма пригодна и в ситуации, когда поверхностная концентрация больше объемной (при относительно низких температурах). В любом случае вместе с начальными условиями краевые задачи приобретают замкнутый вид. Если накопительный эффект пренебрежимо мал, то полагаем левую часть, равной нулю, и получаем граничное условие (п — 1)0дхс (Ь, Ь) = Ьс (Ь, Ь): сколько молекул Н2 диффузионно подошло к выходной поверхности, столько и десорбировалось за единицу времени. Можно подразумевать, что в эффективный коэффициент десорбции Ь «включены» как параметр Н, так и шероховатость выходной поверхности.
Кратко остановимся на начальных условиях. В режиме насыщения образца следует учесть наличие внешнего давления газа (Н2):
1с?(Ь) = ув (Т)р (Ь) — - (1 — П)0(Т)дXc (Ь, Т) — Ь (Т)с?(Ь).
При начальной температуре Т = Т эксперимента установилось равновесие, поэтому все производные равны нулю: увр = Ьс/. В объеме получаем: сд = с/,
*к+ш [1-сстах =к-с[1-шш-ах.
Модель позволяет однозначно определить начальные данные по цепочке р ^ с/ ^ с ^ Ш. Как оценить значения стах, штах? Рассмотрим единичный кубик с максимальным насыщением. Количество молекул Н2 в нем равно, А = пштах + [1 — п] стах. Учтем соотношения П = прУр = пр4пт3/?з = ПрБрт/3 = Бт/3, где пр — количество пор в единичном кубике, т — радиус большой поры, Б — площадь ее поверхности, Ур — объем, Б = прБр — объемная плотность общей поверхности больших пор. Тогда
А = штахБт/3 + стах стахБт/3 =
= стах + & gt-Б[Штах стах]т/3-
Зафиксируем параметр ББ и будем менять т в физически реальном диапазоне. Если уменьшить т, то, чтобы сохранить значение Б, нужно соответствующим образом увеличить количество пор. Кроме того, в более мелких порах молекуле легче удержаться, так что имеется тенденция увеличения штах с уменьшением т. Но тогда коэффициент т [штах — стах] при Б
имеет тенденцию к стабилизации, слабой зависимости от г. Практически линейная зависимость A = kS + k2 (hi = const) при варьировании S обнаружена экспериментально [3]. Тогда k2 = cmax, (ki, A) O (wmax, r). Можно ограничиться некоторыми средними значениями cmax, wmax или грубыми оценками по результатам полной дегазации с учетом того, что из геометрических соображений количество малых шаров (молекул H2) по объему в большой поре радиуса г не должно превысить 4пг3/3. Если в условиях эксперимента c ^ cmax, w ^ wmax, что согласуется с D = D©, то исключаем параметры cmax, wmax из модели: c/cmax «0, w/wmax «0.
Теплопередача в металлоорганических каркасах происходит значительно быстрее диффузии. Поэтому вместо распределенного уравнения теплопроводности целесообразно оперировать квазистационарами. Если в уравнении теплопроводности vdtT = XdxxT считать изменения по времени медленными, то в относительном масштабе получаем dtT ~ 0 и линейную зависимость T = ах + в .С течением времени квазистационар перестраивается: T (t, x) = a (t)x + e (t). На границе слоев два линейных квазистационара стыкуются с изломом из-за различных коэффициентов теплопроводности. Температура эксперимента низка, так что тепловым излучением в вакуум, по-видимому, можно пренебречь. Дальнейшая детализация зависит от набора теплофизических измерений в эксперименте. Для достаточно тонкого образца для начала можно принять, что прогрев равномерный.
Изложенное следует воспринимать лишь как один из вариантов модели. Различных MOF синтезировано уже более тысячи и говорить о достаточно общей модели не приходится. Приведем одну из модификаций представленной модели с учетом необходимости решать обратную задачу параметрической идентификации. Это приводит к необходимости учета лимитирующих факторов минимальным набором неизвестных априори параметров. Коэффициенты имеют смысл обобщенных (интегральных, эффективных) показателей переноса водорода в слое MOF. В частности, диффузию не разделяем на молекулярную и кнудсеновскую, считая D = Deff. Рассмотрим MOF со следующими свойствами. Имеется две системы пор: большие (~ 10А)
0
и малые (~ 5A). Грубо диаметр молекулы водорода ~ 3A. Система больших пор взаимосвязана и является каналом диффузионного переноса (диффузионного в формальном «фиковском» смысле: поток пропорционален градиенту концентрации). Большие поры имеют кубическую форму, малые являются «треугольными карманами» в вершинах куба. Геометрическая и химическая структура каркаса малых пор позволяет считать их ловушками водорода. В «щели» при наличии ионов металла молекула H2 связана с каркасом значительно сильнее, чем в большой поре. Пор огромное число, так что для анализа макропереноса следует ввести подходящие макропараметры. Ограничимся лишь двумя геометрическим характеристиками: пористость n и удельная поверхность S ([•] = cm2/cm3 = 1/cm). Общий баланс подкачки водорода из ловушек в процессе равномерного монотонного (обычно линейного) нагрева T = T (t) моделируем соотношением J+ - J- = a (T)w (t) ([a] = cm/s), где w (t) — концентрация в объеме малых пор. Из материального баланса получаем: c 2c (1 — n) g-t =(1 — n) D (T)2 + a (T)w (t)S, ndw/dt = -a (T)w (t)S.
После интегрирования второго линейного обыкновенного уравнения получаем одно дифференциальное уравнение для c (t, x). Начальное насыщение по постановке ТДС-равномерно: c (0,x) = a, w (0) = w, x? [0, L]. Пренебрегая накоплением на выходе, прини-
СВЕДЕНИЯ ОБ АВТОРЕ:
Заика Юрий Васильевич
зав. лаб. моделирования природнотехнических систем, д. ф. -м. н.
Институт прикладных математических исследований Карельского научного центра РАН ул. Пушкинская, 11, Петрозаводск,
Республика Карелия, Россия, 185 910 эл. почта: zaika@krc. karelia. ru
маем граничные условия
dxc (t, 0) =0, (n — 1) Ddxc (t, L) = bc (t, L).
Ловушки, обращенные в вакуум (x = L), опустошились при t & lt- to ^ 1. Неизвестные значения Do, Ed, a0, Ea, b0, Eb оцениваются по плотности десорбции J (t) = b (T (t))c (t, L), часто известной лишь в относительных единицах (J/Jmax). Для первого приближения (оценка значений Do, Ed, С) можно принять
a = 0, T & lt- Tcrit, a = a & gt- 0, T ^ Tcrit, c (t, L) = 0, J = (n — 1) D xc (t, L).
При необходимости нетрудно учесть в форме ловушки определенной емкости накопление водорода на границе с подложкой.
ЛИТЕРАТУРА
1. Исаева В. И., Кустов Л. М. Металлоорганические каркасы — новые материалы для хранения водорода // Российский химический журнал. 2006. № 6. С. 56−72.
2. Писарев А. А., Цветков И. В., Марен-ков Е. Д., Ярко С. С. Проницаемость водорода через металлы. М.: МИФИ, 2008. 144 с.
3. Hirscher M., Panella B., Schmitz B. Metal-organic frameworks for hydrogen storage // Microporous and mesoporous materials. 2010. Vol. 129. P. 335−339.
4. Varin R.A., Czujko T., Wronski Z.S. Nanomaterials for solid state hydrogen storage. Springer, 2009. P. 338.
Zaika, Yury
Institute of Applied Mathematical Research,
Karelian Research Centre, Russian Academy of Sciences 11 Pushkinskaya St., 185 910 Petrozavodsk, Karelia, Russia
e-mail: zaika@krc. karelia. ru tel.: (8142) 766 312
тел.: (8142) 766 312

ПоказатьСвернуть
Заполнить форму текущей работой