О методике определения кинетических параметров окисления углеродных измельченных материалов

Тип работы:
Реферат
Предмет:
Механика


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

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

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

Серия «Химическое машиностроение и инженерная экология» Литература
1. Носков С. А., Баранов Д. А., Скопинцев И. В., Шибанов А. В. Определение краевого угла смачивания композиций полимерных составов для оросителей градирен // Известия МГТУ «МАМИ», № 3(17), 2013, т. 2, с. 19 — 22.
2. Скопинцев И. В., Мелешкина А. М., Камшад Ф. Новое применение вторичных полимерных композиционных материалов // Известия МГТУ «МАМИ», 2012, № 2(14), т. 4, с. 197 — 201.
3. ГОСТ 11 645–73 — Пластмассы. Метод определения показателя текучести расплава термопластов.
О методике определения кинетических параметров окисления углеродных измельченных материалов
К.ф. -м.н. Обвинцева Н. Ю., д.х.н. проф. Каминский В. А., д.т.н. Эпштейн С. А.
Национальный исследовательский технологический университет «МИСиС»
obvint@gmail. com
Аннотация. Рассмотрены особенности термических методов определения кинетических параметров процесса окисления углеродных материалов. Предложен модифицированный метод определения кинетических параметров на основе величин, характеризующих стационарное состояние системы.
Ключевые слова: процесс окисления угля, кинетические параметры, параметр Франк-Каменецкого, критические условия
В связи с возможностью самовозгорания углеродных сыпучих материалов (угля, сельскохозяйственной продукции и др.) при их хранении и использовании представляет интерес изучение процесса окисления углерода при низких температурах. Для определения параметров, соответствующих условиям самовозгорания, разработаны методы, которые являются общими при изучении возгорания разных материалов. В данной работе эти методы будут рассмотрены в связи с проблемой самовозгорания угля. Для прогнозирования поведения различных типов углей необходимо знать кинетику протекающих химических реакций и физико-химических процессов и условия теплообмена системы с окружающей средой в естественных условиях, поскольку известно, что самовозгоранию углей способствует аэрация штабеля, наличие внешних источников тепла (солнечная радиация), наличие посторонних примесей, смешение разных марок и сортов углей. Учет такого большого числа перечисленных факторов при определении кинетики реакций является достаточно сложной задачей.
Исследование самовозгорания углей, основанное на кинетическом подходе, направлено на описание характерных реакций и оценку вклада каждой из них при разных температурах. В литературе приводятся отдельные данные о влиянии перечисленных факторов на химическую активность угля (см., например, обзоры [1, 2]). Однако в настоящее время нет полной теории, описывающей процесс саморазогрева угля.
При изучении разогрева углей на основе термических методов с целью упрощения весь процесс описывается одной экзотермической брутто-реакцией [3 — 5]. Если определить эффективные кинетические параметры этой реакции, протекающей в естественных условиях, то этого вполне достаточно, чтобы прогнозировать поведение углей.
Классическими термическими методами определения кинетических констант, характеризующих интенсивность экзотермических процессов, приводящих к разогреву, являются методы, основанные на стационарной теории Франк-Каменецкого [6]. С помощью этих методов [3 — 5] при заданной геометрии системы и известных теплофизических параметрах на основе данных термического эксперимента можно определить кинетические константы брутто-реакции. Для наблюдения разогрева в лабораторных условиях используются небольшие объемы угля, при этом его предварительно измельчают до размера частиц 0. 25 — 1 мм.
Серия «Химическое машиностроение и инженерная экология»
Чтобы наблюдать разогрев за счет экзотермической реакции, реактор с образцом угля помещают в печь при температуре выше 100 °C, температура печи поддерживается постоянной. В эксперименте проводятся измерения температуры в некоторых точках образца. В стационарном методе определяется критическая температура на поверхности реактора, при которой происходит резкое повышение температуры в его центральной точке [3]. Этот метод требует проведения серии экспериментов для реакторов разных размеров. Уменьшить количество экспериментов предлагается с помощью нестационарного метода Чена [4]. В этом методе измеряется скорость изменения температуры в центральной точке образца. Наиболее перспективным термическим методом является нестационарный метод, представленный в работе [5]. На основе установленной зависимости между величиной безразмерной температуры & quot-crossing point& quot- и параметром Франк-Каменецкого предложена процедура определения кинетических параметров.
В данной работе проводится сравнительный анализ методов определения кинетических параметров на основе данных термического эксперимента и предлагается модифицированный метод определения кинетических параметров по температуре стационарного разогрева в центральной точке реактора.
Процесс разогрева угля описывается нестационарным уравнением теплопроводности с внутренним источником тепловыделения:
РСР? = kv2T+QW (T), (1)
dt
где k — коэффициент теплопроводности, Ср — удельная теплоемкость. В уравнении первый член в правой части описывает процесс теплопроводности, второй — тепловыделение за счет экзотермической реакции окисления угля. Принимается, что скорость реакции зависит экспоненциально от температуры по закону Аррениуса:
__E_
W (Т) = Ae RgT, (2)
где A — предэкспоненциальный множитель, в который включены концентрации реагентов, Rg — газовая постоянная, Q — теплота процесса, выделяющаяся в результате брутто-реакции, E — энергия активации. Расходом исходных реагентов пренебрегается.
В зависимости от условий нагрева, теплофизических характеристик и размера реактора с течением времени в реакторе устанавливается либо стационарный температурный режим, либо происходит резкое повышение температуры, которое в итоге приводит к самовоспламенению угля. Условия, при которых стационарное распределение температуры в реакторе становится невозможным, называются критическими условиями. Определить критические условия можно с помощью метода Франк-Каменецкого — разложения экспоненты в ряд [6]. Уравнение (1) приводится к безразмерному виду, задавая в качестве масштаба длины характерный размер образца R, масштаба времени — R2/a (а — коэффициент температуропроводности), безразмерную температуру — в = E (T-Ts)/RgTs2:
дв дт
где Ag — оператор Лапласа, записанный в безразмерных переменных,
— = Axe + dexp (e), (3)
? ER 2B (E Л о =-- exp--
RgT2a Ч RgT 0
— безразмерный параметр Франк-Каменецкого,
В = QA|Cp р — параметр, определяющий скорость возрастания температуры за счет экзотермической реакции окисления угля.
Граничные условия к уравнению (3): на поверхности образца 08 = 0.
Видно, что распределение безразмерной температуры в определяется единственным
Серия «Химическое машиностроение и инженерная экология» параметром ё. Параметр ё отражает совокупность всех свойств системы: скорость и тепловой эффект реакции, теплопроводность, геометрические размеры реактора. При некотором значении 8с достигается критическое состояние, при котором стационарное распределение температуры становится невозможным. При значениях ё& lt-ёс — в системе устанавливается стационарное распределение температуры. При ё& gt-ёс — теплота не успевает отводиться в окружающую среду и температура в зоне реакции начинает повышаться, что приводит к самовоспламенению.
В последнее время метод Франк-Каменецкого используется для определения эффективных кинетических параметров — энергии активации и предэкспоненциального множителя. Для сравнения методов определения этих параметров в данной работе выполнен численный эксперимент по разогреву образца угля с известными кинетическими константами. На основе полученных данных проводится описание процедуры восстановления кинетических параметров в соответствии с разными методами. Численные расчеты уравнения (1) выполнены для цилиндрического реактора при равных значениях его высоты и диаметра И = 0,1 м, ё = 0,1 м. Теплофизические параметры взяты из [4]: р = 1350кг/ж3, к = 0,23 Вт/(мК), С = 1132 ДжЦкг • К), а = 1,5−10−7 м2/с, & lt- = 22. 9−106 Дж/кг, Е=101,7103 Дж/моль.
Рисунок 1. Определение критической температуры в стационарном методе: 1 — 7& gt-110 °С, 2 — 7& gt-113 °С, 3 — 7& gt-114 °С, 4 — 7& gt-120 °С
Традиционным методом определения кинетических параметров является стационарный метод. На основе результатов серии экспериментов определяется температура поверхности образца, при которой достигаются критические условия Т5сг. Для этого проводятся измерения температуры в центральной точке реактора от времени Т0. Такие зависимости, полученные на основе численных расчетов при разных значениях температуры реактора, приведены на рисунке 1. Внешняя температура задавалась равной 110 °C, 113 °C, 114 °C, 120 °C. Начальный прямой участок кривых соответствует прогреву угля за счет процесса теплопроводности. При возрастании температуры в образце добавляется тепловыделение за счет химической реакции. При внешней температуре 114 °C и выше наблюдается резкое возрастание температуры в центральной точке реактора (кривые 3, 4). На основании этих данных значение критической температуры для данной системы равно 113,5±0,5 °С. В стационарном методе таким же образом необходимо определить критические температуры для нескольких реакторов разных размеров. Для каждого из них предварительно рассчитывается критический параметр 8С и на основании полученных данных устанавливается соотношение между
Т и параметром 8С. Определить энергию активации и предэкспоненциальный множитель
Серия «Химическое машиностроение и инженерная экология»
можно из графика линейной зависимости 1п (?^Т^/Я2) от уТ8сг. Использование стационарного метода требует многократных измерений температуры при определении критической температуры стенки реактора для реакторов разного размера.
Уменьшить количество измерений для оценки эффективных кинетических параметров окисления угля было предложено с помощью нестационарного метода (метод Чена), описанного в [4]. В этом методе используется факт, что в процессе разогрева угля наступает момент времени, когда член, описывающий теплопроводность в уравнении (1), обращается в нуль. Уравнение теплопроводности в цилиндрической системе координат имеет вид:
рСп — = к
р дг
1 д
г дг
дТ
дг
+ -
д2ТЛ
д22
+л& lt-2в
(4)
где г — координата, направленная по радиусу цилиндра, 2 — по высоте. На рисунке 2 показаны профили температуры вдоль радиальной составляющей на высоте половины реактора в разные моменты времени, полученные на основе численных расчетов. Внешняя температура задавалась 120 °C. Видно, что наступает момент времени, когда профиль температуры в реакторе вблизи его центра становится практически & quot-плоским"-. Следовательно, радиальная составляющая члена теплопроводности в уравнении (4) в этот момент времени обращается в нуль. Величиной вертикальной составляющей в этот момент пренебрегают.
Рисунок 2. Профили температуры в образце при разогреве (7^=120 °С)
Из уравнения (4) для момента времени, при котором образуется & quot-плоский"- профиль,
можно записать:
1п I —
дг
ал
= 1п
Т=Т0 рС
Е
р

