Сетевое издание
Современные проблемы науки и образования
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

ОБРАБОТКА ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ ТЕРМОГРАВИМЕТРИИ НА ОСНОВЕ ИНТЕГРАЛЬНЫХ ОЦЕНОК ИЗМЕНЕНИЯ СКОРОСТЕЙ РЕАКЦИИ С РОСТОМ ТЕМПЕРАТУРЫ

Лощилов С.А. 1 Коробейничев О.П. 2 Масленников Д.А. 1 Котова Ю.В. 1 Катаева Л.Ю. 1 Палецкий А.А. 2 Гончикжапов М.Б. 2
1 ГОУ ВПО «Нижегородский государственный технический университет им. Р.Е. Алексеева»
2 ФГБУН Институт химической кинетики и горения им. В.В. Воеводского Сибирского
В статье рассматриваются вопросы, связанные с обработкой результатов термогравиметрического эксперимента по нагреву хвои. В связи с тем, что при нагреве образца происходит ряд реакций, моделирование каждой из которых слишком сложно, поэтому целесообразно моделировать все процессы в виде одной или двух стадий. В данной работе стадии считаются последовательными, то есть продукт первой стадии является реагирующем веществом второй. Термокинетические постоянные этих стадий подбираются на основе решения обратной задачи, тогда как доля нереагирующего остатка на каждой стадии и границы температурного отрезка для анализа определяются экспертно. Предложенный в статье алгоритм основан на анализе изменения скорости расхода вещества с ростом температуры. Для корректной работы алгоритма необходимо, чтобы на заданном отрезке температур доминировала ключевая реакция, термокинетические параметры которой необходимо найти.
термокинетические постоянные
термогравиметрия
обратная задача
химическая кинетика
1. Катаева Л.Ю. Анализ динамических процессов аварийных ситуаций природного и техногенного характера: дис... д-ра физ.-мат. наук. – Нижний Новгород, 2009. – 328 с.
2. Леванов А. В. Введение в химическую кинетику: монография / А. В. Леванов, Э. Е. Антипенко. – М.: Высш. шк., 2006. – 52 с.
3. Фильков А.И. Кинетические исследования процесса пиролиза торфа // Горение твердого топлива: сб. ст. Всерос. конф. с междунар. уч. (Томск, 13-16 ноября 2012 г.). – Томск, 2012. С. 99.1-99.8.
4. Chen H. X., Liu N. A., Fan W. C. Two-Step Consecutive Reaction Model and Kinetic Parameters Relevant to the Decomposition of Chinese Forest Fuels // Journal of Applied Polymer Science, 2006, Sci 102. P. 571-576.
5. Huang X., Rein G. Smouldering Combustion of Soil Organic Matter: Inverse Modelling of the Thermal and Oxidative Degradation Kinetics // Proceedings of the ЕСМ 2013 (Lund, Sweden, 25– 28 June, 2013). P.1-6 (CD version).
6. Korobeinichev O.P. Combustion chemistry and decomposition kinetics of forest fuels // Procedia Engineering. 2013. Vol. 62. P. 182-193.
7. Rein G. Application of Genetic Algorithms and Thermogravimetry to Determine the Kinetics of Polyurethane Foam in Smoldering Combustion // Combustion and Flame. 2006. Vol. 146 P. 95-108.

Моделирование лесных пожаров является важной составляющей стратегии борьбы с ними. Современные подходы к моделированию лесных пожаров заключается в использовании математических постановок, основанных на физических законах. Химико-физические процессы являются одной из ключевых составляющих лесного пожара. Для их моделирования используется, как правило, закон Аррениуса, использование которого требует наличие информации о термокинетических постоянных этих процессов. Значения этих параметров оказывают существенное влияние на ход моделируемого пожара.

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

Другой сложностью нахождения термокинетических констант является существенное отличие убыли массы в экспериментах от модельных расчётов, проводимых в предположении существования одной реакции, скорость которой определяется законом Аррениуса. Именно эти факторы обуславливают разброс в термокинетических константах. Немаловажным также является и существенные различия этих значений в зависимости от вида материала.

