Макрокинетика гетерогенных каталитических реакций. 
Устойчивость экзотермических реакций. 
Выбор типа реактора и условий реализации промышленного процесса

Федеральное агентство по образованию Государственное образовательное учреждение высшего профессионального образования

" ТОМСКИЙ ПОЛИТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ"

Химико-технологический факультет Процессы и аппараты химических производств Кафедра химической технологии топлива

" Макрокинетика гетерогенных каталитических реакций. Устойчивость экзотермических реакций. Выбор типа реактора и условий реализации промышленного процесса"

Реферат по дисциплине «Дополнительные главы физической химии» .

Исполнитель:

маг. гр.5М550 Афанасьева Ю.И.

Руководитель:

АссистентЧеканцев Н.В.

Томск — 2010

• Макрокинетика гетерогенных каталитических реакций
• Массоперенос и теплопередача в порах
• Внешняя массои теплопередача
• Устойчивость экзотермических реакций
• Выбор типа реакции и условий реализации промышленного процесса. Оптимальный температурный профиль
• Оптимизация стадийных cхем
• Список используемой литературы

Макрокинетика гетерогенных каталитических реакций

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

Массоперенос и теплопередача в порах

Важное значение имеет перенос вещества и тепла внутри пористой частицы катализатора.

Перенос вещества в порах осуществляется исключительно путем молекулярной диффузии.

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

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

Средняя скорость движения молекул зависит от температуры газа и в равновесном газе определяется законом Максвелла:

где k — константа Больцмана; Т - абсолютная температура; т - масса молекулы. Коэффициент молекулярной диффузии в газах может быть определен:

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

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

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

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

Эту область (л «d_п) называют кнудсеновской, а величину D_кн — кнудсеновским коэффициентом диффузии.

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

При исследовании макрокинетики химических реакций в пористом зерне нерационально рассматривать процесс в отдельной поре. Поры реальной частицы катализатора неодинаковы по размеру и образуют запутанную сеть. Поэтому пористое зерно рационально рассматривать как квазигомогенную среду, характеризуя скорость диффузии реагентов эффективным коэффициентом диффузии D, а скорость химической реакции — эффективной кинетической функцией r (С, Т).

Эффективный коэффициент диффузии D<D_м коэффициента молекулярной диффузии. Это объясняется, тем, что в пористом зерне доступно для диффузии не все сечение, а лишь часть, занятая порами, равная доле свободного объема и удлинением пути диффузии в извилистых порах. В общем случае эффективный коэффициент диффузии можно представить в виде:

D = аD_м

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

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

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

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

В пористом катализаторе перенос тепла осуществляется как с помощью молекулярного теплопереноса в порах, так и за счет теплопроводности самого катализатора.

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

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

Внешняя массои теплопередача

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

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

где и_x и_у, и_z - компоненты вектора скорости и, являющиеся функциями координат X, Y, Z.

Величину и направление скорости в каждой точке определяют решением уравнений гидродинамики. В правой части уравнения (7) оставлена производная только по координате X, нормальной к поверхности, так как по всем другим направлениям перенос вещества молекулярной диффузией пренебрежимо мал. Граничные условия для уравнения (7) определяются тем, что диффузионный поток на твердую поверхность катализатора равен скорости химической реакции, а на достаточном удалении от поверхности концентрация равна С_?.

При больших скоростях движения практически весь перепад скорости сосредоточен в тонком гидродинамическом пограничном слое толщиной д₀, а перепад концентрации — в диффузионном пограничном слое толщиной д. Отношение д₀/д тем выше, чем больше отношение кинематической вязкости вещества х к коэффициенту молекулярной диффузии D_м. В жидкостях, где х/D_м>1. В этом случае при решении уравнения (7) используют простые выражения для скорости потока вблизи твердой поверхности, что позволяет найти аналитическое решение при протекании быстрой реакции или реакции первого порядка на поверхности частиц простой геометрической формы (пластина или шар). В газах толщины диффузионного и гидродинамического пограничных слоев одного порядка и уравнение (7) может быть решено только приближенными методами. Аналогичным способом — путем решения уравнения конвективной теплопроводности — определяют и тепловой поток на твердую поверхность

