Динамическая модель реактивной периодической ректификации требует связанных балансов компонентов, общего материального баланса и энергетических балансов для каждого объёма накопления, а также уравнений химической кинетики и фазового равновесия пар–жидкость. Как минимум, это включает конденсатор, каждый объём накопления на тарелке или в насадочной секции, кипятильник и сборник дистиллята. Балансы должны учитывать потоки пара и жидкости, отбор дистиллята, флегму, фазовый переход, а также образование или расходование веществ в результате химической реакции.

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

Определите динамические состояния и контрольные объёмы

Объём накопления жидкости и состав

Для каждого контрольного объёма жидкости (j) определите общий молярный объём накопления (H_j) и мольную долю компонента (i) в жидкости (x_{j,i}). Количество компонента (i) в этом объёме накопления равно

[ n_{j,i}=H_jx_{j,i} ]

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

Состав и поток пара

Для каждой ступени определите поток пара (V_j), поток жидкости (L_j), состав пара (y_{j,i}) и состав жидкости (x_{j,i}). Эти потоки могут рассчитываться на основе тепловой нагрузки кипятильника, режима работы конденсатора, гидравлики тарелок или заданных рабочих политик, таких как постоянная производительность по пару или постоянный отбор дистиллята.

Переменные, специфичные для периодического процесса

Модель также должна отслеживать накопленный дистиллят (D_{\mathrm{cum}}) и, при необходимости, накопленный продукт или степень протекания реакции. Эти переменные необходимы для расчёта общего извлечения, выхода продукта и продолжительности периода.

Применяйте балансы компонентов по всей колонне

Общий баланс компонента для ступени

Для объёма накопления жидкости на ступени (j) динамический баланс компонента имеет вид

[ \frac{d(H_jx_{j,i})}{dt}

L_{j-1}x_{j-1,i} + V_{j+1}y_{j+1,i}

L_jx_{j,i}

V_jy_{j,i} + r_{j,i}H_j ]

где (r_{j,i}) — суммарная скорость образования компонента (i) на моль объёма накопления жидкости.

Баланс включает жидкость, поступающую с расположенной ниже ступени, пар, поступающий с расположенной выше ступени, жидкость и пар, покидающие ступень, а также образование или расходование веществ в результате реакции.

Нереакционные ступени

Для тарелок или насадочных секций, в которых реакция не протекает, примите

[ r_{j,i}=0 ]

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

Баланс компонента для кипятильника

Для объёма накопления в кипятильнике (H_N) репрезентативный баланс компонента имеет вид

[ \frac{d(H_Nx_{N,i})}{dt}

L_{N-1}x_{N-1,i}

L_Bx_{N,i}

V_Ny_{N,i} + r_{N,i}H_N ]

Здесь (L_B) — скорость отбора кубового остатка, а (V_N) — поток пара, образующегося в кипятильнике.

Для периодического куба без отбора кубового остатка (L_B=0). Если кипятильник является основным реактором, его кинетический член может определять реакционное поведение системы.

Баланс компонента для конденсатора

Для полного или частичного конденсатора с объёмом накопления жидкости (H_c)

[ \frac{d(H_cx_{c,i})}{dt}

V_2y_{2,i}

L_1x_{c,i}

D x_{c,i} + r_{c,i}H_c ]

где (V_2) — пар, поступающий в конденсатор, (L_1) — флегма, возвращаемая в колонну, а (D) — скорость отбора дистиллята.

Для полного конденсатора с пренебрежимо малой скоростью реакции (r_{c,i}=0), а состав сконденсированной жидкости обычно непосредственно связан с составом поступающего пара:

[ x_{c,i}\approx y_{2,i} ]

Если конденсатор имеет значительный объём накопления и в нём протекает реакция, его кинетический член необходимо сохранить.

Баланс сборника дистиллята

Для сборника дистиллята с объёмом накопления (H_D)

[ \frac{d(H_Dx_{D,i})}{dt}

D_{\mathrm{in}}x_{c,i}

D_{\mathrm{out}}x_{D,i} + r_{D,i}H_D ]

В обычном режиме периодической ректификации сборник часто опорожняется непрерывно или периодически. Если он используется только как ёмкость для сбора, примите (r_{D,i}=0). Баланс накопленного дистиллята имеет вид

[ \frac{dD_{\mathrm{cum}}}{dt}=D_{\mathrm{out}} ]

а накопленное количество собранного компонента (i) равно

[ \frac{dN_{D,i}}{dt}=D_{\mathrm{out}}x_{D,i} ]

Включите балансы общего количества вещества и объёмов накопления