Так, в работах [1, 4, 6] авторами приводятся константы скорости и значения энергий активации для лесных горючих материалов (ЛГМ) Сибирских и Китайских лесов – хвои сосны, ели, кедра, а также веток и коры сосны. В работах [3, 5] авторы приводят константы скорости и значения энергий активации для различных образцов торфа.

Существуют множество методов определения кинетических констант реакций пиролиза конденсировнных систем (ЛГМ, полимеров) на основе экспериментальных термогравиметрических данных. Одним из известных является метод, основанный на логарифмировании уравнения изменения массы в форме Аррениуса, что позволяет привести его к шаблону линейной регрессии [2]. В результате этого минимизируемая ошибка определяется не абсолютной, а относительной погрешностью, что вполне приемлемо, если исследуется только одна реакция и погрешность в определении её скорости мала. Однако при нагреве органических материалов происходит множество реакций, в результате чего по мере уменьшения остаточной массы образца побочные реакции вносят существенную погрешность. В работе [3] на основе модифицированного закона Герца – Кнудсена авторы исследуют характеристики пиролиза на основе экспериментальных данных по потере массы от времени, полученных при различных значениях начальной температуры в изотермических условиях, и сопоставляют их с термогравиметрическими данными, полученными в неизотермических условиях. Для нахождения констант скорости двухстадийных последовательно-параллельных реакций окислительного пиролиза торфа и полимеров при тлеющем их горении в работах [5, 7] применялся генетический алгоритм, представляющий собой эвристический метод поиска, имитирующий принципы биологической адаптации на основе дарвинистский теории естественного отбора. В работе [5] определено 16 параметров. Определение кинетических параметров было выполнено с помощью минимизации параметра Φ, содержащего массу и скорость ее изменения (данные ТГ и ДТГ)

, (1)

где γ = 0,5. В работах [4, 6] кинетический анализ пиролиза лесных горючих материалов в окислительной среде был проведен с использованием механизма для многокомпонентных веществ. Процесс потери массы лесных горючих материалов моделируется четырьмя параллельными реакциями, в соответствии с разложением трех основных компонентов биомассы (гемицеллюлоза, целлюлоза и лигнин) и окисления углеродного остатка. Определение кинетических параметров для каждой реакции было выполнено с помощью минимизации SDTG , коэффициента экспериментальных данных

(2)

где m, T – масса, температура, индексы “exp” and “simu” означают экспериментальное и рассчитанное значения, число экспериментальных кривых и число точек для каждого эксперимента. Минимизация SDTG выполнена нелинейным алгоритмом подгонки с помощью программы Matlab. Для оценки соответствия экспериментальных и моделируемых данных используется параметр Dev

(3)

Величина этого параметра была равна 2 %.

В данной работе представлен анализ экспериментальных термогравиметрических данных пиролиза иголок сосны, полученных в ИХКГ СО РАН, с помощью алгоритма, разработанного в НГТУ им. Р.Е. Алексеева.

Предполагается, что существуют две последовательные реакции, причём продукт первой является реагентом для второй.

, (4)

, (5)

. (6)

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

(7)

В качестве начальных условий используются соотношения.

(8)

В уравнениях (4–8) приняты обозначения: i – порядковый номер реагирующего вещества; mi(t) – масса i-й компоненты сухого органического вещества в процессе реакции на момент t; kj, Ej – предэкспонента и энергия активации j-й стадии реакции пиролиза, подбираемые с помощью решения обратной задачи; αj – отношение массы продукта j-й стадии реакции к массе исходного вещества на данной стадии. Температура T(t) изменяется по заданному закону. Во всех рассмотренных в данной работе экспериментах эта зависимость линейная.

Для решения обратной задачи в данной работе авторами предлагается алгоритм, вычисляющий термокинетические параметры стадий на основе экспериментальных данных ТГА (термогравианалитической кривой), полученных в результате нагрева образца с заданной скоростью (10 К/с).

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

