Зависимость константы скорости от температуры обычно описывают с помощью уравнения Аррениуса (см. обсуждение и варианты этого уравнения в разделе "Теория"). Здесь мы обсудим, как задавать параметры кинетической схемы для моделирования температурных зависимостей констант скорости процессов.
Возьмём простой процесс
2B → C // k
Пусть k =
Здесь A – предэкспонента, Eₐ – энергия активации по Аррениусу, T – температура в градусах Кельвина.
В соответствием с C-синтаксисом в схеме мы это запишем так: A*exp(-Eₐ/8.3/T).
Нам понадобится объявить следующие параметры схемы:
<schema id="sample_arrhenius"> <parameter name="A" value=""/> <parameter name="Eₐ" value=""/> <parameter name="T"/> <parameter name="k" value="A*exp(-Eₐ/8.3/T)"/> </schema>
Обобщённое выражение для константы скорости реакции мы выносим на уровень схемы процесса. Поскольку A и Eₐ в этом примере постоянны, их значения мы задаём здесь же.
Рассмотрим кинетики нескольких вариантов, отличающихся температурой процесса (100℃, 150℃ и 200℃ переводим в градусы Кельвина):
<instance schema="sample_arrhenius" id="var₁"> <parameter name="T" value="373"/> </instance> <instance schema="sample_arrhenius" id="var₂"> <parameter name="T" value="423"/> </instance> <instance schema="sample_arrhenius" id="var₃"> <parameter name="T" value="473"/> </instance>
Расчёты запускаем командой compute_kinetics:
<compute_kinetics instance="var₁"/>