При обтекании тел сложной геометрической формы уравнения (7) становятся практически неразрешимыми. В этом случае эффективная толщина диффузионного слоя или коэффициент массопередачи в = D_M/д определяются экспериментально. Опыты по измерению коэффициента массопередачи проводят в условиях, когда гетерогенная реакция протекает практически мгновенно и скорость исследуемого процесса лимитируется исключительно переносом вещества к активной поверхности.

При исследовании теплообмена между движущимся веществом и активной поверхностью, тепловой поток q_т может быть выражен формулой

где б — коэффициент теплопередачи; Т_? - Т₀ - перепад температур между веществом и твёрдой поверхностью.

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

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

Устойчивость экзотермических реакций

Система уравнений материального и теплового балансов экзотермической реакции на пористом зерне катализатора:

С учетом сопротивления переносу вещества и тепла на внешней поверхности зерна Г:

(11)

где r — скорость реакции в единице объема пористого катализатора; D_i и ч — эффективные коэффициенты диффузии и теплопроводности в пористом зерне; х_i - стехиометрический коэффициент; h - теплота реакции; ² — оператор Лапласа; C_i0 и Т₀ — концентрации реагентов и температура на внешней поверхности зерна; С_i?, Т_? — концентрация реагентов и температура в ядре потока; в_i, б — коэффициенты массо — и теплопередачи из ядра потока к внешней поверхности зерна; п - направление внешней нормали к поверхности зерна катализатора.

Заметим, что в_i = D_iM/д и D_{i = aDiM (где Diм — коэффициент молекулярной диффузии 1-го вещества, д — толщина диффузионного пограничного слоя и а — коэффициент, зависящий от структуры пористого катализатора). Величины д и a можно с достаточной степенью точности считать одинаковыми для всех веществ, участвующих в реакции. Система уравнений (90) может быть сведена к единственному уравнению для концентрации одного из реагирующих веществ, которое принимают за ключевое. Введем с этой целью вспомогательную величину}

(23)

Из уравнений (19) и граничных условий (21) следует, что величина w_ij определяется решением задачи

(24)

Очевидно, решение этой задачи имеет вид w_ij = const = w_ij^? Из формулы (23) теперь следует, что концентрации всех веществ, участвующих в реакции, связаны между собой линейными соотношениями

(25)

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

(26), (27), (28), (29)

где х - безразмерная координата (с масштабом l) направленная вдоль внешней нормали к поверхности, и Bi_м = вl/D, диффузионное и тепловое числа Био.

Выбор безразмерных переменных в уравнениях (26), (27) несколько необычен. За масштаб концентрации ключевого вещества внутри пористого зерна мы приняли концентрацию ключевого вещества на внешней поверхности зерна, так что с = С/С₀. Аналогично, безразмерная температура внутри зерна определяется как и = (Т - То) /Т₀^* (где Т₀^* - масштаб внутренней температуры, зависящий от температуры внешней поверхности зерна). В соответствии с этим параметры Ш₀ = l [r (С₀, Т₀) /DС₀] ^1/г, И₀=hC_oD/чT₀^* (параметр максимального разогрева), а также безразмерная скорость реакции f (с, и) = r (С, Т) /r (С₀, Т ₀) зависят от значений переменных на внешней поверхности зерна. Поскольку эти значения неизвестны заранее, величины Ш₀ и И₀, как и другие параметры, которые могут входить в выражение для f (с, и), не являются характерными параметрами рассматриваемого процесса. Преимущество такого выбора безразмерных переменных заключается, однако, в том, что в дальнейшем можно воспользоваться результатами исследования хорошо изученной задачи расчета экзотермической реакции на пористом зерне катализатора при фиксированных значениях концентрации и температуры на внешней поверхности зерна. Нетрудно видеть, что Ш₀ и И₀ являются характерными параметрами этой задачи. В граничных условиях, (29) участвуют безразмерная концентрация ключевого вещества на внешней поверхности частицы c_П ⁼C₀/C_? и безразмерная температура внешней поверхности и_П= (T₀-T_?) /T*, масштабы которых связаны уже только со значениями переменных в ядре потока.

Пользуясь уравнениями (26), (27) и граничными условиями (28), (29), можно получить оценку условий, при которых существует заметный перепад концентрации и температуры между поверхностью катализатора и внешней средой. Внешнее сопротивление массопередаче начинает сказываться только, когда реакция локализуется в тонком слое, толщина которого сравнима с толщиной диффузионного пограничного слоя д. Действительно, поскольку величина Ш₀ является мерой проникновения реакции в глубь пористого катализатора, так что dc/dx₀, из граничного условия (28) следует:

(30)

Так как всегда, внешний перепад концентрации сравним с внутренним только при .

Для оценки перепада температуры воспользуемся линейной связью между концентрацией и температурой внутри пористой частицы, которая следует из уравнений (25), (26)

(31)

откуда

(32)

Принимая во внимание, что для газовой среды (где ч_м — коэффициент молекулярной теплопроводности газа) и используя граничное условие (29), находим:

(33)

Так как перенос тепла в пористом катализаторе идет в основном по твердой фазе, то. Для сильно экзотермических реакций величина И₀ может принимать значения порядка нескольких единиц. Поэтому величина и_П, характеризующая влияние внешнего сопротивления теплоотводу, может быть существенной () при сравнительно небольших значениях Ш₀, когда внешнедиффузионное сопротивление еще не оказывает заметного влияния на процесс.

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

(34)

Используя значения фактора эффективности з, найденные путем решения задачи с фиксированными С₀, T₀, можно вычислить правую часть уравнения (34) ц (и_П); после этого значения внешнего разогрева и_П определяются как точки пересечения графика функции ц (и_П) с прямыми с наклоном Вi_Т, проведенными через начало координат.

Поскольку величины Ш₀ и И₀ зависят от поверхностных значений концентрации и температуры и потому не являются характерными параметрами процесса, они должны быть выражены через параметры

,

отнесенные к условиям в ядре потока. Установлено, что при не слишком больших внешних разогревах и_П С₀? С_? и, следовательно, (Т₀*/Т*) И₀? И Ш₀ = Шу (и_П) (где у (и_п) — сильно возрастающая функция и_П).

В случае, когда зависимость з (Шо) неоднозначна, график зависимости правой части ц (и_П) уравнения (34) от и_П имеет форму, показанную на рис, ШЛО. При увеличении и_П функция ц (и_П) резко возрастает за счет сильной зависимости Ш₀ от и_П. При дальнейшем увеличении и_П. когда начинает сказываться внешне-диффузионное торможение процесса, функция ц (и_П) проходит через максимум и резко спадает до нуля. Решения уравнения (111.115) определяются как точки пересечения кривой ц (и_П) и прямой с наклоном Вi_Т. Из рис. 1 видно, что уравнение (34) может иметь до пяти решений, причем при значениях числа Вi_Т меньше некоторых критических, решения, соответствующие малым внешним разогревай, исчезают. При увеличении температуры потока, а следовательно, параметра Ш кривая ц (и_П) сдвигается вверх, и при значениях Ш, больших некоторых критических, также исчезают низкотемпературные решения. Участку кривой ц (и_П) слева от точки 1 соответствует внутрикинетический, а между точками 2 и 3 — внутридиффузиошшй режим протекания реакции. В задаче с фиксированными значениями поверхностных концентраций и температуры внутрикинетический режим существует левее точки 4, а внутридиффузионный — правее точки 5. Наличие сопротивления теплоотводу на внешней поверхности зерна приводит к расширению области промежуточных неустойчивых состояний. Кроме того, при наличии внешнего сопротивления теплоотводу множественные режимы возникают и в том случае, когда зависимость з (Ш₀) однозначна.

Рис. 1. Графическое решение уравнения (34)