(5)
На основе этого соотношения можно определить искомые параметры, А и Е из графика линейной зависимости 1п (дТ/ дг) тт от (1/То). Для этого требуется определить скорость изменения температуры (дТ/дгг в центральной точке реактора при & quot-плоском"- профиле температуры при разных значениях температуры, поддерживаемой на стенках реактора. Такие зависимости в [4] были получены для реакторов разных размеров.
В работе [5] для определения кинетических параметров в работе установлена зависи-
Е
Серия «Химическое машиностроение и инженерная экология» мость между величиной безразмерной температуры и параметром Франк-Каменецкого S. Такая зависимость получена для плоского слоя. Экспериментально измеряются значения температуры в центре реактора и в точке, расположенной рядом, в той же горизонтальной плоскости. Таким способом определяется температура соответствующая & quot-плоскому"- профилю, которая называется & quot-crossing point temperature& quot-. С помощью полученной зависимости можно определить эффективные кинетические параметры как минимум на основе двух экспериментальных значений температуры при & quot-плоском"- профиле.
Для данного метода требуется небольшое количество экспериментов по измерению температуры, что достаточно удобно. Однако, величина температуры & quot-crossing point& quot- мало отличается от температуры на поверхности образца и более удобным в качестве характерной температуры представляется использовать температуру в центральной точке реактора при выходе на стационарный режим. Это видно по приведенным на рисунке 3 результатам расчетов, где показано распределение температуры в образце, соответствующее & quot-плоскому"- профилю и установившееся в стационарном режиме. Внешняя температура задавалась Ts=110 °C.
Рисунок 3. & quot-Плоский"- профиль температуры и профиль в стационарном режиме
(7& gt-110 °С)
В связи с этим можно предложить модифицированный метод определения кинетических параметров реакции окисления угля по температуре стационарного разогрева на основе данных термического эксперимента — достаточно определить температуру в центральной точке реактора (Т0) при установлении стационарного режима. Измерения требуется провести
несколько раз при разных значениях температуры на поверхности реактора (Т).
Восстановление кинетических параметров Е и В по экспериментальным данным сводится к следующей процедуре.
Для системы с заданной геометрией необходимо установить зависимость между безразмерной температурой в центральной точке образца в стационарном режиме (00) и параметром Франк-Каменецкого 5. Для этого проводится численное решение уравнения (3) при граничных условиях на поверхности образца = 0. Полученные значения температуры в0 в
зависимости от параметра 3 для цилиндрического реактора с соотношением высоты к диаметру И/ё=1 приведены в таблице 1. По данным этой таблицы строится график зависимости 5 = / (00). Полученная зависимость Г (5,в0) = 0 хорошо описывается полиномом четвертой
степени 5 = 5. 7192•Ю-4 + 4. 96 576>-0 -3. 3852<-902 +1. 11 692<-903 -0. 16 556<-904. Чтобы определить искомые параметры Е и В последнее выражение требуется записать в размерном виде и
Серия «Химическое машиностроение и инженерная экология»
подставить значения температуры в центральной точке, полученные из эксперимента, и на поверхности образца. Таким образом, составляется система алгебраических уравнений относительно искомых параметров Е и В.
Таблица 1
8 во 8 во
0 0 2.1 0. 664
0. 25 0. 052 2.2 0. 7242
0.5 0. 1082 2.3 0. 792
0. 75 0. 1693 2.4 0. 8704
1 0. 2365 2.5 0. 9644
1. 25 0. 3112 2.6 1. 0842
1.5 0. 3956 2. 65 1. 1617
1. 75 0. 4932 2.7 1. 2615
2 0. 6098 2. 75 1. 4143
Рисунок 4. Зависимость 8 от 00
Применим данную процедуру восстановления кинетических параметров для цилиндрического реактора. Для этого необходимы данные хотя бы двух термических экспериментов. Пусть в эксперименте на поверхности образца заданы 7^=100 °С и Т2 =110 °С. Соответствующие значения температуры в центре реактора в условиях стационара — Т01 и Т02 необходимо определить экспериментально. В данной работе эти значения получены на основе численного решения уравнения (1). При расчетах для реактора определенных размеров задавались кинетические и теплофизические параметры, эти значения приведены ранее. Рассчитанные значения температуры равны Т01 =102,9 °С и Т02 =119,7 °С, соответственно. На основе этих результатов определим искомые параметры Е и В. Чтобы определить функцию Е (8,в0) = 0, необходимо представить в явном виде и подставить данные эксперимента (Тл, Т01), (Т2, Т02). Таким образом, записывается система двух уравнений 0(Т81, Т01) = 0 и G (TS2,Т02) = 0 с двумя неизвестными Е и В.
Решение системы уравнений удобно провести графическим способом, записав эти уравнения в виде В=/(Е). Обозначив эту функцию при значениях (ТбЛ, Т01) и (Т2, Т02) как В] и В2, соответственно, можно графически представить функцию В]-В2 от Е/Яг (см. рисунок 5). Тогда величина энергии активации Е определяется по точке пересечения функции с осью
Серия «Химическое машиностроение и инженерная экология» абсцисс, а затем вычисляется параметр В. С помощью рассмотренной процедуры для указанной системы получены следующие значения кинетических параметров Е=99,7 103 Дж/моль, В=6,81 010 К/с. Полученное значение энергии активации отличается от использованного для расчетов при определении температуры приблизительно на 1%.
Рисунок 5. Зависимость B1-B2 от E/Rg для определения энергии активации
Таким образом, предложенный метод температуры центральной точки позволяет определить кинетические параметры окисления угля на основе величин, характеризующих стационарное состояние, без достижения критических условий. Для определения искомых параметров требуется экспериментально определить значение температуры только в центральной точке образца, хотя бы при двух заданных значениях температуры на его поверхности, а также знать коэффициент температуропроводности. Предложенный подход значительно упрощает проведение эксперимента и сокращает объем требуемых экспериментальных данных по сравнению с другими термическими методами.
Работа выполнена при финансовой поддержке ФЦП & quot-Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России на 2014 -2020 годы& quot-, мероприятие 1. 2, Соглашение № 14. 575. 21. 0062.
Литература
1. Carras J.N., Young B.C. Self-heating of coal and related materials: models, application and test methods//Prog Energy Combust Sci. 1994- 20:1 — 15.
2. Haihui Wang, Bogdan Z. Dlugogorski, Eric M. Kennedy. Progress in Energy and Combustion Science 29 (2003), 487−513.
3. Parr S., Coons C. Carbon dioxide as an index of the critical oxidation temperature for coal in storage // Industrial and Engineering Chemistry, 17 (1925), 118 — 120.
4. Sujanti W., Zhang D., Chen X. D. Low-temperature oxidation of coal studied using wire-mesh reactors with both steady-state and transient methods // Combustion and Flame. 1999. V. 117. P. 646−651.
5. Chen X.D., Sidhu H., Nelson M. A linear relationship between dimensionless crossing-point-temperature and Frank-Kamenetskii reactivity parameter in self-heating test at infinite Biot number for slab geometry //Fire Safety Journal. 61 (2013), 138 — 143.
6. Франк-Каменецкий Д. А. Диффузия и теплопередача в химической кинетике. -М.: Наука, 1987.

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