Главная
База знаний "Allbest"
Химия
Макрокинетика гетерогенных каталитических реакций. Устойчивость экзотермических реакций. Выбор типа реактора и условий реализации промышленного процесса
Макрокинетика гетерогенных каталитических реакций. Устойчивость экзотермических реакций. Выбор типа реактора и условий реализации промышленного процесса
Кинетические закономерности каталитического процесса, их определение истинной кинетикой реакции на активной поверхности и условиями массопереноса и теплопереноса. Определение оптимальной температуры в каждом сечении реактора идеального вытеснения.
| Рубрика | Химия |
| Вид | реферат |
| Язык | русский |
| Дата добавления | 23.10.2010 |
| Размер файла | 693,0 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
35
Федеральное агентство по образованию
Государственное образовательное учреждение
высшего профессионального образования
"ТОМСКИЙ ПОЛИТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ"
Химико-технологический факультет
Процессы и аппараты химических производств
Кафедра химической технологии топлива
"Макрокинетика гетерогенных каталитических реакций. Устойчивость экзотермических реакций. Выбор типа реактора и условий реализации промышленного процесса"
Реферат по дисциплине "Дополнительные главы физической химии".
Исполнитель:
маг. гр.5М550 Афанасьева Ю.И.
Руководитель:
АссистентЧеканцев Н.В.
Томск - 2010
Содержание
- Макрокинетика гетерогенных каталитических реакций
- Массоперенос и теплопередача в порах
- Внешняя массо- и теплопередача
- Устойчивость экзотермических реакций
- Выбор типа реакции и условий реализации промышленного процесса. Оптимальный температурный профиль
- Оптимизация стадийных cхем
- Список используемой литературы
Макрокинетика гетерогенных каталитических реакций
Реальные кинетические закономерности каталитического процесса определяются как истинной кинетикой реакции на активной поверхности, так и условиями массопереноса и теплопереноса; их изучение составляет предмет макрокинетики химических процессов.
Массоперенос и теплопередача в порах
Важное значение имеет перенос вещества и тепла внутри пористой частицы катализатора.
Перенос вещества в порах осуществляется исключительно путем молекулярной диффузии.
Если диаметр поры значительно превышает среднюю длину свободного пробега молекулы, то диффузия в порах протекает так же, как в объеме неподвижной жидкости или газа. В этом случае скорость переноса вещества вдоль поры определяется законом Фика.
Движущей силой диффузии является градиент концентрации диффундирующего вещества С вдоль направления поры X, а поток вещества направлен в сторону уменьшения концентрации.
Средняя скорость движения молекул зависит от температуры газа и в равновесном газе определяется законом Максвелла:
где k - константа Больцмана; Т - абсолютная температура; т - масса молекулы. Коэффициент молекулярной диффузии в газах может быть определен:
Величина Dм возрастает с повышением температуры, чем больше масса и диаметр молекулы, тем труднее, она диффундирует. Определение коэффициентов молекулярной диффузии в многокомпонентных смесях представляет собой чрезвычайно сложную задачу. При расчете химических процессов зависимостью коэффициентов диффузии от состава газовой смеси и температуры можно пренебречь.
В жидкостях величины коэффициентов молекулярной диффузии на несколько порядков ниже, чем в газах. Медленность молекулярной диффузии в жидкостях сильно ограничивает возможность применения пористых катализаторов для проведения жидкофазных реакций.
При уменьшении диаметра пор столкновения молекул со стенками поры начинают играть все большую роль в процессе молекулярной диффузии.
Когда диаметр поры намного меньше, чем средняя длина свободного пробега, коэффициент диффузии будет пропорционален уже не средней длине свободного пробега, а диаметру поры dП:
Эту область (л " dп) называют кнудсеновской, а величину Dкн - кнудсеновским коэффициентом диффузии.
При фиксированном диаметре поры переход в кнудсеновскую область диффузии происходит при уменьшении давления; он является плавным, поэтому в довольно широком интервале давлений для определения коэффициента диффузии пользуются различными приближенными переходными формулами, например, формулой гармонического среднего:
При исследовании макрокинетики химических реакций в пористом зерне нерационально рассматривать процесс в отдельной поре. Поры реальной частицы катализатора неодинаковы по размеру и образуют запутанную сеть. Поэтому пористое зерно рационально рассматривать как квазигомогенную среду, характеризуя скорость диффузии реагентов эффективным коэффициентом диффузии D, а скорость химической реакции - эффективной кинетической функцией r (С, Т).
Эффективный коэффициент диффузии D<Dм коэффициента молекулярной диффузии. Это объясняется, тем, что в пористом зерне доступно для диффузии не все сечение, а лишь часть, занятая порами, равная доле свободного объема и удлинением пути диффузии в извилистых порах. В общем случае эффективный коэффициент диффузии можно представить в виде:
D = аDм
где а - множитель, зависящий от структуры зерна, определяемый экспериментально.
В кнудсеновской области а падает с уменьшением размера пор. При определенных условиях в мелких порах проявляется, однако, дополнительный механизм переноса вещества за счет поверхностной диффузии молекул, адсорбированных на стенках пор. Этот эффект, в принципе, может приводить к величинам а, превышающим единицу.
Существенное влияние на величину D в катализаторах, содержащих узкие поры, оказывает распределение пор по размерам. При резко неоднородном распределении размеров пор само понятие эффективного коэффициента диффузии теряет определенность. Представим себе частицу, свободный объем которой состоит из сети широких транспортных макропор и множества отходящих от них узких капилляров, работающих в кнудсеновской области. Зерна такой структуры находят себе широкое применение, поскольку они сочетают хорошо развитую внутреннюю поверхность с относительно высокой скоростью диффузии, обеспечиваемой системой транспортных макропор. Измерение величины D в подобном зерне даст лишь величину D в макропорах. Между тем, химическая реакция, протекающая в основном в капиллярах, на которые приходится преобладающая часть внутренней поверхности катализатора, может лимитироваться гораздо более медленной диффузией в кнудсеновских микропорах.
При описании макрокинетики каталитической реакции на составных зернах применяют "двойную" диффузионную модель, вводя отдельные эффективные коэффициенты диффузии для системы транспортных макропор и для микропор в мелких гранулах. При этом сначала определяют зависимость скорости реакции в мелких гранулах от концентраций реагентов в транспортных макропорах, а затем вычисляют макроскопическую скорость реакции в зерне в целом с учетом диффузионного торможения в макропорах. Описывать составное зерно как квазигомогенную среду с эффективным коэффициентом диффузии, найденным в отсутствие химической реакции, можно только в предельных случаях, когда реакция либо не тормозится диффузией в микропорах, либо протекает настолько быстро, что локализуется на внешней поверхности малых гранул.
В пористом катализаторе перенос тепла осуществляется как с помощью молекулярного теплопереноса в порах, так и за счет теплопроводности самого катализатора.
В газах коэффициент молекулярной теплопроводности чм примерно равен коэффициенту молекулярной диффузии Dм, умноженному на теплоемкость единицы объема газа г. Эффективный коэффициент теплопроводности пористой частицы можно представить формулой
где ч0 - теплопроводность данной частицы в вакууме, которая определяется переносом тепла по твердой фазе.
Внешняя массо- и теплопередача
Существенное влияние на макроскопическую скорость каталитической реакции может оказывать массо- и теплообмен между внешней поверхностью частицы и омывающим ее потоком. Вблизи поверхности перенос вещества будет определяться молекулярной диффузией реагентов. .
Скорость потока плавно спадает по мере приближения к твердой поверхности. Чтобы найти поток вещества, диффундирующего на твердую поверхность, необходимо решить уравнение конвективной диффузии с граничными условиями, заданными на этой поверхности. В случае ламинарного движения распределение концентрации вещества:
где иx иу, иz - компоненты вектора скорости и, являющиеся функциями координат X, Y, Z.
Величину и направление скорости в каждой точке определяют решением уравнений гидродинамики. В правой части уравнения (7) оставлена производная только по координате X, нормальной к поверхности, так как по всем другим направлениям перенос вещества молекулярной диффузией пренебрежимо мал. Граничные условия для уравнения (7) определяются тем, что диффузионный поток на твердую поверхность катализатора равен скорости химической реакции, а на достаточном удалении от поверхности концентрация равна С?.
При больших скоростях движения практически весь перепад скорости сосредоточен в тонком гидродинамическом пограничном слое толщиной д0, а перепад концентрации - в диффузионном пограничном слое толщиной д. Отношение д0/д тем выше, чем больше отношение кинематической вязкости вещества х к коэффициенту молекулярной диффузии Dм. В жидкостях, где х/Dм>1. В этом случае при решении уравнения (7) используют простые выражения для скорости потока вблизи твердой поверхности, что позволяет найти аналитическое решение при протекании быстрой реакции или реакции первого порядка на поверхности частиц простой геометрической формы (пластина или шар). В газах толщины диффузионного и гидродинамического пограничных слоев одного порядка и уравнение (7) может быть решено только приближенными методами. Аналогичным способом - путем решения уравнения конвективной теплопроводности - определяют и тепловой поток на твердую поверхность
При обтекании тел сложной геометрической формы уравнения (7) становятся практически неразрешимыми. В этом случае эффективная толщина диффузионного слоя или коэффициент массопередачи в = DM/д определяются экспериментально. Опыты по измерению коэффициента массопередачи проводят в условиях, когда гетерогенная реакция протекает практически мгновенно и скорость исследуемого процесса лимитируется исключительно переносом вещества к активной поверхности.
При исследовании теплообмена между движущимся веществом и активной поверхностью, тепловой поток qт может быть выражен формулой
где б - коэффициент теплопередачи; Т? - Т0 - перепад температур между веществом и твёрдой поверхностью.
Функциональные зависимости, описывающие процессы массо - и теплопередачи, почти не отличаются друг от друга и для зернистого слоя сферических частиц выражаются формулами:
Как уже указывалось, процессы массопередачи и химической реакции не являются полностью независимыми. Поэтому нельзя вычислять макроскопическую скорость гетерогенного процесса, просто приравнивая величину q (количество вещества, диффундирующего через единицу поверхности в единицу времени) к скорости реакции на активной поверхности.
Устойчивость экзотермических реакций
Система уравнений материального и теплового балансов экзотермической реакции на пористом зерне катализатора:
С учетом сопротивления переносу вещества и тепла на внешней поверхности зерна Г:
(11)
где r - скорость реакции в единице объема пористого катализатора; Di и ч - эффективные коэффициенты диффузии и теплопроводности в пористом зерне; хi - стехиометрический коэффициент; h - теплота реакции; 2 - оператор Лапласа; Ci0 и Т0 - концентрации реагентов и температура на внешней поверхности зерна; Сi?, Т? - концентрация реагентов и температура в ядре потока; вi, б - коэффициенты массо - и теплопередачи из ядра потока к внешней поверхности зерна; п - направление внешней нормали к поверхности зерна катализатора.
Заметим, что вi = DiM/д и Di = aDiM (где Diм - коэффициент молекулярной диффузии 1-го вещества, д - толщина диффузионного пограничного слоя и а - коэффициент, зависящий от структуры пористого катализатора). Величины д и a можно с достаточной степенью точности считать одинаковыми для всех веществ, участвующих в реакции. Система уравнений (90) может быть сведена к единственному уравнению для концентрации одного из реагирующих веществ, которое принимают за ключевое. Введем с этой целью вспомогательную величину
(23)
Из уравнений (19) и граничных условий (21) следует, что величина wij определяется решением задачи
(24)
Очевидно, решение этой задачи имеет вид wij = const = wij? Из формулы (23) теперь следует, что концентрации всех веществ, участвующих в реакции, связаны между собой линейными соотношениями
(25)
и, таким образом, концентрации всех реагентов можно выразить через концентрацию одного ключевого вещества. Уравнения и граничные условия, определяющие концентрацию ключевого вещества и температуру в пористой частице катализатора, удобно записать в следующей безразмерной форме:
(26), (27), (28), (29)
где х - безразмерная координата (с масштабом l) направленная вдоль внешней нормали к поверхности, и Biм = вl/D, диффузионное и тепловое числа Био.
Выбор безразмерных переменных в уравнениях (26), (27) несколько необычен. За масштаб концентрации ключевого вещества внутри пористого зерна мы приняли концентрацию ключевого вещества на внешней поверхности зерна, так что с = С/С0. Аналогично, безразмерная температура внутри зерна определяется как и = (Т - То) /Т0* (где Т0* - масштаб внутренней температуры, зависящий от температуры внешней поверхности зерна). В соответствии с этим параметры Ш0 = l [r (С0, Т0) /DС0] 1/г, И0=hCoD/чT0* (параметр максимального разогрева), а также безразмерная скорость реакции f (с, и) = r (С, Т) /r (С0, Т 0) зависят от значений переменных на внешней поверхности зерна. Поскольку эти значения неизвестны заранее, величины Ш0 и И0, как и другие параметры, которые могут входить в выражение для f (с, и), не являются характерными параметрами рассматриваемого процесса. Преимущество такого выбора безразмерных переменных заключается, однако, в том, что в дальнейшем можно воспользоваться результатами исследования хорошо изученной задачи расчета экзотермической реакции на пористом зерне катализатора при фиксированных значениях концентрации и температуры на внешней поверхности зерна. Нетрудно видеть, что Ш0 и И0 являются характерными параметрами этой задачи. В граничных условиях, (29) участвуют безразмерная концентрация ключевого вещества на внешней поверхности частицы cП =C0/C? и безразмерная температура внешней поверхности иП= (T0-T?) /T*, масштабы которых связаны уже только со значениями переменных в ядре потока.
Пользуясь уравнениями (26), (27) и граничными условиями (28), (29), можно получить оценку условий, при которых существует заметный перепад концентрации и температуры между поверхностью катализатора и внешней средой. Внешнее сопротивление массопередаче начинает сказываться только, когда реакция локализуется в тонком слое, толщина которого сравнима с толщиной диффузионного пограничного слоя д. Действительно, поскольку величина Ш0 является мерой проникновения реакции в глубь пористого катализатора, так что dc/dx0, из граничного условия (28) следует:
(30)
Так как всегда , внешний перепад концентрации сравним с внутренним только при .
Для оценки перепада температуры воспользуемся линейной связью между концентрацией и температурой внутри пористой частицы, которая следует из уравнений (25), (26)
(31)
откуда
(32)
Принимая во внимание, что для газовой среды (где чм - коэффициент молекулярной теплопроводности газа) и используя граничное условие (29), находим:
(33)
Так как перенос тепла в пористом катализаторе идет в основном по твердой фазе, то . Для сильно экзотермических реакций величина И0 может принимать значения порядка нескольких единиц. Поэтому величина иП, характеризующая влияние внешнего сопротивления теплоотводу, может быть существенной () при сравнительно небольших значениях Ш0, когда внешнедиффузионное сопротивление еще не оказывает заметного влияния на процесс.
Учитывая, что тепловой поток на внешнюю поверхность катализатора равен скорости тепловыделения в результате реакции, граничное условие (22) или (20) можно преобразовать, введя фактор эффективности з, вычисленный для процесса с фиксированными значениями концентрации и температуры на внешней поверхности катализатора:
(34)
Используя значения фактора эффективности з, найденные путем решения задачи с фиксированными С0, T0, можно вычислить правую часть уравнения (34) ц (иП); после этого значения внешнего разогрева иП определяются как точки пересечения графика функции ц (иП) с прямыми с наклоном ВiТ, проведенными через начало координат.
Поскольку величины Ш0 и И0 зависят от поверхностных значений концентрации и температуры и потому не являются характерными параметрами процесса, они должны быть выражены через параметры
,
отнесенные к условиям в ядре потока. Установлено, что при не слишком больших внешних разогревах иП С0 ? С? и, следовательно, (Т0*/Т*) И0 ? И Ш0 = Шу (иП) (где у (ип) - сильно возрастающая функция иП).
В случае, когда зависимость з (Шо) неоднозначна, график зависимости правой части ц (иП) уравнения (34) от иП имеет форму, показанную на рис, ШЛО. При увеличении иП функция ц (иП) резко возрастает за счет сильной зависимости Ш0 от иП. При дальнейшем увеличении иП. когда начинает сказываться внешне-диффузионное торможение процесса, функция ц (иП) проходит через максимум и резко спадает до нуля. Решения уравнения (111.115) определяются как точки пересечения кривой ц (иП) и прямой с наклоном ВiТ. Из рис.1 видно, что уравнение (34) может иметь до пяти решений, причем при значениях числа ВiТ меньше некоторых критических, решения, соответствующие малым внешним разогревай, исчезают. При увеличении температуры потока, а следовательно, параметра Ш кривая ц (иП) сдвигается вверх, и при значениях Ш, больших некоторых критических, также исчезают низкотемпературные решения. Участку кривой ц (иП) слева от точки 1 соответствует внутрикинетический, а между точками 2 и 3 - внутридиффузиошшй режим протекания реакции. В задаче с фиксированными значениями поверхностных концентраций и температуры внутрикинетический режим существует левее точки 4, а внутридиффузионный - правее точки 5. Наличие сопротивления теплоотводу на внешней поверхности зерна приводит к расширению области промежуточных неустойчивых состояний. Кроме того, при наличии внешнего сопротивления теплоотводу множественные режимы возникают и в том случае, когда зависимость з (Ш0) однозначна.
Рис.1. Графическое решение уравнения (34)
Промежуточных неустойчивых состояний. Кроме того, при наличии внешнего сопротивления теплоотводу множественные режимы возникают и в том случае, когда зависимость т) (Ч'о) однозначна.
Реакция первого порядка. Рассмотрим более подробно необратимую экзотермическую реакцию первого порядка на пористом катализаторе, имеющем форму пластины толщиной 2l с "запечатанными" боковыми гранями. Это простейший пример, на котором видны все характерные особенности рассматриваемой задачи. Пусть . Если принять за масштаб температуры Т величину , то
(35)
Величины Ш0 иИ0 связаны с характерными параметрами задачи Ш, И=чR следующими соотношениями, вытекающими из формул (28), (29)
(36), (37)
При записи параметра в за масштаб внешней температуры принята величина . При не слишком больших иП (иП1) величины и Т* практически совпадают и, так как . Заметим, что из формулы (37) следует существование максимально возможного разогрева ип. макс=Ип=mИ.
Используя соотношение (31) между концентрацией и температурой внутри катализатора, сведем систему (26), (27) к единственному уравнению относительно безразмерной температуры и
(38)
где x - безразмерная координата, принимающая значения x=0 в центре пластины и x=±1 - на ее поверхности.
Переходя к интегрированию уравнения (38), будем считать величину иП, а следовательно, и значения параметров Ш0 и И0 известными. Далее можно найти величину иП, подставляя полученное решение в граничное условие (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ШИ0, а минимум выражения (46) достигается при иП = 2. Из выражения (46) следует, что во внутридиффузионном режиме правая часть уравнения (42) не зависит от параметра Ш и, следовательно, значение иП в этой области определяется только параметрами М и И0.
Зависимость иП (М) при известных значениях Ш и И0 может быть построена с помощью графика функции ц (и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 (соответствующие значениям Ш0, при которых в задаче с заданными значениями поверхностных концентраций С0 и температуры Т0 происходил перескок между режимами) попадают в область неустойчивых режимов, разделяющих внутрикинетическую. и внутри-диффузионную области. Таким образом, наличие сопротивления теплоотводу на внешней поверхности катализатора приводит к расширению области неустойчивых режимов. Это проявляется также в возникновении неустойчивой области при И < 4,5, когда задача без внешнего сопротивления вообще не имеет неустойчивых режимов. Скачкообразные переходы во внутренних режимах могут наблюдаться, как это видно из рис.3, до И = 2.
При значениях параметра М, меньших абсолютного минимума функции М (иП) (этот минимум может соответствовать точкам 1 и 5, см. рис.2,3), в рассматриваемой области режимов с небольшими внешними разогревами не существует и процесс переходит во внешнедиффузионный режим (иП ~ Ип). Во внешнедиффузионном режиме существенным становится перепад концентрации у поверхности катализатора, и поэтому для получения зависимости М (иП) в этой области значений иП необходимо учитывать связь И0 с иП, определяемую формулой (37). Поскольку при этом реакция локализована в узком слое близ поверхности катализатора из-за большого внутридиффузионного торможения, ход кривой М (иП) во внешнедиффузионной области можно получить из формулы (46) с учетом зависимости (37). При больших иП можно также заменить приближенную кинетическую функцию на ехр ип/2 (1 + bиП) (где b == RТ?/Е) и учесть различие между величинами T0* и Т*.
При увеличении иП от значений порядка единицы до проходит через максимум, высота которого тем больше, чем больше иП или m, а при иП>Ип функция М (иП) спадает до нуля. При уменьшении т максимум М (иП) снижается, и, начиная с некоторых значений т, максимум исчезает и переход во внешнедиффузионную область при изменении параметра М происходит плавно.
Таким образом, исследование зависимости иП (М) показывает, что при заданном значении параметра М в рассматриваемой системе возможны одно, три или пять соответствующих ему значений иП и соответственно одно, три или пять стационарных решений. Этим решениям отвечают внутрикинетический режим I, внутридиффузионный II и внешнедиффузионный III режимы и чередующиеся с ними два промежуточных неустойчивых режима. Выпишем возможные режимы процесса при И > 4,5 и т 1, а также переходы между ними при изменении гидродинамических условий:
Характер переходов между режимами зависит от того, попадает ли величина Ш в интервал Шмин < Ш < Шмакс, где задача с фиксированными значениями С0, Т0 имеет множественные решения. Возможные режимы при 4,5 > И >2 таковы же, как и в случае 3.
Для исследования переходов между режимами, связанных с изменением температуры потока и, следовательно, изменением параметра Ш, удобно пользоваться кривыми, изображающими зависимость параметра Ш от значения температуры в центре катализатора иl при фиксированных прочих параметрах задачи ВiТ, И, m.
Рассмотрим область небольших внешних разогревов иП 2, где несущественно внешнедиффузионное сопротивление и И0 И.
Из уравнения (41) получаем:
(48)
Подставляя это соотношение в формулу (42), определим иП через температуру в центре катализатора иl и параметры ВiТ, И0, Ш:
(49)
Выразив в уравнении (41) Ш0 через Ш и воспользовавшись соотношением (49), находим искомую связь между иl и Ш:
(50)
Из формулы (50) следует, что при отсутствии сопротивления теплопередаче на поверхности катализатора (ВiT> ?) зависимость Ш от иl совпадает с (41). При наличии внешнего теплового сопротивления начинает сказываться экспоненциальный фактор в формуле (50) и кривая Ш проходит ниже соответствующей кривой при Biт = ?. Во внутридиффузионной области, где Ш (иl) принимает большие значения, функция Ш (иl) проходит через максимум, после чего спадает до малых значений (см. рис.4). Значение функции в максимуме нетрудно определить, если учесть, что этот максимум находится в области значений иl, весьма близких к И0. При этом с высокой точностью можно считать иl - И0 под знаком радикала в уравнении (50):
(51)
Тогда для значения Ш (иl) в максимуме получаем соотношение:
(52)
Величина Ш в максимуме может быть получена подстановкой соотношения (52) в (51):
(54)
Чтобы выяснить дальнейший ход кривой Ш (иl) в области, соответствующей значительным внешним перепадам температуры щ, необходимо учесть, что максимально возможный внутренний разогрев И0 с ростом иП падает, и согласно формуле (48), И0>0 при иП> ИП. При этих значениях ип иl ?И0 и, пользуясь соотношениями (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, BiT=? (a); 20 (б), 10 (в). Стрелки указывают направление скачкообразных переходов между режимами.
Переходы между режимами при изменении параметра Ш (изменении внешней температуры) зависят от относительной высоты двух максимумов на кривой Ш (иl). Для каждого И > 4,5 существует такое значение параметра Biт = Вiт*, что при Biт < Вiт* значение функции в правом максимуме, которое линейно зависит от ВiT, станет меньше значения функции в левом максимуме (см. рис.4). Оценку этого значения ВiТ можно получить, считая, что вследствие слабого влияния экспоненциального фактора в формуле (50) в области левого максимума, значение Ш в нем не зависит от Вiт и равно Шмакс. Тогда, используя (53), получаем:
(56)
Возможны следующие схемы переходов между стационарными режимами процесса при изменении температуры среды и, следовательно, параметра Ш:
Интересной особенностью переходов между режимами в рассмотренной системе, является то, что при некоторых условиях внутри-диффузионный режим, хотя и существует, по не проявляется ни при увеличении, ни при уменьшении температуры и скорости потока.
Выбор типа реакции и условий реализации промышленного процесса. Оптимальный температурный профиль
Основные положения. В любое из расчетных уравнений химического процесса входит ряд переменных: время контакта, температура потока на входе в реактор и температура теплоносителя, скорость потока, диаметр зерна катализатора и т.д., значения которых можно изменять в более или менее широких пределах. Приступая к проектированию химического реактора, необходимо выбрать значения этих переменных так, чтобы добиться наилучшего результата процесса. Число и номенклатура варьируемых переменных определяются принятым типом реактора и его схемой. Последняя также должна быть выбрана оптимальной, а этого в большинстве случаев можно добиться только путем сравнения наилучших результатов процесса, достижимых в реакторах различных типов.
Понятие "наилучшего" должно обладать количественной мерой, называемой критерием оптимальности. Любой применяемый критерий оптимальности имеет экономическую природу и определяется, во-первых, изменением состава, а следовательно, и стоимости реагирующего потока в результате процесса и, во-вторых, затратами на ведение процесса. Не все составляющие критерия оптимальности имеют одинаковое значение. Некоторые из них могут быть настолько малы, что их разумно не принимать во внимание, и в каждом конкретном случае надо решать вопрос о том, каким упрощенным вариантом критерия оптимальности надо пользоваться.
Во многих случаях упрощение критерия оптимальности компенсируется введением дополнительных условий и ограничений, которым должно удовлетворять оптимальное решение. Простейшим, чисто химическим критерием оптимальности может быть выход целевого продукта.
Задача оптимального проектирования реактора определенного типа сводится, таким образом, к разысканию максимального значения критерия оптимальности путем варьирования ряда независимых переменных, допустимые значения которых обычно ограничены технологическими пределами. Проведение процесса в рассчитанном режиме даст наилучший результат, достижимый (в реакторе данного типа) на данном катализаторе при принятых условиях и ограничениях. Сравнивая максимальные значения критерия оптимальности для реакторов различных типов, можно определить, какой тип реактора предпочтителен для осуществления данного процесса.
Рассмотрим одну из простейших и в то же время теоретически наиболее важную задачу оптимизации - определение оптимальной температуры в каждом сечении реактора идеального вытеснения. Состав смеси продуктов реакции на выходе трубчатого реактора зависит от профиля температуры по длине аппарата. Очевидно, должен существовать такой продольный температурный профиль, при котором выход целевого продукта или, в общем случае, значение принятого критерия оптимальности будет максимальным. Этот температурный профиль будет для данного процесса оптимальным. Оптимальный температурный профиль (ОТП), как правило, практически не может быть реализован, однако к нему можно приблизиться, применяя различного типа секционированные реакторы. Теоретическая роль ОТП в реакторе идеального вытеснения состоит в том, что он дает наилучшие возможные показатели, достижимые в процессе с данной кинетикой - тот идеал, к которому следует стремиться при проектировании промышленного процесса. Из аналогии между процессом в реакторе идеального вытеснения и периодическим процессом следует, что в точности тот же результат будет достигнут в периодическом реакторе с температурой, оптимальным образом изменяющейся со временем. Подчеркивая эту аналогию, будем использовать вместо продольной координаты X текущее время контакта t = Х/и (где и - линейная скорость потока в реакторе идеального вытеснения).
Верхний предел температуры. В некоторых случаях можно сразу указать характер ОТП. Если рост температуры, ускоряя процесс, увеличивает также, независимо от состава реагирующей смеси, его избирательность (или, по крайней мере, не влияет на последнюю), оптимальная температура должна быть как можно более высокой.
Так, в случае единственной необратимой реакции повышение температуры только увеличивает ее скорость, а в случае обратимой эндотермической реакции - к тому же и смещает равновесие в сторону образования целевого продукта. Если, помимо основной реакции образования целевого продукта, имеется параллельная или (и) последовательная побочная реакция с энергией активации, меньшей, чем у основной, то повышение температуры увеличивает и скорость, и избирательность процесса. Во всех этих случаях температуру процесса следует поддерживать на верхнем допустимом пределе Т*. Эта предельная температура может определяться, например, условиями скачкообразного перехода процесса в диффузионный режим, при котором, вследствие сильного разогрева активной поверхности плавится или дезактивируется катализатор или начинают идти незаметные при низкой температуре побочные реакции. Другим фактором, ограничивающим допустимую температуру процесса, может быть возникновение при повышенных температурах нежелательных реакций, идущих в объеме (вне поверхности катализатора) по цепному механизму. Предельная температура Т*-зависит от состава реагирующей смеси и поэтому может изменяться по длине реактора. Так, если необратимая сильно экзотермическая реакция первого порядка, идущая на внешней поверхности непористых частиц, должна по технологическим условиям проводиться в кинетическом режиме, то верхний предел температуры определяется критическим соотношением между параметрами м и И в точке "зажигания" реакции. Сами параметры м и И определены формулами (56), где под С? и Т? следует теперь понимать текущие значения концентрации исходного вещества и температуры потока в данном сечении реактора. Из этих формул и вида кривой зажигания следует, что максимально допустимая температура процесса повышается с уменьшением концентрации исходного вещества. Если переход в диффузионную область допустим, то после этого перехода скорость реакции практически не зависит от температуры, так что дальнейшее повышение температуры становится ненужным.
Обратимая экзотермическая реакция. В случае обратимой экзотермической реакции повышение температуры ускоряет обратную реакцию сильнее, чем прямую, так что равновесие процесса смещается в нежелательную сторону. При некоторой температуре скорость образования целевого продукта r (С, Т) проходит через максимум; эта оптимальная температура зависит от состава реагирующей смеси, уменьшаясь по мере уменьшения концентрации исходного вещества. Таким образом, вначале, пока еще не накопилось значительное количество конечного продукта, процесс выгодно вести при высокой температуре, чтобы увеличить скорость прямой реакции, а затем температуру следует снижать, чтобы, сместив равновесие в нужную сторону, добиться максимального выхода целевого продукта. Температура в каждом сечении реактора должна быть выбрана так, чтобы скорость образования целевого вещества в этом сечении была максимальной, т.е., чтобы выполнялось равенство:
?r/?T=0 (57)
Если температурная зависимость константы скорости реакции определяется уравнением Аррениуса, где Е1, Е2 - энергия активации; Z1, Z3 - предэкспоненциальные множители соответственно прямой и обратной реакции; R - газовая постоянная; f (С) и g (С) - кинетические функции прямой и обратной реакций; в случае экзотермической реакции Е2 > Е1. Концентрации исходных веществ и продуктов реакции связаны между собой линейными соотношениями; поэтому функции, выражающие зависимость скорости прямой и обратной реакций от концентраций реагентов, можно выразить через концентрацию С одного из исходных веществ, которое принимается за ключевое. Дифференцируя функцию. (58) по температуре и приравнивая производную нулю, находим оптимальную температуру процесса как функцию концентрации ключевого вещества, где h = E2 - E1 - теплота реакции.
Условие для оптимальной температуры можно привести к форме, не зависящей от кинетических функций f (С), g (С). Обозначим через Тр температуру, при которой смесь того же состава, что и в рассматриваемом сечении, находилась бы в равновесии. Так как r (Тp) = 0.
Сравнение формул (IX.3) и (IX.4) приводит к соотношению.
Конечно, переход от формулы (59) к (61) возможен только в том случае, если равновесная температура, соответствующая данному составу, существует и кинетические функции f (С), g (С) правильно описывают кинетику процесса в окрестности равновесия. Формулы (59), (61) показывают, что оптимальная температура, так же, как и температура равновесия, снижается с увеличением степени превращения. Так как при малых степенях превращения обе величины неограниченно возрастают, в сечениях, близких к входному, оптимальная температура не может быть достигнута и температуру реакции следует фиксировать на верхнем пределе Т*. В частности, в случае реакции первого порядка равновесная температура обращается в бесконечность при степени превращения з = = 1 - С/С0, равной и оптимальная - при степени превращения. Поэтому при степенях превращения з < з* повышение температуры всегда приводит к ускорению реакции.
Оптимизация стадийных cхем
Методы оптимизации стадийных схем. При проектировании реальных технологических процессов оптимальному выбору подлежит лишь ограниченное число параметров. В этом состоит существенное отличие практических задач оптимизации от задачи определения ОТП. В гибкой и эффективной реакторной схеме число варьируемых параметров Ф может, однако, быть весьма большим. Примером такой схемы является цепочка последовательно соединенных реакторов, размеры и режим работы которых должны быть выбраны оптимальным способом.
Классический метод поиска максимума функции Ф переменных состоит, как известно, в следующем. Определяются и приравниваются нулю частные производные функции по всем независимым переменным; в результате получается Ф уравнений, совместное решение которых дает искомое положение максимума. Этот метод чрезвычайно громоздок при большом Ф, а, кроме того, часто неосуществим по той причине, что аналитический вывод уравнений, определяющих точку оптимума, невозможен. Другой причиной непригодности классического метода является наличие технологических пределов варьирования независимых переменных. Может оказаться, что критерий оптимальности вовсе не имеет максимума в аналитическом смысле, а его наивысшее значение достигается на одной из границ разрешенной области, т, е. когда одна или несколько независимых переменных фиксированы на предельных значениях.
Наиболее общим, но и самым трудоемким методом расчетного поиска оптимума является "эксперимент" на математической модели. Задавшись некоторой совокупностью значений независимых переменных, всегда можно путем решения системы расчетных уравнений вычислить соответствующее значение критерия оптимальности. Чтобы найти оптимум, не обязательно испытывать все возможные сочетания значений варьируемых переменных (для этого понадобился бы фантастический объем вычислительной работы); как и при экспериментальном поиске, здесь должен быть применен' один из методов направленного движения к оптимуму типа метода крутого восхождения. Поисковые методы редко бывают эффективны. Кроме того, немаловажно, что ни один поисковый метод не может дать информации об общей структуре оптимального решения для ряда сходных задач.
Единый подход к решению широкого класса задач на разыскание экстремума функции большого конечного числа переменных дает теория динамического программирования Беллмана. Сущность этой теории покажем на примере типичной задачи оптимизации, возникающей в химической технологии. Требуется найти оптимальный режим для последовательности N реакторов (или N-стадийного аппарата), причем на каждой стадии варьируется М независимых переменных. Пронумеруем реакторы в обратном порядке, так что первый номер присваивается последнему, а N-й - первому по ходу потока реактору. Состояние потока на выходе n-го реактора обозначим индексом n; в соответствии с этим исходное состояние потока обозначается индексом N + 1 (рис.5). Состояние реагирующего потока в общем случае описывается некоторым вектором X. Вектор X часто совпадает с вектором состава С; в более сложных случаях, однако, компонентами вектора Смогут быть, помимо концентраций ключевых веществ, также и температура потока, давление и пр. Если осуществить оптимальный выбор значений всех МN варьируемых переменных, получаем максимальное значение критерия оптимальности, - зависящее только от исходного состояния потока:
maxPN=цN (XN+1) (64)
Выделим из N-стадийной последовательности первый по ходу потока реактор и в соответствии с этим разобьем критерий для N-стадийного процесса PN на два слагаемых: критерий для первого по ходу потока реактора pN и критерий для оставшейся (N - 1) - стадийной последовательности РN-1.
PN=PN+PN-1 (65)
Очевидно, что режим N-стадийной последовательности в целом может быть оптимальным только в том случае, если (N - 1) - стадийная последовательность работает в режиме, оптимальном относительно состояния потока, выходящего из N-го реактора. Действительно, сумма (65) не может стать максимальной, пока ее слагаемое РN-1 не достигнет максимального значения, зависящего только от состояния потока XN на выходе из N-го реактора. Подставляя (65) в (64) и заменяя, согласно сказанному, PN-1 на цN-1 (XN), приходим к основному функциональному уравнению
ЦN (XN+1) =max [pN+цN-1 (XN)] (66)
которое представляет собой математическое выражение принципа оптимальности, лежащего в основе теории динамического программирования и сформулированного Беллманом: "Оптимальная стратегия обладает тем свойством, что, каковы бы ни были исходное состояние и решение в начальный момент, последующие решения должны составлять оптимальную стратегию относительно состояния, получающегося в результате первого решения".
Рис 5 Цепочка последовательно соединенных реакторов
В уравнении (66) максимум достигается варьированием только параметров, управляющих процессом на первой по ходу потока стадии. Принцип оптимальности позволяет, таким образом, заменить задачу одновременного выбора оптимальных значений MN независимых переменных гораздо более простой задачей N-стадийного выбора, на каждой стадии которого оптимум достигается варьированием М переменных. Другой отличительной чертой поиска оптимума методом динамического программирования является то, что задача решается не для единственного процесса с какими-то определенными параметрами исходного состояния, как это делается при использовании метода крутого восхождения, а для совокупности процессов с различными исходными состояниями. Действительно, как результат решения получаем зависимость максимального значения критерия оптимальности от параметров исходного состояния цN (XN+1).
Подобные документы
Определение температуры газового потока на входе в реакторе, обеспечивающей максимальную производительность реактора. Программа для расчета, составляется в приложении REAC. График зависимости производительности реактора от температуры газового потока.
контрольная работа [36,0 K], добавлен 14.06.2011Основные параметры реакторов идеального вытеснения и полного смешения. Расчет необходимого времени пребывания реагентов в реакционной зоне. Параметры химико-технологического процесса в потоке полного смешения при изотермическом температурном режиме.
контрольная работа [171,6 K], добавлен 14.06.2011Характеристика процесса проектирования реактора. Описание материальных моделей химических реакторов: идеального вытеснения, полного смешения. Технологическое оформление процесса синтеза аммиака. Основные требования, предъявляемые к промышленным реакторам.
курсовая работа [620,7 K], добавлен 16.05.2012Стадии цепных разветвленных реакций. Стационарный и нестационарный режимы быстрого самоускорения. Зависимость пределов воспламенения от давления, температуры и критических размеров реактора. Кинетика цепных реакций с вырожденным разветвлением цепей.
реферат [182,5 K], добавлен 09.03.2015Технология синтеза аммиака. Материальный и тепловой балансы РИВ и РПС. Выбор адиабатического реактора для синтеза NH3. Расчет адиабатического коэффициента. Анализ зависимости объема реактора от начальной температуры, давления и степени превращения.
курсовая работа [523,3 K], добавлен 22.04.2012Основные требования к промышленным реакторам. Термодинамика и кинетика окисления диоксида серы. Математические модели химических реакторов. Модель реактора идеального вытеснения и полного смешения. Получение максимальной степени окисления диоксида серы.
курсовая работа [284,2 K], добавлен 17.06.2010Физико-химические основы процесса метанирования, применение катализаторов и промышленные схемы. Программа расчета адиабатического коэффициента для выбора типа реактора, определение зависимости производительности от давления и начальной концентрации.
курсовая работа [1008,4 K], добавлен 09.06.2011Последовательность расчета материального баланса реактора синтеза аммиака. Мольные потоки компонентов. Работа реакторов идеального вытеснения и полного смешения. Определение зависимости производительности реактора от давления и начальной концентрации.
контрольная работа [197,0 K], добавлен 06.10.2014Хлороводород: производство, применение. Выбор адиабатического реактора для синтеза HCl. Программа расчета адиабатического коэффициента. Программа и анализ зависимости объема реактора от начальной температуры, степени превращения, начальной концентрации.
курсовая работа [80,2 K], добавлен 17.05.2012Технологические особенности процесса сернокислотного алкилирования изопарафинов олефинами. Выбор типа и конструкции реактора. Механизм пиролиза пентана. Катализаторы риформинга и уравнения протекающих реакций. Этерификация спиртов карбоновыми кислотами.
реферат [1,0 M], добавлен 28.02.2009
