Математическое различие заключается в расположении и единицах измерения источникового члена реакции. Для гомогенной реакции скорость непосредственно входит в дифференциальное уравнение сохранения вещества как объёмный источник, (R_{Vi}), с единицами измерения (\mathrm{mol,m^{-3},s^{-1}}). Для гетерогенной реакции, протекающей на поверхности или границе раздела фаз, скорость задаётся как граничный поток, (R_{Si}), с единицами измерения (\mathrm{mol,m^{-2},s^{-1}}), либо преобразуется в эквивалентный объёмный источник при использовании усреднённой модели реактора.

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

Как баланс веществ учитывает реакцию

Общее уравнение сохранения

Для вещества (i) молярный баланс можно записать в виде

[ \frac{\partial C_i}{\partial t}

-\nabla \cdot \mathbf{N}i + R{Vi} ]

где (C_i) — концентрация вещества, (\mathbf{N}i) — его молярный поток, а (R{Vi}) — суммарная скорость образования в расчёте на единицу объёма.

Дивергентный член учитывает перенос в дифференциальный контрольный объём или из него. Член реакции учитывает химическое образование или потребление вещества внутри этого объёма.

Скорости гомогенных реакций

Гомогенная реакция протекает во всей одной жидкой фазе, например в жидком растворе или газовой смеси. Её скорость непосредственно входит в основное дифференциальное уравнение:

[ \frac{\partial C_i}{\partial t}

-\nabla \cdot \mathbf{N}i + R{Vi} ]

Для механизма, включающего реакции (j=1,\ldots,N_r),

[ R_{Vi}

\sum_{j=1}^{N_r} \nu_{ij} r_{Vj} ]

где (\nu_{ij}) — стехиометрический коэффициент вещества (i) в реакции (j), а (r_{Vj}) — скорость реакции в расчёте на единицу объёма жидкости.

Типичные единицы измерения:

[ r_{Vj}, R_{Vi} ;:; \mathrm{mol,m^{-3},s^{-1}} ]

Отрицательное значение означает суммарное потребление, а положительное — суммарное образование.

Как гетерогенные реакции входят в баланс

Поверхностная реакция как граничное условие

Если реакция протекает только на поверхности катализатора или на границе раздела фаз, она не действует во всём прилегающем объёме жидкости. Вместо этого она изменяет поток вещества на данной поверхности.

На границе раздела (\Gamma) граничное условие может быть выражено как

[ -\mathbf{n}\cdot\mathbf{N}_i

R_{Si} ]

где (\mathbf{n}) — внешняя единичная нормаль, а (R_{Si}) — скорость образования или потребления вещества на поверхности.

Для нескольких поверхностных реакций

[ R_{Si}

\sum_{j=1}^{N_s} \nu_{ij} r_{Sj} ]

с типичными единицами измерения

[ r_{Sj}, R_{Si} ;:; \mathrm{mol,m^{-2},s^{-1}} ]

Таким образом, реакция появляется в граничном условии, а не как объёмный источник в уравнении для основной жидкости.

Скачок потока на границе раздела

Для реакции на границе между двумя фазами более общим условием является разрыв потока:

[ \mathbf{n}\cdot \left( \mathbf{N}{i,2}-\mathbf{N}{i,1} \right)

R_{Si} ]

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

Такая форма полезна для пилотных газожидкостных, газотвёрдых или жидкость-твёрдое реакторов, в которых перенос в обеих фазах описывается отдельными уравнениями.

Преобразование скоростей реакций на катализаторе для моделей реакторов

Скорость в расчёте на единицу массы катализатора

Кинетические параметры гетерогенного катализа часто приводятся в расчёте на единицу массы катализатора:

[ r_{i,\mathrm{cat}} ;:; \mathrm{mol,kg_{cat}^{-1},s^{-1}} ]

Для модели реактора с неподвижным слоем может потребоваться источниковый член в расчёте на единицу объёма реактора. Преобразование выполняется по формуле

[ R_{Vi,\mathrm{eff}}

\rho_{b,\mathrm{cat}},r_{i,\mathrm{cat}} ]

где (\rho_{b,\mathrm{cat}}) — объёмная загрузка катализатора в расчёте на единицу объёма реактора.

Получаемый эффективный источник имеет единицы измерения

[ R_{Vi,\mathrm{eff}} ;:; \mathrm{mol,m^{-3}_{reactor},s^{-1}} ]

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

Скорость в расчёте на единицу площади поверхности катализатора

Если скорость выражена в расчёте на единицу площади поверхности катализатора,

[ r_{i,\mathrm{surf}} ;:; \mathrm{mol,m^{-2}_{cat},s^{-1}} ]