Промежуточных неустойчивых состояний. Кроме того, при наличии внешнего сопротивления теплоотводу множественные режимы возникают и в том случае, когда зависимость т) (Ч'о) однозначна.

Реакция первого порядка. Рассмотрим более подробно необратимую экзотермическую реакцию первого порядка на пористом катализаторе, имеющем форму пластины толщиной 2l с «запечатанными» боковыми гранями. Это простейший пример, на котором видны все характерные особенности рассматриваемой задачи. Пусть. Если принять за масштаб температуры Т величину, то

(35)

Величины Ш₀ иИ₀ связаны с характерными параметрами задачи Ш, И=чR следующими соотношениями, вытекающими из формул (28), (29)

(36), (37)

При записи параметра в за масштаб внешней температуры принята величина . При не слишком больших и_П (и_П1) величины и Т* практически совпадают и, так как. Заметим, что из формулы (37) следует существование максимально возможного разогрева и_{п. макс}=И_п=mИ.

Используя соотношение (31) между концентрацией и температурой внутри катализатора, сведем систему (26), (27) к единственному уравнению относительно безразмерной температуры и

(38)

где x — безразмерная координата, принимающая значения x=0 в центре пластины и x=±1 — на ее поверхности.

Переходя к интегрированию уравнения (38), будем считать величину и_П, а следовательно, и значения параметров Ш₀ и И₀ известными. Далее можно найти величину и_П, подставляя полученное решение в граничное условие (20). Таким образом, будем решать уравнение (38) с граничными условиями:

(39)

Метод решения задачи (38), (39) выражается формулами (88) и (89):

(40), (41)

Метод решения задачи (38), (39) изложен в предыдущем разделе, ее решение выражается формулами (88) и (89):

(42)

где

(43)

безразмерный параметр теплопередачи.

Исследуем сначала область, в которой можно пренебречь внешнедиффузионным торможением. Тогда и задача сводится к совместному решению уравнений (42) и (41) относительно неизвестных величин и_l и и_П.

В двух предельных случаях из этой системы уравнений можно исключить величину и_l. В кинетической области, считая, получим из формулы (41):

(44)

Подставляя и_l, найденное из уравнения (44), в (42) получим в том же приближении:

(45)

Во внутридиффузионной области можно положить в уравнении (42). При этом получаем:

(46)

Правая часть формулы (45) имеет минимум, равный при eШИ₀, а минимум выражения (46) достигается при и_П = 2. Из выражения (46) следует, что во внутридиффузионном режиме правая часть уравнения (42) не зависит от параметра Ш и, следовательно, значение и_П в этой области определяется только параметрами М и И₀.

Зависимость и_П (М) при известных значениях Ш и И₀ может быть построена с помощью графика функции ц (и_l). Действительно, задавая значение и_l, можно по графику определить и, подставляя полученные значения и_п и и_l в формулу (42), построить зависимость и_П (М).

Рис 2. Зависимость внешнего разогрева и_П от параметра М: И=5; Ш=0,37 (а); 0,30 (б); 0, 20 (в).

Рис. 3. Зависимость внешнего разогрева и_п от параметра М: И=4 (а-в); Ш=0,47 (а), 0,45 (б), 0,40 (в); И=3 (г, д); Ш=0,50 (г), 0,45 (д); И=2; Ш0,53 (е).

На рис. 2 и 3 приведен ход кривых и_П (М) в области и_П ~ 1 для нескольких значений И и Ш. Пересекая полученные кривые прямыми М = const, можно найти соответствующие этим значениям параметра теплопередачи безразмерные перепады температуры между поверхностью катализатора и внешним потоком. Обратимся к изучению зависимости и_П (М) при различных комбинациях параметров И и Ш. Представленное на рис. 2 семейство кривых соответствует И = 5 и нескольким значениям Ш. В этом случае, как уже указывалось, зависимость и_l от может быть неоднозначной, поэтому кривая М (и_П) тоже может быть неоднозначной в некоторой области изменения и_П. Границы области неоднозначности кривой М (и_П), отмеченные точками 4 и 5 на рис. 2, в которых |dM/dи_п|=?, соответствуют точкам максимума и минимума функции ц (и_l). Указанная неоднозначность М (и_П) при И >4,5 проявляется в том случае, когда Ш < Ш_макс, где Ш_макс значение ц (и_l) в максимуме. При этом точкам 4 на рис. 2 соответствуют и_п=2 ln (Ш_макс/Ш), а точкам 5 и_П =2 ln (Ш_мин/Ш) (где Ш_мин — значение ц (и_l) в минимуме). Если Ш < Ш_мин, то на графике М (и_П) присутствуют точки типа 5 и 4; при Ш_мин < Ш < Ш_макс точки типа 5 отсутствуют, и при Ш > Ш_макс зависимость М (и_П) является однозначной.

Кривые М (и_П)) соответствующие разным Ш, могут иметь в области и_П ~ 1 минимумы (точки 1 и 3) и максимумы (точки 2) и выходят во внутридиффузионной области на единую асимптотическую кривую, определяемую формулой (46). Наличие на кривой М (и_П) экстремумов означает существование в рассматриваемой задаче множественных режимов, причем значения и_П, при которых dM/dи_П =0 (точки 1-3 на рис. 2 и 3), разделяют чередующиеся области устойчивых и неустойчивых режимов. Заметим, что точки 4 и 5 (соответствующие значениям Ш₀, при которых в задаче с заданными значениями поверхностных концентраций С₀ и температуры Т₀ происходил перескок между режимами) попадают в область неустойчивых режимов, разделяющих внутрикинетическую. и внутри-диффузионную области. Таким образом, наличие сопротивления теплоотводу на внешней поверхности катализатора приводит к расширению области неустойчивых режимов. Это проявляется также в возникновении неустойчивой области при И < 4,5, когда задача без внешнего сопротивления вообще не имеет неустойчивых режимов. Скачкообразные переходы во внутренних режимах могут наблюдаться, как это видно из рис. 3, до И = 2.