Баланс общего объёма накопления

Общий баланс для любого объёма накопления жидкости получают суммированием балансов отдельных компонентов:

[ \frac{dH_j}{dt}

L_{j-1}+V_{j+1}-L_j-V_j + H_j\sum_i r_{j,i} ]

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

Для набора реакций (k)

[ \sum_i r_{j,i}

\sum_k \Delta \nu_k r_{j,k} ]

где (\Delta \nu_k) — суммарное изменение стехиометрического количества молей в реакции (k).

Общий баланс кипятильника

Общий баланс кипятильника имеет вид

[ \frac{dH_N}{dt}

L_{N-1}-L_B-V_N + H_N\sum_i r_{N,i} ]

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

Общий баланс конденсатора

Для конденсатора

[ \frac{dH_c}{dt}

V_2-L_1-D + H_c\sum_i r_{c,i} ]

Баланс без члена отбора дистиллята справедлив только при (D=0), например при полном орошении или в период запуска без отбора.

Общий баланс периодического процесса

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

[ \frac{d}{dt} \left( \sum_j\sum_i H_jx_{j,i} + D_{\mathrm{cum},i} + B_{\mathrm{cum},i} \right)

\text{суммарное образование компонента }i\text{ в результате реакции} ]

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

Свяжите балансы с кинетикой реакций

Выражения для скорости реакции

Для каждой реакции (k) задайте соответствующий закон скорости. Для реакции в жидкой фазе общий вид имеет следующий вид:

[ r_k=f_k(T_j,\mathbf{C}_j) ]

где (T_j) — локальная температура, а (\mathbf{C}_j) — вектор концентраций в жидкой фазе.

Для элементарной реакции пример выражения для скорости имеет вид

[ r_k=k_k(T_j)\prod_i C_{j,i}^{\alpha_{i,k}} ]

с температурной зависимостью Аррениуса:

[ k_k(T_j)=k_{0,k}\exp\left(-\frac{E_{a,k}}{RT_j}\right) ]

Скорости реакций для отдельных компонентов получают из стехиометрии:

[ r_{j,i}=\sum_k \nu_{i,k}r_{j,k} ]

где (\nu_{i,k}) имеет отрицательное значение для реагентов и положительное — для продуктов.

Фаза и расположение реакции

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

Основа расчёта концентрации и активности

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

[ a_{j,i}=\gamma_{j,i}x_{j,i} ]

где (\gamma_{j,i}) — коэффициент активности. Выбранная кинетическая основа должна соответствовать экспериментальным кинетическим данным, используемым для пилотной установки.

Добавьте энергетические балансы

Общий энергетический баланс ступени

Для каждой ступени со значительным тепловым объёмом накопления энергетический баланс можно записать следующим образом:

[ \frac{dU_j}{dt}

L_{j-1}h^L_{j-1} + V_{j+1}h^V_{j+1}

L_jh^L_j

V_jh^V_j + Q_j

Q_{j,\mathrm{loss}} + H_j\sum_k(-\Delta H_k)r_{j,k} ]

где (U_j) — внутренняя энергия, (h^L) и (h^V) — молярные энтальпии жидкости и пара, а (Q_j) — внешний тепловой поток.

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

Энергетический баланс кипятильника

Упрощённый баланс кипятильника имеет вид

[ \frac{dU_N}{dt}

L_{N-1}h^L_{N-1}

L_Bh^L_N

V_Nh^V_N + Q_R + H_N\sum_k(-\Delta H_k)r_{N,k} ]

Тепловая нагрузка кипятильника (Q_R) может быть задана, изменяться регулятором или рассчитываться на основе модели теплоносителя.

Энергетический баланс конденсатора

Для конденсатора с объёмом накопления жидкости

[ \frac{dU_c}{dt}

V_2h^V_2

L_1h^L_c

Dh^L_c + Q_c + H_c\sum_k(-\Delta H_k)r_{c,k} ]

Тепловая нагрузка конденсатора (Q_c) обычно отрицательна при соглашении о знаках, согласно которому положительное значение (Q) означает подвод тепла к системе.

Упрощения энергетического баланса

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

Замкните модель термодинамическими и гидравлическими соотношениями

Равновесие пар–жидкость

Для каждой равновесной ступени требуется соотношение между (x_{j,i}), (y_{j,i}), температурой и давлением. Распространённая форма для низких давлений имеет вид

[ y_{j,i}P

x_{j,i}\gamma_{j,i}P_i^{\mathrm{sat}}(T_j) ]

или, эквивалентно