то соответствующий источник в расчёте на объём реактора равен

[ R_{Vi,\mathrm{eff}}

a_s r_{i,\mathrm{surf}} ]

где (a_s) — площадь поверхности катализатора в расчёте на единицу объёма реактора, с единицами измерения (\mathrm{m^2_{cat},m^{-3}_{reactor}}).

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

Варианты моделирования на пилотном масштабе

Явная гетерогенная модель

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

Например, баланс для жидкой фазы может быть схематически представлен как

[ \frac{\partial C_{i,f}}{\partial t}

-\nabla\cdot\mathbf{N}{i,f} + a_s k{i,m}(C_{i,f}-C_{i,s}) ]

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

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

Псевдогомогенная модель

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

[ \frac{\partial C_i}{\partial t}

-\nabla\cdot\mathbf{N}i + a_s r{Si} ]

или, для основы в виде массы катализатора,

[ \frac{\partial C_i}{\partial t}

-\nabla\cdot\mathbf{N}i + \rho{b,\mathrm{cat}}r_{i,\mathrm{cat}} ]

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

Выбор представления

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

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

Понимание компромиссов

Не смешивайте единицы площади и объёма

Распространённая ошибка — подставлять (\mathrm{mol,m^{-2},s^{-1}}) непосредственно в уравнение, источниковый член которого требует (\mathrm{mol,m^{-3},s^{-1}}). Сначала скорость необходимо умножить на площадь границы раздела в расчёте на объём реактора либо применить как граничный поток.

Та же проблема возникает при использовании скорости в расчёте на килограмм катализатора. Перед включением в объёмный баланс вещества её необходимо умножить на массу катализатора в расчёте на объём реактора.

Не путайте гетерогенную химию с гетерогенным моделированием

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

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

Учитывайте внутренние поверхности катализатора

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

Представление внутренней каталитической активности как внешнего граничного условия без учёта доступной внутренней площади может существенно исказить скорость реакции.

Сохраняйте стехиометрическую согласованность

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

[ R_i=\sum_j \nu_{ij}r_j ]

Изменение основы нормирования изменяет единицы измерения и масштабирование, но не исходную стехиометрию.

Как применить это к вашему проекту

Корректная реализация зависит от физического расположения реакции и основы, на которой представлены доступные кинетические данные.

Если основное внимание уделяется проектированию гомогенного реактора: Используйте (R_{Vi}) непосредственно в уравнении сохранения вещества с единицами измерения (\mathrm{mol,m^{-3},s^{-1}}).
Если основное внимание уделяется разрешённой поверхностной кинетике: Применяйте (R_{Si}) как межфазный поток или граничное условие с единицами измерения (\mathrm{mol,m^{-2},s^{-1}}).
Если основное внимание уделяется масштабированию реактора с неподвижным слоем: Преобразуйте скорости на основе массы катализатора или площади поверхности в эффективный объёмный источник с использованием (\rho_{b,\mathrm{cat}}) или (a_s).
Если основное внимание уделяется ограничениям катализатора или переноса: Используйте отдельные уравнения для жидкости и твёрдой фазы, чтобы учитывать перенос через плёнку, температурные различия и внутрин частичную диффузию.
Если основное внимание уделяется быстрому моделированию пилотной установки: Используйте псевдогомогенную модель с тщательно определённой эффективной объёмной скоростью реакции.

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

Сводная таблица:

Тип реакции или модели	Математическое представление	Типичные единицы измерения	Применение на пилотном масштабе
Гомогенная реакция	Объёмный источник, `Rᵥᵢ = Σνᵢⱼrᵥⱼ`	mol·m⁻³·s⁻¹	Реакции в объёме жидкости или газа
Разрешённая гетерогенная реакция	Поверхностный поток или скачок межфазного потока, `Rₛᵢ`	mol·m⁻²·s⁻¹	Моделирование поверхностей и границ раздела фаз
Основа в виде массы катализатора	`Rᵥᵢ,eff = ρᵦ,cat rᵢ,cat`	mol·m⁻³ reactor·s⁻¹	Модели реакторов с неподвижным слоем
Основа в виде поверхности катализатора	`Rᵥᵢ,eff = aₛ rᵢ,surf`	mol·m⁻³ reactor·s⁻¹	Пористые катализаторы и усреднённые модели
Псевдогомогенная модель	Эффективный объёмный источник в балансе	mol·m⁻³·s⁻¹	Ускоренное моделирование и масштабирование реакторов
Явная многофазная модель	Отдельные балансы фаз с межфазным переносом и поверхностными условиями	Зависит от основы	Анализ ограничений, связанных с плёнкой, порами и тепловыми эффектами

Разработайте подходящую модель реактора пилотного масштаба

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

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