При значениях параметра М, меньших абсолютного минимума функции М (и_П) (этот минимум может соответствовать точкам 1 и 5, см. рис. 2,3), в рассматриваемой области режимов с небольшими внешними разогревами не существует и процесс переходит во внешнедиффузионный режим (и_П ~ И_п). Во внешнедиффузионном режиме существенным становится перепад концентрации у поверхности катализатора, и поэтому для получения зависимости М (и_П) в этой области значений и_П необходимо учитывать связь И₀ с и_П, определяемую формулой (37). Поскольку при этом реакция локализована в узком слое близ поверхности катализатора из-за большого внутридиффузионного торможения, ход кривой М (и_П) во внешнедиффузионной области можно получить из формулы (46) с учетом зависимости (37). При больших и_П можно также заменить приближенную кинетическую функцию на ехр и_п/2 (1 + bи_П) (где b == RТ_?/Е) и учесть различие между величинами T₀* и Т*.

При увеличении и_П от значений порядка единицы до проходит через максимум, высота которого тем больше, чем больше и_П или m, а при и_П>И_п функция М (и_П) спадает до нуля. При уменьшении т максимум М (и_П) снижается, и, начиная с некоторых значений т, максимум исчезает и переход во внешнедиффузионную область при изменении параметра М происходит плавно.

Таким образом, исследование зависимости и_П (М) показывает, что при заданном значении параметра М в рассматриваемой системе возможны одно, три или пять соответствующих ему значений и_П и соответственно одно, три или пять стационарных решений. Этим решениям отвечают внутрикинетический режим I, внутридиффузионный II и внешнедиффузионный III режимы и чередующиеся с ними два промежуточных неустойчивых режима. Выпишем возможные режимы процесса при И > 4,5 и т 1, а также переходы между ними при изменении гидродинамических условий:

Характер переходов между режимами зависит от того, попадает ли величина Ш в интервал Ш_мин < Ш < Ш_макс, где задача с фиксированными значениями С₀, Т₀ имеет множественные решения. Возможные режимы при 4,5 > И >2 таковы же, как и в случае 3.

Для исследования переходов между режимами, связанных с изменением температуры потока и, следовательно, изменением параметра Ш, удобно пользоваться кривыми, изображающими зависимость параметра Ш от значения температуры в центре катализатора и_l при фиксированных прочих параметрах задачи Вi_Т, И, m.

Рассмотрим область небольших внешних разогревов и_П 2, где несущественно внешнедиффузионное сопротивление и И₀ И.

Из уравнения (41) получаем:

(48)

Подставляя это соотношение в формулу (42), определим и_П через температуру в центре катализатора и_l и параметры Вi_Т, И₀, Ш:

(49)

Выразив в уравнении (41) Ш₀ через Ш и воспользовавшись соотношением (49), находим искомую связь между и_l и Ш:

(50)

Из формулы (50) следует, что при отсутствии сопротивления теплопередаче на поверхности катализатора (Вi_T> ?) зависимость Ш от и_l совпадает с (41). При наличии внешнего теплового сопротивления начинает сказываться экспоненциальный фактор в формуле (50) и кривая Ш проходит ниже соответствующей кривой при Bi_т = ?. Во внутридиффузионной области, где Ш (и_l) принимает большие значения, функция Ш (и_l) проходит через максимум, после чего спадает до малых значений (см. рис.4). Значение функции в максимуме нетрудно определить, если учесть, что этот максимум находится в области значений и_l, весьма близких к И₀. При этом с высокой точностью можно считать и_l — И₀ под знаком радикала в уравнении (50):

(51)