Далее алгоритм приводит очищенные от повторов экспериментальные данные к заданной равномерной сетке температур с использованием сплайн-интерполяции.

На следующем этапе вычисляется скорость относительной потери массы с помощью соотношения:

, (9)

где Mi – скорость относительной убыли массы i-го компонента, ti,– время на соответствующем j-м шаге сетки (tj+1), t0 – начальный момент времени в анализе данной стадии, αi – величина условно нереагирующего остатка после i-й стадии. Соотношение (9) учитывает также условно нереагирующий остаток.

Перед запуском итерационного процесса определяются достаточно грубые верхняя и нижняя оценки энергии активации Emin и Emax. На каждой итерации происходит следующий процесс:

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

2) Выполняется оценка скорости расхода вещества на основе предполагаемого значения энергии активации

, (10)

3) Вычисляется предэкспонента, при которой убыль массы по эксперименту будет совпадать с модельным расчётом при энергии активации Etest.

(11)

где jmin – индекс, соответствующий минимальной температуре анализируемого отрезка, а jmax – максимальной.

4) Вычисляется расчётная относительная скорость потери массы

(12)

5) Оценивается квадрат ошибки (с учётом знака) отклонения экспериментальных данных от расчётных значений по соотношению (13). Положительный знак указывает на то, что расчётная величина скорости относительной потери массы больше экспериментальной.

(13)

6) Определяется знак интеграла

(14)

7) Если интеграл (14) положителен, значит, расчётная скорость относительной потери массы принимает большие значения, чем экспериментальная, при более низкой температуре, но меньшие – при более высокой. Это в свою очередь означает, что рост скорости реакции с температурой недооценён, и, следовательно, предполагаемая энергия активации является нижней оценкой Emin. В противном случае она является верхней оценкой Emax.

8) Если заданная точность не достигнута, то выполняется следующая итерация, иначе полученные значения Etest и Ktest считаются окончательным решением.

На основе описанного алгоритма был разработан программный комплекс. С его помощью были выполнены расчёты термокинетических постоянных по данным ТГА при нагреве хвоинки в окислительной среде (20 % O2). На рис. 1 показано сопоставление экспериментальных данных и решения прямой задачи при найденных константах для данного эксперимента.

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

При нагреве образца в окислительной среде (20 % O2) было выявлено две стадии: первая стадия доминирует при температуре от 400 до 600К, а вторая от 660 до 730К. При этом термокинетические параметры первой стадии: k1=1.082·106 сек-1, Е1=96696 Дж/моль, второй стадии – k2=9.489·107 сек-1; Е=140654 Дж/моль. Условно нереагирующие остатки α1=0.37, α2=0.216. При температурах ниже выбранного диапазона могут происходить другие реакции, которые не влияют на результат, что видно из рис. 1.

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

Работа выполнена при финансовой поддержке 13-03-91164-ГФЕН_а «Экспериментальное исследование кинетики и механизма термического разложения лесных горючих материалов и процессов распространения пламени по их слою».

Рецензенты:

Петрухин Н.С., д.ф.-м.н., ординарный профессор Национального исследовательского университета «Высшая школа экономики», г. Нижний Новгород.

Карпухин В.Б., д.ф.-м.н., профессор кафедры «Высшая и прикладная математика» Российской открытой академии транспорта Московского государственного университета путей сообщения, г. Москва.


Библиографическая ссылка

Лощилов С.А., Коробейничев О.П., Масленников Д.А., Котова Ю.В., Катаева Л.Ю., Палецкий А.А., Гончикжапов М.Б. ОБРАБОТКА ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ ТЕРМОГРАВИМЕТРИИ НА ОСНОВЕ ИНТЕГРАЛЬНЫХ ОЦЕНОК ИЗМЕНЕНИЯ СКОРОСТЕЙ РЕАКЦИИ С РОСТОМ ТЕМПЕРАТУРЫ // Современные проблемы науки и образования. – 2013. – № 6. ;
URL: https://science-education.ru/ru/article/view?id=10792 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674