|
Макрокинетика гетерогенных каталитических реакций. Устойчивость экзотермических реакций. Выбор типа реактора и условий реализации промышленного процессаМакрокинетика гетерогенных каталитических реакций. Устойчивость экзотермических реакций. Выбор типа реактора и условий реализации промышленного процесса35 Федеральное агентство по образованию Государственное образовательное учреждение высшего профессионального образования "ТОМСКИЙ ПОЛИТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ" Химико-технологический факультет Процессы и аппараты химических производств Кафедра химической технологии топлива "Макрокинетика гетерогенных каталитических реакций. Устойчивость экзотермических реакций. Выбор типа реактора и условий реализации промышленного процесса" Реферат по дисциплине "Дополнительные главы физической химии". Исполнитель: маг. гр.5М550 Афанасьева Ю.И. Руководитель: АссистентЧеканцев Н.В. Томск - 2010 Содержание
откуда (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).Уравнение (66) представляет собой соотношение, позволяющее найти максимальное значение критерия оптимальности для N-стадийной последовательности цN, если максимальное значение критерия для (N - 1) - стадийной последовательности цN-1 уже известно. Этот факт и определяет процедуру поиска оптимума методом динамического программирования.Так как выбор оптимальных условий для реакторов, стоящих в начале последовательности, должен осуществляться с учетом влияния работы этих реакторов на процесс в последующих аппаратах, а обратное влияние отсутствует, расчет N-стадийной последовательности следует начинать с расчета реактора, последнего походу потока, и вести, оптимизируя цепочки реакторов со все большим числом членов, подключая к этой цепочке реакторы, все более удаленные от конца последовательности. Начинаем с "последовательности", содержащей ноль реакторов. В такой "последовательности", разумеется, никаких превращений не происходит и максимальное значение критерия оптимальности ц0 (X1) тождественно равно нулю. Далее рассматриваем последовательность, состоящую из одного реактора. Подставляя в (66) N = 1 и ц0 = 0, получаем:ц1 (X2) =max p1 (67)Вектор Х2 описывает состояние потока, выходящего из реактора 2, и потому пока неизвестен; поэтому приходится определять оптимальные значения варьируемых параметров и соответствующие им максимальные значения критерия оптимальности ц1 (X2) для некоторой более или менее широкой совокупности исходных составов и табулировать полученные результаты. После этого можно перейти к расчету двухстадийной последовательности. Пользуясь принципом оптимальности (66), находим оптимальные значения параметров, управляющих процессом в реакторе 2. Варьируя эти параметры, меняем состояние потока Х2 на выходе из реактора 2; при этом изменяются как величина р2 (критерий для реактора 2), так и уже вычисленное максимальное значение критерия оптимальности ц1 (Х2) для реактора 1. Максимизируется сумма этих величин; для реактора 2 оптимальный режим определяется, таким образом, с учетом не только "локальной пользы", но и влияния работы этого реактора на дальнейший ход процесса. После того как вычислена функция ц2 (X3) (в определенной области значений исходных состояний X3), можно приступать к расчету трехстадийной последовательности и т.д., вплоть до любого N.Описанная общая схема метода динамического программирования, впервые примененная к расчету химических реакторов Арисом, является универсальной, но содержит два неприятных момента. Первый из них - разыскание максимума в уравнении (66), выполняемое на каждой стадии. Здесь можно воспользоваться поисковым методом крутого восхождения, что не вызовет больших затруднений, так как поиск ведется только по малому числу М варьируемых переменных. Второй, еще более неприятный момент - необходимость вычислять на каждой стадии: максимальные значения критерия оптимальности как функции параметров исходного состояния и табулировать вычисленные величины. Эта чрезвычайно трудоемкая задача сильно ограничивает возможности практического применения общей схемы метода динамического программирования и делает ее подчас менее эффективной, чем поисковый метод крутого восхождения. Оба метода являются но своему характеру типично машинными, так как связаны с многократным использованием одних и тех же расчетных формул и уравнений, т.е. многократным повторением относительно простой программы. Вопрос о том, какой из методов рационально применить при расчете, должен решаться в зависимости от конкретных условий задачи. Заранее можно сказать лишь то, что преимущества метода динамического программирования будут сказываться при увеличении числа стадий и при большей неопределенности исходного состояния потока. Рост числа компонентов вектора состояния X, затрудняя задачу табулирования, делает предпочтительным применением метода крутого восхождения.От недостатков общей схемы метода динамического программирования можно, однако, в значительной мере избавиться, используя аналитический метод поиска оптимума на каждой стадии. Именно этот способ будет применен к решению задач оптимизации цепочек реакторов, рассматриваемых ниже. Отметим, что основные расчетные формулы, которые получим, могут быть выведены не только с помощью метода динамического программирования, но и на основе дискретного варианта принципа максимума Понтрягина или классических вариационных методов. Список используемой литературы1. Иоффе И.И., Письмен Л.М. Инженерная химия гетерогенного катализа. Изд-во "Химия", Л., 1972, стр.464, табл.8, рис 102. |
 |
||||||||||||||||||
| ©2011 |