Тогда для значения Ш (и_l) в максимуме получаем соотношение:

(52)

Величина Ш в максимуме может быть получена подстановкой соотношения (52) в (51):

(54)

Чтобы выяснить дальнейший ход кривой Ш (и_l) в области, соответствующей значительным внешним перепадам температуры щ, необходимо учесть, что максимально возможный внутренний разогрев И₀ с ростом и_П падает, и согласно формуле (48), И₀>0 при и_П> И_П. При этих значениях и_п и_l ?И₀ и, пользуясь соотношениями (51) и (37), можно найти следующую зависимость Ш (и_l):

(55)

В случае, когда величиной bИ_П нельзя пренебречь по сравнению с единицей, показатель экспоненты в формуле (54) следует заменить на m (и_l — И) /2 [1 + bm (И — и_l)]. Из формулы (54) следует, что при т 1 функция Ш (и_l) резко спадает до весьма малых значений при уменьшении и_l и, переходя через минимум, поднимается вдоль оси ординат при и_l 1.

Таким образом, задавая значения параметра Ш, можно определить стационарные температуры в центре катализатора, пересекая график Ш (и_l) прямыми Ш = const. На рис. 4 изображена зависимость Ш от и_l при И = 5 и Вi_т = ?, 20 и 10. При различных значениях параметра Ш может существовать одно, три или пять соответствующих ему значений и_l и соответственно одно, три или пять стационарных решений, которым отвечают внутрикинетический режим I, внутридиффузионный II и внешнедиффузионный III режимы и два промежуточных неустойчивых решения.

Рис. 4 Зависимость разогрева в центре пластины от модуля Тиле: И=5, Bi_T=? (a); 20 (б), 10 (в). Стрелки указывают направление скачкообразных переходов между режимами.

Переходы между режимами при изменении параметра Ш (изменении внешней температуры) зависят от относительной высоты двух максимумов на кривой Ш (и_l). Для каждого И > 4,5 существует такое значение параметра Bi_т = Вi_т^*, что при Bi_т < Вi_т^* значение функции в правом максимуме, которое линейно зависит от Вi_T, станет меньше значения функции в левом максимуме (см. рис.4). Оценку этого значения Вi_Т можно получить, считая, что вследствие слабого влияния экспоненциального фактора в формуле (50) в области левого максимума, значение Ш в нем не зависит от Вi_т и равно Ш_макс. Тогда, используя (53), получаем:

(56)

Возможны следующие схемы переходов между стационарными режимами процесса при изменении температуры среды и, следовательно, параметра Ш:

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

Выбор типа реакции и условий реализации промышленного процесса. Оптимальный температурный профиль

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

Понятие «наилучшего» должно обладать количественной мерой, называемой критерием оптимальности. Любой применяемый критерий оптимальности имеет экономическую природу и определяется, во-первых, изменением состава, а следовательно, и стоимости реагирующего потока в результате процесса и, во-вторых, затратами на ведение процесса. Не все составляющие критерия оптимальности имеют одинаковое значение. Некоторые из них могут быть настолько малы, что их разумно не принимать во внимание, и в каждом конкретном случае надо решать вопрос о том, каким упрощенным вариантом критерия оптимальности надо пользоваться.

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

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

Рассмотрим одну из простейших и в то же время теоретически наиболее важную задачу оптимизации — определение оптимальной температуры в каждом сечении реактора идеального вытеснения. Состав смеси продуктов реакции на выходе трубчатого реактора зависит от профиля температуры по длине аппарата. Очевидно, должен существовать такой продольный температурный профиль, при котором выход целевого продукта или, в общем случае, значение принятого критерия оптимальности будет максимальным. Этот температурный профиль будет для данного процесса оптимальным. Оптимальный температурный профиль (ОТП), как правило, практически не может быть реализован, однако к нему можно приблизиться, применяя различного типа секционированные реакторы. Теоретическая роль ОТП в реакторе идеального вытеснения состоит в том, что он дает наилучшие возможные показатели, достижимые в процессе с данной кинетикой — тот идеал, к которому следует стремиться при проектировании промышленного процесса. Из аналогии между процессом в реакторе идеального вытеснения и периодическим процессом следует, что в точности тот же результат будет достигнут в периодическом реакторе с температурой, оптимальным образом изменяющейся со временем. Подчеркивая эту аналогию, будем использовать вместо продольной координаты X текущее время контакта t = Х/и (где и - линейная скорость потока в реакторе идеального вытеснения).

