Слайд 1
Непрерывные реакторы для гомогенных гомофазных процессов.
Слайд 2Непрерывные реакторы для гомогенных гомофазных процессов.
Существуют реакции, проводимые в пламени и
требующие реакторов специального типа. К ним относятся, например, реакторы для парциального окисления углеводородов.
Слайд 3 Конструкционные особенности непрерывных реакторов для гомогенных гомофазных процессов, в наибольшей мере
зависят от способа организации теплообмена. Это, в первую очередь, касается газофазных фракций.
Относительно низкотемпературные жидкофазные реакции чаще всего реализуются в политропических реакторах, имеющих поверхность теплопередачи с хладоагентом (теплоноситель). В реакторах вытеснения обмен теплот осуществляется через поверхность трубы или пучка труб, внутри которых двигается реакционная смесь, а с внешней стороны – хладоагент (теплоноситель). Жидкофазные реакторы смешения снабжены рубашками, встроенными или выносными теплообменниками.
Слайд 4 Для газофазных гомогенных реакций широко распространены и политропический и адиабатический режимы.
В
адиабатических условиях тепло в реактор может вводиться:
С водяным паром
С горючими газами
С предварительно нагретым твердым мелкозернистым инертным теплоносителем, который движется в реакционном пространстве вместе с потоком реагентов.
Слайд 5 В случае экзотермических процессов реактор обычно имеет керамическую насадку, которая аккумулирует
тепло, благодаря чему при подаче холодной реакционной смеси реакция не затухает.
Слайд 6 Иногда реактор представляет собой последовательность отдельных адиабатических секций, соединенных последовательно. В
промежутках между этими секциями возможно организация теплообмена через поверхность или путем введения в реакционную смесь инертных разбавителей или исходных реагентов. Если в первом случае происходит только изменение температуры, то во втором – также и состава, что должно учитываться при соответствующих расчетах.
Слайд 7 Политропические режимы для гомогенных газофазных реакций представляют собой цилиндрический или трубчатый
аппараты, обогреваемые топочными газами (хлораторы) или пламенем (печи пиролиза).
1-корпус; 2-радиантная секция; 3-перегородки; 4-горелки; 5-конвективная секция.
Слайд 8
Расчет непрерывных идеальных, изотермических гомогенных гомофазных реакторов.
Слайд 9Расчет непрерывных реакторов, как и периодических осуществляется с использованием соответствующих характеристических
уравнений.
Исходными данными для расчета служат:
Производительность по целевому продукту
Степень превращения реагента ХА
Начальные концентрации реагентов, инертных разбавителей, гомогенных катализаторов. СА0
Слайд 10 В изотермических условиях для определения объема реактора достаточно решения характеристического уравнения,
полученного при составлении материального баланса только по определенному реагенту А.
Слайд 11 Если возникает задача нахождения состава реакционной смеси на выходе из реактора,
то в этом случае число уравнений возрастает до числа ключевых реагентов в системе.
FB =
Слайд 12 Переход от FB к FA0 проводится по известным уравнениям с использованием
ХА для простых или ХА и ФB для сложных реакций. После определения FA0 объем реактора будет являться уже только функцией параметра, характеризующей глубину протекания реакций, например ХА .
Следует отметить, что это справедливо для всех жидкофазных, а так же газофазных реакций, если перепадом давления в реакторе от его входа до выхода можно пренебречь.
Слайд 13 В общем случае прохождение газового потока через реактор может сопровождаться заметными
потерями давления.
Методика расчета потери давления в трубчатых пустотелых реакторах та же, что и в трубопроводах. В расчетах принимают, что потеря давления равно потерям при трении, пренебрегая остальными составляющими уравнения баланса механической энергии.
Слайд 14Где: P – общее давление в реакторе
λ- коэффициент трения
-плотность реакционной массы
- линейная скорость потока
- эквивалентный диаметр канала, по которому движется поток (для труб крупного сечения ).
Для ламинарного потока
Тогда для элементарной длинны реактора можно записать:
Слайд 15
При турбулентном движении находятся пути обработки экспериментальных данных с использованием теории
подобия и является функцией числа Re и шероховатости стенок труб.
Слайд 16 Дифференциальное уравнение падения давления можно проинтегрировать, если известно, каким образом будет
меняться Л, w и p по длине реактора, в котором протекает химическая реакция. Факторами, влияющие на величину Л, w и p является температура и состав реакционной смеси, последний из которых будет определяющим.
Слайд 17 Представим линейную скорость потока и его плотность как функцию параметра, характеризующего
глубину протекания реакции, в качестве которого выберем степень превращения реагента А:
ε- коэффициент изменения объема.
Слайд 18 Если раздельно W и p изменяются при изменении ХА, то их
произведение остается постоянной величиной, равной Wopo .
Очевидно, что если считать m не зависит от состава реакционной смеси, что достаточно справедливо для газовой фазы, то число Re, а также λ не изменяются при увеличении ХА.
Слайд 19 Теперь можно записать:
Для оценочных расчетов падение давления как функции длинны реактора
можно заменить 1+εХА на её среднеарифметическое значение при ХА=0, ХА=ХАLТогда:
Где ΔР– падение давления в реакторе длиной
Объем раствора:
n-число труб круглого сечения в реакторе
Слайд 20
При точных расчетах уравнение падения давления необходимо рассматривать совместно с характеристическим
уравнением реактора идеального вытеснения.
n-число труб
Слайд 22
Зависимость Р от L находится при совместном решении характеристического уравнения и
уравнения учитывающего падение давления, как функцию ХА. Одновременно с этим определяется объем реактора.
Слайд 24
Основные модели реальных изотермических гомогенных реакторов
Слайд 25Основные модели реальных изотермических гомогенных реакторов
Для реальных моделей введено понятие: «распределение
времени пребывания элементов потока».
Часто проводя изучение времени пребывания элементов потока на «холодных», имеющие основные размеры моделируемых аппаратов. Для этого применяют метод, основанный на искусственном нанесении возмущения на систему и анализе отклика на него системы. Под системой подразумевается непрерывно работающий реактор (или его «холодная» модель), а под возмущением – введенный тем или иным способом в поток индикатор (трассер).
Слайд 26
Отклик на возмущение представляется в виде кривых зависимости концентрации трассера от
времени, которая называется кривыми или функциями отклика. При соответствующей обработке этих кривых получают функцию распределения времени пребывания элементов потока в реакторе.
Слайд 27
Наиболее распространены два следующих способа введения трассера в реактор. В первом
случае в поток, поступающий в аппарат и не содержащий трассирующего вещества, вводится некоторое количество индикатора и затем его концентрация постоянно поддерживается на входе в реактор на одном и там же уровне. Во втором случае индикатор вводится в поток одноразово в виде импульса, как например, проба в газожидкостный хроматограф. Последний способ наиболее распространен. В этом случае выходную кривую называют С-кривой, а импульсный выходной сигнал – дельта функцией или функцией Дарка.
Слайд 28
Часто выходные кривые представляют в безразмерных координатах. По оси ординат откладывается
безразмерная концентрация трассера С=с/С0, c - текущая концентрация трассера на выходе из реактора, С0 - начальная концентрация трассера в реакторе, определяемая в виде: М0/V, где – М0 начальное количество молей трассера.
Слайд 29
По оси абсцисс откладывается безразмерное время
, где –
текущее время,
– среднее время пребывания элементов потока в реакторе, определяемое из соотношения
где W– объемная скорость потока неизменной плотности.
Слайд 30 Зависимость c от называется дифференциальной,
а зависимость
от – интегральной функциями
Слайд 31 Рассмотрим вид кривых отклика реакторов идеального вытеснения и смешения при нанесении
импульсного возмущения. Введенный трассер в виде тонкой пленки равномерно распределяется на входе по всему сечению реактора идеального вытеснения и в «поршневом» режиме, не размываясь, движется к выходу.
Слайд 32 Свойство С-кривой таково, что
Теоретически импульсный ввод определенного количества трассера осуществляется
за бесконечно малый промежуток времени. Учитывая эти два обстоятельства, кривая отклика для реактора идеального вытеснения выражается прямой, уходящей в бесконечность.
Слайд 33
Когда время пребывания частиц трассера ( ) совпадает со средним
временем пребывания элементов потока ( ), то все его качества будет обнаружено на выходе из реактора. (рис1)
Слайд 34
В реакторе идеального смешения количество трассера, введенного импульсно (мгновенно), т.е. за
время много меньше, чем среднее время пребывания, перемешивается с содержимым реакционного пространства. Затем происходит «вымывание» индикатора уходящим из реактора потоком, согласно закону (рис2).
Слайд 35Это уравнение получатся при составлении материального баланса по индикатору.
Слайд 36
В реальных аппаратах наблюдается промежуточный вид кривой отклика – с более
или менее размытым максимумом.
Слайд 37
Очевидно, что время пребывания элементов потока в реакторе является случайной
величиной. Её функцию распределения (кривую отклика) относительно просто охарактеризовать некоторыми величинами, которые называются моментами распределения случайной величины. Или можно пользоваться для сопоставления распределений без сравнения соответствующих кривых. Наибольшее значение и применимость имеют первые два момента: среднее значение и дисперсия распределения.
Слайд 38
Кривые отклика чаще всего строят по точкам, в которых бывает известна
концентрация трассера нВ выходе из реактора в определенный момент времени – время отбора пробы. В этом случае среднее время пребывания определяется по формуле:
среднее время пребывания
-концентрация трассера в –пробе
-время отбора пробы.
Слайд 39
Дисперсия распределения времени пребывания может быть рассчитана по уравнению:
Слайд 40Причины отклонения от идеальности проточных реакторов
Слайд 41Причины отклонения от идеальности проточных реакторов
Время пребывания отдельных элементов потока в
проточном реакторе – это случайная величина которая меняется значительно от 0 до ∞.
Застойная зона. (часть потока не принимают участия в реакции)
Слайд 42 Причины отклонения от идеальности проточных реакторов
2) Зоны циркуляции.
1) Наличие
внутренних байпасов
Слайд 43
В реальных реакторах, даже близких к реактору идеального вытеснения может происходить
частичное перемешивание элементов потока в основном направлении. Это может возникнуть в результате молекулярной диффузии: в двух точках потока концентрации разные, в результате – движущая сила диффузии.
Наряду с молекулярной диффузией – турбулентная диффузия – хаотичное пульсации скорости относительно среднего значения в продольном направлении. Диффузия в осевом направлении происходит в следствии неравномерности поля скоростей (при ламинарном течении жидкости).
Слайд 44Модели, учитывающие неидеальность потока, упрощающие расчет реактора
Ячеечная модель – реальный аппарат
расчленяют на N последовательно соединенных ячеек идеального смешения. Суммарный объем всех ячеек равен полному объему реактора.
Число ячеек N – единственный параметр ячеечной модели.
Диффузионная модель – некоторый промежуточный случай между режимами идеального смешения и идеального вытеснения. Учет диффузии в соответствии с первым законом Фика.
Слайд 45 Коэффициент продольной диффузии - DL - единственный параметр однопараметрической дифференциальной модели.
В
теории реакторов разработаны модели позволяющие учесть неидеальность потока. Модели эти являются до некоторой степени приближенными, однако они более точно описывают реальный процесс.
Слайд 46 Модели реальных реакторов могут быть построены на основе двух подходов:
1. Первый
основан на мысленной замене реактора реального той или иной комбинацией идеальных аппаратов.
2. Второй подход имеет большее физическое обоснование – при составлении материального описания процесса стремятся учесть все реальные физические явления, происходящие в аппарате и внести их в уравнение с помощью соответствующих математических операторов.
Слайд 47
При разработке модели необходимо иметь в виду, что, как правило, теория
дает общий вид уравнения, а численные коэффициенты этих уравнений должны быть найдены экспериментально. Эти численные коэффициенты носят название параметров математической модели. Стремятся, чтобы число параметров было минимальным.
Наиболее распространены 2 модели однопараметрические: ячеечная и диффузионная.
Слайд 48Ячеечная модель.
Модуль потока с продольным перемешиванием можно представить в виде КРПС.
В ячеечной модели использован первый подход к описанию реальных реакторов, а именно: реальный аппарат мысленно расщепляют на N последовательно соединенных ячеек идеального смешения.
Суммарный объем всех ячеек равен полному объему реактора. Правомерность такой замены вытекает из сравнения КРПС с РИВ и РПС.
N=1 РИВ N=∞ для реального реактора.
Слайд 49Однопараметрическая диффузионная модель
Слайд 50 Диффузионная модель, как и ячеечная, описывает реальную гидродинамическую обстановку в проточном
реакторе как некоторый промежуточный случай между режимами идеального смешения и идеального вытеснения. При построении диффузионной модели идеального смешения учитывается неравномерность распределения параметров процесса С/V (в частности, концентрации) по объему аппарата.
Слайд 51
Но неравномерным является и распределение концентрации по длине реактора идеального вытеснения.
В отличие от модели идеального вытеснения в диффузионной модели учитывается наличие перемешивания реакционной среды в осевом направлении, вызванное различными видами диффузии. Последнее условие и легло в основу названия модели – диффузионная.
Слайд 52
Учет диффузии связан с появлением в математическом уравнении дифференциальных уравнений, т.к.
с первым законом Фика перенос вещества (изменение его количества в элементарном объеме) за счет диффузии пропорционален градиенту концентрации в направлении переноса.
Слайд 53Однопараметрическая диффузионная модель
В реальном аппарате в связи с неравномерностью распределения концентрации
в пространстве имеет место диффузионный перенос как в продольном, так и радиальном направлении. Чтобы чрезмерно не усложнять математическое описание в первом приближении считают, что в радиальном направлении распределение концентрации равномерное и диффузия имеет место только в продольном (осевом направлении).
Слайд 54
Перенос вещества за счет турбулентной диффузии Dтруб. Но разделить экспериментально различные
виды диффузии сложно, поэтому целесообразно все их объединить одним уравнением с коэффициентом продольной диффузии DL - это единственный параметр однопараметрической диффузионной модели.
В первом приближении, если длинна реактора больше его диаметра в 50÷60 раз, то отклонением от режима идеального вытеснения, вызванным продольной турбулентной диффузией, можно принебречь.
Слайд 55
Диффузионные модели, их характеристики.
Слайд 56
Под диффузие будет пониматься перенос части жидкости или газа под влиянием
градиента концентрации, независисмо от механизма этого переноса. Следовательно, в общем случае это может быть молекулярная и турбулентная диффузия; перемешивание, вызванное пульсацией потока и насадкой, если такая имеется, конвективный перенос причиной которого могут быть градиенты температуры и т.д.
Слайд 57
Рассмотрим вначале модель, в которой учитывается влияние продольной диффузии на
режим идеального вытеснения. Такая модель используется для описания работы трубчатых аппаратов, подразумевая отсутствие а них застойных зон и струйных потоков.
Слайд 58
Как и процесс молекулярной диффузии, осевое смешение в силу его
статистической природы описывается уравнением, анологичным по виду закона Фика. Однако в нем вместо коэффициента молекулярной диффузии используется эффективный коэффициент продольного переноса DL, полностью характеризующий степень осевого смешения в потоке. Коэффициент DL является функцией длинны реактора, но в технических расчетах используется его среднее значение DL.
Слайд 59
Рассмотрим реактор длинной L, в котором движется установившийся (стационарный во времени)
поток с осевым перемешиванием веществ. Запишем для элементарного объема реактора уравнение материального баланса. (повторить РИВ).
Слайд 60
Если в реакторе идеального вытеснения существует только один конвективный поток, то
теперь добавляется ещё и осевое диффузионный поток. Очевидно, в этом случае характеристическое уравнение будет отличаться от полученного для реактора идеального вытеснения на некоторое слагаемое, учитывающее вклад продольного перемешивания:
DL – эффективный коэффициент продольного переноса.
Слайд 61
Проведем некоторые преобразования этого уравнения, облегчающее его дальнейшее решение:
Окончательно
имеем:
w - линейная скорость потока.
Слайд 62
В случае жидкофазных реакций линейная скорость потока практически не меняется
по длине реактора и поэтому может быть вынесена из-под знака дифференциала. Для газофазных реакций, протекающих с изменением объема в технологических расчетах часто используют среднее значение w, легко определяемого для известной конечной степени превращения ключевого реагента А. Тогда:
Слайд 63
Решение уравнения возможно, если известны начальные условия. При l=0, учитывая
непрерывность переноса А, получим:
Для конечной величины DL это условие влечет за собой прерывистое уменьшение начальной концентрации А при L=0. Другим условием является
=0 при l=L
Слайд 64
Введем следующие обозначения:
Теперь можно записать
Слайд 65
Аналитическое решение этой системы уравнений возможно только в отдельных случаях,
например, для реакции первого порядка. В этом случае решение имеет вид:
Полученный результат справедлив для любых значений критерия Пекле ( ).
Слайд 66
В предельных случаях (
или )
разлагая в ряд этот результат решения, получаем известные расчетные формулы реактора идеального вытеснения
и реактора идеального смешения
Слайд 67 Для более сложных кинетических закономерностей, чем реакция первого порядка, применяют
численные методы решения записанной ранее системы уравнений.
Вообще её любое решение предполагает, что численное значение критерия Пекле известно. Рассмотрим метода его нахождения.
Слайд 68 В случае импульсной формы возмущения при подаче трассера в поток вытеснения
с продольной диффузией, получают семейство С – кривых, вид которых зависит от интенсивности осевого перемешивания.
Функция отклика на импульсный ввод трассера в реактор диффузионной модели.
Слайд 69 а)
Найдена однозначная связь между дисперсией распределения и критерием PeL в
зависимости от условий на границе аппаратов.
Различают «закрытые», «открытые», и «частично закрытые» типы аппаратов.
В «закрытых» реакторах ввод трассера осуществляется на входе в аппарат, а вывод - на выходе. В «открытых» реакторах ввод и вывод трассера осуществляются на участках реакционной трубки, удаленном на некотором расстоянии от входа и выхода аппарата. «Полузакрытый» типы аппаратов являются комбинацией первых двух.
Слайд 70 Во всех рассмотренных типах реакторов вытеснения зависимость
Вид второго слагаемого
меняется в зависимости от используемого типа реактора. Однако, чаще всего его вклад по сравнению с первым слагаемым мал, поэтому в технических расчетах пользуются соотношением:
Слайд 71 Величина критерия Пекле может быть найдена по критериальным зависимостям, которые
имеются в специальной литературе. Критериальные зависимости получены при экспериментах на лабораторных установках менее надежны, чем определение критерия Пекле на «холодных» моделях.
Слайд 72 Нахождение объема реактора для заданной производительности по целевом продукту В
В случае простых реакций, когда ФВ=1, учитывая стехиометрические коэффициенты реакции и конечную степень превращения реагента А, переходят от FB и FA0 затем, зная СА о к W. Выбирают сечение трубrи реактора S и затем принимают такое их количество n, чтобы в трубках создавался нужный гидродинамический режим. Затем находят w=W/nS и PeL критерий Пекле при этой линейной скорости потока. Задаются некоторой длиной реактора L и определяют r=L/w.
Полученное значение w, PeL, r используется для расчета САL при численном решений рассмотренной ранее системы уравнений. Полученное значение САL сравнивают с заданной величиной СА, и если они совпадают, то задача решена. В противном случае решение повторяют при другом выбранном значения длинны реактора.
Слайд 73 Обратное перемешивание уменьшает удельную производительность реактора. Для сложных реакций от
распределения времени пребывания зависит не только удельная производительность, но и селективность. Это особенно существенно в случае последовательных реакций, когда один из промежуточных продуктов – целевой. Величина максимальной селективности снижается при возникновении продольного перемешивания тем более, чем оно интенсивнее. Поэтому в случае сложных реакций при наличие в трубчатом реакторе осевого перемешивания интегральная селективность даже для изотермических условий не может быть рассчитаны, исходя из фВ и ХА.
Слайд 74 В связи с этим методику расчета объема реактора с продольным
перемешиванием для сложных реакций при заданной производительности по целевому продукту можно представить следующим образом.
Задаемся некоторым значением ФВ. Это позволяет как и для простых реакций найти w, PeL, r . Затем решатся такое число характеристических уравнений, которое позволяет определить текущие концентрации веществ A и B в реакторе.
Слайд 75 Например. Пусть имеется реакция
для которой
и
В данном примере необходимо решить совместно два характеристических уравнения с учетом их граничных условий.
Слайд 77 Правильность выбора необходимой длинны реактора осуществляется так же, как и для
простых реакций. В результате окончательного расчета находят соответствующее значение. Затем проверяют по формуле:
Насколько совпадают принятые и рассчитанные значения интегральной селективности по целевому продукту В. Если совпадение ниже требуемой точности, то расчет полностью повторяют при новом выбранном значении ФВ.