[ y_{j,i}=K_{j,i}x_{j,i} ]

где (K_{j,i}) — коэффициент фазового равновесия.

Для неидеальных смесей может потребоваться модель коэффициентов активности, например NRTL, UNIQUAC или UNIFAC. При более высоких давлениях более подходящим может быть подход на основе уравнения состояния, например SRK.

Ограничения суммирования для ступени

Составы фаз должны удовлетворять условиям

[ \sum_i x_{j,i}=1 ]

[ \sum_i y_{j,i}=1 ]

Эти алгебраические ограничения следует соблюдать одновременно с балансами компонентов.

Соотношения для потоков и объёмов накопления

Модели также необходимы уравнения для гидравлики тарелок и скоростей потоков, например

[ H_j=g_j(L_j) ]

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

Рабочие ограничения

Уравнения должны отражать фактическую рабочую политику, включая:

Постоянную или переменную тепловую нагрузку кипятильника (Q_R(t))
Постоянную или переменную производительность по пару (V_N(t))
Политику орошения (L_1(t)) или флегмовое число (R(t))
Постоянный, регулируемый или нулевой отбор дистиллята (D(t))
Рабочее давление (P(t))
Критерии завершения периода на основе состава, конверсии или накопленного дистиллята

Понимание компромиссов

Сложность модели и идентифицируемость

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

Допущения о постоянных потоках

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

Допущения о равновесных ступенях

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

Связь реакции и разделения

Реакция изменяет состав, летучесть, плотность жидкости и, потенциально, объём накопления. Разделение изменяет концентрации, определяющие скорости реакции, поэтому независимый расчёт реакции и ректификации обычно не отражает ключевое поведение реактивной периодической ректификации.

Отсутствующая динамика сборника

Пренебрежение объёмом накопления в конденсаторе или сборнике может привести к неправильному расчёту состава дистиллята при запуске, изменениях орошения и остановке периода. Эти ёмкости следует моделировать динамически, если их объём накопления значителен по сравнению с общей загрузкой.

Выберите подходящий вариант для своей цели

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

Если ваша основная цель — предварительное моделирование: Используйте балансы компонентов и общего количества вещества для кипятильника, конденсатора и ступеней при заданных давлении, производительности по пару, орошении и политике отбора дистиллята.
Если ваша основная цель — конверсия и выход реакции: Добавьте проверенные кинетические выражения, стехиометрические скорости реакций, расположение реакционных зон и члены, учитывающие энтальпию реакции.
Если ваша основная цель — прогнозирование состава: Включите поэтапное равновесие пар–жидкость и, при необходимости, неидеальные термодинамические модели свойств.
Если ваша основная цель — разработка системы управления: Сохраните динамические объёмы накопления, энергетические балансы, ограничения исполнительных механизмов, измерения температуры и регулируемые переменные, такие как орошение и тепловая нагрузка кипятильника.
Если ваша основная цель — оптимизация продолжительности периода или экономики процесса: Отслеживайте накопленные количества дистиллята и продукта, позволяя изменять во времени орошение, тепловую нагрузку, скорость отбора и условия завершения процесса.

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

Итоговая таблица:

Элемент модели	Необходимые уравнения или соотношения	Назначение
Объёмы накопления на ступенях	Балансы компонентов и общего количества вещества	Отслеживание изменения составов и запаса жидкости
Кипятильник	Балансы компонентов, общего количества вещества и энергии	Представление образования пара, реакции, подвода тепла и отбора кубового остатка
Конденсатор	Балансы компонентов, общего количества вещества и энергии	Моделирование конденсации, флегмы, отбора дистиллята и теплового объёма накопления
Сборник дистиллята	Динамический баланс компонентов и баланс накопленного продукта	Прогнозирование состава собранного продукта, извлечения и хода периода
Кинетика реакций	Законы скорости, температурная зависимость Аррениуса и стехиометрия	Расчёт образования и расходования компонентов
Термодинамика	Равновесие пар–жидкость, коэффициенты активности, давление насыщенного пара и ограничения для фаз	Определение составов пара и фазового поведения
Гидравлика	Соотношения для потоков, объёма накопления на тарелках или в насадке и времени пребывания	Связывание внутренних потоков с запасом вещества в оборудовании
Рабочая политика	Орошение, производительность по пару, давление, тепловая нагрузка, скорость отбора и условия завершения	Представление фактической работы пилотной установки
Общая валидация	Элементные и общие материальные балансы	Обнаружение пропущенных потоков, неправильных знаков или несогласованностей модели

Уверенно моделируйте свою пилотную установку

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

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

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