Верхний предел температуры. В некоторых случаях можно сразу указать характер ОТП. Если рост температуры, ускоряя процесс, увеличивает также, независимо от состава реагирующей смеси, его избирательность (или, по крайней мере, не влияет на последнюю), оптимальная температура должна быть как можно более высокой.

Так, в случае единственной необратимой реакции повышение температуры только увеличивает ее скорость, а в случае обратимой эндотермической реакции — к тому же и смещает равновесие в сторону образования целевого продукта. Если, помимо основной реакции образования целевого продукта, имеется параллельная или (и) последовательная побочная реакция с энергией активации, меньшей, чем у основной, то повышение температуры увеличивает и скорость, и избирательность процесса. Во всех этих случаях температуру процесса следует поддерживать на верхнем допустимом пределе Т*. Эта предельная температура может определяться, например, условиями скачкообразного перехода процесса в диффузионный режим, при котором, вследствие сильного разогрева активной поверхности плавится или дезактивируется катализатор или начинают идти незаметные при низкой температуре побочные реакции. Другим фактором, ограничивающим допустимую температуру процесса, может быть возникновение при повышенных температурах нежелательных реакций, идущих в объеме (вне поверхности катализатора) по цепному механизму. Предельная температура Т*-зависит от состава реагирующей смеси и поэтому может изменяться по длине реактора. Так, если необратимая сильно экзотермическая реакция первого порядка, идущая на внешней поверхности непористых частиц, должна по технологическим условиям проводиться в кинетическом режиме, то верхний предел температуры определяется критическим соотношением между параметрами м и И в точке «зажигания» реакции. Сами параметры м и И определены формулами (56), где под С_? и Т_? следует теперь понимать текущие значения концентрации исходного вещества и температуры потока в данном сечении реактора. Из этих формул и вида кривой зажигания следует, что максимально допустимая температура процесса повышается с уменьшением концентрации исходного вещества. Если переход в диффузионную область допустим, то после этого перехода скорость реакции практически не зависит от температуры, так что дальнейшее повышение температуры становится ненужным.

Обратимая экзотермическая реакция. В случае обратимой экзотермической реакции повышение температуры ускоряет обратную реакцию сильнее, чем прямую, так что равновесие процесса смещается в нежелательную сторону. При некоторой температуре скорость образования целевого продукта r (С, Т) проходит через максимум; эта оптимальная температура зависит от состава реагирующей смеси, уменьшаясь по мере уменьшения концентрации исходного вещества. Таким образом, вначале, пока еще не накопилось значительное количество конечного продукта, процесс выгодно вести при высокой температуре, чтобы увеличить скорость прямой реакции, а затем температуру следует снижать, чтобы, сместив равновесие в нужную сторону, добиться максимального выхода целевого продукта. Температура в каждом сечении реактора должна быть выбрана так, чтобы скорость образования целевого вещества в этом сечении была максимальной, т. е., чтобы выполнялось равенство:

?r/?T=0 (57)

Если температурная зависимость константы скорости реакции определяется уравнением Аррениуса, где Е₁, Е₂ - энергия активации; Z₁, Z₃ — предэкспоненциальные множители соответственно прямой и обратной реакции; R — газовая постоянная; f (С) и g (С) — кинетические функции прямой и обратной реакций; в случае экзотермической реакции Е₂ > Е₁. Концентрации исходных веществ и продуктов реакции связаны между собой линейными соотношениями; поэтому функции, выражающие зависимость скорости прямой и обратной реакций от концентраций реагентов, можно выразить через концентрацию С одного из исходных веществ, которое принимается за ключевое. Дифференцируя функцию. (58) по температуре и приравнивая производную нулю, находим оптимальную температуру процесса как функцию концентрации ключевого вещества, где h = E₂ — E₁ - теплота реакции.

Условие для оптимальной температуры можно привести к форме, не зависящей от кинетических функций f ©, g (С). Обозначим через Т_р температуру, при которой смесь того же состава, что и в рассматриваемом сечении, находилась бы в равновесии. Так как r (Т_p) = 0.

Сравнение формул (IX.3) и (IX.4) приводит к соотношению.

