
- •Раздел 2. Основы физического моделирования. 7
- •Раздел 3.2 .Модель идеального вытеснения. 26
- •Раздел 3.4. Диффузионная модель 29
- •Раздел 5. Синтез моделей технологических объектов на базе их гидродинамических моделей и уравнений химической кинетики. 64
- •Раздел 6. Методы численной реализации математических моделей. 68
- •1. Введение. Основные понятия систем
- •1.1.Очень большая система
- •Раздел 2. Основы физического моделирования. Виды испытаний.
- •2.1.Геометрическое и физическое подобие.
- •2.2. Способы получения критериальных уравнений.
- •1.2.2.Математическое описание.
- •1.2.3.Моделирующий алгоритм.
- •2. Общие принципы и этапы построения математических моделей систем.
- •2.1. Структурный анализ и структурный синтез сложных технологических систем
- •3. Модели структуры потоков для технологических объектов.
- •3.1 Модель идеального перемешивания.
- •3.3. Ячеечная модель аппарата
- •Раздел 3.4. Диффузионная модель
- •Стационарный метод определения критерия Пекле.
- •3.5.4.Гидродинамические модели многофазных потоков.
- •4. Кинетические модели для описания химических превращений.
- •4.1.Основные закономерности химической кинетики
- •4.2. Методы определения параметров кинетических моделей.
- •4.2.1.Определение констант скорости параллельных реакций:
- •Раздел. Кинетика гетерогенных процессов
- •2.2 Методы определения лимитирующей стадии гетерогенного процесса.
- •4.3.Определение кинетических констант сложных реакций методами нелинейного программирования.
- •Раздел 5. Синтез моделей технологических объектов на базе их гидродинамических моделей и уравнений химической кинетики.
- •Модель идеального вытеснения:
- •Где: сi – концентрация I-го компонента, wij- скорость расхода (или образования) I-го компонента в j-й реакции. Нj – тепловой эффект j-й реакции
- •Раздел 6. Методы численной реализации математических моделей.
- •6 Решение уравнений ячеечной модели.
- •6.1 Моделирование процесса нейтрализации сточных вод в каскаде реакторов идеального смешения.
- •6.1.1 Разработка модели процесса нейтрализации сточных вод в каскаде реакторов идеального смешения.
- •6.1.2 Моделирование процесса нейтрализации сточных вод в каскаде реакторов идеального смешения.
2.2 Методы определения лимитирующей стадии гетерогенного процесса.
1. По зависимости общей скорости процесса от скорости обтекания частиц твердого вещества потоком жидкости (сжимаемой – газа, или капельной –несжимаемой).
Если процесс протекает в диффузионной области, то его скорость, как и скорость массопередачи будут зависеть от скорости движения жидкости. Если режим движения жидкости при обтекании является турбулентным, то наблюдаемая скорость процесса будет пропорциональна скорости движения жидкости в степени 0.8. Если режим движения ламинарный, то скорость процесса будет зависеть от скорости потока в степени 0.5.
Т
аким
образом,
Для турбулентного режима:
Д
ля
ламинарного режима
скорость
процесса при отсутствии обтекания
скорость
обтекания твердой частицы жидкостью
2. По температурной зависимости наблюдаемой константы скорости.
Если так называемая эффективная энергия активации гетерогенного процесса составляет величину более 40 кДж/моль, то процесс лимитируется скоростью химической реакции и протекает в кинетической области.
Если эффективная энергия активации равна 10-20 кДж/ моль, то процесс проходит в диффузионной области. Если эффективная энергия активации лежит в пределах от 20 до 40 кДж/моль, то процесс происходит в смешанной, диффузионно-кинетической области.
Если исследовать температурную зависимость наблюдаемой константы скорости гетерогенного процесса в широких температурных пределах, то можно наблюдать переход из кинетической области в диффузионную, вследствие того, что при дальнейшем росте скорости реакции транспорт реагента к поверхности реакции становится меньше скорости его расхода и процесс переходит в диффузионную область.
Рис.2.4.Температурная зависимость наблюдаемой константы скорости гетерогенного процесса в координатах уравнения Аррениуса
Р
ис.2.4а.
Температурная зависимость константы
скорости в кинетической области
Р
ис.2.4б.
Температурная зависимость константы
скорости в диффузионной области.
Учет внутренней диффузии в гетерогенных каталитических
процессах.
Катализаторы для каталитических процессов изготавливают обычно в виде частиц различной формы чаще всего в виде таблеток цилиндрической или сферической формы. Таблетки обычно имеют развитую пористую структуру, что позволяет создавать катализаторы, имеющие внутреннюю поверхность пор, достигающих нескольких сотен м2 на грамм. При такой структуре частиц большая часть химических превращений протекает внутри частицы катализатора. Рассмотрим вывод уравнения математической модели каталитического процесса, протекающего внутри пористой частицы катализатора. Для анализа процессов и вывода уравнений математической модели чаще всего используют так называемую квазигомогенную модель частицы катализатора. В соответствии с этой моделью частицу рассматривают как некую изотропную среду, обладающую проницаемостью для реагентов и продуктов реакции и обладающую каталитическими свойствами к превращению исходного реагента в продукты реакции.
При этих допущениях можно считать, что транспорт вещества внутри пористой частицы катализатора подчинятся закону диффузии Фика, а скорость гетерогенной реакции на внутренней поверхности катализатора может быть отнесена к единице объема катализатора при введении так называемой удельной поверхности, т.е. поверхности, отнесенной к единице объема катализатора.
В этом случае, частица катализатора может характеризоваться эффективным коэффициентом диффузии, Deff и константой скоростью, отнесенной к единице объема катализатора kv.
С учетом введенных понятий, можно записать:
Где - коэффициент молекулярной диффузии, м2/с
- пористость частицы, отношение свободного объема пор Vпор к объему
катализатора Vкат, м3/ м3.
- коэффициент
извилистости поры, отношение средней
длины
транспортной поры к радиусу частицы катализатора , м/м
константа скорости гетерогенного процесса, отнесенная к единице
объема пористой частицы катализатора, 1/сек (для реакции первого
порядка)
khet - константа скорости гетерогенного процесса, м/сек (для реакции
первого порядка).
Sуд- удельная поверхность катализатора, отнесенная к единице объема
твердой фазы, м2/м3.
Рассмотрим сферическую частицу катализатора, окруженную сплошной средой с концентрацией реагента сv, кмоль/м3. Концентрацию на поверхности частицы будем считать равной концентрации вдали от поверхности вследствие интенсивного массообмена. На поверхности катализатора происходит реакция 1-го порядка, в результате которой исходное вещество превращается в продукт реакции по схеме:
Рис. 2.5. Схема диффузии реагента
внутри частицы катализатора
В качестве начальных условий обычно задают профиль концентрации по радиусу в начальный момент времени, а в качестве граничных условий задают значения концентраций и ее производных на границах частицы и внутри.
Р
ассмотрим
решение уравнения (2.19) с граничными
условиями (2.21) при стационарном режиме
работы катализатора, т.е. при условии,
что
Уравнение
для стационарных условий будет иметь
вид:
(2.22)
Разделим обе части
уравнения (2.22) на коэффициент диффузии
D
и введем новую переменную – безразмерный
радиус
, где
-радиус
частицы
У
равнение
примет вид:
(2.23)
Уравнение (2.23) имеет общее решение в виде:
(2.24)
Где А1 и А2 постоянные интегрирования определяемые из граничных условий, учитывающих массобмен на границе и изотропность свойств катализатора.
-
модуль Тиле, учитывающий соотношение
интенсивности потока внутрь частицы
за счет диффузии и расход вещества за
счет химической реакции. Это важнейший
массообменный параметр, учитывающей
доступность внутренней поверхности
катализатора для химической реакции.
Доступность внутренней поверхности катализатора для химической реакции принято характеризовать коэффициентом использования внутренней поверхности катализатораю
Коэффициент эффективности использования внутренней поверхности катализатора, или фактор эффективности - отношение количества вещества реагирующего в объеме частицы катализатора с учетом распределения концентрации по радиусу в соответствии с уравнением (2.30) к предельно возможному количеству при отсутствии внутри диффузионного торможения. Это количество будет равно произведению объема частицы на скорость реакции в единице объема при концентрации внутри частицы, равной концентрации в объеме:
Таким образом, можно определить размер частиц катализатора, обеспечивающий достаточно эффективное использование внутренней поверхности катализатора.
На рисунке показана зависимость фактора эффективности от величины модуля Тиле
Рисунок.3 Зависимость фактора эффективности от модуля Тиле.
Из рисунка видно, что наиболее приемлемые величины модуля Тиле, обеспечивающие эффективное использование внутренней поверхности катализатора находятся в интервале =0.11.0. Поскольку величина Модуля Тиле состоит из двух сомножителей, то очевидно, что для уменьшения этой величины при постоянной константе скорости, определяемой химической природой реагентов и катализатора, нужно либо уменьшать размер частиц, либо увеличивать коэффициент диффузии реагентов.
Использование квазигомогенной модели для моделирования гетерогенных процессов.