Конечно, переход от формулы (59) к (61) возможен только в том случае, если равновесная температура, соответствующая данному составу, существует и кинетические функции f ©, g (С) правильно описывают кинетику процесса в окрестности равновесия. Формулы (59), (61) показывают, что оптимальная температура, так же, как и температура равновесия, снижается с увеличением степени превращения. Так как при малых степенях превращения обе величины неограниченно возрастают, в сечениях, близких к входному, оптимальная температура не может быть достигнута и температуру реакции следует фиксировать на верхнем пределе Т*. В частности, в случае реакции первого порядка равновесная температура обращается в бесконечность при степени превращения з = = 1 — С/С₀, равной и оптимальная — при степени превращения. Поэтому при степенях превращения з < з* повышение температуры всегда приводит к ускорению реакции.

Оптимизация стадийных cхем

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

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

Наиболее общим, но и самым трудоемким методом расчетного поиска оптимума является «эксперимент» на математической модели. Задавшись некоторой совокупностью значений независимых переменных, всегда можно путем решения системы расчетных уравнений вычислить соответствующее значение критерия оптимальности. Чтобы найти оптимум, не обязательно испытывать все возможные сочетания значений варьируемых переменных (для этого понадобился бы фантастический объем вычислительной работы); как и при экспериментальном поиске, здесь должен быть применен' один из методов направленного движения к оптимуму типа метода крутого восхождения. Поисковые методы редко бывают эффективны. Кроме того, немаловажно, что ни один поисковый метод не может дать информации об общей структуре оптимального решения для ряда сходных задач.

Единый подход к решению широкого класса задач на разыскание экстремума функции большого конечного числа переменных дает теория динамического программирования Беллмана. Сущность этой теории покажем на примере типичной задачи оптимизации, возникающей в химической технологии. Требуется найти оптимальный режим для последовательности N реакторов (или N-стадийного аппарата), причем на каждой стадии варьируется М независимых переменных. Пронумеруем реакторы в обратном порядке, так что первый номер присваивается последнему, а N-й — первому по ходу потока реактору. Состояние потока на выходе n-го реактора обозначим индексом n; в соответствии с этим исходное состояние потока обозначается индексом N + 1 (рис.5). Состояние реагирующего потока в общем случае описывается некоторым вектором X. Вектор X часто совпадает с вектором состава С; в более сложных случаях, однако, компонентами вектора Смогут быть, помимо концентраций ключевых веществ, также и температура потока, давление и пр. Если осуществить оптимальный выбор значений всех МN варьируемых переменных, получаем максимальное значение критерия оптимальности, — зависящее только от исходного состояния потока:

maxP_N=ц_N (X_N+1) (64)

Выделим из N-стадийной последовательности первый по ходу потока реактор и в соответствии с этим разобьем критерий для N-стадийного процесса P_N на два слагаемых: критерий для первого по ходу потока реактора p_N и критерий для оставшейся (N — 1) — стадийной последовательности Р_N-1.

P_N=P_N+P_N-1 (65)

Очевидно, что режим N-стадийной последовательности в целом может быть оптимальным только в том случае, если (N - 1) — стадийная последовательность работает в режиме, оптимальном относительно состояния потока, выходящего из N-го реактора. Действительно, сумма (65) не может стать максимальной, пока ее слагаемое Р_N-1 не достигнет максимального значения, зависящего только от состояния потока X_N на выходе из N-го реактора. Подставляя (65) в (64) и заменяя, согласно сказанному, P_N-1 на ц_N-1 (X_N), приходим к основному функциональному уравнению

Ц_N (X_N+1) =max [p_N+ц_N-1 (X_N)] (66)

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

Рис 5 Цепочка последовательно соединенных реакторов

В уравнении (66) максимум достигается варьированием только параметров, управляющих процессом на первой по ходу потока стадии. Принцип оптимальности позволяет, таким образом, заменить задачу одновременного выбора оптимальных значений MN независимых переменных гораздо более простой задачей N-стадийного выбора, на каждой стадии которого оптимум достигается варьированием М переменных. Другой отличительной чертой поиска оптимума методом динамического программирования является то, что задача решается не для единственного процесса с какими-то определенными параметрами исходного состояния, как это делается при использовании метода крутого восхождения, а для совокупности процессов с различными исходными состояниями. Действительно, как результат решения получаем зависимость максимального значения критерия оптимальности от параметров исходного состояния ц_N (X_N+1).